应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 工科数学 » 大型稀疏矩阵的预条件
71/1返回列表
查看: 1596  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

正版马甲

银虫 (正式写手)

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

请问各位大侠，
对于非对称非正定的大型稀疏矩阵有什么好的预条件处理方法吗？
矩阵非常病态，例子中的条件数达到1600。
我用广义极小残余法能够求得近似解，但是效率太低了。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-11-22 16:28:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
http://muchong.com/html/201108/3504902.html
赞一下 回复此楼
2楼2014-11-22 20:59:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

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

正版马甲

银虫 (正式写手)
引用回帖:
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的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

引用回帖:
5楼: Originally posted by 正版马甲 at 2014-11-23 17:20:48
这个还真不太好讲，因为微分方程很简单。
可以直接解析得解。
最后形成的矩阵是由大量的类似方程组成的。
如果非要安个名字的话，应该叫输运方程吧。
网格很不规则，属于二维问题。
其中非零元的分布如图所示 ...

输运方程应该是椭圆2阶方程。但如果输运系数很小，流速很大，需要网格很小，是会有点病态。但二维问题，迭代收敛一般不会太差，虽然我没认真算过条件数（好像没听说过实际方程谁真算条件数的）。
赞一下 回复此楼
6楼2014-11-23 19:06:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

正版马甲

银虫 (正式写手)
引用回帖:
6楼: Originally posted by pippi6 at 2014-11-23 19:06:20
输运方程应该是椭圆2阶方程。但如果输运系数很小，流速很大，需要网格很小，是会有点病态。但二维问题，迭代收敛一般不会太差，虽然我没认真算过条件数（好像没听说过实际方程谁真算条件数的）。...

这个与严格意义上的输运方程有些区别，
条件数是我用matlab计算的，
也用matlab求过特征值了，
条件数确实很大。
赞一下 回复此楼
7楼2014-11-23 20:06:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 正版马甲 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 2JOx3r2CYEgw 2026-08-21 9/450 2026-08-23 18:46 by lXgU0RvNTDC9
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 h4CP7TrQR8Lg 2026-08-22 3/150 2026-08-23 16:30 by lXgU0RvNTDC9
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 4/200 2026-08-23 15:30 by KM5EcsNQRBPn
[公派出国] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 6/300 2026-08-23 14:45 by KM5EcsNQRBPn
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +3 process2012 2026-08-23 4/200 2026-08-23 12:54 by huixian257
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 h4CP7TrQR8Lg 2026-08-22 5/250 2026-08-23 12:09 by LR9qGULyN2ew
[基金申请] 什么时候开奖？ +9 CrisMessi 2026-08-18 10/500 2026-08-23 12:03 by 丶昵称占用
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 9/450 2026-08-23 10:44 by IXZuIJ2Q7OVy
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 03:40 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 11/550 2026-08-23 02:55 by OEbVnUOu01ol
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今天基金会出结果吗？20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 今天系统多次维护，明天很可能放榜！ +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 时间戳变了，能看出什么问题？ +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
关闭
请填处理意见
关闭