We describe our three-dimensional numerical relativity code for the ev
olution of inhomogeneous cosmologies. During the evolution, the constr
aint 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.