Arpack を使用して、複雑なシフトを伴う逆シフト法で問題 \lambda Mx = Ax の固有値を計算しています。次に、znaupd e zneupd を使用します。すべて問題ないようです。Arpack で計算された固有値と固有ベクトルを取り、それらを元の問題に入れる関数を書きました。残差のノルムは、ほぼゼロまたは非常に小さいと想定されています。とにかく、かなり高い 10^-3 の誤差が得られます。問題がどこから来るのか、何か考えがありますか?
Arpack を使用して、複雑なシフトを伴う逆シフト法で問題 \lambda Mx = Ax の固有値を計算しています。次に、znaupd e zneupd を使用します。すべて問題ないようです。Arpack で計算された固有値と固有ベクトルを取り、それらを元の問題に入れる関数を書きました。残差のノルムは、ほぼゼロまたは非常に小さいと想定されています。とにかく、かなり高い 10^-3 の誤差が得られます。問題がどこから来るのか、何か考えがありますか?