24小时热门版块排行榜    

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

yanrenyi

铜虫 (初入文坛)

[交流] 【求助】多重网格预处理子如何构造? 已有2人参与

对于最简单的椭圆方程
-a u_{xx}-u_{yy}=f,
其中a是一常数。

我想用Krylov子空间法求解上述方程离散后的线性系统,在求解之前用多重网格进行预处理,请问这个预处理子怎么构造?
非常感谢!
回复此楼

» 收录本帖的淘帖专辑推荐

机械科研经验专辑

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yanrenyi

铜虫 (初入文坛)

引用回帖:
Originally posted by saladin983 at 2010-05-31 20:34:08:


所谓多重网格的preconditioner,实际上就是在粗网格上求方程的近似解(比较典型的是用几步Gauss-Seidel迭代求解),然后用来作为preconditioner,实现和分析都比显式的矩阵形式的preconditioner要复杂一些。这 ...

我对多重网格的了解:读过王烈衡和许学军老师编的有限元方法的数学基础那本书,听过许老师关于该方法的讲课,知道方法的基本原理,会编写最基本的多重网格程序。再复杂的就不清楚了。

我构造的 preconditioner 类似于一种块形式的 hierarchical basis 方法,可以将我上面写的那个方程离散后的系数矩阵的条件数(系数epsilon非常小时)降低到 O(1)。
5楼2010-06-02 16:13:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

saladin983

铁杆木虫 (正式写手)

yanrenyi(金币+2):感谢! 2010-05-31 17:05:47
预处理子?你是说preconditioner?建议把线性方程贴上来。

[ Last edited by saladin983 on 2010-5-28 at 18:34 ]
2楼2010-05-28 18:09:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yanrenyi

铜虫 (初入文坛)

引用回帖:
Originally posted by saladin983 at 2010-05-28 18:09:15:
预处理子?你是说preconditioner?建议把线性方程贴上来。

[ Last edited by saladin983 on 2010-5-28 at 18:34 ]

就是 preconditioner, 考虑的是如下的各项异性方程

我自己设计了一个preconditioner, 想和多重网格preconditioner的预处理效果做一个比较。即将该方程离散(我用的有限差分(中心差分)方法),得到一线性系统Ax=b, 然后对该系统分别用我设计的方法和多重网格进行预处理,比较预处理后的线性系统系数矩阵的条件数。我现在不会构造多重网格preconditioner, 希望您能给予帮助,非常感谢!!!
3楼2010-05-31 17:04:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

saladin983

铁杆木虫 (正式写手)


yanrenyi(金币+10): 2010-06-02 16:13:53
javeey(金币+1):谢谢解答 2010-06-02 20:40:39
引用回帖:
Originally posted by yanrenyi at 2010-05-31 11:04:24:



就是 preconditioner, 考虑的是如下的各项异性方程

我自己设计了一个preconditioner, 想和多重网格preconditioner的预处理效果做一 ...

所谓多重网格的preconditioner,实际上就是在粗网格上求方程的近似解(比较典型的是用几步Gauss-Seidel迭代求解),然后用来作为preconditioner,实现和分析都比显式的矩阵形式的preconditioner要复杂一些。这一块我不敢说太熟悉,也不知道你了解多少,不好展开说。推荐一本参考书吧,P. Wesseling的 an introduction to multigrid methods,不厚,我手头有本1999年版的,貌似2004年又有一版。多说一句,用多重网格的方法直接构造solver或者用它构造preconditioner本质上区别不大。

此外,我不知道你自己的算法是什么样子的,你也可以大致说一下。如果你只关心条件数,那多重网格或许是很好的选择,但是通常情况下我们还会考虑preconditioner应用的计算代价。当然,如果你的问题尺度不大,那就无所谓了。

[ Last edited by saladin983 on 2010-5-31 at 15:03 ]
4楼2010-05-31 20:34:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 366求调剂 +9 不知名的小卅 2026-04-11 9/450 2026-04-13 01:19 by 幸免 ..
[考研] 材料工程281还有调剂机会吗 +29 xaw. 2026-04-11 30/1500 2026-04-13 00:20 by 1005715100
[考研] 材料考研调剂 +27 云木达达 2026-04-11 29/1450 2026-04-12 22:48 by kyle12138
[考研] 211本科材料化工求调剂 +16 YHLAH 2026-04-11 19/950 2026-04-12 22:43 by kyle12138
[考研] 2本,初试303,0860求调剂 +6 floriea 2026-04-12 8/400 2026-04-12 18:13 by zhouxiaoyu
[基金申请] 山东省基金2026 +5 jerry681 2026-04-08 6/300 2026-04-12 08:33 by kudofaye
[考研] 280求调剂 +13 wzzz王 2026-04-09 13/650 2026-04-12 00:31 by 勇攀高峰0126
[考研] 材料工程日语考生求调剂 +7 0856?调剂 2026-04-10 7/350 2026-04-11 21:33 by 蓝云思雨
[考研] 求调剂 +10 璃茉一定上岸 2026-04-10 10/500 2026-04-11 13:31 by 1005715100
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +17 努力奋斗112 2026-04-06 20/1000 2026-04-11 00:31 by wangjihu
[考研] 中科院总分315求调剂 +8 lallalh 2026-04-09 8/400 2026-04-10 19:30 by dick_runner
[考研] 085404 298分求调剂 +10 呼啦呼啦呼呼呼 2026-04-10 11/550 2026-04-10 16:44 by wangy0907
[考研] 生物学调剂,一志愿西南大学348,Top期刊一区二作、二区三作,三等奖学金三次 +4 candyyyi 2026-04-09 4/200 2026-04-09 18:39 by l_paradox
[考研] 材料307分求大佬组收留 +17 Hll胡 2026-04-07 17/850 2026-04-09 10:53 by liuhuiying09
[考研] 331求调剂 +5 luoxin0706. 2026-04-08 5/250 2026-04-08 22:15 by zhouyuwinner
[考研] 一志愿吉大化学327求调剂 +12 王王白石 2026-04-06 13/650 2026-04-08 16:05 by luoyongfeng
[考研] 化工学硕 285求调剂 +26 Wisjxn 2026-04-07 26/1300 2026-04-08 14:42 by screening
[考研] 一志愿西电085401求调剂 +4 sunw1306 2026-04-07 4/200 2026-04-07 16:40 by 啵啵啵0119
[考研] 307求调剂 +3 Youth@@ 2026-04-07 3/150 2026-04-07 09:25 by 小黑不怕难
[考研] 362求调剂一志愿中国石油大学 +4 我要考大 2026-04-06 6/300 2026-04-06 14:11 by 无际的草原
信息提示
请填处理意见