Eigen可以在其线性代数Ax = b中使用超过Double精度

Question

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

Answer 1

Eigen可以使用任意精度算术。有an MPFR module随时可用。通常，这不应该是必要的（除非您的原始值已经以高精度存储）。您应该首先尝试使用更稳定的分解。看看at the available decompositions here。哪一个最好，在很大程度上取决于您的用例。