标准化 Jacobi 预处理共轭梯度中的右侧

计算科学 线性求解器 预处理 共轭梯度
2021-12-07 21:47:13

我一直在阅读以下论文:CG 与 MINRES:实证比较其中,共轭梯度求解器应用于系统矩阵AJacobi 对两边进行预处理。具体来说,假设A是 spd,我们要解决x

Ax=b.

矩阵是 Jacobi 预处理的d=diag(A),D12=diag(1/d),但右手边也被归一化:

D12AD12y=D12bD12b.

我不明白为什么作者将右侧标准化 - 这不会改变原始系统吗?隐含考虑(这部分不在论文中)? 这种标准化有什么目的(例如更快的收敛)?x=D12bD12y

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