?? gmres.m
字號:
function [x,k]=gmres(A,b,m,toll,x0)r0=b-A*x0; nr0=norm(r0,2);if nr0 ~= 0 v1=r0/nr0; V=[v1]; H=[]; k=0; residual=1; while k <= m-1 & residual > toll, [k,V,H] = GSarnoldi(A,m,k,V,H); [nr,nc]=size(H); y=(H'*H) \ (H'*nr0*[1;zeros(nr-1,1)]); x=x0+V(:,1:nc)*y; residual = norm(b-A*x,2)/nr0; endelse x=x0;endreturn
?? 快捷鍵說明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -