如何使用ctypes将矩阵从C ++函数返回到Python

Question

我想从C ++函数返回一个矩阵到Python函数。 我已经检查了这个解决方案，这是一个返回数组的例子。

作为一个例子，我想返回一个10x10填充的阵列10的。

function.cpp：

extern "C" int* function(){
int** information = new int*[10];
for(int k=0;k<10;k++) {
    information[k] = new int[10];
    }
for(int k=0;k<10;k++) {
    for(int l=0;l<10;l++) {
        information[k][l] = 10;
        }
}
return *information;
}

Python代码是： wrapper.py

import ctypes
from numpy.ctypeslib import ndpointer

lib = ctypes.CDLL('./library.so')
lib.function.restype = ndpointer(dtype=ctypes.c_int, shape=(10,))

res = lib.function()
print res

为了编译这个我使用：

g++ -c -fPIC function.cpp -o function.o
g++ -shared -Wl,-soname,library.so -o library.so function.o

如果soname不起作用，请使用install_name ：

g++ -c -fPIC function.cpp -o function.o
g++ -shared -Wl,-install_name,library.so -o library.so function.o

运行python程序之后， python wrapper.py这是输出我有：
[10 10 10 10 10 10 10 10 10 10]

只有一行10个元素。 我想要10x10矩阵。 我做错了什么？ 提前致谢。

Answer 1

function.cpp ：

extern "C" int* function(){
    int* result = new int[100];
    for(int k=0;k<100;k++) {
        result[k] = 10;
    }
    return result;
}

在wrapper.py

lib.function.restype = ndpointer(dtype=ctypes.c_int, shape=(10,)) //incorrect shape
lib.function.restype = ndpointer(dtype=ctypes.c_int, ndim=2, shape=(10,10)) // Should be two-dimensional