6

通过一些黑客攻击,我使用 python c-types 在 Python 和 C 之间来回传输数据。看起来有点乱,所以我希望有更多经验的人能告诉我我做错了什么,或者让这件事变得比需要的更困难。

蟒蛇代码:

from ctypes import cdll, POINTER, byref, c_double
lib = cdll.LoadLibrary('./test.so')

def alloc_arrays(length, python_list):
    # Create array of float*
    data = ((POINTER(c_double)) * length)()
    out_data = ((POINTER(c_double)) * length)()
    for i in xrange(length):
        # Create arrays of float
        data[i] = (c_double * 3)()
        out_data[i] = (c_double * 3)()
        # Set Values
        for j in xrange(3):
            data[i][j] = python_list[i][j]
            out_data[i][j] = 0.0
    return data, out_data

if __name__ == "__main__":
    a = [[1.0, 11.0, 21.0],
         [2.0, 12.0, 22.0],
         [3.0, 13.0, 23.0],
         [4.0, 14.0, 24.0],
         [5.0, 15.0, 25.0],
         [6.0, 16.0, 26.0],
         [7.0, 17.0, 27.0]]
    in_data, out_data = alloc_arrays(len(a), a)
    out_len = lib.smain(byref(in_data), len(a), out_data)

    print "out_len", out_len
    clean_out = [None, None, None]
    clean_out = [clean_out[:] for i in xrange(out_len)]
    for i in xrange(out_len):
        for j in xrange(3):
            clean_out[i][j] = out_data[i][j]
    print "Out:", clean_out

C代码(test.c):

int smain(double* in_points[3], int in_len, double* out_points[3]){
    int i;
    int j;
    printf("%s", "\n");
    for(i = 0; i < in_len; i++){
        for(j = 0; j < 3; j++){
            printf("%g, ", *in_points[i][j]);
        }
        printf("%s", "\n");
    }
    printf("%s", "\n");

    //*out_points = malloc(len*sizeof(float[3]));
    int out_len = in_len-2; //Randomly chosen shorter length
    for(i = 0; i < out_len; i++){
        for(j = 0; j < 3; j++){
            //Random function just to see I can do this
            *out_points[i][j] = i*j;
        }
    }
    return out_len;
}

要构建我使用:

gcc -c -fPIC test.c
gcc -shared -o test.so test.o

我想知道的主要事情是是否有更好的方法来创建空二维向量以传递给 C,以及是否有更好的方法将变量转换为/从 ctypes 转换为 python 类型。我不喜欢 alloc_arrays 函数必须进行转换,而且我不喜欢让我的 for 循环读取 c_types 2D 列表(刚刚意识到我这样做也有点乱)。

编辑:只是想清楚这段代码按预期工作,我只是认为必须有一些方法来改进它。

4

1 回答 1

0

这不是关于您的确切代码和 ctypes 的答案,但我想贡献...

我最近刚接触到cython。我什至还不是中级,但它似乎对编写扩展更加友好,并且有很多关于使用 numpy 数组的示例。

到目前为止,我发现 cython 非常棒,因为它可以让你保持类似于 python 的语法,并且可以自由地混合 python。在帮助将代码转换为包装器中的 C 方面,它为您处理了很多繁重的工作,这是在引擎盖下完成的。并且在produced中没有中间依赖.so。只是一个纯粹的扩展。

如果您要在 Python 代码中使用 numpy,那么将它们传递到您的扩展中将毫不费力,因为它不必进行任何转换。numpy 本身在 python 端已经是一个很大的加速,但是做额外的自定义东西可以让你将 numpy 指针传递到 cython 扩展端。

于 2012-08-13T16:56:32.413 回答