Inexact GMRES iterations and relaxation strategies with fast-multipole boundary element method

DOI10.1007/S10444-022-09932-8MaRDI QIDQ2673503zbMATH OpenOpenAlexFDO

Authors Tingyu Wang, Simon K. Layton, Lorena A. Barba

Publication date 10 June 2022

Published in Advances in Computational Mathematics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1506.05957

zbMATH Keywords

Stokes flowboundary integral equationcollocation methodboundary element methoditerative solversfast multipole methodKrylov methods

Mathematics Subject Classification ID

Numerical quadrature and cubature formulas (65D32) Iterative numerical methods for linear systems (65F10) Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation (35J05) PDEs in connection with fluid mechanics (35Q35) Electro- and magnetostatics (78A30) PDEs in connection with optics and electromagnetic theory (35Q60) Spectral, collocation and related methods for boundary value problems involving PDEs (65N35) Boundary element methods for boundary value problems involving PDEs (65N38) Stokes and related (Oseen, etc.) flows (76D07) Physiological flows (76Z05) Physiological flow (92C35) Fredholm integral equations (45B05) Multipole methods applied to problems in optics and electromagnetic theory (78M16)

Abstract: Boundary element methods produce dense linear systems that can be accelerated via multipole expansions. Solved with Krylov methods, this implies computing the matrix-vector products within each iteration with some error, at an accuracy controlled by the order of the expansion, $p$. We take advantage of a unique property of Krylov iterations that allow lower accuracy of the matrix-vector products as convergence proceeds, and propose a relaxation strategy based on progressively decreasing $p$. In extensive numerical tests of the relaxed Krylov iterations, we obtained speed-ups of between $2.1imes$ and $3.3imes$ for Laplace problems and between $1.7imes$ and $4.0imes$ for Stokes problems. We include an application to Stokes flow around red blood cells, computing with up to 64 cells and problem size up to 131k boundary elements and nearly 400k unknowns. The study was done with an in-house multi-threaded C++ code, on a hexa-core CPU. The code is available on its version-control repository, href{https://github.com/barbagroup/fmm-bem-relaxed}{https://github.com/barbagroup/fmm-bem-relaxed}.

Recommendations

• Inexact Matrix-Vector Products in Krylov Methods for Solving Linear Systems: A Relaxation Strategy
• Inexact Krylov Subspace Methods for Linear Systems
• Convergence in Backward Error of Relaxed GMRES
• Relaxation strategies for nested Krylov methods
• Relaxed Krylov subspace approximation

Cites work

• scientific article; zbMATH DE number 46317 (Why is no real title available?)
• scientific article; zbMATH DE number 479949 (Why is no real title available?)
• scientific article; zbMATH DE number 6327408 (Why is no real title available?)
• scientific article; zbMATH DE number 3324542 (Why is no real title available?)
• A fast adaptive multipole algorithm in three dimensions
• A fast algorithm for particle simulations
• Biomolecular electrostatics using a fast multipole BEM on up to 512 GPUs and a billion unknowns
• Explicit expressions for 3D boundary integrals in potential theory
• Explicit expressions for three-dimensional boundary integrals in linear elasticity
• GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems
• Inexact Krylov Subspace Methods for Linear Systems
• Inexact Matrix-Vector Products in Krylov Methods for Solving Linear Systems: A Relaxation Strategy
• The fast multipole boundary element method for potential problems: a tutorial
• Theory of Inexact Krylov Subspace Methods and Applications to Scientific Computing

Cited in

(2)

• A GMRES iterative solution of FEM‐BEM global systems in skin effect problems
• fmm-bem-relaxed

Describes a project that uses

Uses Software

• exafmm

This page was built for publication: Inexact GMRES iterations and relaxation strategies with fast-multipole boundary element method

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2673503)