Application of the inverse fast multipole method as a preconditioner in a 3D Helmholtz boundary element method
Abstract
Here we investigate an efficient preconditioning of iterative methods (such as GMRES) for solving dense linear systems Ax = b that follow from a boundary element method (BEM) for the 3D Helmholtz equation, focusing on the low-frequency regime. While matrix–vector products in GMRES can be accelerated through the low-frequency fast multipole method (LFFMM), the BEM often remains computationally expensive due to the large number of GMRES iterations. We propose the application of the inverse fast multipole method (IFMM) as a preconditioner to accelerate the convergence of GMRES. The IFMM is in essence an approximate direct solver that uses a multilevel hierarchical decomposition and low-rank approximations. The proposed IFMM-based preconditioning has a tunable parameter ε that balances the cost to construct a preconditioner M, which is an approximation of A-1, and the cost to perform the iterative process by means of M. Namely, using a small (respectively, large) value of ε takes a long (respectively, short) time to construct M, while the number of iterations can be small (respectively, large). A comprehensive set of numerical examples involving various boundary value problems with complicated geometries and mixed boundary conditions is presented to validate the efficiency of the proposed method. We show thatmore »
- Authors:
-
- Nagoya University (Japan)
- Katholieke University Leuven, Heverlee (Belgium); Stanford University, CA (United States)
- Stanford University, CA (United States)
- Publication Date:
- Research Org.:
- Stanford Univ., CA (United States)
- Sponsoring Org.:
- USDOE National Nuclear Security Administration (NNSA); US Army Research Laboratory (USARL); Japan Society for the Promotion of Science (JSPS)
- OSTI Identifier:
- 1533958
- Alternate Identifier(s):
- OSTI ID: 1397802
- Grant/Contract Number:
- NA0002373; 15K06683; W911NF-07-0027; NA0002373-1
- Resource Type:
- Accepted Manuscript
- Journal Name:
- Journal of Computational Physics
- Additional Journal Information:
- Journal Volume: 341; Journal Issue: C; Journal ID: ISSN 0021-9991
- Publisher:
- Elsevier
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 97 MATHEMATICS AND COMPUTING; 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; fast multipole method; boundary element method; inverse fast multipole method; iterative solver; preconditioning; low rank compression
Citation Formats
Takahashi, Toru, Coulier, Pieter, and Darve, Eric. Application of the inverse fast multipole method as a preconditioner in a 3D Helmholtz boundary element method. United States: N. p., 2017.
Web. doi:10.1016/j.jcp.2017.04.016.
Takahashi, Toru, Coulier, Pieter, & Darve, Eric. Application of the inverse fast multipole method as a preconditioner in a 3D Helmholtz boundary element method. United States. https://doi.org/10.1016/j.jcp.2017.04.016
Takahashi, Toru, Coulier, Pieter, and Darve, Eric. Fri .
"Application of the inverse fast multipole method as a preconditioner in a 3D Helmholtz boundary element method". United States. https://doi.org/10.1016/j.jcp.2017.04.016. https://www.osti.gov/servlets/purl/1533958.
@article{osti_1533958,
title = {Application of the inverse fast multipole method as a preconditioner in a 3D Helmholtz boundary element method},
author = {Takahashi, Toru and Coulier, Pieter and Darve, Eric},
abstractNote = {Here we investigate an efficient preconditioning of iterative methods (such as GMRES) for solving dense linear systems Ax = b that follow from a boundary element method (BEM) for the 3D Helmholtz equation, focusing on the low-frequency regime. While matrix–vector products in GMRES can be accelerated through the low-frequency fast multipole method (LFFMM), the BEM often remains computationally expensive due to the large number of GMRES iterations. We propose the application of the inverse fast multipole method (IFMM) as a preconditioner to accelerate the convergence of GMRES. The IFMM is in essence an approximate direct solver that uses a multilevel hierarchical decomposition and low-rank approximations. The proposed IFMM-based preconditioning has a tunable parameter ε that balances the cost to construct a preconditioner M, which is an approximation of A-1, and the cost to perform the iterative process by means of M. Namely, using a small (respectively, large) value of ε takes a long (respectively, short) time to construct M, while the number of iterations can be small (respectively, large). A comprehensive set of numerical examples involving various boundary value problems with complicated geometries and mixed boundary conditions is presented to validate the efficiency of the proposed method. We show that the IFMM preconditioner (with a nearly optimal ε of 10-2) clearly outperforms some common preconditioners for the BEM, achieving 1.2–10.8 times speed-up of the computations, in particular when the scale of the underlying scatterer is about five wavelengths or more. In addition, the IFMM preconditioner is capable of solving complicated problems (in a reasonable amount of time) that BD preconditioner can not.},
doi = {10.1016/j.jcp.2017.04.016},
journal = {Journal of Computational Physics},
number = C,
volume = 341,
place = {United States},
year = {Fri Apr 07 00:00:00 EDT 2017},
month = {Fri Apr 07 00:00:00 EDT 2017}
}
Web of Science
Works referenced in this record:
Parallel Hierarchical Solvers and Preconditioners for Boundary Element Methods
journal, January 1998
- Grama, Ananth; Kumar, Vipin; Sameh, Ahmed
- SIAM Journal on Scientific Computing, Vol. 20, Issue 1
An application of fast multipole method to isogeometric boundary element method for Laplace equation in two dimensions
journal, December 2012
- Takahashi, Toru; Matsumoto, Toshiro
- Engineering Analysis with Boundary Elements, Vol. 36, Issue 12
A fast boundary element method for the scattering analysis of high-intensity focused ultrasound
journal, November 2015
- van 't Wout, Elwin; Gélat, Pierre; Betcke, Timo
- The Journal of the Acoustical Society of America, Vol. 138, Issue 5
The Fast Multipole Method: Numerical Implementation
journal, May 2000
- Darve, Eric
- Journal of Computational Physics, Vol. 160, Issue 1
Rapid solution of integral equations of classical potential theory
journal, September 1985
- Rokhlin, V.
- Journal of Computational Physics, Vol. 60, Issue 2
Electromagnetic scattering by surfaces of arbitrary shape
journal, May 1982
- Rao, S.; Wilton, D.; Glisson, A.
- IEEE Transactions on Antennas and Propagation, Vol. 30, Issue 3
Geometry Based Preconditioner for Radiation Problems Involving wire and Surface Basis Functions
journal, January 2009
- Araujo, Marta Gomez; Bertolo, Jose Manuel; Obelleiro, Fernando
- Progress In Electromagnetics Research, Vol. 93
Diagonal Forms of Translation Operators for the Helmholtz Equation in Three Dimensions
journal, December 1993
- Rokhlin, V.
- Applied and Computational Harmonic Analysis, Vol. 1, Issue 1
Recent Advances and Emerging Applications of the Boundary Element Method
journal, May 2011
- Liu, Y. J.; Mukherjee, S.; Nishimura, N.
- Applied Mechanics Reviews, Vol. 64, Issue 3
The fast multipole method for the wave equation: a pedestrian prescription
journal, June 1993
- Coifman, R.; Rokhlin, V.; Wandzura, S.
- IEEE Antennas and Propagation Magazine, Vol. 35, Issue 3
Adaptive Low-Rank Approximation of Collocation Matrices
journal, February 2003
- Bebendorf, M.; Rjasanow, S.
- Computing, Vol. 70, Issue 1
A fast multipole boundary element method based on the improved Burton–Miller formulation for three-dimensional acoustic problems
journal, May 2011
- Li, Shande; Huang, Qibai
- Engineering Analysis with Boundary Elements, Vol. 35, Issue 5
GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems
journal, July 1986
- Saad, Youcef; Schultz, Martin H.
- SIAM Journal on Scientific and Statistical Computing, Vol. 7, Issue 3
The Fast Multipole Method I: Error Analysis and Asymptotic Complexity
journal, January 2000
- Darve, Eric
- SIAM Journal on Numerical Analysis, Vol. 38, Issue 1
Hierarchical matrix techniques for low- and high-frequency Helmholtz problems
journal, March 2007
- Banjai, L.; Hackbusch, W.
- IMA Journal of Numerical Analysis, Vol. 28, Issue 1
Advances in Iterative Methods and Preconditioners for the Helmholtz Equation
journal, December 2007
- Erlangga, Yogi A.
- Archives of Computational Methods in Engineering, Vol. 15, Issue 1
A wideband fast multipole accelerated boundary integral equation method for time-harmonic elastodynamics in two dimensions: A WIDEBAND FAST MULTIPOLE ACCELERATED BIEM FOR ELASTODYNAMICS
journal, March 2012
- Takahashi, Toru
- International Journal for Numerical Methods in Engineering, Vol. 91, Issue 5
A Preconditioned 3-D Multi-Region Fast Multipole Solver for Seismic Wave Propagation in Complex Geometries
journal, February 2012
- Chaillat, S.; Semblat, J. F.; Bonnet, M.
- Communications in Computational Physics, Vol. 11, Issue 2
Minimizing the Condition Number of Boundary Integral Operators in Acoustic and Electromagnetic Scattering
journal, January 1985
- Kress, R.
- The Quarterly Journal of Mechanics and Applied Mathematics, Vol. 38, Issue 2
A Fast Solver for HSS Representations via Sparse Matrices
journal, January 2007
- Chandrasekaran, S.; Dewilde, P.; Gu, M.
- SIAM Journal on Matrix Analysis and Applications, Vol. 29, Issue 1
Controlling sound with acoustic metamaterials
journal, February 2016
- Cummer, Steven A.; Christensen, Johan; Alù, Andrea
- Nature Reviews Materials, Vol. 1, Issue 3
The black-box fast multipole method
journal, December 2009
- Fong, William; Darve, Eric
- Journal of Computational Physics, Vol. 228, Issue 23
Finding Structure with Randomness: Probabilistic Algorithms for Constructing Approximate Matrix Decompositions
journal, January 2011
- Halko, N.; Martinsson, P. G.; Tropp, J. A.
- SIAM Review, Vol. 53, Issue 2
A low frequency elastodynamic fast multipole boundary element method in three dimensions
journal, September 2015
- Wilkes, D. R.; Duncan, A. J.
- Computational Mechanics, Vol. 56, Issue 5
Incomplete LU preconditioning for large scale dense complex linear systems from electromagnetic wave scattering problems
journal, February 2003
- Lee, Jeonghwa; Zhang, Jun; Lu, Cai-Cheng
- Journal of Computational Physics, Vol. 185, Issue 1
A Fast Direct Solver for Structured Linear Systems by Recursive Skeletonization
journal, January 2012
- Ho, Kenneth L.; Greengard, Leslie
- SIAM Journal on Scientific Computing, Vol. 34, Issue 5
Task‐based FMM for heterogeneous architectures
journal, December 2015
- Agullo, Emmanuel; Bramas, Berenger; Coulaud, Olivier
- Concurrency and Computation: Practice and Experience, Vol. 28, Issue 9
Preconditioning based on Calderon’s formulae for periodic fast multipole methods for Helmholtz’ equation
journal, January 2012
- Niino, Kazuki; Nishimura, Naoshi
- Journal of Computational Physics, Vol. 231, Issue 1
A broadband fast multipole accelerated boundary element method for the three dimensional Helmholtz equation
journal, January 2009
- Gumerov, Nail A.; Duraiswami, Ramani
- The Journal of the Acoustical Society of America, Vol. 125, Issue 1
A high-order accurate accelerated direct solver for acoustic scattering from surfaces
journal, July 2014
- Bremer, James; Gillman, Adrianna; Martinsson, Per-Gunnar
- BIT Numerical Mathematics, Vol. 55, Issue 2
The fast multipole method on parallel clusters, multicore processors, and graphics processing units
journal, February 2011
- Darve, Eric; Cecka, Cris; Takahashi, Toru
- Comptes Rendus Mécanique, Vol. 339, Issue 2-3
Task-Based FMM for Multicore Architectures
journal, January 2014
- Agullo, Emmanuel; Bramas, Bérenger; Coulaud, Olivier
- SIAM Journal on Scientific Computing, Vol. 36, Issue 1
A fast multipole method for Maxwell equations stable at all frequencies
journal, January 2004
- Darve, Eric; Havé, Pascal
- Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences, Vol. 362, Issue 1816
An direct solver for integral equations on the plane
journal, March 2015
- Corona, Eduardo; Martinsson, Per-Gunnar; Zorin, Denis
- Applied and Computational Harmonic Analysis, Vol. 38, Issue 2
Fast multipole accelerated boundary integral equation methods
journal, July 2002
- Nishimura, N.
- Applied Mechanics Reviews, Vol. 55, Issue 4
A fast direct solver for boundary integral equations in two dimensions
journal, May 2005
- Martinsson, P. G.; Rokhlin, V.
- Journal of Computational Physics, Vol. 205, Issue 1
Hierarchical Interpolative Factorization for Elliptic Operators: Differential Equations
journal, May 2015
- Ho, Kenneth L.; Ying, Lexing
- Communications on Pure and Applied Mathematics, Vol. 69, Issue 8
Efficient fast multipole method for low-frequency scattering
journal, June 2004
- Darve, Eric; Havé, Pascal
- Journal of Computational Physics, Vol. 197, Issue 1
A Flexible Inner-Outer Preconditioned GMRES Algorithm
journal, March 1993
- Saad, Youcef
- SIAM Journal on Scientific Computing, Vol. 14, Issue 2
Fast directional multilevel summation for oscillatory kernels based on Chebyshev interpolation
journal, February 2012
- Messner, Matthias; Schanz, Martin; Darve, Eric
- Journal of Computational Physics, Vol. 231, Issue 4
Fourier-Based Fast Multipole Method for the Helmholtz Equation
journal, January 2013
- Cecka, Cris; Darve, Eric
- SIAM Journal on Scientific Computing, Vol. 35, Issue 1
Regularized integral equations and fast high-order solvers for sound-hard acoustic scattering problems: REGULARIZED EQUATIONS AND EFFICIENT SOLVERS FOR SOUND-HARD SCATTERING
journal, June 2012
- Bruno, Oscar; Elling, Tim; Turc, Catalin
- International Journal for Numerical Methods in Engineering, Vol. 91, Issue 10
On efficient preconditioners for iterative solution of a Galerkin boundary element equation for the three-dimensional exterior Helmholtz problem
journal, July 2003
- Harris, Paul J.; Chen, Ke
- Journal of Computational and Applied Mathematics, Vol. 156, Issue 2
Efficient mesh deformation based on radial basis function interpolation by means of the inverse fast multipole method
journal, August 2016
- Coulier, Pieter; Darve, Eric
- Computer Methods in Applied Mechanics and Engineering, Vol. 308
Combining Fast Multipole Techniques and an Approximate Inverse Preconditioner for Large Electromagnetism Calculations
journal, January 2005
- Carpentieri, B.; Duff, I. S.; Giraud, L.
- SIAM Journal on Scientific Computing, Vol. 27, Issue 3
Multipole Translation Theory for the Three-Dimensional Laplace and Helmholtz Equations
journal, July 1995
- Epton, Michael A.; Dembart, Benjamin
- SIAM Journal on Scientific Computing, Vol. 16, Issue 4
An efficient blocking M2L translation for low-frequency fast multipole method in three dimensions
journal, May 2016
- Takahashi, Toru; Shimba, Yuta; Isakari, Hiroshi
- Computer Physics Communications, Vol. 202
A fast algorithm for particle simulations
journal, December 1987
- Greengard, L.; Rokhlin, V.
- Journal of Computational Physics, Vol. 73, Issue 2
Combining analytic preconditioner and Fast Multipole Method for the 3-D Helmholtz equation
journal, March 2013
- Darbas, M.; Darrigrand, E.; Lafranche, Y.
- Journal of Computational Physics, Vol. 236