演示 BiCGStab 的示例问题

计算科学 数值分析 迭代法 非线性规划
2021-12-20 06:59:22

所以我们的团队已经能够为一个类项目编写一个 BiCGStab 实现,我们想要一个潜在的示例问题来尝试它。

到目前为止,我们已经讨论了具有 Neumann 边界的一维拉普拉斯算子(因此矩阵变得非对称),但我强烈倾向于证明至少可以使用 BiCGStab 显示的更奇特的问题。

考虑到时间和技能的限制,一个稳态的问题描述将是一个想法,因为团队中没有人对时间传播有一个像样的想法,而且我们对 BiCGStab 的实现相当幼稚(缓慢)。

对于我们的求解器的演示问题的任何建议将不胜感激。

1个回答

查看矩阵市场它收集了来自不同领域(流体动力学、弹性、声学等)的大量矩阵,因此您无需自己组装系统即可提供体面的测试用例。

矩阵以 ASCII 格式提供(Matrix Market 和 Harwell-Boeing)。我不知道您的团队使用什么语言,但您可以在网站上轻松找到 C/C++、Matlab、Python 和 FORTRAN 的 I/O 例程。Mathematica 内置了这些例程(MTXHB)。