Secondly, I do not know if you know this or not, but if you interept the answer (even if its jogging on cluster), you should be capable to get the NR residuals for the last 3 (I generally use 3 however, you can enter all the more) iterations. The solution https://a-solution.info/