Hessian Calculations
We can define a search direction in the Newton Optimisation Method as:
where is the Hessian; a matrix of partial second order derivatives, and the is the Jacobian matrix stacked into a column vector form.
The partial derivatives of the Hessian matrix can be calculated with the grape algorithm (Gradient Ascent Pulse Engineering)