Secondly, I don't know if you recognize this or not, but if you interept the solution (even though its operating on cluster), you ought to be able to obtain the NR residuals for the last three (I generally use 3 however you can enter all the more) iterations. The answer https://a-solution.info/