右预处理和定点线性迭代

计算科学 线性求解器 预处理 固定点
2021-11-30 19:09:12

给定一个线性系统,我们可以将其表达为更容易求解的右预条件形式:Ax=b

AM1y=b,y=Mx

另一方面,左预条件系统是:

M1Ax=M1b

从文献中我们知道很容易证明上面的左预处理线性系统等价于定点线性迭代:

x(k+1)=(IM1A)x(k)+M1b

或等效地:

x(k+1)=x(k)+M1r(k)

我的问题是:我怎样才能得到一个像从正确的预处理公式开始的更新方案?

1个回答

从...开始

AM1y=b,其中yMx

我们可以操纵

y(IAM1)y=b

的一个实例,\替换另一个实例,得到更新方程yykyk+1

yk+1=(IAM1)yk+b

也可以写成

yk+1=yk+(bAM1yk)