| 查看: 3177 | 回复: 50 | |||
| 本帖产生 1 个 数学EPI ，点击这里进行查看 | |||
[交流]
【讨论】预调件共轭梯度法（PCG） 已有4人参与
|
|||
|
有限元计算经常碰到大型稀疏矩阵，由于此类线性方程组通常条件数是比较大的，方程组的性态不好，所以最好用迭代方法求解，比方说是预调件共轭梯度法，但此方法在选择预调件矩阵时似乎没有一个同一的标准，大多推荐的是采用incomplete LU decomposition做为预调件矩阵。incomplete LU decomposition的计算方法似乎又有很多种。 1. incomplete LU decomposition 的计算时间应该比 LU decomposition要快速的多吧，不然直接用LU decomposition不就解出来了吗，又何必再来PCG迭代呢？ 2. 采用PCG方法的前提应该是系数矩阵对称、正定吧，因为其原理是一个相当于势函数的东西取极小值。那对于非正定的系数矩阵能求解吗，我构造了几个非正定的，有的似乎是能够收敛到正确结果的。 希望各位虫用解答和讨论。 |
回复此楼» 猜你喜欢
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有72人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
-
招聘 金属平磨液，抛光液研发工程师
已经有3人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于渤海大学应用化学研究生调剂
已经有13人回复
- 不知道选择调剂还是再考一年
已经有47人回复
- 面上基金项目预算里"劳务费"自动生成的 "博士后"能删掉吗?
已经有5人回复
- [2012-04-15]第12届多尺度、多功能梯度材料国际会议（12th International Symposium o
已经有56人回复
- 如果用高级氧化方法处理生活污水，都需要测什么指标呢？
已经有28人回复
- 溶胀平衡法测定交联聚合物的交联密度
已经有433人回复
- 关于 共轭梯度法 求助～～～好心人帮忙 有金币的哦
已经有11人回复
- 请教：Anglient 1200 数据分析方法
已经有12人回复
- 【求助】求四氢呋喃的除水干燥方法
已经有20人回复
- 求助N,N-二甲基乙酰胺缩二甲醇的合成方法，金币可加，还有没有毒性小一点的方法。
已经有8人回复
- 【求助】求共轭梯度法中的Armijo线性搜索程序
已经有5人回复
匿名
用户注销 (知名作家)
- 应助: 1 (幼儿园)
- 金币: 3228.2
- 散金: 24
- 红花: 1
- 帖子: 6711
- 在线: 384.6小时
- 虫号: 0
- 注册: 2009-11-17
- 专业: 大气环境与全球气候变化
27楼2011-11-14 10:10:48
匿名
用户注销 (知名作家)
- 应助: 1 (幼儿园)
- 金币: 3228.2
- 散金: 24
- 红花: 1
- 帖子: 6711
- 在线: 384.6小时
- 虫号: 0
- 注册: 2009-11-17
- 专业: 大气环境与全球气候变化
49楼2012-09-02 14:29:27
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
匿名
用户注销 (知名作家)
- 应助: 1 (幼儿园)
- 金币: 3228.2
- 散金: 24
- 红花: 1
- 帖子: 6711
- 在线: 384.6小时
- 虫号: 0
- 注册: 2009-11-17
- 专业: 大气环境与全球气候变化
5楼2011-06-14 10:37:26
saladin983
铁杆木虫 (正式写手)
- 数学EPI: 4
- 应助: 6 (幼儿园)
- 贵宾: 0.2
- 金币: 9197.9
- 红花: 3
- 帖子: 687
- 在线: 115.7小时
- 虫号: 448168
- 注册: 2007-11-01
- 专业: 计算数学与科学工程计算
6楼2011-06-15 03:32:09
匿名
用户注销 (知名作家)
- 应助: 1 (幼儿园)
- 金币: 3228.2
- 散金: 24
- 红花: 1
- 帖子: 6711
- 在线: 384.6小时
- 虫号: 0
- 注册: 2009-11-17
- 专业: 大气环境与全球气候变化
7楼2011-06-15 08:57:22
saladin983
铁杆木虫 (正式写手)
- 数学EPI: 4
- 应助: 6 (幼儿园)
- 贵宾: 0.2
- 金币: 9197.9
- 红花: 3
- 帖子: 687
- 在线: 115.7小时
- 虫号: 448168
- 注册: 2007-11-01
- 专业: 计算数学与科学工程计算
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
nono2009(金币+2): 鼓励应助 2011-07-01 07:10:32
小木虫(金币+0.5):给个红包，谢谢回帖
nono2009(金币+2): 鼓励应助 2011-07-01 07:10:32
|
通常意义下的PCG，preconditioner只要求对称正定，这里的CC(i)只要都是正数，显然P可以用作preconditioner。但是需要知道的是，对于不同的线性方程组，preconditioner只能依照系数矩阵来做具体的选择。这里为什么这么用，针对这个代码要解决的问题，你可以看看P的应用有没有达到这样的效果： (1) 使系数矩阵的特征值聚集； (2) 改善条件数。 因为P是个简单的对角矩阵，除非系数矩阵也近似于相关的对角矩阵，否则，我更偏向于(2)。 另外，你可以让P=I（单位矩阵），这样PCG就变成标准CG了，可以对比一下，这个P有没有改善你的计算结果。 |
8楼2011-06-18 00:29:42
匿名
用户注销 (知名作家)
- 应助: 1 (幼儿园)
- 金币: 3228.2
- 散金: 24
- 红花: 1
- 帖子: 6711
- 在线: 384.6小时
- 虫号: 0
- 注册: 2009-11-17
- 专业: 大气环境与全球气候变化
9楼2011-06-18 08:58:08
saladin983
铁杆木虫 (正式写手)
- 数学EPI: 4
- 应助: 6 (幼儿园)
- 贵宾: 0.2
- 金币: 9197.9
- 红花: 3
- 帖子: 687
- 在线: 115.7小时
- 虫号: 448168
- 注册: 2007-11-01
- 专业: 计算数学与科学工程计算
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
nono2009(金币+1, 数学EPI+1): 鼓励应助 2011-07-01 07:11:01
小木虫(金币+0.5):给个红包，谢谢回帖
nono2009(金币+1, 数学EPI+1): 鼓励应助 2011-07-01 07:11:01
|
注意一点，通常选择preconditioner的时候一个常规的做法就是选择一个系数矩阵的近似，这样的话，条件预优后得到的矩阵就会接近于单位矩阵。 假定使用LU分解，然后用L和U的乘积作为P，那么实际上就是就是用系数矩阵自身作为P，亦即条件预优后的矩阵严格地等于单位矩阵。这是最理想的情况，但是，事实上这么做没有任何意义，因为已经有LU分解，那么可以直接用这个分解求解原方程组，只需要简单的向前向后代入就足够，无需使用CG这样的迭代方法。 这里使用ILU（不完全的LU分解）是因为LU分解代价太大，尤其是对于大规模线性方程组而言。ILU分解得到的L和U作乘积之后可以充当系数矩阵的近似，按照之前所说的，也就是P的一种理想选择。 你可以看看ILU分解的具体做法，Yousef Saad的Iterative methods for sparse linear systems一书里面有专门的章节介绍，有空可以看看。这是本很好的书，在作者的个人主页可以直接下载（http://www-users.cs.umn.edu/~saad/PS/all_pdf.zip）。 |
10楼2011-06-18 15:21:12
