LeastSquares.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2008 by Willem van Straten
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
// psrchive/More/MEAL/MEAL/LeastSquares.h
 
#ifndef __MEAL_LeastSquares_H
#define __MEAL_LeastSquares_H
 
#include "Reference.h"
 
namespace MEAL {
 
  class LeastSquares : public Reference::Able
  {
 
  public:
 
    typedef std::vector< std::vector<double> > matrix;
 
    LeastSquares ();
 
    virtual std::string get_name () const = 0;
 
    void set_maximum_iterations (unsigned maximum_iterations);
 
    unsigned get_maximum_iterations () const { return maximum_iterations; }
 
    void set_convergence_chisq (float chisq);
 
    void set_convergence_delta (float delta_parameters);
 
    void set_maximum_reduced_chisq (float maximum_reduced_chisq);
 
    void set_debug (bool flag = true) { debug = flag; }
 
    void set_solved (bool val) { solved = val; }
 
    bool get_solved () const { return solved; }
 
    void set_singular (bool val) { singular = val; }
 
    bool get_singular () const { return singular; }
 
    unsigned get_iterations () const { return iterations; }
 
    float get_chisq () const { return best_chisq; }
 
    unsigned get_nfree () const { return nfree; }
 
    unsigned get_nparam_infit () const { return nparam_infit; }
 
    unsigned get_ndat_constraint () const { return ndat_constraint; }
 
    void get_covariance (matrix& c) const { c = covariance; }
 
  protected:
 
    unsigned maximum_iterations;
 
    float convergence_chisq;
 
    float convergence_delta;
 
    float maximum_reduced;
 
    unsigned iterations;
 
    float best_chisq;
 
    unsigned nfree;
 
    unsigned nparam_infit;
 
    unsigned ndat_constraint;
 
    bool debug;
 
    bool solved;
 
    bool singular;
 
    matrix covariance;
 
  };
 
}
 
#endif