23.12.2012 Views

MATLAB Function Reference Volume 1: A - E - Bad Request

MATLAB Function Reference Volume 1: A - E - Bad Request

MATLAB Function Reference Volume 1: A - E - Bad Request

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

icgstab<br />

flag2 is 0 because bicgstab converges to the tolerance of 3.1757e-016 (the<br />

value of relres2) at the sixth iteration (the value of iter2) when<br />

preconditioned by the incomplete LU factorization with a drop tolerance of<br />

1e-6. resvec2(1) = norm(b) and resvec2(13) = norm(b-A*x2). You can<br />

follow the progress of bicgstab by plotting the relative residuals at the halfway<br />

point and end of each iteration starting from the initial estimate (iterate<br />

number 0).<br />

semilogy(0:0.5:iter2,resvec2/norm(b),'-o')<br />

xlabel('iteration number')<br />

ylabel('relative residual')<br />

relative residual<br />

10 0<br />

10 −2<br />

10 −4<br />

10 −6<br />

10 −8<br />

10 −10<br />

10 −12<br />

10 −14<br />

10<br />

0 1 2 3 4 5 6<br />

−16<br />

iteration number<br />

See Also bicg, cgs, gmres, lsqr, luinc, minres, pcg, qmr, symmlq<br />

@ (function handle), \ (backslash)<br />

2-157

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!