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