我的问题是周期性泊松矩阵的残差的 L2 范数最初是减少的,但在一定次数的迭代后开始爆炸。矩阵越大,爆炸发生得越早。
另一方面,泊松矩阵具有齐次狄利克雷边界条件的多次迭代以指数方式收敛,没有任何明显的振荡。
作为参考,我使用的是 Wikipedia 提供的 CG算法。
我遇到的唯一收敛标准似乎是系统矩阵的正定性。什么可能导致矩阵爆炸?
我的问题是周期性泊松矩阵的残差的 L2 范数最初是减少的,但在一定次数的迭代后开始爆炸。矩阵越大,爆炸发生得越早。
另一方面,泊松矩阵具有齐次狄利克雷边界条件的多次迭代以指数方式收敛,没有任何明显的振荡。
作为参考,我使用的是 Wikipedia 提供的 CG算法。
我遇到的唯一收敛标准似乎是系统矩阵的正定性。什么可能导致矩阵爆炸?
当您的矩阵条件不佳时,通常会发生这种行为。原因:
建议的行动: