Radiation trapping simulations using the propagator function method: complete and partial frequency redistribution