我有一个方阵S(160 x 160)和一个巨大的矩阵X(160 x 250000).两者都是密集的numpy数组. 我的目标：找到Q,使得Q = inv(chol(S))* X,其中chol(S)是S的较低的cholesky因式分解. 当然,一个简单的解决方案是 cholS = scipy.linalg.cholesky( S, lower=True)scipy.li