在 C 中通过 void 类型为 GSL 传递双精度数组

Question

我正在尝试使用 GSL 库来解决 ODE，但使用 void 指针时遇到一些困难

我需要在应该包含数组的数组上发送一个参数：

double k1[2][4];

发送到此

gsl_odeiv_system sys = {func, jac, 2, &k1};

的参数会传递给两者func 和 jac 作为

int func (double t, const double y[], double f[], void *params)

func 中的 *params，我试图通过以下方式提取 k1：

double k1[2][4];
k1 = *(double[][])params;

或

k1 = (double[][])params;

或...

k1 = *(double *)params;

等等

我想问题是: 有没有单行解决方案？

Answer 1

我不认为你可以像这样转换为数组类型（多维数组）。您可能需要声明一个临时变量来保存指向数组第一个元素的指针。

当然，您需要指定每行的元素数量才能正常工作。否则，编译器不知道如何访问结果数组中的元素（请记住，x[i][j] 在内部转换为*(x + i*n + j) code> 其中 n 是每行中的元素数量）。

即顺便说

 double x[5][2];

 int main()
 {
     double (*y)[5][2];

     void *z = x;

     y = z;


     /* after you extract the pointer from 'z' you 
        can access the elements using (*y)[][] */
     (*y)[1][1] = 1.0;
 }

一句，将数组传递给函数时不必使用 &k1 。数组的名称可以用作其地址（指向第一个元素的指针）。