
![]() | The wrapper code must allocate extra data:
double **app; double *ap;
ap = (double *) a->data; /* a is a PyArrayObject* pointer */
app = (double **) malloc(nx*sizeof(double*));
for (i = 0; i < nx; i++) {
app[i] = &(ap[i*ny]); /* point row no. i in a->data */
}
/* clean up when app is no longer needed: */ free(app);
|