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

Question

我正在寻找迭代线性系统求解器来计算连续变化的场。为了使模拟正常工作，我需要为每个时间步重新计算字段（可能是几次）。幸运的是，我对每个时间步都有一个很好的初始猜测，所以我可以将它提供给迭代求解器。系数矩阵非常密集。

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

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

非常感谢！

他

Answer 1

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