Addlbfgs

I would merge the computation of rho into the line 161, it is not used anywhere else.