SingleAxisSolver.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2004 by Willem van Straten
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
// psrchive/More/Polarimetry/Pulsar/SingleAxisSolver.h
 
#ifndef __Calibration_SingleAxisSolver_H
#define __Calibration_SingleAxisSolver_H
 
#include "Quaternion.h"
#include "Stokes.h"
#include "Estimate.h"
#include "MEAL/ScalarMath.h"
 
namespace Calibration {
 
  class SingleAxis;

  class SingleAxisSolver : public Reference::Able {
 
  public:

    SingleAxisSolver ();

    SingleAxisSolver (const SingleAxisSolver&);

    ~SingleAxisSolver ();

    SingleAxisSolver& operator = (const SingleAxisSolver&);

    void set_input (const Stokes< Estimate<double> >& stokes);

    void set_output (const Stokes< Estimate<double> >& stokes);

    void solve (SingleAxis* model);
 
  protected:

    Quaternion<MEAL::ScalarMath,Hermitian> input;

    Quaternion<MEAL::ScalarMath,Hermitian> output;

    Vector<3, MEAL::ScalarMath> axis;

    MEAL::ScalarMath gain;

    MEAL::ScalarMath diff_gain;

    MEAL::ScalarMath diff_phase;
 
  private:

    void init ();
 
  };
 
}
 
#endif