発生した問題
pythonのライブラリnumpyを用いて作業をしているのですが、内積を求めるのにnp.dotを用いると以下のように行列の形で返ってきてしまいます。どうすればいいのでしょうか?
np.matrixを用いて書き直すべきでしょうか?
[[ 0.16666667 0. 0. ] [ 0. 0.91666667 -0.27777778] [ 0. 0.625 -0.16666667]] C:\Users\Test\Anaconda3\lib\site-packages\ipykernel_launcher.py:21: RuntimeWarning: invalid value encountered in true_divide --------------------------------------------------------------------------- ValueError Traceback (most recent call last) <ipython-input-20-9acb8a5e22ef> in <module>() 25 26 I = np.array((1,1,1)) ---> 27 power_method(I.T, A1) <ipython-input-20-9acb8a5e22ef> in power_method(x_0, A) 21 l = np.dot(x,x)/ np.dot(x,pre_x) 22 H.append(v) ---> 23 if np.abs(v - (pre_v)) < eps: 24 return H 25 ValueError: The truth value of an array with more than one element is ambiguous. Use a.any() or a.all()
該当のソースコード
Python
1import numpy as np 2A1 = np.array([[1,0,0], [0,5,-1], [0,3,0]]) 3A2 = np.array([[1,0,0], [0,0,-1], [0,3,0]]) 4def power_method(x_0,A): 5 eps =1e-4 6 iteration_max = 100 7 x = x_0 8 y = np.array(A * x) 9 pre_x = x 10 x = y / np.linalg.norm(y) 11 12 v = np.dot(x , x, out = None) / np.dot(x , pre_x, out = None) 13 print (v) 14 H = [] 15 16 for i in range(iteration_max): 17 pre_x = x 18 y = x * A 19 x = y / np.linalg.norm(y) 20 pre_v = v 21 l = np.dot(x,x)/ np.dot(x,pre_x) 22 H.append(v) 23 if np.abs(v - (pre_v)) < eps: 24 return H 25 26I = np.array((1,1,1)) 27power_method(I.T, A1)
試したこと
Aをかけたりする過程で型がかわってしまったりしているのかなとy = np.array(~)としてみるなどはためしてみたのですがうまくいきません…
コードのどこの部分に関して言っていて、何を期待していて期待とどう違うのかがわかりづらいので、どの部分なのか、np.dotに渡っている引数、期待しているnp.dotの返り値、実際に返っている値などを示すと回答しやすいです。