Polynomial.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2004 by Willem van Straten
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
// psrchive/More/MEAL/MEAL/Polynomial.h
 
#ifndef __Polynomial_H
#define __Polynomial_H
 
#include "MEAL/Univariate.h"
#include "MEAL/Scalar.h"
#include "MEAL/Parameters.h"
 
namespace MEAL {

  class Polynomial : public Univariate<Scalar> {
 
  public:

    Polynomial (unsigned ncoef = 0);

    Polynomial (const Polynomial&);

    Polynomial& operator = (const Polynomial&);

    Polynomial* clone () const;

    void resize (unsigned ncoef);
 
    unsigned get_ncoef () const { return get_nparam(); }
    void set_ncoef (unsigned nc) { resize (nc); }

    template<class T>
    Polynomial (const std::vector<T>& coefs) : parameters (this, coefs.size())
    {
      init ();
      for (unsigned ic=0; ic < coefs.size(); ic++)
        set_param (ic, coefs[ic]);
    }

    void set_abscissa_offset(double x0)
    { x_0 = x0; }
 
    // ///////////////////////////////////////////////////////////////////
    //
    // Function implementation
    //
    // ///////////////////////////////////////////////////////////////////

    std::string get_name () const;

    void set_param_name (unsigned icoeff, const std::string& name);

    void parse (const std::string& text);

    class Interface;

    TextInterface::Parser* get_interface ();
 
  protected:

    void calculate (double& x, std::vector<double>* grad=0);

    void print_parameters (std::string& text, const std::string& sep) const;
 
  private:

    double x_0;

    Parameters parameters;

    void init ();
 
  };
 
}
 
#endif