numerical-methods - 为什么在 eigen 中使用 ILUT perconditioner 执行 BICGSTAB 线性求解器时会出现“估计错误”-1.#IND

Question

我想问一个问题，当我使用 eigen（C++ 库，进行数值线性计算）来求解线性方程组时，我使用双共轭梯度 BICGSTAB 和 Incomplete LU preconditioner，但是，结果 bg.error()("estimated错误”）变为-1.#IND，1次迭代后，线性方程组的大小为20000，有人可以帮我解决这个问题吗？我也尝试了其他小尺寸，例如200，没关系。