24小时热门版块排行榜    

查看: 1446  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

正版马甲

银虫 (正式写手)

[求助] 大型稀疏矩阵的预条件 已有1人参与

请问各位大侠,
对于非对称非正定的大型稀疏矩阵有什么好的预条件处理方法吗?
矩阵非常病态,例子中的条件数达到1600。
我用广义极小残余法能够求得近似解,但是效率太低了。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

正版马甲

银虫 (正式写手)

引用回帖:
2楼: Originally posted by feixiaolin at 2014-11-22 20:59:53
http://muchong.com/html/201108/3504902.html

谢谢关注,方法我已经确定下来了,
但是由于A矩阵太病态了,
所以直接求解不经济。
并且还是一个多右边项MRHS问题。
所以想请教下有没有好的预条件方法。
4楼2014-11-23 17:11:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

feixiaolin

荣誉版主 (文坛精英)

优秀版主

2楼2014-11-22 20:59:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

感谢参与,应助指数 +1
泛泛的谈pre-conditioning 难以得到好的效果,因为恐怕没有放之四海而皆准的办法。实用中大型稀疏矩阵经常由离散偏微分方程而来,谈谈你的应用背景也许会有帮助。比如,有限体积离散SOR经常就很有效了,规则网格可以考虑multigrid 类型的pre-conditioning,  有限元常用 ICCG。关键的问题有两个(i) 什么方程 (ii)怎么离散
3楼2014-11-23 09:12:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

正版马甲

银虫 (正式写手)

引用回帖:
3楼: Originally posted by pippi6 at 2014-11-23 09:12:01
泛泛的谈pre-conditioning 难以得到好的效果,因为恐怕没有放之四海而皆准的办法。实用中大型稀疏矩阵经常由离散偏微分方程而来,谈谈你的应用背景也许会有帮助。比如,有限体积离散SOR经常就很有效了,规则网格可以 ...

这个还真不太好讲,因为微分方程很简单。
可以直接解析得解。
最后形成的矩阵是由大量的类似方程组成的。
如果非要安个名字的话,应该叫输运方程吧。
网格很不规则,属于二维问题。
其中非零元的分布如图所示。
也许看起来形式对称,但是数据不对称。
大型稀疏矩阵的预条件
untitled.gif

5楼2014-11-23 17:20:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见