/Users/jmonk/Physics/BSpline/Spline/BSpline.hh

#ifndef SPLINE_BSPLINE_HH
#define SPLINE_BSPLINE_HH

#include "Spline/DataPoint.hh"
#include "Spline/Matrix.hh"

#include <gsl/gsl_bspline.h>
#include <gsl/gsl_multifit.h>

#include <ostream>
#include <vector>

namespace Spline{
    
  using std::vector;
  
  class BSpline{
    
  public:
    
    BSpline(size_t degree = 4);
    
    void setKnotVector(const vector<double> &knots);
    
    void setNCoefficientsUniform(size_t nCoeffs, double min, double max);
    
    size_t nCoefficients()const;
    
    double basisFunction(size_t index, double affine)const;
    
    vector<double> basisFunctions(double affine)const;
    
    void fitLeastSquares(const DataSet &data);
    
    void fitLeastSquares(const DataSet &data, const DataSet &signalTemplate);
    
    double chi2()const;
    
    double chi2_dof()const;
    
    DataPoint evaluate(double affine)const;
    
    
    const vector<double> &knots()const;
    
    const vector<double> &coefficients()const;
    
    double signalCoefficientSize()const;
    
    const Matrix<double> &covarianceMatrix()const;
    
    const PointSet &controlPoints()const;
    
  private:
    
    size_t m_nCoeffs;
    size_t m_nBreakPoints;
    size_t m_degree;
    
    bool m_canEvaluate;
    
    vector<double> m_knots;
    
    gsl_vector *m_gslCoefficients;
    gsl_matrix *m_covarianceMatrix;
    
    gsl_vector *_basisFuncCoeffs()const;
    mutable gsl_vector *m_basisFuncCoeffs;
    
    bool m_haveSignalTemplate;
    
    // The same coefficients in a more usable std::vector
    vector<double> m_coefficients;
    // The same covariance matrix in a more friendly std::vector<std::vector>
    Matrix<double> m_covariance;
    
    mutable PointSet m_controlPoints;
    
    double m_chi2;
    double m_chi2_dof;
    
    gsl_bspline_workspace *m_bSplineWorkspace;
    
  };
}

#endif