numerical-methods - 固有値で ILUT パーコンディショナーを使用して BICGSTAB 線形ソルバーを実行すると、「推定誤差」-1.#IND が発生するのはなぜですか?

Question

質問したいのですが、固有値 (C++ ライブラリ、数値線形計算を行う) を使用して線形方程式を解く場合、双共役勾配 BICGSTAB を不完全 LU 前提条件で使用しますが、結果 bg.error()("estimated error") は -1 になります。#IND、1 回の反復の後、線形方程式のサイズは 20000 です。この問題について誰か助けてくれませんか? 他の小さなサイズ、たとえば 200 も試してみましたが、それで問題ありません。