重叠加法施瓦茨

计算科学 优化 域分解
2021-12-21 19:54:58

我正在解决用有限差分(3 点模板)离散化的一维拉普拉斯问题。我想以经典形式使用加法施瓦茨方法:

Uk+1=Uk+M1rk,

在哪里rk=FAUkM1=i=1NRiT(RiARiT)1Ri.

限制运算符Ri被定义为二进制矩阵,其中标记域的现有部分。

对于非重叠情况(Jacobi 块),方法很好地收敛。但是,如果有任何重叠,我无法获得收敛。

如果我使用受限加法施瓦茨,方法再次收敛。这种行为的原因是什么?

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