A Robust and Accurate Solver of Laplace's Equation with General Boundary Conditions on General Domains in the Plane

Abstract

A robust and general solver for Laplace's equation on the interior of a simply connected domain in the plane is described and tested. The solver handles general piecewise smooth domains and Dirichlet, Neumann, and Robin boundary conditions. It is based on an integral equation formulation of the problem. Difficulties due to changes in boundary conditions and corners, cusps, or other examples of non-smoothness of the boundary are handled using a recent technique called recursive compressed inverse preconditioning. The result is a rapid and very accurate solver which is general in scope, and its performance is demonstrated via some challenging numerical tests.