The l_bfgs_b includes another feature which I would like to take advantage of -- this is the combined evaluation of a function and a derivative.  But I am putting this off until a later PR.  Maybe this later PR could also include optional bounds?