DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: 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 » 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.« less

Authors:
 [1];  [2];  [3]
  1. Nagoya University (Japan)
  2. Katholieke University Leuven, Heverlee (Belgium); Stanford University, CA (United States)
  3. 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}
}

Journal Article:

Citation Metrics:
Cited by: 9 works
Citation information provided by
Web of Science

Save / Share:

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
  • DOI: 10.1137/S1064827596313322

An application of fast multipole method to isogeometric boundary element method for Laplace equation in two dimensions
journal, December 2012


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
  • DOI: 10.1121/1.4932166

The Fast Multipole Method: Numerical Implementation
journal, May 2000


Rapid solution of integral equations of classical potential theory
journal, September 1985


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
  • DOI: 10.1109/TAP.1982.1142818

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
  • DOI: 10.2528/PIER09042104

Diagonal Forms of Translation Operators for the Helmholtz Equation in Three Dimensions
journal, December 1993


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
  • DOI: 10.1115/1.4005491

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
  • DOI: 10.1109/74.250128

Adaptive Low-Rank Approximation of Collocation Matrices
journal, February 2003


A fast multipole boundary element method based on the improved Burton–Miller formulation for three-dimensional acoustic problems
journal, May 2011


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
  • DOI: 10.1137/0907058

The Fast Multipole Method I: Error Analysis and Asymptotic Complexity
journal, January 2000


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
  • DOI: 10.1093/imanum/drm001

Advances in Iterative Methods and Preconditioners for the Helmholtz Equation
journal, December 2007


A Preconditioned 3-D Multi-Region Fast Multipole Solver for Seismic Wave Propagation in Complex Geometries
journal, February 2012


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
  • DOI: 10.1093/qjmam/38.2.323

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
  • DOI: 10.1137/050639028

Controlling sound with acoustic metamaterials
journal, February 2016

  • Cummer, Steven A.; Christensen, Johan; Alù, Andrea
  • Nature Reviews Materials, Vol. 1, Issue 3
  • DOI: 10.1038/natrevmats.2016.1

The black-box fast multipole method
journal, December 2009


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
  • DOI: 10.1137/090771806

A low frequency elastodynamic fast multipole boundary element method in three dimensions
journal, September 2015


Incomplete LU preconditioning for large scale dense complex linear systems from electromagnetic wave scattering problems
journal, February 2003


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
  • DOI: 10.1137/120866683

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
  • DOI: 10.1002/cpe.3723

Preconditioning based on Calderon’s formulae for periodic fast multipole methods for Helmholtz’ equation
journal, January 2012


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
  • DOI: 10.1121/1.3021297

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
  • DOI: 10.1007/s10543-014-0508-y

The fast multipole method on parallel clusters, multicore processors, and graphics processing units
journal, February 2011


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
  • DOI: 10.1137/130915662

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
  • DOI: 10.1098/rsta.2003.1337

An O ( N ) 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
  • DOI: 10.1016/j.acha.2014.04.002

Fast multipole accelerated boundary integral equation methods
journal, July 2002


A fast direct solver for boundary integral equations in two dimensions
journal, May 2005


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
  • DOI: 10.1002/cpa.21582

Efficient fast multipole method for low-frequency scattering
journal, June 2004


A Flexible Inner-Outer Preconditioned GMRES Algorithm
journal, March 1993

  • Saad, Youcef
  • SIAM Journal on Scientific Computing, Vol. 14, Issue 2
  • DOI: 10.1137/0914028

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
  • DOI: 10.1016/j.jcp.2011.09.027

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
  • DOI: 10.1137/11085774X

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
  • DOI: 10.1002/nme.4302

On efficient preconditioners for iterative solution of a Galerkin boundary element equation for the three-dimensional exterior Helmholtz problem
journal, July 2003


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
  • DOI: 10.1016/j.cma.2016.05.029

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
  • DOI: 10.1137/040603917

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
  • DOI: 10.1137/0916051

An efficient blocking M2L translation for low-frequency fast multipole method in three dimensions
journal, May 2016


A fast algorithm for particle simulations
journal, December 1987


Combining analytic preconditioner and Fast Multipole Method for the 3-D Helmholtz equation
journal, March 2013