Secondly, I Don't know if you are aware of this or not, but for those who interept the answer (even though its running on cluster), you should be able to obtain the NR residuals for the last three (I usually use 3 but you can enter more) iterations. The solution https://a-solution.info/