A Method to Accelerate Convergence and to Preserve Radiative Energy Balance in Solving the P1 Equation by Iterative Methods