Next, I don't know if you know this or not, but for those who interept the answer (even though its running on cluster), you should be able to get the NR residuals for the last three (I typically use three however you can enter a lot more) iterations. The solution https://a-solution.info/