Secondly, I do not know if you realize this or not, but for those who interept the answer (even when its jogging on cluster), you should be able to obtain the NR residuals for the final three (I typically use 3 however , you can enter a lot more) iterations. https://a-solution.info/