skip to main content
OSTI.GOV 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

Authors:
; ;
Publication Date:
Sponsoring Org.:
USDOE
OSTI Identifier:
1397802
Grant/Contract Number:
NA0002373-1
Resource Type:
Journal Article: Publisher's Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 341; Journal Issue: C; Related Information: CHORUS Timestamp: 2017-10-04 22:14:01; Journal ID: ISSN 0021-9991
Publisher:
Elsevier
Country of Publication:
United States
Language:
English

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. doi:10.1016/j.jcp.2017.04.016.
Takahashi, Toru, Coulier, Pieter, and Darve, Eric. Sat . "Application of the inverse fast multipole method as a preconditioner in a 3D Helmholtz boundary element method". United States. doi:10.1016/j.jcp.2017.04.016.
@article{osti_1397802,
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 = {},
doi = {10.1016/j.jcp.2017.04.016},
journal = {Journal of Computational Physics},
number = C,
volume = 341,
place = {United States},
year = {Sat Jul 01 00:00:00 EDT 2017},
month = {Sat Jul 01 00:00:00 EDT 2017}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record at 10.1016/j.jcp.2017.04.016

Save / Share:
  • The numerical solution of the Poisson Boltzmann (PB) equation is a useful but a computationally demanding tool for studying electrostatic solvation effects in chemical and biomolecular systems. Recently, we have described a boundary integral equation-based PB solver accelerated by a new version of the fast multipole method (FMM). The overall algorithm shows an order N complexity in both the computational cost and memory usage. Here, we present an updated version of the solver by using an adaptive FMM for accelerating the convolution type matrix-vector multiplications. The adaptive algorithm, when compared to our previous nonadaptive one, not only significantly improves themore » performance of the overall memory usage but also remarkably speeds the calculation because of an improved load balancing between the local- and far-field calculations. We have also implemented a node-patch discretization scheme that leads to a reduction of unknowns by a factor of 2 relative to the constant element method without sacrificing accuracy. As a result of these improvements, the new solver makes the PB calculation truly feasible for large-scale biomolecular systems such as a 30S ribosome molecule even on a typical 2008 desktop computer.« less
  • We investigate the parallel performance of an iterative solver for 3D heterogeneous Helmholtz problems related to applications in seismic wave propagation. For large 3D problems, the computation is no longer feasible on a single processor, and the memory requirements increase rapidly. Therefore, parallelization of the solver is needed. We employ a complex shifted-Laplace preconditioner combined with the Bi-CGSTAB iterative method and use a multigrid method to approximate the inverse of the resulting preconditioning operator. A 3D multigrid method with 2D semi-coarsening is employed. We show numerical results for large problems arising in geophysical applications.
  • Cited by 5
  • This work introduces new functions based on the spherical harmonic and the solid harmonic which be used to construct a multipole development for the 3D Stokes problem, in order to reduce the operations costs in the BEM method. We show that the major properties of those functions are inherited from the solid harmonics. The contribution of this paper is the introduction of new formulas that serves to calculate the multipole moments and the transfer functions that are necessar for the schemes of order O(NlogN). Moreover, new translations formulas were introduced to obtain O(N) scheme. The error truncation of the resultingmore » scheme is discussed.« less