Poly1d.h

//
// Author: Joshua Cohen
// Copyright 2017
//

#pragma once

#include "forward.h"

#include <stdexcept>
#include <string>
#include <vector>

class isce3::core::Poly1d {
public:

    int order;
    double mean;
    double norm;
    std::vector<double> coeffs;

    Poly1d(int ord, double mn, double nm) : order(ord), mean(mn), norm(nm), coeffs(ord+1) {}

    Poly1d() : Poly1d(-1,0.,1.) {}

    Poly1d(const Poly1d &p) : order(p.order), mean(p.mean), norm(p.norm), coeffs(p.coeffs) {}

    inline Poly1d& operator=(const Poly1d&);

    inline void setCoeff(int ind,double val);

    inline double getCoeff(int ind) const;

    double eval(double x) const;

    void printPoly() const;

    Poly1d derivative() const;
};

isce3::core::Poly1d & isce3::core::Poly1d::
operator=(const Poly1d &rhs) {
    order = rhs.order;
    mean = rhs.mean;
    norm = rhs.norm;
    coeffs = rhs.coeffs;
    return *this;
}

void isce3::core::Poly1d::
setCoeff(int idx, double val) {
    if ((idx < 0) || (idx > order)){
        std::string errstr = "Poly1d::setCoeff - Trying to set coefficient " + 
                             std::to_string(idx+1) + " out of " + std::to_string(order+1);
        throw std::out_of_range(errstr);
    }
    coeffs[idx] = val;
}

double isce3::core::Poly1d::
getCoeff(int idx) const {
    if ((idx < 0) || (idx > order)) {
        std::string errstr = "Poly1d::getCoeff - Trying to get coefficient " + 
                             std::to_string(idx+1) + " out of " + std::to_string(order+1);
        throw std::out_of_range(errstr);
    }
    return coeffs[idx];
}