Secondly, I Don't know if you understand this or not, but for those who interept the answer (even though its managing on cluster), you should be capable to obtain the NR residuals for the final three (I commonly use three however, you can enter a lot more) iterations. The solution https://a-solution.info/