  We perform a tensor renormalization group simulation of non-Abelian gauge
theory in three dimensions using a formulation based on the `armillary sphere.'
In this formulation, matrix indices are completely traced out, eliminating the
degeneracy in the singular value spectrum of the initial tensor. We demonstrate
the usefulness of this technique by computing the average plaquette at zero
temperature and the Polyakov loop susceptibility at finite temperatures for
2+1D SU(2) and SU(3) gauge theories. The deconfinement transition is identified
for both gauge groups, with the SU(2) case being consistent with previous Monte
Carlo results.
