Polynomial.h

Go to the documentation of this file.

#ifndef _utl_Polynomial_h_
#define _utl_Polynomial_h_

#include <utl/AugerException.h>
#include <utl/Math.h>
#include <vector>
#include <iostream>


namespace utl {

  class Polynomial {

  public:
    Polynomial(const double coeff = 0) : fCoeff(1, coeff) { }

    Polynomial(const std::vector<double>& coeff) : fCoeff(coeff) { Check(); }

    template<typename T, unsigned int n>
    Polynomial(const T (&coeff)[n]) : fCoeff(coeff, coeff + n) { Check(); }

    double operator()(const double x) const { return EvalPoly(fCoeff, x); }

    const std::vector<double>& GetCoefficients() const { return fCoeff; }

    bool
    IsZero()
      const
    {
      for (const auto& c : fCoeff)
        if (c)
          return false;
      return true;
    }

    explicit operator bool() const { return !IsZero(); }

  private:
    void
    Check()
      const
    {
      if (fCoeff.empty())
        throw DoesNotComputeException("Polynomial should have at least one coefficient!");
    }

    std::vector<double> fCoeff;

  };


  std::istream& operator>>(std::istream& is, Polynomial& p);
  std::ostream& operator<<(std::ostream& is, const Polynomial& p);

}


#endif