我处理的是N = 10 ^ 7 x N = 10 ^ 7个矩阵;非零元素的数量约为6 x N.(这些元素围绕对角线分组。)我的RAM有16 Gbt大小;所以我显然需要稀疏矩阵求解器。我运行Ubuntu LINUX,并使用fortran90(gfortran),或者说正好使用ratfor90。
我有LAPACK,但它似乎不支持稀疏矩阵求解。 (我错了吗?)MATLAB一定很好,但我不想花太多时间熟悉它;时间紧迫。我安装了旧的/金色slatec并将其用于规格。功能;它有稀疏矩阵例程吗?
我听说过ARPACK,但它可以用作普通解算器吗?它可以从gfortran调用吗?
还有其他建议吗?
谢谢, - Alex
答案 0 :(得分:0)
你是对的。 Lapack不适用于此问题。
直接稀疏求解器由MUMPS,UMFPACK,SuperLU库提供。 此外,PETSc是一个图书馆馆藏,您可以在其中找到大量信息
您可以找到适用于所有这些库的Ubuntu软件包。
答案 1 :(得分:0)
ARPACK是一个解决特征值问题的包,但它本身并不是解算器。 我不确定你能否在16 Gb上解决你的问题。我建议看一下freefem ++