#ifndef _RESIDFUNC_H_
#define _RESIDFUNC_H_

class residFunc {
 public:
  virtual ~residFunc();
  virtual void resid(const vector<double> *x, void *data, vector<double> *f) = 0;
  virtual void deriv(const vector<double> *x, void *data, vector< vector<double> > *J) = 0;
  virtual void residAndDeriv(const vector<double> *x, void *data, vector<double> *f, vector< vector<double> > *J) = 0;
  virtual unsigned int getn() = 0;
  virtual unsigned int getp() = 0;
};

extern "C" {
#include <gsl/gsl_vector.h> 
#include <gsl/gsl_blas.h>
#include <gsl/gsl_multifit_nlin.h>

  int fGsl(const gsl_vector *x, void *data, gsl_vector *f);
  int dfGsl(const gsl_vector *x, void *data, gsl_matrix *J);
  int fdfGsl(const gsl_vector *x, void *data, gsl_vector *f, gsl_matrix *J);
}

#endif // _RESIDFUNC_H_