SplineFit.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2010 by Paul Demorest
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
#ifndef __Pulsar_SplineFit_h
#define __Pulsar_SplineFit_h
 
#include "Estimate.h"
#include "Reference.h"
 
#include <gsl/gsl_bspline.h>
#include <gsl/gsl_matrix.h>

class SplineFit : public Reference::Able {
 
public:

  SplineFit ();

  virtual ~SplineFit ();

  void reset();

  void set_order(int n) { order=n; calculated=false; }

  int get_order() { return order; }

  void set_uniform_breaks(int nint);

  void add_data(double x, Estimate<double> y);

  void compute();

  double evaluate(double x);

  double evaluate_deriv(double x);

  double get_rchi2();
 
protected:

  std::vector<double> x;

  std::vector< Estimate<double> > y;

  std::vector<double> bp;

  int order;

  double chi2;

  int ndof;

  bool calculated;

  gsl_vector *coeffs;

  gsl_matrix *cov;

  bool check_range(double x);

  void interval_check(bool fix=false);

  void free_workspaces();

  gsl_bspline_workspace *bwork;
 
private:
 
};
 
#endif