Numerical study of discrete models in the class of the nonlinear molecular beam epitaxy equation