陰的 Runge-Kutta 法における線形方程式のソルバーのベンチマーク

陰的 Runge-Kutta 法における線形方程式のソルバーのベンチマーク#

陰的 Runge-Kutta 法において,Jacobian の計算が困難な場合に方程式を解く際に使用する 線形方程式のソルバーについて,以下の解法の計算時間を測定した.

  • GMRES 法 [9] [10] (numerical-collection-cpp リポジトリ [1] の実装)

  • BiCGstab 法 [9] (numerical-collection-cpp リポジトリ [1] の実装)

結果#

Hide code cell source

from num_anal_plots.show_plot_in_jupyter import show_plot_in_jupyter

show_plot_in_jupyter("ode-runge-kutta-sparse-solvers", version=5)

環境#

  • CPU: Intel(R) Core(TM) i7-9750H CPU @ 2.60GHz

  • コンパイラ:Clang 21.1.8

ソースコード#

numerical-collection-cpp リポジトリ [1] のコミット aada5a2ffd442f3ce90121ac78fac102a17b602d 時点のものを使用した.