我正在使用 GSL 示例页面上的代码来尝试求解一个由四个微分方程组成的系统。我一直在努力将五个参数传递给 ODE 系统,并且遇到了一个最终的(我希望如此!)编译时错误。下面是一个片段,给了我错误
114:57: error: invalid conversion from ‘int (*)(double, const double*,
double*, double**, void*)’ to ‘int (*)(double, const double*, double*,
double*, void*)’ [-fpermissive]
这对应于以 开头的行gsl_odeiv2_system
:
int main()
{
double t = 0.0;
double y[4] = { 0.0, 0.0, 1.0, 1.0 };
int i, s;
struct pendula_params * info;
info->m2 = 1.0;
info->m1 = 1.0;
info->l1 = 1.0;
info->l2 = 1.0;
info->g = 1.0;
gsl_odeiv2_system sys = { pendula, jacobian, 4, &info };
gsl_odeiv2_driver *d =
gsl_odeiv2_driver_alloc_y_new (&sys, gsl_odeiv2_step_msadams,
1e-3, 1e-8, 1e-8);
对可能发生的事情有任何想法吗?
非常感谢,
马克 C。