4

C++ 中有没有一种方法可以有效地创建一个作为函数指针的闭包？我正在使用 Gnu 科学图书馆，我必须创建一个gsl_function。这个函数需要在我创建它时有效地“关闭”几个可用的参数。是否有一个很好的技巧来创建一个闭包，这样我就不必将所有这些作为参数传递给 gsl_function 结构？如果不是，我应该只传递一个指向包含这些参数的数组的指针吗？

编辑 我曾尝试像这样使用 boost::bind：

#include <gsl/gsl_integration.h>
#include <boost/bind.hpp>

#include "bondpricecalculator.h"
#include "functions.h"

double integrand (double xi, double t, double x, void * p) {
        Functions *functions = (Functions *) p;
        double vx = functions->v(x);
        return functions->rho0(x)*exp(vx * xi - 0.5 * vx * vx * t);
     }

double BondPriceCalculator::value(double t, double T, double xi)
{
    gsl_integration_workspace * w
         = gsl_integration_workspace_alloc (10000);

    gsl_function F;

    F.function = &boost::bind(integrand, xi, t, _1, _2);
    F.params = &functions;

    double integral_t;
    double integral_T;
    double error;

    int res = gsl_integration_qags(&F, T, 1e+14, 0, 1e-7, 10000, w, &integral_T, &error);
    if(res)
    {
        throw "Error intgrating";
    }

    int res = gsl_integration_qags(&F, T, 1e+14, 0, 1e-7, 10000, w, &integral_t, &error);
    if(res)
    {
        throw "Error intgrating";
    }

    return integral_T/integral_t;
}

但我收到以下错误消息：

/home/ga/svn/PhD/inflation/cpp/ioi/bondpricecalculator.cpp:20: error: cannot convert ‘boost::_bi::bind_t<double, double (*)(double, double, double, void*), boost::_bi::list4<boost::_bi::value<double>, boost::_bi::value<double>, boost::arg<1>, boost::arg<2> > >*’ to ‘double (*)(double, void*)’ in assignment
4

4 回答 4

2

看看这个结合 boost::bind 和 boost::function 的简单例子。

于 2010-07-07T11:15:23.957 回答
2

我在下面找到了代码。

http://bytes.com/topic/c/answers/657124-interface-problem

// Use in combination with boost::bind.
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;
}

并像这样使用

gslFunction gslF = convertToGslFunction( boost::bind( &Sde::drift, &sde, _1 ) );
于 2010-07-07T12:48:52.213 回答
2

我从所有这些“gsl_”前缀中猜测该库不是 C++，而是普通 C。这意味着它不会理解 C++ 闭包（函子）。您不能将 C++ 仿函数传递给 C 函数。您必须传递 void 指针，交叉手指并将它们重新解释为 C 遗忘。

于 2010-07-07T16:11:40.747 回答
1

尽管 bradgonesurfing 给出了一个很好的答案，可以将闭包转换为gsl_functions 而无需进一步考虑，但我想与您分享从 C++ 直接翻译成 C 的习语。

假设你有闭包：

double a;
[&a](double x){return a+x;}

您可以将其转换为等效的函数指针习语，如下所示：

struct paramsAPlusX{
    double* a;
    paramsAPlusX(double & a_):a(&a_){}
}
double funcAPlusX(double x, void* params){
   paramsAPlusX* p= (paramsAPlusX*)params;
   return *(p->a) + x;
}

//calling code:
double a;
paramsAPlusX params(a);
gsl_function f;
f.function=funcAPlusX;
f.params=&paramsAPlusX;
//use f here.

许多 C 库都使用这种习惯用法，并且它们并不都使用结构体（它们经常将其作为两个单独的参数传递给函数），因此并不总是可以进行自动转换。

于 2010-07-07T22:23:43.043 回答