you might need to tweak the surface semi-diameter/radius a bit to make rays traceable towards the image plane. your calculation based on PW method is an approximation I would assume, so make a little change to start with for few cycles of optimization then you should be able to get around the problem.