2

初期条件が異なる ODE の複数の解を平均するので、すべての解が同時に値を持つことが重要です。たとえば、増分は 0.01 です。

数値レシピ 3 (nr3) の ODE ルーチンを使用しています。それらは適応サイズステップを行い、計算された値を使用して同じ順序の補間を行います。ブーストと競合するため、使用できません。他の同様のルーチンはありますか？

私はGSLを見ました、それはとてもいいですが、補間が組み込まれていません。私ができる1つの方法は、適応サイズでODEを解き、Akima補間を実行することです。しかし、nr3 ソリューションの方が高速で正確なようです。

4

2 に答える 2

1

http://www.unige.ch/~hairer/software.htmlのDOPRI5を密な出力 = 補間で使用しました。信頼できると思いました。元のバージョン (Fortran) を使用しました。同じ Web ページに C バージョンもありますが、私自身は使用していませんが、人々が満足していたことを覚えているようです。

于 2011-10-05T08:39:11.817 に答える
1

odeintを使用できます。Dopri5、Rosenbrock4、Burlish-Stoer があり、高密度の出力が得られます。

于 2011-10-31T21:27:47.527 に答える