本征在其线性代数Ax=b中能否使用双精度以上?

Question

我正在用Eigen做一些线性代数。即在Ax=b中求出x。其中A为(n,n)矩阵，x和b为(n,1)大小向量。我想知道是否有可能以超过两倍的精度计算它们

Answer 1

Eigen可以使用任意精度的算法。有一个MPFR模块随时可用。通常，这应该是不必要的(除非您的原始值已经以高精度存储)。您应该首先尝试使用更稳定的分解。在这里查看可用的分解。哪一个是最好的，很大程度上取决于你的用例。