Figured out boundary bug for Krylov.  When doing final free iterations to pick up the angular and boundary fluxes, I was forgetting to set the boundary condition again.  Since vacuum and reflect do nothing, I never caught it before.