CSR矩阵的LU分解

Question

我以压缩稀疏行格式(CSR)存储稀疏矩阵A。我想计算A的LU分解。常用算法效率不高，因为它们必须循环矩阵的所有系数。是否有一种有效的算法可以利用CSR格式来计算LU分解？谢谢你的帮助!

Answer 1

我是la4j (线性代数for Java)库的作者。la4j支持CRS格式，也支持LU分解。因此，您可以尝试使用它。但是，在执行分解时，la4j (0.4.0)并不真正处理具体格式的细节。但这就是我下一步要做的事情。你可以拉取更新，或者尝试自己实现这样的算法，然后向la4j发送一个subscribe -request。

如何使用LU分解程序：

Matrix a = new CRSMatrix(new double[][]{ 
  {1.0, 2.0}, 
  {3.0, 4.0}
});

// lu[0] - L, lu[1] - U
Matrix[][] lu = a.decompose(Matrices.LU_DECOMPOSITOR);