Nonlinear gyrokinetic simulation of ion temperature gradient turbulence based on a numerical Lie-transform perturbation method