An direct solver for integral equations on the plane

• Published in: Applied and computational harmonic analysis Vol. 38; no. 2; pp. 284 - 317
• Main Authors: Corona, Eduardo, Martinsson, Per-Gunnar, Zorin, Denis
• Format: Journal Article
• Language: English
• Published: 01.03.2015
• Subjects: Asymptotic properties; Complexity; Compressing; Computation; Integral equations; Mathematical models; Matlab; Solvers
• DOI: 10.1016/j.acha.2014.04.002
• ISSN: 1063-5203

An efficient direct solver for volume integral equations with complexity for a broad range of problems is presented. The solver relies on hierarchical compression of the discretized integral operator, and exploits that off-diagonal blocks of certain dense matrices have numerically low rank. Technically, the solver is inspired by previously developed direct solvers for integral equations based on "recursive skeletonization" and "Hierarchically Semi-Separable" (HSS) matrices, but it improves on the asymptotic complexity of existing solvers by incorporating an additional level of compression. The resulting solver has optimal complexity for all stages of the computation, as demonstrated by both theoretical analysis and numerical examples. The computational examples further display good practical performance in terms of both speed and memory usage. In particular, it is demonstrated that even problems involving 10 super(7) unknowns can be solved to precision using a simple Matlab implementation of the algorithm executed on a single core.