  The quenched spectrum of glueballs with positive charge parity is calculated
from two-body bound state equations. As input, a self-contained solution for
the primitively divergent correlation functions from Dyson-Schwinger equations
is used. It only has one parameter to be set which is the physical scale. An
important feature of this setup is the consistent construction of the bound
state kernels along the same lines as the equations from which the input was
obtained. Keeping only the one-particle exchanges, already good agreement with
lattice results is obtained. For the tensor glueball, we present first results
including two-loop contributions, elevating its calculation to the same level
of truncation as for the spin zero glueballs for which such calculations have
been done previously.
