/**
 *  \file QuasiNewton.cpp
 *  \brief XXXX.
 *
 *  Copyright 2007-2013 IMP Inventors. All rights reserved.
 *
 */

#include "IMP/gsl/QuasiNewton.h"
#include <gsl/gsl_multimin.h>

IMPGSL_BEGIN_NAMESPACE

QuasiNewton::QuasiNewton(Model *m): GSLOptimizer(m) {
  initial_step_=.01;
  line_step_=.01;
  min_gradient_=.001;
}

void QuasiNewton::do_show(std::ostream &) const {
}

Float QuasiNewton::do_optimize(unsigned int nsteps) {
  /* The GnuWin32 Windows build of GSL is only 1.8, which doesn't have bfgs2;
     use the older less efficient bfgs optimizer instead. */
  const gsl_multimin_fdfminimizer_type *t
#if defined _MSC_VER
    =gsl_multimin_fdfminimizer_vector_bfgs;
#else
    =gsl_multimin_fdfminimizer_vector_bfgs2;
#endif
  return GSLOptimizer::optimize(nsteps, t, initial_step_,
                                line_step_, min_gradient_);
}

IMPGSL_END_NAMESPACE