令是来自拉普拉斯算子的一维有限差分离散化的众所周知的三对角矩阵,模板为。
系统非常大,因此必须使用迭代方法。另外,它是对称正定的,所以我可以使用 PCG 来解决它。最大的问题是:我应该设置什么预处理器?我在网上搜索过,但令人惊讶的是,我找不到任何东西。
任何参考、MatLab/C++ 代码或介于两者之间的任何东西都将受到高度赞赏。
最好的祝福
令是来自拉普拉斯算子的一维有限差分离散化的众所周知的三对角矩阵,模板为。
系统非常大,因此必须使用迭代方法。另外,它是对称正定的,所以我可以使用 PCG 来解决它。最大的问题是:我应该设置什么预处理器?我在网上搜索过,但令人惊讶的是,我找不到任何东西。
任何参考、MatLab/C++ 代码或介于两者之间的任何东西都将受到高度赞赏。
最好的祝福