PythonのNumpyのNdarrayで、関数の引数で行列を渡して、その関数の中で代入すると値が変わってこないことにハマった。def f(x,y): y=x みたいな。y[...]=xとかy[:,:]=xとかy[:]=xで予期した効果になった。
Pythonは本当に電卓的とかグラフ表示とかScikit-RFで高周波の簡単な解析しかしてこなかった初心者なのですが、、、
ちょっとプログラムを組んでみようかと思ってハマった話です。
Numpyで2次元の行列、例えばx,yがあって、ある関数fを定義してx,yを引数にとってfの中でy=xの何か複雑な式
をしようとしておかしくなった。簡単な例として
def f(x,y):
y=x
とかしておいて
x=np.zeros((3,3))
y=np.eye(3)
f(x,y)
とするとyが0行列になってほしいじゃないですか。参照渡しなんだから。
でも結果はそのまま渡す前とidも同じで
[[1. 0. 0.]
[0. 1. 0.]
[0. 0. 1.]]
あれ?と思って調べたら(というかおそらく簡単な話なので逆に調べ方がわからなくて苦労、、、)、
関数の中で配列そのもの(要素ではなくて)に代入してしまうと、別の配列になってしまって
引数が変更されていない様子。
こういう時は参照渡しで受け取った行列をそのまま返すというのが必要で、
(要素ごとに変えるようなもの)
def f(x,y):
y[...]=x
とか、
※三点リーダーはhttps://qiita.com/yubessy/items/cc1ca4dbc3161f84285e Ellipsistというらしい。
配列の全要素を表す。
def f(x,y):
y[:,:]=x
か、後ろの:が省略できるので
def f(x,y):
y[:]=x
とすると、idも変わらず
[[0. 0. 0.]
[0. 0. 0.]
[0. 0. 0.]]
とできた。
まあこれはいいのだが、、、逆パターンの
y=x[:]
はリストとNumpyの配列で振る舞いが違うとか、、、Pythonはなかなかよくわからん、、、
« マスクがもう無くなるので楽天で100枚入りのを昨日の夜注文したら今日の午前中に届いた。 | トップページ | 吉野家でスタミナ超特盛丼を食す。確かに肉(牛、豚、鶏)だらけででかい、のとマヨネーズとにんにくがジャンキーな味にしている。 »
「日記・コラム・つぶやき」カテゴリの記事
- 高周波・RFニュース 2025年2月13日 iFixitがSamsung Galaxy S25 ultraを分解、ぐりぐり回せる3D CTモデルに、Spirentの5G Outlook レポート、QualcommのオンデバイスAIホワイトペーパー、InfineonがCoolSiC™ MOSFETs 650 Vの Q-DPAKとTOLLパッケージで発表(2025.02.13)
- 高周波・RFニュース 2025年2月12日 GSAが5G standaloneウェビナー開催、CTIAと5G americasがケーブルに対抗してSpectrum for Broadband Competitionで連携、NXPがエッジAIのKinaraを買収、PBKreviewsがSamsung Galaxy S25/plus/ultra全部分解していて5Gミリ波AiPは全部1個(2025.02.12)
« マスクがもう無くなるので楽天で100枚入りのを昨日の夜注文したら今日の午前中に届いた。 | トップページ | 吉野家でスタミナ超特盛丼を食す。確かに肉(牛、豚、鶏)だらけででかい、のとマヨネーズとにんにくがジャンキーな味にしている。 »
コメント