1

我想编写一个包含用于求解 ODE 系统的所有函数（包括 GSL 函数）和参数的结构。从主函数开始，我只想调用结构中定义的更新函数，以将系统推进一个时间步长。但是，当我尝试此操作时，出现错误：

Line 27, ERROR:  cannot convert ‘ODE::funcS’ from type ‘int (ODE::)(double, const double*, double*, void*)’ to type ‘int (*)(double, const double*, double*, void*)’ Below is a minimal code. \

这是我的代码的最小版本：

#include <iostream>
#include <gsl/gsl_errno.h>
#include <gsl/gsl_matrix.h>
#include <gsl/gsl_odeiv.h>

struct ODE
{
    void update(double dt)
    {
        // code to advance ODE solution by one time-step dt
    }

    int
    funcS (double t, const double y[], double f[],
          void *params)
    {
        return GSL_SUCCESS;
    }

    double mu = 10;

    gsl_odeiv_system sysS;

    void
    initializeSys()
    {
        sysS.function = funcS; //Line 27
    }
};

int
func (double t, const double y[], double f[],
          void *params)
{
    return GSL_SUCCESS;
}

int main()
{
    // GIVES ERROR
    ODE mySys;
    mySys.update(0.01);


    // WORKS
    double mu = 10;
    gsl_odeiv_system sys;
    sys.function = func;

    return 0;
}
4

1 回答 1

2

您不需要直接使用静态函数。相反，您可以编写一个非常通用的包装器。

我相信这是一个重复的问题。我对刚刚链接的问题的回答基于此处提供的包装器。但是，我使用模板对其进行了概括，以避免由于 std::function 持有的函子的堆分配而导致std::function 的性能损失（原始答案仅警告读者有关由涉及的多重间接引起的损失） std::function 实现，这与堆分配引起的问题相比可以忽略不计）。

编辑1：这里也讨论了这个问题

编辑 2（回答您在对我的回答的第一条评论中提出的问题）。第一个警告是，您必须确保std::function在 GSL 完成计算之前不会删除任何持有的内容。此外，@Managu 指出，当 GSL 工作时，包装器本身不能超出范围。如果您仔细编码，这并不难执行。错误代码示例：

 // BAD PROGRAM - EXAMPLE OF WHAT YOU MUST NOT DO. DO NOT COPY THIS CODE
 // HERE THE WRAPPER GETS PREMATURELY OUT OF SCOPE => CRASH
 gsl_function *F 
 auto ptr = [](double x)->double{return 2*x;};
 std::function<double(double)> FF1(std::cref(ptr))

 {      
   gsl_function_pp Fp(FF1);
   F = static_cast<gsl_function*>(&Fp);  
 }
 (...)
 // CALL GSL
于 2013-10-04T06:31:47.507 回答