ReceptionModelSolver.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2004 by Willem van Straten
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
// psrchive/More/Polarimetry/Pulsar/ReceptionModelSolver.h
 
#ifndef __ReceptionModelSolver_H
#define __ReceptionModelSolver_H
 
#include "Pulsar/ReceptionModel.h"
#include "MEAL/LeastSquares.h"
#include <fstream>
 
namespace Calibration
{
  class ReceptionModel::Solver : public MEAL::LeastSquares
  {
 
  public:
 
    Solver () { verbose = 0; }
    
    static bool report_chisq;
 
    virtual Solver* clone () const = 0;
 
    void set_equation (ReceptionModel*);
 
    void add_convergence_condition( Functor< bool(ReceptionModel*) > );
 
    void add_acceptance_condition( Functor< bool(ReceptionModel*) > );
 
    virtual void solve ();
 
    void set_verbosity (unsigned level) { verbose = level; }
    unsigned get_verbosity () const { return verbose; }
    
  protected:
 
    virtual void fit () = 0;
 
    std::vector< Functor< bool(ReceptionModel*) > > convergence_condition;
 
    std::vector< Functor< bool(ReceptionModel*) > > acceptance_condition;
 
    std::string prefit_report_filename;
 
    void count_infit ();
 
    void count_constraint ();
 
    void check_constraints ();
    void set_variances ();
    void check_solution ();
 
    Reference::To<ReceptionModel, false> equation;
 
    std::vector<CoherencyMeasurementSet>& get_data ();
 
    unsigned verbose;
    
 private:
 
    std::vector<bool> state_observed;
 
    std::vector<bool> path_observed;
  };
 
}
 
#endif