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