Abstract
A description and explanation of a simple multigrid algorithm for solving finite element systems is given. Numerical results from an implementation are reported for a number of elliptic equations, including cases with singular coefficients and indefinite equations. The method shows the high efficiency, essentially independent of the grid spacing, predicted by the theory.