对于一个密集的矩阵, 是块雅可比预条件子,包括对角块的逆最优块对角预处理器?是否存在另一个矩阵具有相同的稀疏模式,其乘积在控制小特征值方面做得更好?
块雅可比预条件子的最优性
计算科学
线性代数
迭代法
预处理
2021-12-12 01:22:40
1个回答
考虑矩阵 a矩阵:
的奇异值是:
产生条件数(2-范数).
如果我们考虑一个块对角雅可比预条件子块大小为 1(有效地使其成为对角线或普通 Jacobi 预处理器):
这里,表示大小为 2 的单位矩阵。
现在,如果我们能够找到一个矩阵具有相同的结构(对于我们的例子是对角线),这样,那么我们可以说 不是最优块对角预条件子。
让我们考虑一下, 如下:
并通过一个优化器运行它,目标是最小化. 我能够达到和和.
因此,通过这种蛮力方法,我们能够找到具有相同的结构,这导致更好的 2 范数条件数。
您还可以查看数学 SE 的讨论,其中传达了这一点
最优对角缩放没有简单的关系,它可以最小化矩阵的光谱条件数,除了几个特殊情况
由于这对于对角缩放(Jacobi)是正确的,因此对于块对角线也应该是正确的。该答案还谈到了导致凸优化问题的最佳缩放的公式不确定对于一般块对角线缩放是否有类似的工作(块大小1)。
其它你可能感兴趣的问题