具有大而稀疏矩阵的线性代数

Question

我想为未知矩阵x求解线性方程Ax = b。 A和b既大又稀疏，形状分别为30,000 x 25和30,000 x 100,000（转换为密实形状）。

我尝试同时使用scipy.sparse.linalg.lsqr和scipy.sparse.linalg.lsmr，但是它们都要求b是致密的，这在计算上非常昂贵且禁止。 我该怎么办？

Answer 1

您可以使用numpy.linalg.pinv查找“ x”值