c++ - 三重対角行列で一次方程式系を解くためのライブラリ?

Question

私は熱伝導を伴う物理システムをモデル化しています。数値計算を行うには、三重対角行列を使用して線形方程式系を解く必要があります。このアルゴリズムを使用して結果を取得しています: http://en.wikipedia.org/wiki/Tridiagonal_matrix_algorithmしかし、私の方法は簡単で最適ではないのではないかと心配しています。そのシステムを最速で解決するには、どの C++ ライブラリを使用する必要がありますか? また、マトリックスは頻繁に変更されないことにも言及する必要があります (方程式の右側の部分のみが変更されます)。ありがとう！

Answer 1

エイゲンをチェックしてください。

Answer 2

このアルゴリズムのパフォーマンスは、浮動小数点除算によって支配される可能性があります。SSE2 を使用して (ci_とd _iの) 2 つの除算を一度に実行すると、最適なパフォーマンスに近づくことができます。