gslpp_function_adapter.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2017 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef GSLPP_FUNCTION_ADAPTER_H
#define GSLPP_FUNCTION_ADAPTER_H
 
#include <gsl/gsl_math.h>
#include <assert.h>
 
template<class F>
static double gslFunctionAdapter( double x, void* p)
{
    // Here I do recover the "right" pointer, safer to use static_cast
    // than reinterpret_cast.
    F* function = static_cast<F*>( p );
    return (*function)( x );
}
 
template<class F>
gsl_function convertToGslFunction( const F& f )
{
    gsl_function gslFunction;
    
    const void* p = &f;
    assert (p != 0);
    
    gslFunction.function = &gslFunctionAdapter<F>;
    // Just to eliminate the const.
    gslFunction.params = const_cast<void*>( p );
    
    return gslFunction;
}
 
#endif /* GSLPP_FUNCTION_ADAPTER_H */