我在求解稀疏线性方程 Ax = b 时遇到了一些问题

我的矩阵 A 是稀疏的，尺寸为 500 万乘 500 万。实际上，它是两个矩阵的组合。一个是三对角矩阵，另一个矩阵具有更大的频带。例如$29450^{th}$ 行具有非零对角线条目。加上非零元素$14725^{th}$,$44175^{th}$,$29451^{th}$, $29449^{th}$,$29605^{th}$,$29295^{th}$，柱子。所以每一行都有几乎 6 个非零元素加上非零对角线条目。

我尝试在 matlab 中使用几个迭代求解器，但它们似乎并没有给出很好的准确性。此外，如果我长时间运行迭代，我将有四舍五入的错误和所有。我该怎么做才能得到最好的估计？