我有一段代码类型：
nnt = np.real(np.einsum('xa,xb,yc,yd,abcde->exy',evec,evec,evec,evec,quartic))
其中evec是（比如说）一个 L x Lnp.float32阵列，并且quartic是一个 L x L x L x L x Tnp.complex64阵列。
我发现这个例程相当慢。
我认为由于所有evec的都是相同的，可能有更快的方法吗？
提前致谢。