Interpolator.h

#ifndef CARMA_SERVICES_INTERPOLATOR_H
#define CARMA_SERVICES_INTERPOLATOR_H

#include "carma/services/Types.h"
#include <vector>
#include <cstddef>

extern "C" {
#include <gsl/gsl_errno.h>
#include <gsl/gsl_interp.h>
#include <gsl/gsl_spline.h>
}


namespace carma {
  namespace services {

  class Interpolator {
   public:

      Interpolator ( 
          const ::std::vector<double> & xvalue,
          const ::std::vector<double> & yvalue,
          const interpolationType type );

      virtual ~Interpolator ( );

      double evaluate( double x );


   private:
      void freeAllocations();
      void freeGslAllocations();
      void freeDoubleAllocations();
      void initialize();
      void setInterpType( interpolationType type );


      ::gsl_interp_accel  * acc_;
      ::gsl_interp        * interp_;
      ::gsl_spline        * spline_;
      // eventually do away with these if not allowing
      // interp type to be changed after instantiation.
      ::std::vector<double> xval_;
      ::std::vector<double> yval_;
      size_t size_;
      double * x_;
      double * y_;
      double answer_;
      interpolationType     interpType_;
  }; // class Interpolator

} // namespace services
} // namespace carma


#endif //CARMA_SERVICES_INTERPOLATOR_H