  We show that any polyhomogeneous asymptotically hyperbolic
constant-mean-curvature solution to the vacuum Einstein constraint equations
can be approximated, arbitrarily closely in H\"older norms determined by the
physical metric, by shear-free smoothly conformally compact vacuum initial
data.
