Gerard L.G. Sleijpen, Peter Sonneveld and Martin B. van Gijzen, Bi ...
Gerard L.G. Sleijpen, Peter Sonneveld and Martin B. van Gijzen, Bi ...
Gerard L.G. Sleijpen, Peter Sonneveld and Martin B. van Gijzen, Bi ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Author's personal copy<br />
1112 G.L.G. <strong>Sleijpen</strong> et al. / Applied Numerical Mathematics 60 (2010) 1100–1114<br />
Fig. 1. Convergence of <strong>Bi</strong>-CGSTAB(4) <strong>and</strong> IDR(4) for SHERMAN4.<br />
Table 1<br />
Number of matrix–vector multiplications to solve the SHERMAN4 system<br />
such that the scaled norm of the updated residual is less than 10 −8 .<br />
Method Number of MVs ‖b − Ax‖/‖b‖<br />
GMRES 120 9.8×10 −9<br />
IDR(1) 179 6.6×10 −9<br />
IDR(2) 161 8.3×10 −9<br />
IDR(3) 153 8.7×10 −9<br />
IDR(4) 146 3.0×10 −9<br />
IDR(5) 150 1.2×10 −9<br />
<strong>Bi</strong>-CGSTAB(1) 180 9.0×10 −9<br />
<strong>Bi</strong>-CGSTAB(2) 162 8.6×10 −9<br />
<strong>Bi</strong>-CGSTAB(3) 156 4.6×10 −9<br />
<strong>Bi</strong>-CGSTAB(4) 150 1.8×10 −7<br />
<strong>Bi</strong>-CGSTAB(5) 144 7.1×10 −9<br />
Table 2<br />
Number of matrix–vector multiplications to solve the ADD20 system such<br />
that the scaled norm of the updated residual is less than 10 −8 .<br />
Method Number of MVs ‖b − Ax‖/‖b‖<br />
GMRES 295 1.0×10 −8<br />
IDR(1) 672 9.0×10 −9<br />
IDR(2) 581 9.9×10 −9<br />
IDR(3) 588 4.8×10 −9<br />
IDR(4) 480 5.3×10 −9<br />
IDR(5) 444 9.4×10 −9<br />
<strong>Bi</strong>-CGSTAB(1) 728 9.0×10 −9<br />
<strong>Bi</strong>-CGSTAB(2) 648 9.8×10 −9<br />
<strong>Bi</strong>-CGSTAB(3) 528 2.0×10 −7<br />
<strong>Bi</strong>-CGSTAB(4) 545 7.1×10 −6<br />
<strong>Bi</strong>-CGSTAB(5) 702 0.026