Next: Algorithm for the new Up: Crack problems for three-dimensional Previous: New FM-BIEM

## Rotation of coefficients

The discussion in the previous section has been restricted to the case where Ct is in +x3 direction of Cs. In this section we shall remove this assumption to generalise the discussion in the same way as in 4.1. If the target cell is included in lists except the uplist of Cs (See Fig.4.2) we rotate the coordinate system so that the target cell is in the positive direction viewed from the source cell, where denotes the new axis. In general the multipole moments in the new coordinate system are obtained as follows:

 = (4.34) = (4.35)

where is the coefficient of rotation (See (4.12)), is a unit vector parallel to the rotation axis, is a rotation angle and is a rotation matrix.

We next describe the generalised M2L translation process in the new FMM.

1.
Rotation:

First we rotate the multipole moments via (4.34) and (4.35) so as to make the procedure presented in 4.2.1 applicable. The specific forms of (4.34) and (4.35) depend on the location of Ct and are described as follows:

• uplist
 = (4.36) = M2n,m(O) (4.37)

• downlist
 = (4.38) = (4.39)

• northlist
 = (4.40) = (4.41)

• southlist
 = (4.42) = (4.43)

• eastlist
 = (4.44) = (4.45)

• westlist
 = (4.46) = (4.47)

where is the base vector for the Cartesian coordinates and superposed indices {U, D, N, S, E, W} correspond to the initial letters of {uplist, downlist, northlist, southlist, eastlist, westlist}, respectively.
2.
Compute the coefficients of the exponential expansion:

Compute the coefficients of the exponential expansion via (4.28) and (4.29) as follows:

 = (4.48) = (4.49)

where is an element of .
3.
Translation of the coefficients of the exponential expansion:

As the centre of the exponential expansion is shifted from the centroid of Cs (O) to the centroid of Ct ( ), the coefficients of the exponential expansion are translated according to (4.30) and (4.31) as follows:

 = (4.50) = (4.51) =

where is an element of .
4.
Compute the coefficients of the local expansion:

Compute the coefficients of the local expansion from the exponential expansion according to (4.32) and (4.33) as follows:

 = (4.52) = (4.53)

where is an element of . Then rotate as follows:
• uplist
 = (4.54) = (4.55)

• downlist
 = (4.56) = (4.57)

• northlist
 = (4.58) = (4.59)

• southlist
 = (4.60) = (4.61)

• eastlist
 = (4.62) = (4.63)

• westlist
 = (4.64) = (4.65)

Finally add LU,LD,LN,LS,LE and LW together
 = (4.66) = (4.67)

to obtain the coefficients of the local expansion.

Next: Algorithm for the new Up: Crack problems for three-dimensional Previous: New FM-BIEM
Ken-ichi Yoshida
2001-07-28