块雅可比预条件子的最优性

计算科学 线性代数 迭代法 预处理
2021-12-12 01:22:40

对于一个密集的N×N矩阵A, 是块雅可比预条件子,包括对角块的逆A最优块对角预处理器?是否存在另一个矩阵P具有相同的稀疏模式,其乘积PA在控制小特征值方面做得更好A?

1个回答

考虑矩阵 a2×2矩阵A

A=(1021)

的奇异值A是:

σ1=2+1,σ2=21

产生条件数(2-范数)κ2(A)=3+22.

如果我们考虑一个块对角雅可比预条件子J=I2块大小为 1(有效地使其成为对角线或普通 Jacobi 预处理器):

JA=I2Aκ2(JA)=κ2(A)=3+225.828

这里,I2表示大小为 2 的单位矩阵。

现在,如果我们能够找到一个矩阵P具有相同的结构J(对于我们的例子是对角线),这样κ2(PA)<κ2(JA),那么我们可以说J 不是最优块对角预条件子。

让我们考虑一下P, 如下:

P(x,y)=(x00y)
并通过一个优化器运行它,目标是最小化κ2(PA). 我能够达到κ2(PA)4.236x0.401y0.179.

因此,通过这种蛮力方法,我们能够找到P具有相同的结构J,这导致更好的 2 范数条件数。


您还可以查看数学 SE 的讨论,其中传达了这一点

最优对角缩放没有简单的关系,它可以最小化矩阵的光谱条件数,除了几个特殊情况

由于这对于对角缩放(Jacobi)是正确的,因此对于块对角线也应该是正确的。该答案还谈到了导致凸优化问题的最佳缩放的公式不确定对于一般块对角线缩放是否有类似的工作(块大小1)。