|
|
[求助]
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]()
cg.png
![gmres和cg程序求助-3]()
主程序.png |
» 猜你喜欢
售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急
已经有6人回复
售一区SCI文章T0P,我:8O.551.O54,科目全,可十急
已经有7人回复
售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急
已经有10人回复
售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急
已经有8人回复
售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急
已经有9人回复
售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急
已经有8人回复
售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急
已经有13人回复
|