主题:求解大型稀疏方程组
szr1215
[专家分:110] 发布于 2010-06-10 18:04:00
现有一大型稀疏方程组,系数疏矩阵主对角不为零,不正定、主对角也不占优,有谁能快速求解,我用的消去法求解的,慢死了。。。谢谢
回复列表 (共2个回复)
沙发
pasuka [专家分:670] 发布于 2010-06-10 19:47:00
1、大型稀疏矩阵?具体多大?10w×10w以下的现在不能算做大型稀疏矩阵;
2、快速求解的话,迭代法如果选择适当,速度是最快的,但是不能保证所有问题都速度快;
3、直接分解是比较稳妥的办法,一般学术界公认的开源且速度比较快的稀疏矩阵求解器是PARDISO
安装ivf专业版就自带pardiso,具体使用参看帮助文档
此外还有UMFPACK、SuperLU等等开源的求解器可供选择
又,可以导入到matlab中,形成稀疏矩阵,利用UMFPACK求解(已经被matlab集成),速度不会比自己直接调用上述开源包慢一个数量级的。一般4G内存、64位Matlab下,计算10w×10w以下的稀疏矩阵,没啥大问题
板凳
szr1215 [专家分:110] 发布于 2010-06-10 22:46:00
谢谢1楼的回帖,我用的稀疏矩阵的系数压缩到一维向量中(去掉所有0元素),而不是在二维的向量中,不知道pardiso还能不能解决。楼主有pardiso的开源代码吗???我的qq425288059,有空详聊(明天可能出差)。。。再次感谢1楼
我来回复