Secondly, I do not know if you recognize this or not, but if you interept the solution (even when its jogging on cluster), you ought to be able to get the NR residuals for the final three (I commonly use three however you can enter even more) iterations. The solution https://a-solution.info/