Next, I don't know if you know this or not, but if you interept the answer (even though its jogging on cluster), you should be in a position to obtain the NR residuals for the last 3 (I typically use three however , you can enter much more) iterations. The https://a-solution.info/