The line search algorithm is based on a quasi-Newton variable storage method which was implemented by [21].
TO BE CONTINUED...