Change to diffusion solver to update operator for new keff