回 帖 发 新 帖 刷新版面

主题:fortran求解matlab生成的大型稀疏矩阵如何做呢?

我现在想使用fortran  来求解matlab程序产生的大型稀疏矩阵,阶数有几十万阶甚至百万阶。

以前使用matlab求解这种矩阵,内存花费很高,30万阶的矩阵,在matlab要花费 100G左右的内存,

听群里的朋友们说,fortran里的MKL库可以解千万阶梯的矩阵方程组,所以很想 使用 fortran里的库来

求解matlab的方程组。因为我不太熟悉fortran,而时间又比较紧,导师催着完成,所以求助各位大侠

们几个问题:

1、fortran里读取 matlab产生的大型稀疏矩阵 ,并调用MKL或其他库 求解,可行吗?难度大不大?就是

fortran读取matlab的稀疏矩阵数据并调用库求解难度大吗?

2、fortran求解百万阶的矩阵花费时间内存一般多少?

3、superLU或MUMPS直接求解大型稀疏矩阵方程组的库,有人用过吗?对处理的稀疏矩阵有什么要求?


第一次来这个论文,因为搜索superLU 和MUmps 来到这个论坛,希望能得到各位大侠的帮助。

先谢谢了。

回复列表 (共1个回复)

沙发

如果是对称的,你可以一维变带宽存储,调用‘dgetrf’与dgetrs就可以了,储存方式,上面也给出了,仔细的看MKL手册

我来回复

您尚未登录,请登录后再回复。点此登录或注册