これまで、固有値/ベクトルを計算し、常微分方程式 (ODE) を C++ で解くために GSL ルーチンを使用してきました。
ここで、GNU 多倍精度ライブラリ (GMP) を使用する必要がありますが、GSL は double 以外のデータ型をサポートしていません。
固有値を計算し、GMP データ型 (mpf_class) と互換性のある ODE を解くための良い代替手段を知っている人はいますか?
現時点では、Numerical Recipes で提供されているルーチンを使用することを考えています。これが私のタスクに適しているかどうかは誰にもわかりませんか?
よろしくドミニク