If BFGS is used, then by default, the
Broyden [19]-Fletcher [20]-Goldfarb [21]-Shanno [22]
method will be used to optimize geometries. The older
Davidon-Fletcher-Powell method can be invoked by also specifying DFP.
This is intended to be used for comparison of the two methods.