Rapid solution of integral equations of classical potential theory

From MaRDI portal





A rapid solution of integral equations is described which is applicable to Dirichlet and Neumann boundary value problems for the Laplace equation. The method has computational requirements proportional to n, where n is the number of nodes on the boundary. It uses the classical approach of transforming the problem to an integral equation for the single or double layer potential on the boundary. This equation is then discretized using the Nyström algorithm associated with the trapezoidal quadrature rule. The resulting system is solved by the generalized conjugate residual algorithm (GCRA). The decrease of computational requirements is achieved by reducing the number of operations needed for applying a matrix to a vector in the process of solving by the GCRA. This is made possible by approximations based on harmonic expansions. The algorithm is tested on some standard problems which confirm the theoretically predicted properties. It must, however, be kept in mind that the method is superior to fast Poisson solvers only when the solution in a limited number of points outside the boundary is required.




Cited in
(only showing first 100 items - show all)








This page was built for publication: Rapid solution of integral equations of classical potential theory

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1094119)