python - for 循环等效于 einsum 表达式

Question

我有以下 Einstein Sum (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 表达式，但它并没有多大帮助。

有人知道怎么做吗？

Answer 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