GslFunctionAdapter.hpp

Go to the documentation of this file.

#ifndef JULIAN_GSLFUNCTIONADAPTER_HPP
#define JULIAN_GSLFUNCTIONADAPTER_HPP

#include <gsl/gsl_math.h>


namespace julian {
  template< typename F >
  class GslFunctionAdapter : public gsl_function {
  public:

    GslFunctionAdapter(const F& functor) : functor_(functor) {
      function = &GslFunctionAdapter::call;
      params=this;
    }
  private:
    static double call(double x, void *params) {
      return static_cast<GslFunctionAdapter*>(params)->functor_(x);
    }

    const F& functor_; 
  };
}
#endif