我希望反复解决X=B(A^-1)问题。也就是说,求解一个线性系统。对于C++,哪些数值解算器支持128位长的双精度(quads)?
使用C风格的数组是一个主要优势,因为我所有的2D数据都存储在一个std::向量中。
我希望用GCC或ICC在linux上编译代码。
许多C++线性代数库都基于模板,包括NT2、Boost.uBLAS、Eigen(请参阅最广泛使用的C++向量/矩阵数学/线性代数库是什么,以及它们的成本和收益权衡?以获取链接)。因此,如果您的编译器/库可以使用quad进行数学运算,那么它们应该能够支持quad。例如,在Eigen中,类型Eigen::Matrix<long double, Dynamic, Dynamic>
表示包含长二重的任意大小的矩阵,您可以使用标准函数来求解这些矩阵。