The Fast Multipole Method: Numerical Implementation
β Scribed by Eric Darve
- Publisher
- Elsevier Science
- Year
- 2000
- Tongue
- English
- Weight
- 615 KB
- Volume
- 160
- Category
- Article
- ISSN
- 0021-9991
No coin nor oath required. For personal study only.
β¦ Synopsis
We study integral methods applied to the resolution of the Maxwell equations where the linear system is solved using an iterative method which requires only matrix-vector products. The fast multipole method (FMM) is one of the most efficient methods used to perform matrix-vector products and accelerate the resolution of the linear system. A problem involving N degrees of freedom may be solved in CN iter N log N floating operations, where C is a constant depending on the implementation of the method. In this article several techniques allowing one to reduce the constant C are analyzed. This reduction implies a lower total CPU time and a larger range of application of the FMM. In particular, new interpolation and anterpolation schemes are proposed which greatly improve on previous algorithms. Several numerical tests are also described. These confirm the efficiency and the theoretical complexity of the FMM.
π SIMILAR VOLUMES
We apply a new method of reducing the computational effort required for solution of the Electric Field Integral Equation used for modelling microstrip structures. The Fast Multipole Method is used to compute the radiation pattern and input impedance of single layer microstrip antennas.
The numerical solution to N-body problems in gravitation or electrostatics has traditionally been obtained via particle-in-cell methods (PLC) since direct evaluation of all pairwise interparticle forces, requiring t!~( N 2) operations, is too expensive. Recently, hierarchical solvers, which use tree
One def'ming characteristic of a fast multipole calculation is the number of tiers (depth of tree) used to group the particles. For three dimensions, the standard boxing scheme restricts the number of lowest level boxes to be a power of eight. We present a method which through a simple scaling of th