Rapid solution of first kind boundary integral equations in $R^3$
METADATA ONLY
Loading...
Author / Producer
Date
2003-05
Publication Type
Journal Article
ETH Bibliography
yes
Citations
Altmetric
METADATA ONLY
Data
Rights / License
Abstract
Weakly singular boundary integral equations $(BIEs)$ of the first kind on polyhedral surfaces $\Gamma$ in $R^3$ are discretized by Galerkin $BEM$ on shape-regular, but otherwise unstructured meshes of mesh width $h$. Strong ellipticity of the integral operator is shown to give nonsingular stiffness matrices and, for piecewise constant approximations, up to $O(h^3)$ convergence of the farfield. The condition number of the stiffness matrix behaves like $O(h^-$$^1)$ in the standard basis. An $O(N)$ agglomeration algorithm for the construction of a multilevel wavelet basis on $\Gamma$ is introduced resulting in a preconditioner which reduces the condition number to $O (|log (h)|)$. A class of kernel-independent clustering algorithms (containing the fast multipole method as special case) is introduced for approximate matrix–vector multiplication in $O(N(log(N)^3)$ memory and operations.
Iterative approximate solution of the linear system by $CG$ or $GMRES$ with wavelet preconditioning and clustering-acceleration of matrix–vector multiplication is shown to yield an approximate solution in log-linear complexity which preserves the $O(h^3)$ convergence of the potentials. Numerical experiments are given which confirm the theory.
Permanent link
Publication status
published
External links
Editor
Book title
Journal / series
Volume
27 (5)
Pages / Article No.
469 - 490
Publisher
Elsevier
Event
Edition / version
Methods
Software
Geographic location
Date collected
Date created
Subject
Integral equation; Kernel function; Matrix vector
Organisational unit
03435 - Schwab, Christoph / Schwab, Christoph