Secondly, I do not know if you already know this or not, but if you interept the answer (even if its running on cluster), you should be ready to obtain the NR residuals for the last 3 (I typically use three however, you can enter a lot more) iterations. The https://a-solution.info/