  We calculate the mean link in Landau gauge for Wilson and improved SU(3)
anisotropic gauge actions, using two loop perturbation theory and Monte Carlo
simulation employing an accelerated Langevin algorithm. Twisted boundary
conditions are employed, with a twist in all four lattice directions
considerably improving the (Fourier accelerated) convergence to an improved
lattice Landau gauge. Two loop perturbation theory is seen to predict the mean
link extremely well even into the region of commonly simulated gauge couplings
and so can be used remove the need for numerical tuning of self-consistent
tadpole improvement factors. A three loop perturbative coefficient is inferred
from the simulations and is found to be small. We show that finite size effects
are small and argue likewise for (lattice) Gribov copies and double Dirac
sheets.
