Secondly, I Do not know if you know this or not, but should you interept the answer (even if its running on cluster), you ought to be able to get the NR residuals for the final three (I ordinarily use three however you can enter a lot more) iterations. The https://a-solution.info/