3

我正在寻找一个迭代线性系统求解器来计算不断变化的场。为了使模拟正常工作,我需要为每个时间步重新计算场(可能多次)。幸运的是,我对每个时间步都有一个很好的初始猜测,所以我最好将它输入到迭代求解器中。并且系数矩阵非常密集。

问题是我在线检查了几个迭代求解器,如 Gmm++、IML++、ITL、DUNE/ISTL 等。它们要么用于稀疏系统,要么不提供用于输入初始猜测的接口(我可能错了,因为我没有时间浏览所有文档)。

所以我有两个问题: 1 网上有没有这样的 c++ 求解器?2 由于系数矩阵可以大到数千 * 数千,直接求解器是否比具有非常好的初始猜测的迭代求解器更快?

太谢谢了!

4

1 回答 1

2

如果您检查 IML++ 中的共轭梯度 (http://math.nist.gov/iml++/cg.h.txt) 的标题,您会发现您可以很容易地在非常变量中提供解决方案的初始猜测你期望得到解决方案的地方。

于 2012-01-05T23:02:08.403 回答