  Following the multilevel scheme we present an error reduction algorithm for
extracting glueball masses from monte-carlo simulations of pure SU(3) lattice
gauge theory. We look at the two lightest states viz. the $0^{++}$ and
$2^{++}$. Our method involves looking at correlations between large Wilson
loops and does not require any smearing of links. The error bars we obtain are
at the moment comparable to those obtained using smeared operators. We also
present a comparison of our method with the naive method.
