我有以下爱因斯坦和(einsum)表达式
import numpy as np
x = np.random.rand(1,8,2,8,10)
y = np.random.rand(8,10,10)
z = np.einsum('nkctv,kvw->nctw', x, y)
print (z.shape)
输出为:
(1, 2, 8, 10)
我想做的是使用一组嵌套的for循环构建等效功能,我遵循了tutorial,它显示了如何使用嵌套的for循环实现一些einsum表达式,但是没有做到。没有那么大的帮助。
有人知道怎么做吗?
答案 0 :(得分:1)
感谢@hpaulj
解决方案是:
R = np.zeros((1,2,8,13))
for n in range(1):
for c in range(2):
for t in range(8):
for w in range(13):
total = 0
# These are the variables to sum over
for v in range(13):
for k in range(8):
total += x[n,k,c,t,v]*y[k,v,w]
R[n,c,t,w] = total