boost - boost::numeric::odeint をスティッフ システムに使用すると、解析的なヤコビ行列を取得できない場合はどうなるでしょうか。

Question

boost::numeric::odeint rosenbrock4 インテグレーターを使用しています。問題は、解析的ヤコビ行列を取得できないことです。これを修正する解決策はありますか?

Answer 1

数値スキームを使用して近似ヤコビアンを計算できます。ただし、これは自分で実装するか、他のライブラリで実装する必要があります。odeint はこの機能を提供しません。

好奇心から: ヤコビアンが見つからないシステムは正確には何ですか? rhs が微分可能でない場合、数値に関する他の問題が発生する可能性があります。