Addlbfgs

But s_k := x_{k+1}-x_k, the same for y_k. So we don't need l.first and prevGamma is simply gamma, which should be computed after updating l.y and l.s.