278 research outputs found
Fast iterative boundary element methods for high-frequency scattering problems in 3D elastodynamics
International audienceThe fast multipole method is an efficient technique to accelerate the solution of large scale 3D scattering problems with boundary integral equations. However, the fast multipole accelerated boundary element method (FM-BEM) is intrinsically based on an iterative solver. It has been shown that the number of iterations can significantly hinder the overall efficiency of the FM-BEM. The derivation of robust preconditioners for FM-BEM is now inevitable to increase the size of the problems that can be considered. The main constraint in the context of the FM-BEM is that the complete system is not assembled to reduce computational times and memory requirements. Analytic preconditioners offer a very interesting strategy by improving the spectral properties of the boundary integral equations ahead from the discretization. The main contribution of this paper is to combine an approximate adjoint Dirichlet to Neumann (DtN) map as an analytic preconditioner with a FM-BEM solver to treat Dirichlet exterior scattering problems in 3D elasticity. The approximations of the adjoint DtN map are derived using tools proposed in [40]. The resulting boundary integral equations are preconditioned Combined Field Integral Equations (CFIEs). We provide various numerical illustrations of the efficiency of the method for different smooth and non smooth geometries. In particular, the number of iterations is shown to be completely independent of the number of degrees of freedom and of the frequency for convex obstacles
Hierarchical interpolative factorization for elliptic operators: integral equations
This paper introduces the hierarchical interpolative factorization for
integral equations (HIF-IE) associated with elliptic problems in two and three
dimensions. This factorization takes the form of an approximate generalized LU
decomposition that permits the efficient application of the discretized
operator and its inverse. HIF-IE is based on the recursive skeletonization
algorithm but incorporates a novel combination of two key features: (1) a
matrix factorization framework for sparsifying structured dense matrices and
(2) a recursive dimensional reduction strategy to decrease the cost. Thus,
higher-dimensional problems are effectively mapped to one dimension, and we
conjecture that constructing, applying, and inverting the factorization all
have linear or quasilinear complexity. Numerical experiments support this claim
and further demonstrate the performance of our algorithm as a generalized fast
multipole method, direct solver, and preconditioner. HIF-IE is compatible with
geometric adaptivity and can handle both boundary and volume problems. MATLAB
codes are freely available.Comment: 39 pages, 14 figures, 13 tables; to appear, Comm. Pure Appl. Mat
- …