An adaptive fast direct solver for boundary integral equations in two dimensions

Abstract We describe an algorithm for the rapid direct solution of linear algebraic systems arising from the discretization of boundary integral equations of potential theory in two dimensions. The algorithm is combined with a scheme that adaptively rearranges the parameterization of the boundary in order to minimize the ranks of the off-diagonal blocks in the discretized operator, thus obviating the need for the user to supply a parameterization r of the boundary for which the distance ‖ r ( s ) − r ( t ) ‖ between two points on the boundary is related to their corresponding distance | s − t | in the parameter space. The algorithm has an asymptotic complexity of O ( N log 2 N ) , where N is the number of nodes in the discretization. The performance of the algorithm is illustrated with several numerical examples.

[1]  Y. Meyer,et al.  Wavelets: Calderón-Zygmund and Multilinear Operators , 1997 .

[2]  G. Verchota.  Layer potentials and regularity for the Dirichlet problem for Laplace's equation in Lipschitz domains , 1984 .

[3]  GuMing,et al.  Efficient algorithms for computing a strong rank-revealing QR factorization , 1996 .

[4]  Ming Gu,et al.  Efficient Algorithms for Computing a Strong Rank-Revealing QR Factorization , 1996, SIAM J. Sci. Comput..

[5]  G. Verchota,et al.  Layer potentials and boundary value problems for laplace's equation on lipschitz domains : a thesis submitted to the faculty of the graduate school of the University of Minnesota , 1982 .

[6]  J. Stoer,et al.  Introduction to Numerical Analysis , 2002 .

[7]  J. CARRIERt,et al.  A FAST ADAPTIVE MULTIPOLE ALGORITHM FOR PARTICLE SIMULATIONS * , 2022 .

[8]  Wolfgang Hackbusch,et al.  A Sparse Matrix Arithmetic Based on H-Matrices. Part I: Introduction to H-Matrices , 1999, Computing.

[9]  V. Rokhlin,et al.  A fast randomized algorithm for the approximation of matrices ✩ , 2007 .

[10]  O. D. Kellogg.  Foundations of potential theory , 1934 .

[11]  G. Folland.  Introduction to Partial Differential Equations , 1976 .

[12]  Steffen Börm,et al.  Data-sparse Approximation by Adaptive ℋ2-Matrices , 2002, Computing.

[13]  Shivkumar Chandrasekaran,et al.  A Fast ULV Decomposition Solver for Hierarchically Semiseparable Representations , 2006, SIAM J. Matrix Anal. Appl..

[14]  V. Hutson.  Integral Equations , 1967, Nature.

[15]  V. Rokhlin,et al.  A randomized algorithm for the approximation of matrices , 2006 .

[16]  G. Chandler,et al.  Galerkin's method for boundary integral equations on polygonal domains , 1984, The Journal of the Australian Mathematical Society. Series B. Applied Mathematics.

[17]  Gene H. Golub,et al.  Matrix computations , 1983 .

[18]  Vladimir Rokhlin,et al.  On the numerical solution of two-point boundary value problems II , 1994 .

[19]  Leslie Greengard,et al.  A fast algorithm for particle simulations , 1987 .

[20]  J. Navarro-Pedreño.  Numerical Methods for Least Squares Problems , 1996 .

[21]  V. Rokhlin.  Rapid solution of integral equations of classical potential theory , 1985 .

[22]  Carlos E. Kenig,et al.  ELLIPTIC BOUNDARY VALUE PROBLEMS ON LIPSCHITZ DOMAINS , 1987 .

[23]  V. Rokhlin,et al.  A fast direct solver for boundary integral equations in two dimensions , 2003 .

[24]  Vladimir Rokhlin,et al.  An Improved Fast Multipole Algorithm for Potential Fields , 1998, SIAM J. Sci. Comput..

[25]  Johan Helsing,et al.  Corner singularities for elliptic problems: Integral equations, graded meshes, quadrature, and compressed inverse preconditioning , 2008, J. Comput. Phys..

[26]  Weng Cho Chew,et al.  An N2 algorithm for the multiple scattering solution of N scatterers , 1989 .

[27]  H. Keller.  Numerical Solution of Two Point Boundary Value Problems , 1976 .

[28]  Yu Chen,et al.  A Fast, Direct Algorithm for the Lippmann–Schwinger Integral Equation in Two Dimensions , 2002, Adv. Comput. Math..

[29]  W. Hackbusch.  A Sparse Matrix Arithmetic Based on $\Cal H$-Matrices. Part I: Introduction to ${\Cal H}$-Matrices , 1999, Computing.

[30]  Rainer Kress,et al.  A Nyström method for boundary integral equations in domains with corners , 1990 .

[31]  Gregory Beylkin,et al.  LU Factorization of Non-standard Forms and Direct Multiresolution Solvers , 1998 .