Jacobi预处理器不减少条件数?

计算科学 线性代数 预处理 条件数
2021-12-14 15:06:02

假设您有一个通用矩阵A, 对角线条目aii=d>0. (不对非对角元素进行任何假设。)然后 Jacobi 预处理不会改善条件数,因为

κ(A)=|λmax(A)λmin(A)|
变成
κ(M1A)=|λmax(A)/dλmin(A)/d|=κ(A).
在这种情况下,使用 Jacobi 预处理器可以做任何有用的事情吗?

1个回答

不,雅可比只纠正相对比例。它对“平滑”病态处理没有任何作用,例如κ(A)O(h2)二阶椭圆问题的渐近。

如果您使用的是 Krylov 方法,则会自动校正全局比例,但是对于固定迭代,需要以某种方式进行(恒定)缩放(可能只是在固定方法的阻尼参数中)。