Next, I don't know if you recognize this or not, but in case you interept the answer (even if its running on cluster), you should be able to obtain the NR residuals for the final three (I ordinarily use three however , you can enter much more) iterations. The solution https://a-solution.info/