CurveFitVal.h

#pragma once
 
#include <exception>
#include <iostream>
#include <stdexcept>
#include <vector>
class CurveFitVal {
  public:
    CurveFitVal(std::vector<double> xcoord, std::vector<double> ycoord, const std::size_t pdegree,
                const double loadFactor = 0)
        : pdegree(pdegree), xcoord(std::move(xcoord)), ycoord(std::move(ycoord)), loadFactor(loadFactor) {
        if (this->xcoord.size() != this->ycoord.size()) {
            throw std::runtime_error("X and Y coordinate vectors must be the same size");
        }
 
        coeff = Fit_Coefficients();
    }
 
    double calculate() const;
    double calculate(double) const;
 
  public:
    // coeff public so quadratic equation can be used
    std::vector<double> coeff;
 
  private:
    std::size_t pdegree;
    std::vector<double> xcoord, ycoord;
    double loadFactor;
 
    std::vector<double> Fit_Coefficients();
};