从3.2版开始，PETSc 支持 gcc/gfortran 4.6 和更高版本上的稀疏四精度计算。

您将需要四精度 BLAS 和 LAPACK，PETSc 可以通过以下（部分）配置命令为您提供（以及四精度支持）：

./configure --with-precision=__float128 --download-f2cblaslapack

有关更多信息，请参阅常见问题解答。

另外，我同意nOOb，如果可能的话，在切换到四精度之前尝试对系统进行正则化。

病态系统通过正则化比通过增加数值精度更好地解决。搜索“regularisation ill-posed”以获取血腥细节。

我可以试试 Trilinos 库。他们在 Tpetra 下有模板化的稀疏矩阵库（应该替换 Epetra，他们原来的稀疏矩阵库）。您可以模板化 double、complex、quad 等，它们可能拥有仅次于 PETSc 的最大求解器选择（直接和迭代）。

编辑：阅读评论后，Tpetra 的直接用途似乎有点可疑 wrt 四精度...

MATLAB 的多精度计算工具箱支持稀疏矩阵，并专门针对四倍精度计算进行了优化。

以下是四精度稀疏求解器的时序细节： 稀疏矩阵的直接求解器

（我是工具箱的作者）。