计算科学 - 反转一个巨大的稀疏算子； - 吾爱随笔录

方程中反转一个巨大的稀疏（非循环）： $Ax=y$

(λ I + β D + σ C) x = y

$(\lambda I+ \beta D+ \sigma C)x=y$ 其中 I 是单位矩阵，D 是对角矩阵，C 是循环矩阵，而、、是正常数。

λ

$\lambda$

β

$\beta$

σ

$\sigma$

并且都具有相同的 10000x10000 大小（实际上它可能更大，因为我计算处理大小为 1000 x 1000 的图像，所以我的矩阵将是 1000000x1000000）。

我使用了 scipy 求解器，并且反转仅适用于 100x100 的图像大小，这意味着 100000x100000 矩阵。