其他类的成员函数对成员函数的嵌套引用

我正在尝试构建一个由类de_solver编码的通用随机微分方程求解器，该类需要一组由model类给出的微分方程。该模型通过类sde提供给求解器，该类在模型和求解器之间工作接口，将原始模型的方程重新转换为确定性 + 随机方程的两个子集。

目前我的代码有几个问题，我认为必须处理继承和类型转换。特别是，我无法将指定模型方程的成员函数传递给积分器。首先假设我有一个这种类型的model类：

class model{
...
public:
size_t N_eq;
model(...); // Assign N_eq and parameters and other stuff
...
int equations_det(double t, const double y_[], double dy_[]); // Define deterministic RHS of equations
static int equations_gsl_wrapper(double t, const double y_[], double dy_[],void *params); // A wrapper that make equations_det suitable for a GSL solver
void equations_stoch(double t,double dt,const double y_[], double dy_[],); // The stochastic RHS of the model equations
static void ode_s_wrapper(double t,const double y_[], double dy_[], void *params); // A wrapper to pass equations_det to the stochastic integrator (in the class `de_solver`).
int simulate(t,tfin,tstep); // The actual simulator which will invoke the 'de_solver'
};

正如本文所指出的，static规范源于对模型的确定性部分使用 GSL 集成器的需要。

然后，接口类sde如下：

class sde{
...
public:
size_t N_eq;
sde(size_t N_,
int (*dtr)(double, const double*, double*, void*),
void (*stc)(double, double, const double*, double*, void*));
int (*deterministic)(double t, const double* y_[], double* dy_[], void * params);
void (*stochastic)(double t,const double y_[], double dy_[], void *params);
//Then again, akin to the model class, use some `wrappers`
static int deterministic_wrapper(double t, const double y_[], double dy_[], void * params);
static void stochastic_wrapper(double t, const double y_[], double *dy_, void *params);
};

这个想法是在sde类中具有从给定的模型继承的成员函数。至于两个"包装纸"，我将澄清我稍后介绍它们的原因。

最后，de_solver类是这样的：

class de_solver{
sde *sys;
public:
de_solver(sde *system); // Will initialize solver with the system put in the `sde` form
...
void integrate(void *params, double *ts, double **sol);
};  

model在单独的文件(model.H和model.CPP(中声明/定义，涉及sde和de_solver(solvers.h和solvers.cpp(。

这个想法是在model类中有一个simulate成员函数，以便

int model::simulate(double t, double tfin, double dt){
// Prepare solver
// 1. Create the `sde` object from model `sys`
sde recast_sys(NEQ, model::deterministic_wrapper, model::stochastic_wrapper);
// 2. Instantiate solver with recast system
de_solver integrator(&recast_sys);
// Run simulation
double *ts = ...   // Output time instants
double **sol = ... // Output solution
void *params_base = static_cast<void*>(std::addressof(this)); 
integrator.integrate(params_base);
return 1;  // In practice there is some error check on this return condition (omitted here for brevity)
}

总之，sys调用了通过recast_sys提供的模型方程的确定性和随机性部分的integrator。由于积分器的确定性部分依赖于 GSL 求解器，因此我使用附加参数参数将指向求解器的指针传递给实际类成员函数。以这种方式，在我拥有的integrator.integrate成员函数中(见上面提到的帖子(

de_solver::integrate(void *params_base, ...){
...
// I allocate an array of two void pointers: the first to the `model` class (assumed to be passed by `params_base`), and the second to the `sde` class
void **params = (void**)calloc(2,sizeof(void*));
params[0] = params_base;
params[1] = reinterpret_cast<void *>(std::addressof(sys)); // the recast system as private member of the sde class
gsl_odeiv2_driver * d;
gsl_odeiv2_system system = {sys->deterministic_wrapper, nullptr, sys->NEQ, params};
d = gsl_odeiv2_driver_alloc_y_new (&system, gsl_odeiv2_step_bsimp, opts.dt, opts.atol, opts.rtol);
...
}
int sde::deterministic_wrapper(double t, const double y_[], double dy_[], void * params){
assert(params);
return(static_cast<sde*>(params[1])->deterministic(t,y_,dy_,params)); // This will issue an error: ‘void*’ is not a pointer-to-object type
}
int model::equations_gsl_wrapper(double t, const double y_[], double dy_[], void * params){
assert(params);
return(static_cast<model*>(params[0])->ode_gsl(t,y_,dy_)); // This will issue an error: ‘void*’ is not a pointer-to-object type
}

从这篇文章中获取了两个指向空隙的指针数组。但是，似乎一旦在包装器中使用，它就会产生错误，可能是因为 void 数组上的算术不清楚(如此处指出的(？

目前，我无法为上面报告的错误编译代码。同样出于某种原因，编译器告诉我model.simulate成员函数中的this指针生成

error: use of deleted function ‘const _Tp* std::addressof(const _Tp&&) [with _Tp = model*]’

我怀疑我弄乱了静态和非静态成员函数，并且我没有正确传递它们。任何意见将不胜感激。