Secondly, I do not know if you understand this or not, but in case you interept the answer (even though its running on cluster), you should be able to acquire the NR residuals for the final three (I typically use three however, you can enter even more) iterations. The solution https://a-solution.info/