24小时热门版块排行榜    

查看: 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过程吧?
修改了很多回都是发散的,请教各位帮我看看,谢谢

gmres和cg程序求助
迭代法.png


gmres和cg程序求助-1
ihns.png


gmres和cg程序求助-2
主程序.png


gmres和cg程序求助-3
cg.png
回复此楼

» 猜你喜欢

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

辉煌在拼搏之后
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 duinichixin 的主题更新
信息提示
请填处理意见