やりたいことは時間と位置座標の値が入った配列でそれぞれの差分を取って時間ごとの成長率を見たいのですが、具体的には
data1=[t,x,y,z]
data2=[t,x,y,z]
があって、これの差分を取ればいいと思い、
l=data1.shape[0]
li = data1.shape[1]
lj = data1.shape[2]
lk = data1.shape[3]
lists = []
s=0
for t in range (0,l):
for i in range (0,li):
for i in range (0,li):
for j in range (0,lk):
s+=sum(np.power((data1[t,i,yn,j]-data2[t,i,yn,j])*2))
lists.insert(s)
print(lists)
としたのですが、
s+=sum(np.power((data1[t,i,yn,j]-data2[t,i,yn,j])*2))
ValueError: invalid number of arguments
とエラーをはきます。どうすれば良いのでしょうか?
また、初心者でfor文を回す方法しか思いつかなかったのですが、効率の良い方法があれば教えてください。
nump.power()
に与えるパラメータが不適切なためで、np.power(data1[t,i,yn,j]-data2[t,i,yn,j], 2)
とする必要があります。ところで、data1
やdata2
の具体的な中身を一部分で良いので追記していただけませんか?s
を0
に初期化する場所とlists
に挿入する場所はこれであっていますでしょうか? それぞれのt
について距離の二乗を計算したいということかなと思ったのですが、今のコードだと全てのt
についても和を取っているように見えます。