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

Question

我想解决线性方程Ax = b，对于未知矩阵x。 A和b都很大且稀疏，并且分别具有30,000×25和30,000×100,000的形状（当转换成密集时）。

我尝试过使用scipy.sparse.linalg.lsqr和scipy.sparse.linalg.lsmr，但它们都要求b密集，这在计算上非常昂贵且过于昂贵。我怎样才能做到这一点？

Answer 1

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