将Numpy数组传递给C函数进行输入和输出

Question

哦，我的话我是个傻瓜。 调用函数时，我只是省略了第二个和第三个参数。像个傻瓜。因为那就是我。原来愚蠢的问题如下：

这似乎是它必须做一个很平常的事，但是我无法找到相关的教程，我太无知了约Numpy，并ctypes弄明白自己。

我的文件中有一个C函数ctest.c。

#include <stdio.h>

void cfun(const void * indatav, int rowcount, int colcount, void * outdatav) {
    //void cfun(const double * indata, int rowcount, int colcount, double * outdata) {
    const double * indata = (double *) indatav;
    double * outdata = (double *) outdatav;
    int i;
    puts("Here we go!");
    for (i = 0; i < rowcount * colcount; ++i) {
        outdata[i] = indata[i] * 2;
    }
    puts("Done!");
}

（您可能会猜到，我最初的参数为double 而不是void ，但无法弄清楚在Python方面该怎么做。我当然很想将其改回来，但我并不挑剔只要有效。）

我用它制作了一个共享库。gcc -fPIC -shared -o ctest.so ctest.c

然后在Python中，我有几个numpy数组，我想将它们传递给C函数，一个作为输入，另一个作为输出。

indata = numpy.ones((5,6), dtype=numpy.double)
outdata = numpy.zeros((5,6), dtype=numpy.double)
lib = ctypes.cdll.LoadLibrary('./ctest.so')
fun = lib.cfun
# Here comes the fool part.
fun(ctypes.c_void_p(indata.ctypes.data), ctypes.c_void_p(outdata.ctypes.data))

print 'indata: %s' % indata
print 'outdata: %s' % outdata

这不会报告任何错误，但是会打印出来

>>> Here we go!
Done!
indata: [[ 1.  1.  1.  1.  1.  1.]
 [ 1.  1.  1.  1.  1.  1.]
 [ 1.  1.  1.  1.  1.  1.]
 [ 1.  1.  1.  1.  1.  1.]
 [ 1.  1.  1.  1.  1.  1.]]
outdata: [[ 0.  0.  0.  0.  0.  0.]
 [ 0.  0.  0.  0.  0.  0.]
 [ 0.  0.  0.  0.  0.  0.]
 [ 0.  0.  0.  0.  0.  0.]
 [ 0.  0.  0.  0.  0.  0.]]

outdata数组未修改。实际上，如果我再次调用该函数，则会遇到段错误。这并不令我感到惊讶-我真的不知道我在这里做什么。谁能指出我正确的方向？

Answer 1

只需将所有四个参数传递给C函数。从以下位置更改Python代码：

fun(ctypes.c_void_p(indata.ctypes.data), ctypes.c_void_p(outdata.ctypes.data))

至：

fun(ctypes.c_void_p(indata.ctypes.data), ctypes.c_int(5), ctypes.c_int(6),
    ctypes.c_void_p(outdata.ctypes.data))