f90的稀疏矩阵求解器

时间:2013-08-28 18:52:16

标签: matrix gfortran lapack sparse-matrix

我处理的是N = 10 ^ 7 x N = 10 ^ 7个矩阵;非零元素的数量约为6 x N.(这些元素围绕对角线分组。)我的RAM有16 Gbt大小;所以我显然需要稀疏矩阵求解器。我运行Ubuntu LINUX,并使用fortran90(gfortran),或者说正好使用ratfor90。

我有LAPACK,但它似乎不支持稀疏矩阵求解。 (我错了吗?)MATLAB一定很好,但我不想花太多时间熟悉它;时间紧迫。我安装了旧的/金色slatec并将其用于规格。功能;它有稀疏矩阵例程吗?

我听说过ARPACK,但它可以用作普通解算器吗?它可以从gfortran调用吗?

还有其他建议吗?

谢谢, - Alex

2 个答案:

答案 0 :(得分:0)

你是对的。 Lapack不适用于此问题。

直接稀疏求解器由MUMPS,UMFPACK,SuperLU库提供。 此外,PETSc是一个图书馆馆藏,您可以在其中找到大量信息

您可以找到适用于所有这些库的Ubuntu软件包。

答案 1 :(得分:0)

ARPACK是一个解决特征值问题的包,但它本身并不是解算器。 我不确定你能否在16 Gb上解决你的问题。我建议看一下freefem ++