ODE を解くときは、小数点以下 40 桁以上の精度が必要です。おそらくlibquadmathを使用して、gsl ODEソルバーを使用してより高い精度の結果を得ることが可能かどうか誰か教えてください。http://www.gnu.org/software/gsl/manual/gsl-ref.html#Driverにあるgsl_odeiv2 のドキュメントを見ると、doubleのみを使用しているように見えますが、これは本当ですか?
そうでない場合、誰かが代替案を提案できますか? 「odeint」のようなものは、倍精度よりも高い精度をサポートしていますか? ゼロから自分で書かなければなりませんか？
（参考までに、私の特定の問題は次のとおりです：Daniel Fischerのおかげでバグが修正され、gslを使用してRunge-Kutta ODEソルバーを書く）