質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
85.50%
NumPy

NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。

Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

Q&A

解決済

2回答

16879閲覧

PythonのNumpy配列が関数で処理できません

Sunset_Yuhi

総合スコア15

NumPy

NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。

Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

0グッド

0クリップ

投稿2016/07/09 06:08

編集2016/07/09 06:48

###前提・実現したいこと
PythonのNumpy配列を用いた行列を,関数で編集したいのですが,上手くいきません.
関数の使い方が間違っていると思うのですが,調べてもよく分かりませんでした.
(Pythonは初心者です.多分,初歩的なミスだと思うのですが)

###該当のソースコード
零行列A_gloのある行に対し,対角成分は1,それ以外の成分は10にするプログラムです.
零ベクトルB_gloも少し編集します.
関数を使わないNumpy配列の編集は簡単に出来ました.

import numpy as np #NumPyライブラリ def Edit(NodeNum, Dirichlet, Neumann): if (Dirichlet!=9999): #Dirichlet=9999の時は処理しない A_glo[NodeNum][:] = 10.0 #行を全て10にする A_glo[NodeNum][NodeNum] = 1.0 #対角成分は1にする B_glo[NodeNum][0] = Dirichlet #任意の値で固定 if (Neumann!=9999): #Neumann=9999の時は処理しない B_glo[NodeNum][0] += Neumann #任意の傾きで固定 Num = 11 #行列の次元 A_glo = np.zeros((Num,Num), dtype=np.float64) #正方行列 B_glo = np.zeros(Num, dtype=np.float64) #列ベクトル #行列の確認 print('Pre matrix') for i in range(Num): for j in range(Num): print(format(A_glo[i][j], '6.2f'), end='') print(';', format(B_glo[i], '6.2f')) #行列を編集 Edit(0, -10, 9999) #左端は固定境界 Edit(Num-1, 9999, 10) #右端は自由境界 #行列の確認 print('Post matrix') for i in range(Num): for j in range(Num): print(format(A_glo[i][j], '6.2f'), end='') print(';', format(B_glo[i], '6.2f'))

###発生している問題・エラーメッセージ
「Numpyオブジェクトには値を割り当て出来ません」と出てきます.
そんな訳無いだろうと思うのですが……

Pre matrix 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00; 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00; 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00; 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00; 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00; 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00; 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00; 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00; 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00; 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00; 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00; 0.00 Traceback (most recent call last): File "FEM1Da.py", line 28, in <module> Edit(0, -10, 9999) #左端は固定境界 File "FEM1Da.py", line 8, in Edit B_glo[NodeNum][0] = Dirichlet #任意の値で固定 TypeError: 'numpy.float64' object does not support item assignment

###補足情報
Pythonの環境はAnaconda4.1で入れました.
Pythonは3.5.2,Numpyは1.11.1です.

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

回答2

0

B_gloは1次元配列になっているので、
B_glo[NodeNum][0]のように2次元配列としてアクセスしている箇所に問題があります。

投稿2016/07/09 06:57

argius

総合スコア9388

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

Sunset_Yuhi

2016/07/09 07:04

ご回答ありがとうございます.解決いたしました.
guest

0

自己解決

申し訳ありません.自分で解決しました.
関数側でB_glo配列の次元を間違えていました.以下の関数で動きます.

def Edit(NodeNum, Dirichlet, Neumann): if (Dirichlet!=9999): #Dirichlet=9999の時は処理しない A_glo[NodeNum][:] = 10.0 #行を全て10にする A_glo[NodeNum][NodeNum] = 1.0 #対角成分は1にする B_glo[NodeNum] = Dirichlet #任意の値で固定 ← 変更箇所 if (Neumann!=9999): #Neumann=9999の時は処理しない B_glo[NodeNum] += Neumann #任意の傾きで固定 ← 変更箇所

つまらないミスでお騒がせいたしました.

投稿2016/07/09 07:02

編集2016/07/20 20:49
Sunset_Yuhi

総合スコア15

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.50%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問