CSR矩陣的LU分解

Question

我以壓縮稀疏行格式（CSR）存儲稀疏矩陣A.我想計算A的LU分解。通用算法效率不高，因爲它們必須在矩陣的所有係數上循環。有沒有一種有效的算法可以利用CSR格式來計算LU分解？ 謝謝你的幫助！

Answer 1

我是la4j（Linear Algebra for Java）庫的作者。 la4j支持CRS格式以及LU分解。所以，你可以嘗試使用它。但是，la4j（0.4.0）在處理分解時並不真正處理結構體格式細節。但是這正是我計劃在下一次發佈中做的。您可以subscribe更新或嘗試自行實施此類算法，並向la4j發送拉取請求。

如何使用LU decompositor：

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);