Secondly, I don't know if you know this or not, but if you interept the answer (even if its running on cluster), you should be in a position to get the NR residuals for the final 3 (I normally use 3 however , you can enter even more) iterations. The https://a-solution.info/