我正在寻找迭代线性系统求解器来计算连续变化的场。为了使仿真正常工作，我需要为每个时间步骤重新计算该字段(可能是几次)。幸运的是，我对每个时间步都有很好的初始猜测，因此最好将其输入到迭代求解器中。并且系数矩阵非常密集。
问题是我在线检查了多个迭代求解器，例如Gmm ++，IML ++，ITL，DUNE / ISTL等。它们要么用于稀疏系统，要么不提供用于输入初始猜测的界面(由于我没有时间浏览所有文档，所以我可能是错的)。
所以我有两个问题:
1在线有没有这样的c++求解器？
2由于系数矩阵可以大到成千上万，因此直接求解器会比具有真正良好初始猜测的迭代求解器更快吗？
十分感谢!
他
最佳答案
如果您在IML ++(http://math.nist.gov/iml++/cg.h.txt)中检查“共轭渐变” header ，您会发现您可以非常轻松地在非常可变的变量中提供解决方案的初始猜测您期望得到解决方案的地方。