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