数值分析、旋转和不完全 LU 分解

计算科学 线性求解器 稀疏矩阵 预处理 矩阵分解
2021-12-12 13:10:44

在进行LU分解时,如果任何对角元素,算法都会分解xii为零。因此,我们可以在矩阵上使用旋转使得xii不再为零。那是而不是看xii我们看另一个元素x在矩阵中。

LU 分解的一个问题是当我们的矩阵A是稀疏的,我们希望保持稀疏模式A,当我们进行 LU 分解时不会发生这种情况。因此,我们可以使用不完全 LU 分解,如果xijA为零,我们也跳过xijL要么U取决于在哪里xijA.

我的问题是:我们会为不完整的 LU 分解进行旋转吗?如果我们这样做,我们该怎么做?如果我们不这样做,我们如何确保我们的算法不会出错?

0个回答
没有发现任何回复~