Theoretical and Practical Aspects of a Multigrid Method

Abstract

A multigrid method is described. A novel item is the use of incomplete $LU$ decomposition for smoothing. Numerical experiments show that its speed and robustness compare favorably with other multigrid methods. A fairly simple rate of convergence proof is presented.