24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1325  |  回复: 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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0831生医工第一轮调剂失败求助 +12 小熊睿睿_s 2026-04-11 16/800 2026-04-12 16:28 by 钰璞
[考研] 313求调剂 +5 余摆摆是大哥 2026-04-11 5/250 2026-04-12 10:34 by 猪会飞
[考研] 332求调剂 +14 蕉蕉123 2026-04-10 14/700 2026-04-12 00:27 by 蓝云思雨
[考研] 求调剂 +18 MAX怅惘 2026-04-09 20/1000 2026-04-11 23:31 by zhen~
[考研] 求调剂 +11 月@163.com 2026-04-07 13/650 2026-04-11 22:55 by BruceLiu320
[硕博家园] 有没有学校材料专业收跨调(一志愿085410) +6 momo(上岸版) 2026-04-06 9/450 2026-04-11 22:38 by wj165256
[考研] 296求调剂 +6 汪!?! 2026-04-09 6/300 2026-04-11 11:25 by zhq0425
[考研] 生物学调剂 +8 小冉要努力 2026-04-10 9/450 2026-04-11 10:22 by wwj2530616
[考研] 农业管理302分求调剂 +3 xuening1 2026-04-10 3/150 2026-04-11 10:18 by zhq0425
[考研] 一志愿京区985,085401电子信息,本科电子信息 +3 阳光开朗的男孩 2026-04-10 3/150 2026-04-10 16:29 by sophia_93
[考研] 青岛科技大学材料学院,环境学院调剂补录4月10日以前都可以 +3 1青科大。 2026-04-09 5/250 2026-04-10 09:58 by 翩翩一书生
[考研] 347材料专硕求调剂 +19 zj8215216 2026-04-06 19/950 2026-04-10 09:36 by 690616278
[考研] 08工学 309分求调剂 +6 Yin DY 2026-04-08 6/300 2026-04-10 09:18 by Delta2012
[考研] 材料专硕初试分332一志愿西北工业大学, +12 故人?? 2026-04-09 12/600 2026-04-09 18:34 by Ccclqqq
[考研] 085400电子信息类(川大控制工程)求调剂可跨专业 求老师联系 +3 626776879 2026-04-08 3/150 2026-04-09 16:05 by 猪会飞
[考研] 328求调剂 +17 lftmya 2026-04-07 18/900 2026-04-09 08:05 by 5268321
[考研] 296求调剂 +3 汪!?! 2026-04-08 3/150 2026-04-08 22:00 by zhouyuwinner
[考研] 293分求调剂,外语为俄语 +7 加一一九 2026-04-07 10/500 2026-04-08 20:14 by yutian743
[考研] 305求调剂 +4 77Qi 2026-04-06 4/200 2026-04-07 20:06 by shanqishi
[考研] 生物学363调剂求助 +7 fanzhang6666 2026-04-06 9/450 2026-04-07 17:37 by lijunpoly
信息提示
请填处理意见