精确算术中的稀疏非平方线性方程组

计算科学 线性代数
2021-12-21 13:26:12

精确求解大型稀疏线性方程组的最著名算法是什么?我正在研究的系统不是对称的,不是正定整数。唯一的好处是稀疏。我还需要指出,矩阵不是正方形的。维度为 m×n,一般不会低估或高估。

2个回答

Krylov 迭代方法是通常的选择。

如果您碰巧可以访问 Mathematica,它提供了一种测试不同方法的好方法:如果 A 是您的矩阵,请编写 B=SparseArray[A]; 然后使用带有 Method->"Krylov" 的 LinearSolve 函数。您还可以测试以查看保留整数位是否有优势。转换为实数可能会产生更快的结果,但可能会以准确性为代价。