polynomial1d.h

#pragma once
 
#include <cmath>
#include <vector>
 
class Polynomial1d
{
   public:
    // A helper struct used to represent constraints
    typedef struct
    {
        double input;
        double output;
    } Constraint;
 
    explicit Polynomial1d();
 
    explicit Polynomial1d(const std::vector<double> &coeffs);
 
    explicit Polynomial1d(const std::initializer_list<double> &coeffs);
 
    explicit Polynomial1d(const std::vector<Constraint> constraints);
 
    double getCoeff(unsigned int order) const;
 
    void setCoeff(unsigned int order, double coeff);
 
    unsigned int getOrder() const;
 
    double valueAt(double val) const;
 
   private:
    // the coefficients for the polynomial
    // of the form coeffs[0] + coeffs[1]*x^(1)
    // + ... + coeff[n-1]*x^(n-1)
    std::vector<double> coeffs;
};
 
Polynomial1d operator+(const Polynomial1d &p1, const Polynomial1d &p2);
 
Polynomial1d operator-(const Polynomial1d &p1, const Polynomial1d &p2);
 
Polynomial1d operator*(const Polynomial1d &p1, const Polynomial1d &p2);
 
Polynomial1d &operator+=(Polynomial1d &p1, const Polynomial1d &p2);
 
Polynomial1d &operator-=(Polynomial1d &p1, const Polynomial1d &p2);
 
Polynomial1d &operator*=(Polynomial1d &p1, const Polynomial1d &p2);
 
bool operator==(const Polynomial1d &p1, const Polynomial1d &p2);