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