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