Next, I Do not know if you know this or not, but in case you interept the answer (even if its functioning on cluster), you should be equipped to obtain the NR residuals for the last three (I typically use three however you can enter much more) iterations. The solution https://a-solution.info/