24小时热门版块排行榜

返回列表

duinichixin

金虫 (著名写手)

应助: 1 (幼儿园)
金币: 892
散金: 596
红花: 4
帖子: 1009
在线: 338.5小时
虫号: 2475163
注册: 2013-05-18
性别: GG
专业: 计算数学与科学工程计算

[求助] 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