哦,我的话我是个傻瓜。 调用函数时,我只是省略了第二个和第三个参数。像个傻瓜。因为那就是我。原来愚蠢的问题如下:
这似乎是它必须做一个很平常的事,但是我无法找到相关的教程,我太无知了约Numpy,并ctypes弄明白自己。
Numpy
ctypes
我的文件中有一个C函数ctest.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数组未修改。实际上,如果我再次调用该函数,则会遇到段错误。这并不令我感到惊讶-我真的不知道我在这里做什么。谁能指出我正确的方向?
只需将所有四个参数传递给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))