c++ - 支持长 double 的最喜欢/最佳 A.X=B 求解器？

Question

我希望反复锤炼一个 X=B(A^-1) 问题。即求解一个线性系统。对于 C++，哪些数值求解器支持 128 位长 double (四元组)？

使用 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>表示包含 long double 的任意大小的矩阵，您可以使用标准函数来求解此类矩阵。