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

DOI10.1007/S10444-022-09932-8zbMATH Open1491.35350arXiv1506.05957OpenAlexW4229443128MaRDI QIDQ2673503FDOQ2673503

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)

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}.

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

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

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)

Cites Work

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

Cited In (2)

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

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)