I\'m working now on the code improvement namely its performance which realizes automatic differentiation (AD). The core of the algorithm has the following implementation: