Efficient parallelization of the multilevel fast multipole algorithm for the solution of large-scale scattering problems

Download
2008-08-01
We present fast and accurate solutions of large-scale scattering problems involving three-dimensional closed conductors with arbitrary shapes using the multilevel fast multipole algorithm (MLFMA). With an efficient parallelization of MLFMA, scattering problems that are discretized with tens of millions of unknowns are easily solved on a cluster of computers. We extensively investigate the parallelization of MLFMA, identify the bottlenecks, and provide remedial procedures to improve the efficiency of the implementations. The accuracy of the solutions is demonstrated on a scattering problem involving a sphere of radius 110 lambda discretized with 41 883 638 unknowns, the largest integral-equation problem solved to date. In addition to canonical problems, we also present the solution of real-life problems involving complicated targets with large dimensions.
IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION

Suggestions

Efficient solution of the electric-field integral equation using the iterative LSQR algorithm
Ergül, Özgür Salih (Institute of Electrical and Electronics Engineers (IEEE), 2008-01-01)
In this letter, we consider iterative solutions of the three-dimensional electromagnetic scattering problems formulated by surface integral equations. We show that solutions of the electric-field integral equation (EFIE) can be improved by employing an iterative least-squares QR (LSQR) algorithm. Compared to many other Krylov subspace methods, LSQR provides faster convergence and it becomes an alternative choice to the time-efficient no-restart generalized minimal residual (GMRES) algorithm that requires la...
Enhancing the accuracy of the interpolations and anterpolations in MLFMA
Ergül, Özgür Salih (Institute of Electrical and Electronics Engineers (IEEE), 2006-01-01)
We present an efficient technique to reduce the interpolation and anterpolation (transpose interpolation) errors in the aggregation and disaggregation processes of the multilevel fast multipole algorithm (MLFMA), which is based on the sampling of the radiated and incoming fields over all possible solid angles, i.e., all directions on the sphere. The fields sampled on the sphere are subject to various operations, such as interpolation, aggregation, translation, disaggregation, anterpolation, and integration....
Hierarchical parallelisation strategy for multilevel fast multipole algorithm in computational electromagnetics
Ergül, Özgür Salih (Institution of Engineering and Technology (IET), 2008-01-03)
A hierarchical parallelisation of the multilevel fast multipole algorithm (MLFMA) for the efficient solution of large-scale problems in computational electromagnetics is presented. The tree structure of MLFMA is distributed among the processors by partitioning both the clusters and the samples of the fields appropriately for each level. The parallelisation efficiency is significantly improved compared to previous approaches, where only the clusters or only the fields are partitioned in a level.
Fast and accurate solutions of extremely large integral-equation problems discretised with tens of millions of unknowns
Gurel, L.; Ergül, Özgür Salih (Institution of Engineering and Technology (IET), 2007-04-26)
The solution of extremely large scattering problems that are formulated by integral equations and discretised with tens of millions of unknowns is reported. Accurate and efficient solutions are performed by employing a parallel implementation of the multilevel fast multipole algorithm. The effectiveness of the implementation is demonstrated on a sphere problem containing more than 33 million unknowns, which is the largest integral-equation problem ever solved to our knowledge.
Stabilization of integral-equation formulations for the accurate solution of scattering problems involving low-contrast dielectric objects
Ergül, Özgür Salih (Institute of Electrical and Electronics Engineers (IEEE), 2008-03-01)
The solution of scattering problems involving low-contrast dielectric objects with three-dimensional arbitrary shapes is considered. Using the traditional forms of the surface integral equations, scattered fields cannot be calculated accurately if the contrast of the object is low. Therefore, we consider the stabilization of the formulations by extracting the nonradiating parts of the equivalent currents. We also investigate various types of stable formulations and show that accuracy can be improved systema...
Citation Formats
Ö. S. Ergül, “Efficient parallelization of the multilevel fast multipole algorithm for the solution of large-scale scattering problems,” IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, pp. 2335–2345, 2008, Accessed: 00, 2020. [Online]. Available: https://hdl.handle.net/11511/36798.