Next, I do not know if you recognize this or not, but in case you interept the answer (even when its jogging on cluster), you should be ready to obtain the NR residuals for the last 3 (I usually use three but you can enter much more) iterations. The solution https://a-solution.info/