Addlbfgs

See above, but my version of Numerical optimization uses s_{k-1} and y_{k-1}, not the new values of s and y.