我正在尝试求解线性方程组：。为此，我目前正在使用英特尔 MKL Pardiso 求解器。当的阶数在及以下时效果很好。在此之上，求解器给出了错误的答案（也许我给出了在调用求解器之前需要设置的错误参数值）。和的维度需要分别高达 ( ) 和 ( )。我的程序是 Fortran 固定格式，我使用的是 Intel Visual Fortran 编译器。$AX=B$$A$$13500\times13500$$A$$B$$15\cdot 10^6\times15\cdot 10^6$$15\cdot 10^6 \times 6$

使用稀疏线性求解器的原因：

a)矩阵主要是稀疏的。我附上了用于可视化的矩阵的形态，其中黄色像素代表非零值。$A$$A$

b) 我之前使用的是DGELSDLAPACK 子程序，但它花费了很多时间。稀疏求解器比非稀疏求解器快约 180 倍。此外，Pardiso 求解器的内存消耗非常少，因为它适用于矩阵的 CSR 格式。

我的问题是：如何为 Pardiso 求解器设置参数，或者是否有任何其他我可以使用的预设要求较少的求解器。

目前我在 Windows 10 上工作。系统规格：

Intel Xeon E 系列处理器 32 核，64 GB RAM。

图片信息：

$A$矩阵形态，共有元素。$135\times135$

黄色像素：666 个。矩阵尺寸为（最简单的情况）。请注意，随着尺寸的增加，矩阵结构完全保持相同的时尚，A 矩阵也是非对称的。 $135\times135$

谢谢