具有复数的任意精度线性代数 c/c++ 库答案

【问题标题】：arbitrary precision linear algebra c/c++ library with complex numbers具有复数的任意精度线性代数 c/c++ 库
【发布时间】：2013-11-23 21:48:18
【问题描述】：

我正在进行一项涉及在 4 维相空间中具有复系数的线性微分方程的研究。为了能够检查一些关于解的根的假设，我需要能够以任意精度对这些方程进行数值求解。我曾经使用mpmath Python 模块，但它运行缓慢，所以我更喜欢用 C/C++ 重写我的程序以实现最大性能。所以我有一个问题：

是否存在支持任意精度算术和复数的 C/C++ 线性代数库？我需要一些基本功能，如点积等。（其实我也需要矩阵指数，但如果需要我可以自己实现）。

我尝试将Eigen 与MPFR C++ 一起使用，但由于它不支持复数而失败（并且像complex <mpreal> 这样的构造不起作用，因为它假定基本类型是标准的浮动）。

【问题讨论】：

【解决方案1】：

尝试使用任意精度数库（例如 GMP http://gmplib.org/）和支持复数的线性代数数学库（例如 Eigen http://eigen.tuxfamily.org/）

【讨论】：

【解决方案2】：

最后，zkcm 似乎做了我想做的事。我不确定从性能角度来看它是否好（没有做任何基准测试），但至少它可以工作并提供所有必要的功能。

【讨论】：

【解决方案3】：

您可以查看uBLAS from boost。

【讨论】：