Secondly, I do not know if you recognize this or not, but in case you interept the answer (even though its functioning on cluster), you should be capable to acquire the NR residuals for the final 3 (I commonly use three however, you can enter even more) iterations. The solution https://a-solution.info/