BFGS improvements

These commits improve the BFGS-based optimization routine and make BFGS the default algorithm when a point_estimate is requested.