在只有numpy的cython中调用此函数的最佳方法是什么?我不打算使用ctype,memcpy,malloc等等。
职能1)
#include <stdio.h>
extern "C" void cfun(const void * indatav, int rowcount, int colcount,
void * outdatav);
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] * 4;
}
puts("Done!");
}
职能2)
#include <stdio.h>
extern "C" __declspec(dllexport) void cfun(const double ** indata, int
rowcount, int colcount, double ** outdata) {
for (int i = 0; i < rowcount; ++i) {
for (int j = 0; j < colcount; ++j) {
outdata[i][j] = indata[i][j] * 4;
}
}
}
旺君湾
发布于 2012-01-24 18:51:09
您可以直接从Cython‘调用’函数,方法是将它声明为extern。
cdef extern from "mylibraryheader.h":
void cfun1(void* indatav, int rowcount, int colcount, void* outdatav)
void cfun2(double** indata, int rowcount, int colcount, doubke** outdata)
您现在可以像在C/C++中那样调用这些函数。注意,Cython中没有const关键字,可以省略它。不幸的是,我不能给您一个示例,说明如何将NumPy数组转换为double数组。但是,这里有一个从双倍列表中运行它的例子。
cdef extern from "mylibraryheader.h":
void cfun1(void* indatav, int rowcount, int colcount, void* outdatav)
void cfun2(double** indata, int rowcount, int colcount, double** outdata)
cdef extern from "stdlib.h":
ctypedef int size_t
void* malloc(size_t)
void free(void*)
def py_cfunc1(*values):
cdef int i = 0
cdef int size = sizeof(double)*len(values)
cdef double* indatav = <double*> malloc(size)
cdef double* outdatav = <double*> malloc(size)
cdef list outvalues = []
for v in values:
indatav[i] = <double>v
i += 1
cfun1(<void*>indatav, 1, len(values), <void*>outdatav)
for 0 <= i < len(values):
outvalues.append(outdatav[i])
return outvalues
注:未经测试
发布于 2012-01-10 15:51:44
您不能从cython调用这个函数--您必须将它写成一个cython函数--一些好的例子是这里。作为参考,职能1):
cimport numpy as np
def cfun(np.ndarray indata, int rowcount, int colcount, np.ndarray outdata):
cdef int i
print("Here we go!")
for i in range(rowcount * colcount):
outdata[i] = indata[i] * 4
print("Done!")
如果您真的想调用它,您必须使用ctype或编写您自己的包装器。使用大口也是可行的。
https://stackoverflow.com/questions/8797335
复制相似问题