我正在做一个python项目,并使用numpy。我经常需要通过单位矩阵来计算矩阵的Kronecker乘积。这些都是我代码中的一个很大的瓶颈,所以我想对它们进行优化。有两种产品我必须带上。第一个是:
np.kron(np.eye(N), A)
只需使用scipy.linalg.block_diag就可以很容易地优化这一点。该产品相当于:
la.block_diag(*[A]*N)
这大约快了10倍。然而,我不确定如何优化第二种产品:
np.kron(A, np.eye(N))
有没有类似的技巧我可以使用?
https://stackoverflow.com/questions/44461658
复制相似问题