Numerical methods for nuclear mean-field dynamics