c++ - Long Doubles をサポートするお気に入り/最高の AX=B ソルバーは?

Question

X=B(A^-1) の問題を繰り返し打ちたいと思っています。つまり、線形システムを解きます。C++ の場合、128 ビット長の double (quad) をサポートする数値ソルバーはどれですか?

すべての 2D データが単一の std::vector として格納されるため、C スタイルの配列を使用することは大きなメリットです。

LinuxでGCCまたはICCのいずれかを使用してコードをコンパイルしたいと考えていました。

Answer 1

多くの C++ 線形代数ライブラリは、NT2、Boost.uBLAS、Eigen などのテンプレートに基づいています (リンクについては、最も広く使用されている C++ ベクトル/行列数学/線形代数ライブラリとそのコストと利益のトレードオフを参照してください)。したがって、コンパイラ/ライブラリがクワッドを使用して計算できる場合、それらはクワッドをサポートできるはずです。たとえば、Eigen では、型Eigen::Matrix<long double, Dynamic, Dynamic>は long double を含む任意のサイズの行列を表し、標準関数を使用してそのような行列を解くことができます。