GSL ODE解决方案中的指针类型编译时错误

Question

我正在使用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.

Answer 1

您可以通过比较错误消息中出现的两种函数类型来查看错误是什么。

gsl_odeiv2_system结构期望jacobian成员是指向以double*作为第四个参数的函数的指针。 但是您的jacobian函数将double**作为第四个参数，这使其不兼容。