We develop a multilevel method suitable for solving operator equations. Thi
s method combines the multiresolution structure of the spaces used to solve
the operator equation with a Gauss-Seidel strategy to solve the associated
matrix equations. We prove that this multilevel scheme has an optimal orde
r of convergence and provide an application of it to the solution of second
kind integral equations. (C) 2001 Academic Press.