  We address a study of glueball masses in the confining regime of SU(2) in D=3
using an algorithm inspired by the multi-level scheme. Our method, which
exploits the locality of the action to achieve high precision results, is based
on a technique already used for compact QED, and generalises it to the
non-Abelian case. We discuss the main features of this method, in comparison
with other algorithms that have been used in similar studies.
