A fast operator perturbation method for the solution of the special relativistic equation of radiative transfer in spherical symmetry