c++ - c++ 中是否有任何免费的 ITERATIVE 线性系统求解器允许我输入任意初始猜测？

Question

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

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

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

非常感谢!

他

Answer 1

如果您在 IML++ (http://math.nist.gov/iml++/cg.h.txt) 中查看 Conjugate Gradient 的标题，您会发现您可以非常轻松地为解决方案提供初始猜测您希望获得解决方案的变量非常多。