Adaptive BEM with inexact PCG solver yields almost optimal computational costs

DOI10.1007/S00211-018-1011-1MaRDI QIDQ670306zbMATHOpenAlexFDO

Authors Thomas Führer, Alexander Haberl, Dirk Praetorius, Stefan Schimanko

Publication date 18 March 2019

Published in Numerische Mathematik (Search for Journal in Brave)

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

zbMATH Keywords

Laplace operatorboundary element methodcondition numberspreconditioned conjugate gradient methodweakly singular integral equationadditive Schwarz preconditioner

Mathematics Subject Classification ID

Complexity and performance of numerical algorithms (65Y20) Preconditioners for iterative methods (65F08) Numerical methods for integral equations (65R20) Numerical solution of discretized equations for boundary value problems involving PDEs (65N22) Boundary element methods for boundary value problems involving PDEs (65N38) Mesh generation, refinement, and adaptive methods for boundary value problems involving PDEs (65N50) Rate of convergence, degree of approximation (41A25) Integral equations with kernels of Cauchy type (45E05)

Abstract: We consider the preconditioned conjugate gradient method (PCG) with optimal preconditioner in the frame of the boundary element method (BEM) for elliptic first-kind integral equations. Our adaptive algorithm steers the termination of PCG as well as the local mesh-refinement. Besides convergence with optimal algebraic rates, we also prove almost optimal computational complexity. In particular, we provide an additive Schwarz preconditioner which can be computed in linear complexity and which is optimal in the sense that the condition numbers of the preconditioned systems are uniformly bounded. As model problem serves the 2D or 3D Laplace operator and the associated weakly-singular integral equation with energy space $widetilde{H}^{-1/2}(Gamma)$. The main results also hold for the hyper-singular integral equation with energy space $H^{1/2}(Gamma)$.

Recommendations

• Optimal additive Schwarz preconditioning for adaptive 2D IGA boundary element methods
• Quasi-optimal convergence rate for an adaptive boundary element method
• Optimal additive Schwarz methods for the \(hp\)-BEM: the hypersingular integral operator in 3D on locally refined meshes
• Adaptive BEM with optimal convergence rates for the Helmholtz equation
• Adaptive BEM for elliptic PDE systems. I: Abstract framework, for weakly-singular integral equations

Cites work

• scientific article; zbMATH DE number 2113718 (Why is no real title available?)
• scientific article; zbMATH DE number 777426 (Why is no real title available?)
• scientific article; zbMATH DE number 6159604 (Why is no real title available?)
• scientific article; zbMATH DE number 1446717 (Why is no real title available?)
• A Convergent Adaptive Algorithm for Poisson’s Equation
• A Posteriori Error Estimates for Boundary Element Methods
• A posteriori error estimate and h-adaptive algorithm on surfaces from Symm's integral equation
• Adaptive FEM with coarse initial mesh guarantees optimal convergence rates for compactly perturbed elliptic problems
• Adaptive FEM with optimal convergence rates for a certain class of nonsymmetric and possibly nonlinear problems
• Adaptive boundary element methods with convergence rates
• Adaptive finite element methods with convergence rates
• An a posteriori error estimate for a first-kind integral equation
• An adaptive finite element eigenvalue solver of asymptotic quasi-optimal computational complexity
• Averaging Techniques for the Effective Numerical Solution of Symm's Integral Equation of the First Kind
• Axioms of adaptivity
• Data Oscillation and Convergence of Adaptive FEM
• Discrete extension operators for mixed finite element spaces on locally refined meshes
• Efficiency and Optimality of Some Weighted-Residual Error Estimator for Adaptive 2D Boundary Element Methods
• Extension by zero in discrete trace spaces: inverse estimates
• HILBERT -- a MATLAB implementation of adaptive 2D-BEM. \(\underline {\text H}\)ilbert \(\underline {\text I}\)s a \(\underline {\text L}\)ovely \(\underline {\text B}\)oundary \(\underline {\text E}\)lement \(\underline {\text R}\)esearch \(\underline {\text T}\)ool
• Hierarchical matrices: algorithms and analysis
• Interplay between discretization and algebraic computation in adaptive numerical solution of elliptic PDE problems
• Local inverse estimates for non-local boundary integral operators
• Local multigrid in \(\boldsymbol{H} (\mathbf{curl})\)
• On 2D newest vertex bisection: optimality of mesh-closure and \(H ^{1}\)-stability of \(L _{2}\)-projection
• Optimal additive Schwarz methods for the \(hp\)-BEM: the hypersingular integral operator in 3D on locally refined meshes
• Optimal additive Schwarz preconditioning for hypersingular integral equations on locally refined triangulations
• Optimal multilevel methods for \(H\)(grad), \(H\)(curl), and \(H\)(div) systems on graded and unstructured grids
• Optimal preconditioning for the symmetric and nonsymmetric coupling of adaptive finite elements and boundary elements
• Optimality of a standard adaptive finite element method
• Quasi-Optimal Convergence Rate for an Adaptive Finite Element Method
• Quasi-optimal convergence rate for an adaptive boundary element method
• Quasi-optimal convergence rates for adaptive boundary element methods with data approximation. I: Weakly-singular integral equation
• Quasi-optimal convergence rates for adaptive boundary element methods with data approximation. II: Hyper-singular integral equation
• Rate optimal adaptive FEM with inexact solver for nonlinear operators
• Residual-based a posteriori error estimate for hypersingular equation on surfaces
• Solving boundary integral problems with BEM++
• Stable multilevel splittings of boundary edge element spaces
• Stopping criteria for adaptive finite element solvers
• The Conditioning of Boundary Element Equations on Locally Refined Meshes and Preconditioning by Diagonal Scaling
• The completion of locally refined simplicial partitions created by bisection
• Uniform Convergence of Adaptive Multigrid Methods for Elliptic Problems and Maxwell's Equations
• Uniform preconditioners for problems of negative order

Cited in

(12)

• An optimal adaptive Fictitious Domain Method
• Rate optimality of adaptive finite element methods with respect to overall computational costs
• Uniform preconditioners of linear complexity for problems of negative order
• Dörfler marking with minimal cardinality is a linear complexity problem
• Adaptive boundary element methods for the computation of the electrostatic capacity on complex polyhedra
• Optimal additive Schwarz preconditioning for adaptive 2D IGA boundary element methods
• Optimal quasi-diagonal preconditioners for pseudodifferential operators of order minus two
• A short note on plain convergence of adaptive least-squares finite element methods
• Adaptive IGAFEM with optimal convergence rates: T-splines
• Adaptive Uzawa algorithm for the Stokes equation
• Adaptive BEM with optimal convergence rates for the Helmholtz equation
• Multilevel decompositions and norms for negative order Sobolev spaces

Describes a project that uses

Uses Software

• Matlab
• BEM++
• HILBERT

This page was built for publication: Adaptive BEM with inexact PCG solver yields almost optimal computational costs

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