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);
 
    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 parse (const std::string& text);
 
  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