  We describe our 3-dimensional numerical relativity code for the evolution of
inhomogeneous cosmologies. During the evolution the constraint equations are
monitored but not solved. The code has been tested against perturbation theory
with good results. We present some runs of inhomogeneous inflation with
strongly inhomogeneous initial data.
