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