An Adaptive Fast Direct Solver for Boundary Integral Equations in Two Dimensions
YALE UNIV NEW HAVEN CT DEPT OF COMPUTER SCIENCE
Pagination or Media Count:
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 rs -rt between two points on the boundary is related to their corresponding distance r- s in the parameter space. The algorithm has an asymptotic complexity of Onlog2 n, where n is the number of nodes in the discretization. The performance of the algorithm is illustrated with several numerical examples.
- Numerical Mathematics
- Theoretical Mathematics