我正在尝试执行以下操作，并重复直到收敛:其中每个Xi是nxp,还有r他们中的一个rxnxp名为samples的数组.U是nxn,V是pxp.(我得到了matrixnormaldistribution的MLE。)尺寸都可能很大；我期待的事情至少在r=200的顺序上,n=1000,p=1000.我当前的代码可以V=np.einsum('aji,jk,akl->il',samples,np.linalg.inv(U)/(r*n),samples)U=np.einsum('aij,jk,alk->il',samples,np.linalg.inv(V)/(r*p),samples)这没问题，但