Question

我有以下爱因斯坦和（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

等效于einsum表达式的for循环

1 个答案: