是否有四对偶算术稀疏矩阵包？

我正在研究一些病态的大型稀疏线性方程组。我想使用双精度双精度算术或四进制双精度算术来解决它们。我知道有一个由Nakata Maho开发的名为MPACK的程序包，它可以在四对偶算术下执行数值线性代数计算。但是，它是为密集矩阵而不是稀疏矩阵设计的。您是否知道是否有四重对偶算术稀疏矩阵包？

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

您将需要一个四精度BLAS和LAPACK，PETSc可以使用以下（部分）configure命令向您提供这些精度：

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

有关更多信息，请参阅FAQ。

另外，我也同意nOOb，如果可能的话，请尝试使系统规范化，然后再切换到四精度模式。

通过正则化比通过增加数值精度更好地解决病态系统。搜索“正规化不适定”以获取详细信息。

我可以尝试Trilinos库。他们在Tpetra下建立了模板化的稀疏矩阵库（应该替代其最初的稀疏矩阵库Epetra）。您可以模板化double，complex，quad等，它们可能是PETSc旁边最多的求解器选择（直接求解和迭代求解）。

编辑：阅读评论后，Tpetra的即时实用性似乎对wrt quad精度有点可疑...

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

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

（我是工具箱的作者）。