0

I am dealing with up to N=10^7 x N=10^7 matrices; number of nonzero elements is about 6 x N. (Those elements are grouped around diagonal.) My RAM has 16 Gbt size; so I clearly need sparse matrix solver. I run Ubuntu LINUX, and use fortran90 (gfortran), or precisely speaking, ratfor90.

I have LAPACK, but it doesn't seem to support sparse matrix solving. (am I wrong with that?) MATLAB must be good, but I don't want to spend much time to get familiar with it; the time is pressing. I have old/gold slatec installed and use it for spec. functions; does it have sparse matrix routins?

I hear about ARPACK, but can it be used as a plain solver? could it be called from gfortran?

Any other suggestion?

Thanks, -- Alex

4

2 回答 2

0

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

于 2013-10-11T23:58:10.453 回答
0

你说的对。Lapack 不适用于这个问题。

直接稀疏求解器由 MUMPS、UMFPACK、SuperLU 库提供。PETSc 也是一个图书馆收藏，您可以在其中找到很多信息

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

于 2013-08-28T22:59:44.797 回答