Seminar für Angewandte Mathematik
ETH-Zürich
Rämistrasse 101
CH-8092 Zürich, Switzerland
The implementation of a wavelet-based Galerkin discretization of the double layer potential operator on polyhedral surfaces is described. The algorithm generates an approximate stiffness matrix with O(N(log N)2) entries in O(N(log N)3) operations where N is the number of degrees of freedom on the boundary. The condition number of the compressed stiffness matrix is bounded uniformly with respect to N. This work was presented in the 13th GAMM-Seminar on Numerical Treatment of Multi-Scale Problems, January 1997.