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