解決すべき2*m+3 STIFF常微分方程式があります。m=1 に対して MATLAB ode15s を使用して解決しようとしましたが、問題なく動作します。しかし、方程式を解くために Sundials パッケージ CVODE を使用したいと思います。そうしようとしている間、私は「後方微分公式」とニュートン反復を使用しました。私はヤコビアンを提供せず、数値的に計算できるようにします。しかし、それは機能しておらず、エラーが表示されます:
[CVODE WARNING] CVode
Internal t = 0 and h = 0 are such that t + h = t on the next step. The solver will continue anyway.
[CVODE ERROR] CVode
At t = 0 and h = 0, the correction convergence test failed repeatedly or with :h: = hmin.
SUNDIALS_ERROR: CVode<> failed with flag -4
- CVODE は ode15s と同じ後方微分を使用していると思います。では、なぜ機能しないのですか？
- CVODE で事前条件付きの Krylov ソルバーを使用する必要がありますか?
どんな助けも楽しみにしています。ありがとうございました。