c++ - 支持长双打的最爱/最佳 AX=B 求解器？

Question

我希望反复敲击 X=B(A^-1) 的问题。也就是说，求解一个线性系统。对于 C++，哪些数值求解器支持 128 位长双精度（四边形）？

使用 C 样式数组是一个主要优点，因为我所有的 2D 数据都存储为单个 std::vector。

我希望使用 GCC 或 ICC 在 linux 上编译代码。

Answer 1

许多 C++ 线性代数库都基于模板，包括 NT2、Boost.uBLAS、Eigen（请参阅什么是最广泛使用的 C++ 向量/矩阵数学/线性代数库，以及它们的成本和收益权衡？链接）。因此，如果您的编译器/库可以使用四边形进行数学运算，它们应该能够支持四边形。例如，在 Eigen 中，类型Eigen::Matrix<long double, Dynamic, Dynamic>表示包含长双精度的任意大小的矩阵，您可以使用标准函数来求解此类矩阵。