02.11.2014 Views

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 ...

SHOW MORE
SHOW LESS

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

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

Saved successfully!

Ooh no, something went wrong!