NthOrderPolynomial.cpp

 
/* SPDX-License-Identifier: CC0-1.0 */
#include <math.h>
#include <QString>
 
#include "FileName.h"
#include "Constants.h"
#include "IException.h"
#include "NthOrderPolynomial.h"
 
using namespace std;
namespace Isis {

  NthOrderPolynomial::NthOrderPolynomial(int degree) : 
    Isis::BasisFunction("NthOrderPolynomial", 2, degree) {
    p_degree = degree;
  }
  

  void NthOrderPolynomial::Expand(const std::vector<double> &vars) {
 
    if((int) vars.size() != Variables()) {
      QString mess = "Number of variables given (" + QString::number(vars.size())
          + ") does not match expected (" + QString::number(Variables()) + ")!";
      throw IException(IException::Programmer, mess, _FILEINFO_);
    }
    
    double t1 = vars[0];
    double t2 = vars[1];
    p_terms.clear();
    for (int i = p_degree; i >= 1; i--) {
      p_terms.push_back(pow(t1, i) - pow(t2, i));
    }
    return;
  }
} // end namespace isis