chisq_convergence.h

Go to the documentation of this file.

#ifndef CHISQ_CONVERGENCE_H
#define CHISQ_CONVERGENCE_H
#include "convergence_check.h"

namespace FullPhysics {
/****************************************************************/
class ChisqConvergence : public ConvergenceCheck {
public:
  ChisqConvergence(double stopping_criteria = 0.001,
                   double dropf = 0.1, double boostf = 10,
                   double min_chisq = 0.01,
                   int max_iteration = 50);
  virtual ~ChisqConvergence() {}
  virtual void convergence_check(const FitStatistic& fit_stat_last,
                                 FitStatistic& fit_stat,
                                 bool& has_converged,
                                 bool& convergence_failed,
                                 double& gamma,
                                 bool& step_diverged);
  virtual void evaluate_quality(FitStatistic& fit_stat,
                                const blitz::Array<double, 1>& Residual,
                                const blitz::Array<double, 1>& Residual_cov_diag);
  virtual void print(std::ostream& Os) const;
private:
  double stopping_criteria_, dropf_, boostf_, min_chisq_, max_iteration_;
};
}
#endif