| 查看: 2460 | 回复: 50 | |||
| 本帖产生 1 个 数学EPI ，点击这里进行查看 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【讨论】预调件共轭梯度法（PCG） 已有4人参与
|
|||
|
有限元计算经常碰到大型稀疏矩阵，由于此类线性方程组通常条件数是比较大的，方程组的性态不好，所以最好用迭代方法求解，比方说是预调件共轭梯度法，但此方法在选择预调件矩阵时似乎没有一个同一的标准，大多推荐的是采用incomplete LU decomposition做为预调件矩阵。incomplete LU decomposition的计算方法似乎又有很多种。 1. incomplete LU decomposition 的计算时间应该比 LU decomposition要快速的多吧，不然直接用LU decomposition不就解出来了吗，又何必再来PCG迭代呢？ 2. 采用PCG方法的前提应该是系数矩阵对称、正定吧，因为其原理是一个相当于势函数的东西取极小值。那对于非正定的系数矩阵能求解吗，我构造了几个非正定的，有的似乎是能够收敛到正确结果的。 希望各位虫用解答和讨论。 |
回复此楼» 猜你喜欢
AI 太可怕了，写基金时，提出想法，直接生成的文字比自己想得深远，还有科学性
已经有5人回复
-
有院领导为了换新车，用横向课题经费买了俩车
已经有9人回复
-
酰胺脱乙酰基
已经有13人回复
-
博士延得我，科研能力直往上蹿
已经有8人回复
-
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗
已经有4人回复
-
有时候真觉得大城市人没有县城人甚至个体户幸福
已经有10人回复
-
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）
已经有5人回复
-
遇见不省心的家人很难过
已经有22人回复
saladin983
铁杆木虫 (正式写手)
- 数学EPI: 4
- 应助: 6 (幼儿园)
- 贵宾: 0.2
- 金币: 9197.9
- 红花: 3
- 帖子: 687
- 在线: 115.7小时
- 虫号: 448168
- 注册: 2007-11-01
- 专业: 计算数学与科学工程计算
★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小雨萌萌: 金币+5, 3Q 2012-09-03 13:01:48
小木虫: 金币+0.5, 给个红包，谢谢回帖
小雨萌萌: 金币+5, 3Q 2012-09-03 13:01:48
|
就是你所说的先求解(5)式，得到y之后，再求解(4)式，这里因为L和U都是三角矩阵，直接求解的计算量很小，求解y和x的时候，按顺序计算各个分量即可，这里是不需要使用CG法。之前说的两个相关的方程组，就是指可以轻易地直接求解的（5）和（4）式。 事实上，P的应用，就是求解一个形如Py=r的方程组，这里也可以用其它的办法，ILU是一种比较廉价的做法，也经常见用（代数）多重网格法。还有些时候，构造出来的P可能会有特殊结构，使得Py=r非常易解，比如P是三对角矩阵。 你的理解有些偏差。CG法里应用preconditioner所需要的不一定是显式的P，算法里事实上只用到Py=r里的y，这个y总是用尽可能低廉的算法达到合适的精度，这里如何求解y会比较需要经验。CG求解的是Ax=b这个方程组，而如果利用precondioner的话，CG算法的每一步都要求解一个Py=r的线性方程组，这是两个不同的层级。在Py=r这个级别上，对于算法的高效性要求更高，而对精度的要求相对要低点，所以才会容许使用ILU分解这种近似算法。当然，你也可以直接用CG方法求解Py=r到适当的精度，但是这个时候，这里的CG法和求解Ax=b的CG法已经是一个嵌套的关系了。 |
40楼2012-08-06 00:24:28
saladin983
铁杆木虫 (正式写手)
- 数学EPI: 4
- 应助: 6 (幼儿园)
- 贵宾: 0.2
- 金币: 9197.9
- 红花: 3
- 帖子: 687
- 在线: 115.7小时
- 虫号: 448168
- 注册: 2007-11-01
- 专业: 计算数学与科学工程计算
2楼2010-09-06 22:58:01
3楼2010-09-07 14:04:33
saladin983
铁杆木虫 (正式写手)
- 数学EPI: 4
- 应助: 6 (幼儿园)
- 贵宾: 0.2
- 金币: 9197.9
- 红花: 3
- 帖子: 687
- 在线: 115.7小时
- 虫号: 448168
- 注册: 2007-11-01
- 专业: 计算数学与科学工程计算
4楼2010-09-07 16:06:00