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