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 check_transformations ();

    void erase_transformation (unsigned index);

    void check_inputs ();

    void erase_input (unsigned index);
 
    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