我在调试多重网格时遇到问题。该怎么办?

计算科学 有限元 有限差分 matlab 多重网格
2021-12-11 03:33:01

我花了太多时间编码和调试多重网格。虽然我显然不能发布我的所有代码,因为要求某人浏览所有代码是很愚蠢的,但有什么我应该特别注意的吗?

注意:我发布这个是因为我不想让人们浪费我的时间。这可能与计算科学并不完全相关,但我认为它对任何初学者都有帮助(因为我当然仍然是其中之一)。

2个回答

平滑器、限制器和延长器可以独立测试。本质上,取一个周期域,放入特定的傅里叶模式,看看代码的行为是否符合傅里叶分析。此外,在 Briggs、Henson 和 McCormick 的A Multigrid 教程中,您会发现一个名为诊断工具的部分(第 4 章)。这应该可以帮助您进行系统调试。

我们在这里收集了大量关于如何调试有限元代码的经验: https ://code.google.com/p/dealii/wiki/FrequentlyAskedQuestions