** Next:** One penny-shaped crack
** Up:** Crack problems for three-dimensional
** Previous:** Algorithm for the new

The proposed techniques have been implemented in Fortran 77 and
have been tested on a computer having a DEC Alpha 21264(500 MHz) as the
CPU. The integrals in the multipole moments in (3.53) and
(3.54) are computed numerically with Gaussian quadrature. The
sums in the finite series (3.52) and
(3.59) are truncated at 10 terms and the series in
(4.48), (4.49), (4.52) and (4.53) are
computed with the 109 point generalised Gaussian quadrature formula
given in Yarvin and Rokhlin[83]. The maximum number of boundary
elements in a leaf is set to be 100. To solve the resulting matrix
equation we use the preconditioned GMRES and adopt the block diagonal
matrix corresponding to the leaves as the preconditioner according to
Nishida and Hayami[62]. In GMRES the iteration is stopped when
the relative residual norm is below 10^{-5}.

*Ken-ichi Yoshida*

*2001-07-28*