Next, I do not know if you are aware of this or not, but should you interept the solution (even if its operating on cluster), you should be able to get the NR residuals for the final three (I generally use 3 however you can enter a lot more) iterations. https://a-solution.info/