| 查看: 781 | 回复: 1 | ||||
duinichixin金虫 (著名写手)
|
[求助]
gmres和cg程序求助
|
|
我在编写算法1程序时碰到了以下疑问： 1.gmres方法求解（alpha*I+i*W)*z=r时可以直接用z=gmres(B,r,10,1e-8,500）吗，还是需要进行 [L,U,P]=luinc(B,'0') 这样的不完全分解或LU分解？restart、tol和maxit任取？ 2.cg求解（alpha*T+W*W)*z=r时我是自己编的cg程序，有没有系统自带cg程序呢？ 3.算法里的判断eta k 和tau k就是我上面gmres和cg里的tol和ecg吗，还是得另外加？ 4.我对gmres和cg算法本身并不算太了解，还见过预处理的gmres和gmres（20），那它不包含cg过程吧？ 修改了很多回都是发散的，请教各位帮我看看，谢谢  迭代法.png  ihns.png  主程序.png  cg.png |
回复此楼» 猜你喜欢
参与限项
已经有3人回复
-
假如你的研究生提出不合理要求
已经有7人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
不自信的我
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 大型线性方程组的求解方法讨论
已经有8人回复
- 大型非对称稀疏矩阵求解
已经有9人回复
- 【讨论】预调件共轭梯度法（PCG）
已经有50人回复
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3430
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
|
ttp://www.cnblogs.com/kmliang/archive/2012/12/11/2813208.html http://www.pudn.com/downloads137/sourcecode/app/detail583463.html |
2楼2015-04-17 21:08:55