Numerical methods for stochastic parabolic PDEs