GslFunctionFdfAdapter.hpp

Go to the documentation of this file.

#ifndef JULIAN_GSLFUNCTIONFDFADAPTER_HPP
#define JULIAN_GSLFUNCTIONFDFADAPTER_HPP

#include <gsl/gsl_math.h>


namespace julian {
 
  template< typename F, typename dF >
  class GslFunctionFdfAdapter : public gsl_function_fdf {
  public:

    GslFunctionFdfAdapter(const F& kf, const dF& kdf) : f_(kf), df_(kdf) {
      f   = &GslFunctionFdfAdapter::call_f;
      df  = &GslFunctionFdfAdapter::call_df;
      fdf = &GslFunctionFdfAdapter::call_fdf;
      params=this;
    }
  private:
    static double call_f(double x, void *params) {
      return static_cast<GslFunctionFdfAdapter*>(params)->f_(x);
    }
    
    static double call_df(double x, void *params) {
      return static_cast<GslFunctionFdfAdapter*>(params)->df_(x);
    }
    
    static void call_fdf(double x, void * params, double* f, double* df) {
      (*f)  = static_cast<GslFunctionFdfAdapter*>(params)->f_(x);
      (*df) = static_cast<GslFunctionFdfAdapter*>(params)->df_(x);
    }

    const F& f_;   
    const dF& df_; 
  }; 
}
#endif /* GSLFDFFUNCTIONADAPTER_HPP */