nlls_solver_gsl_lmder.h

Go to the documentation of this file.

#ifndef NLLS_SOLVER_GSL_LMDER_H
#define NLLS_SOLVER_GSL_LMDER_H
#include <nlls_solver_gsl.h>

namespace FullPhysics {
/******************************************************************
  This class is the implementation of J. J. More's
  version of the Levenberg-Marquardt NLLS solver with 
  one difference.  The diagonal weighing matrix in
  described in the More's paper is replace with the
  identity matrix.
*******************************************************************/
class NLLSSolverGSLLMDER : 
    public NLLSSolverGSL {

public:

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  NLLSSolverGSLLMDER(const boost::shared_ptr<NLLSProblem>& p, int max_cost_function_calls, 
                     double dx_tol_abs=0.000001, double dx_tol_rel=0.000001, double g_tol=6.0555e-06, 
                     bool vrbs=false)
    : NLLSSolverGSL(p, max_cost_function_calls, dx_tol_abs, dx_tol_rel, g_tol, vrbs)
  {}

  virtual ~NLLSSolverGSLLMDER() {}

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  virtual void print(std::ostream& Os) const 
  { Os << "NLLSSolverGSLLMDER"; }

protected:

  virtual const gsl_multifit_fdfsolver_type* get_gsl_multifit_fdfsolver()
  { return gsl_multifit_fdfsolver_lmder; }

};
}
#endif