polynomial_eval.h

Go to the documentation of this file.

#ifndef POLYNOMIAL_EVAL_H
#define POLYNOMIAL_EVAL_H

#include "array_ad.h"

namespace FullPhysics {
/****************************************************************/
class Poly1d : public Printable<Poly1d> {
public:

  //-----------------------------------------------------------------------
  //-----------------------------------------------------------------------
  Poly1d(const ArrayAd<double, 1>& Coefficients, const bool Decreasing_order = true) : coeffs_(Coefficients), decreasing_order_(Decreasing_order) {}

  //-----------------------------------------------------------------------
  //-----------------------------------------------------------------------
  double operator()(double Value) const;
  AutoDerivative<double> operator()(const AutoDerivative<double>& Value) const;

  //-----------------------------------------------------------------------
  //-----------------------------------------------------------------------
  blitz::Array<double,1> operator()(const blitz::Array<double,1>& Arr) const;
  ArrayAd<double, 1> operator()(const ArrayAd<double, 1>& Arr) const;
  
  virtual void print(std::ostream& Os) const;

private:

  blitz::Array<int,1> eval_indexes() const;

  ArrayAd<double, 1> coeffs_;
  bool decreasing_order_;
};

}
#endif