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