# distutils: libraries = GSL_LIBRARIES
cdef extern from "gsl/gsl_interp.h":

  ctypedef struct gsl_interp_accel

  ctypedef struct gsl_interp_type
  ctypedef struct gsl_interp

  gsl_interp_type * gsl_interp_linear
  gsl_interp_type * gsl_interp_polynomial
  gsl_interp_type * gsl_interp_cspline
  gsl_interp_type * gsl_interp_cspline_periodic
  gsl_interp_type * gsl_interp_akima
  gsl_interp_type * gsl_interp_akima_periodic

  gsl_interp_accel * gsl_interp_accel_alloc()

  size_t gsl_interp_accel_find(gsl_interp_accel * a,  double x_array[], size_t size, double x)

  int gsl_interp_accel_reset (gsl_interp_accel * a)

  void gsl_interp_accel_free(gsl_interp_accel * a)

  gsl_interp * gsl_interp_alloc( gsl_interp_type * T, size_t n)

  int gsl_interp_init(gsl_interp * obj,  double xa[],  double ya[], size_t size)

  char * gsl_interp_name( gsl_interp * interp)
  unsigned int gsl_interp_min_size( gsl_interp * interp)


  int gsl_interp_eval_e( gsl_interp * obj,
                     double xa[],  double ya[], double x,
                    gsl_interp_accel * a, double * y)

  double gsl_interp_eval( gsl_interp * obj,
                   double xa[],  double ya[], double x,
                  gsl_interp_accel * a)

  int gsl_interp_eval_deriv_e( gsl_interp * obj,
                           double xa[],  double ya[], double x,
                          gsl_interp_accel * a,
                          double * d)

  double gsl_interp_eval_deriv( gsl_interp * obj,
                         double xa[],  double ya[], double x,
                        gsl_interp_accel * a)

  int gsl_interp_eval_deriv2_e( gsl_interp * obj,
                            double xa[],  double ya[], double x,
                           gsl_interp_accel * a,
                           double * d2)

  double gsl_interp_eval_deriv2( gsl_interp * obj,
                          double xa[],  double ya[], double x,
                         gsl_interp_accel * a)

  int gsl_interp_eval_integ_e( gsl_interp * obj,
                           double xa[],  double ya[],
                          double a, double b,
                          gsl_interp_accel * acc,
                          double * result)

  double gsl_interp_eval_integ( gsl_interp * obj,
                         double xa[],  double ya[],
                        double a, double b,
                        gsl_interp_accel * acc)

  void gsl_interp_free(gsl_interp * interp)

  size_t gsl_interp_bsearch( double x_array[], double x,
                            size_t index_lo, size_t index_hi)



cdef extern from "gsl/gsl_spline.h":
  ctypedef struct gsl_spline

  gsl_spline * gsl_spline_alloc( gsl_interp_type * T, size_t size)

  int gsl_spline_init(gsl_spline * spline,  double xa[],  double ya[], size_t size)


  int gsl_spline_eval_e( gsl_spline * spline, double x,
                    gsl_interp_accel * a, double * y)

  double gsl_spline_eval( gsl_spline * spline, double x, gsl_interp_accel * a)

  int gsl_spline_eval_deriv_e( gsl_spline * spline, double x,
                          gsl_interp_accel * a, double * y)

  double gsl_spline_eval_deriv( gsl_spline * spline, double x, gsl_interp_accel * a)

  int gsl_spline_eval_deriv2_e( gsl_spline * spline, double x,
                           gsl_interp_accel * a, double * y)

  double gsl_spline_eval_deriv2( gsl_spline * spline, double x,
                         gsl_interp_accel * a)

  int gsl_spline_eval_integ_e( gsl_spline * spline, double a, double b,
                          gsl_interp_accel * acc, double * y)

  double gsl_spline_eval_integ( gsl_spline * spline, double a, double b,
                        gsl_interp_accel * acc)

  void gsl_spline_free(gsl_spline * spline)

