"""
Optimisation routine based on a diagonal approximation of the Hessian matrix,
using a first order Taylor expansion of the change in the gradient vector, with
AR/MA/ARMA smoothing of diagonal terms
"""

def taylor_newton_optimiser():
    raise NotImplementedError