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