Next, I do not know if you know this or not, but in case you interept the solution (even though its functioning on cluster), try to be able to obtain the NR residuals for the final three (I commonly use 3 however , you can enter much more) iterations. The https://a-solution.info/