next up previous contents
Next: One penny-shaped crack Up: Three-dimensional scattering of scalar Previous: Algorithm

Numerical examples

The proposed techniques have been implemented in Fortran 77, and have been tested on a computer having a DEC Alpha 21164(533 MHz) chip as the CPU (except for the results in Fig.3.49). In FMM, we truncate the infinite series via (3.133). Also, we set the maximum number of boundary elements in a leaf to be 100. The integrals in (3.128) are computed numerically with Gaussian quadrature. In GMRES we terminate the iteration when the relative error is less than 10-5. To solve the resulting matrix equation we use the preconditioned GMRES (The issue of the choice of the preconditioners will be touched upon later).


Ken-ichi Yoshida