GslMultiminFunctionAdapter.hpp

Go to the documentation of this file.

#ifndef JULIAN_GSLMULTIMINFUNCTIONADAPTER_HPP
#define JULIAN_GSLMULTIMINFUNCTIONADAPTER_HPP

#include <gsl/gsl_multimin.h>

namespace julian {

  template< typename F >
  class GslMultiminFunctionAdapter : public gsl_multimin_function {
  public:
    GslMultiminFunctionAdapter(const F& func,unsigned int n_) : func_(func) {
      f = &GslMultiminFunctionAdapter::call;
      n = n_;
      params=this;

    }
    private:
    static double call(const gsl_vector* x, void *params) {
      std::vector<double> x_std;

      for(unsigned int i = 0; i < static_cast<GslMultiminFunctionAdapter*>(params)->n; ++i) {
    x_std.push_back( gsl_vector_get(x,i));
      }
      return static_cast<GslMultiminFunctionAdapter*>(params)->func_(x_std);
    }

    const F& func_; 
  };
}
#endif /* GSLMULTIMINFUNCTIONADAPTER_HPP */