最喜欢/最好的A.X =支持长双打的B解算器？

Question

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

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

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

Answer 1

许多C ++线性代数库都基于模板，包括NT2，Boost.uBLAS，Eigen（链接参见What are the most widely used C++ vector/matrix math/linear algebra libraries, and their cost and benefit tradeoffs?）。因此，如果您的编译器/库可以使用四边形进行数学运算，它们应该能够支持四边形。例如，在Eigen中，类型Eigen::Matrix<long double, Dynamic, Dynamic>表示包含长双精度的任意大小的矩阵，您可以使用标准函数来解决这些矩阵。