回 帖 发 新 帖 刷新版面

主题:求解大型稀疏方程组

现有一大型稀疏方程组,系数疏矩阵主对角不为零,不正定、主对角也不占优,有谁能快速求解,我用的消去法求解的,慢死了。。。谢谢

回复列表 (共2个回复)

沙发

1、大型稀疏矩阵?具体多大?10w×10w以下的现在不能算做大型稀疏矩阵;
2、快速求解的话,迭代法如果选择适当,速度是最快的,但是不能保证所有问题都速度快;
3、直接分解是比较稳妥的办法,一般学术界公认的开源且速度比较快的稀疏矩阵求解器是PARDISO
安装ivf专业版就自带pardiso,具体使用参看帮助文档
此外还有UMFPACK、SuperLU等等开源的求解器可供选择
又,可以导入到matlab中,形成稀疏矩阵,利用UMFPACK求解(已经被matlab集成),速度不会比自己直接调用上述开源包慢一个数量级的。一般4G内存、64位Matlab下,计算10w×10w以下的稀疏矩阵,没啥大问题

板凳

谢谢1楼的回帖,我用的稀疏矩阵的系数压缩到一维向量中(去掉所有0元素),而不是在二维的向量中,不知道pardiso还能不能解决。楼主有pardiso的开源代码吗???我的qq425288059,有空详聊(明天可能出差)。。。再次感谢1楼

我来回复

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