Question

我有一个带有以下签名的c函数

__declspec( dllexport ) void* setup(int c_force, int c_stepping, int c_iteration, int c_roots,
                                    struct para* c_userdata, double* c_y0,
                                    double c_reltol, double c_abstol)
{
    ....
    return "a pointer";
}

其中para定义如下。我想用这个函数（和其他人）构建一个DLL并从Python访问它。重要的Python行是

import ctypes as ct
lib = ct.cdll.LoadLibrary('lib_path.dll')
getattr(lib, 'setup')
lib.setup.restype = ct.c_void_p
# pystruct as defined below to avoid clogging code
lib.setup.argtypes = [ct.c_int, ct.c_int, ct.c_int, ct.c_int,
                      ct.POINTER(para), ct.POINTER(ct.c_double),
                      ct.c_double, ct.c_double]
ptr = lib.setup(cf, cs, ci, cr, ct.byref(para), c_y0, c_reltol, c_abstol) 
myobj = ct.c_void_p(ptr)

其中cf，cs，ci，cr是（python）整数，para是结构类型，定义如下，{{1 （应该是长度为2，来自python端和c端）并且c_y0=(ct_c_double * 2)()和c_reltol被转换为ct.c_double为c_abstol。

当我尝试运行我的主应用程序时，我在c_reltol = ct.c_double(reltol) - 函数调用时得到WindowsError: exception: access violation writing xxx，我无法理解为什么...在调用lib.setup之前打印给出以下输出传递给函数的参数的值和类型（按顺序）

lib.setup

我一直在尝试使用this和this问题进行调试，但没有成功。由于c函数的调用签名非常简单，我无法理解为什么它会破坏。

P.S。它在Ubuntu上完美运行，1 2 2 1 <cparam 'P' (0000000003C84EB0)> <cvode_library.c_double_Array_2 object at 0x0000000003F95BC8> c_double(1e-06) c_double(1e-08) <type 'int'> <type 'int'> <type 'int'> <type 'int'> <type 'CArgObject'> <class 'cvode_library.c_double_Array_2'> <class 'ctypes.c_double'> <class 'ctypes.c_double'>替换为__declspec...

c-struct定义为

extern

和相应的pystruct为

typedef struct para PARA;
struct para
{
    double a;
    double b;
    double c;
};

编辑 class para(ct.Structure): _fields_ = [('a', ct.c_double), ('b', ct.c_double), ('c', ct.c_double)]定义为

c_y0

所有＆＃34; cvode＆＃34; -functions和y0 = np.array([0., 0.]) c_y0 = (ct.c_double * 2)() c_y0[0] = y0[0] c_y0[1] = y0[1]是Sundials suite for solving nonlinear equations的一部分

N_Vector

Answer 1

这是一个MCVE。它显示您的声明是正确的，因此问题很可能在函数实现中。如果以下内容对您不起作用，请使用类似的MCVE更新您的问题，以便重现您的失败。

<强> test.c的

#include <stdio.h>

typedef struct para PARA;
struct para
{
    double a;
    double b;
    double c;
};

__declspec(dllexport) void* setup(int c_force, int c_stepping, int c_iteration, int c_roots,
                                  struct para* c_userdata, double* c_y0,
                                  double c_reltol, double c_abstol)
{
    printf("%d %d %d %d %lf %lf %lf %lf %lf %lf %lf\n",c_force,c_stepping,c_iteration,c_roots,c_userdata->a,c_userdata->b,c_userda
ta->c,c_y0[0],c_y0[1],c_reltol,c_abstol);
    return NULL;
}

<强> test.py

import ctypes as ct

class para(ct.Structure):
    _fields_ = [('a', ct.c_double),
               ('b', ct.c_double),
               ('c', ct.c_double)]

lib = ct.CDLL('test')
lib.setup.restype = ct.c_void_p
# pystruct as defined below to avoid clogging code
lib.setup.argtypes = [ct.c_int, ct.c_int, ct.c_int, ct.c_int,
                      ct.POINTER(para), ct.POINTER(ct.c_double),
                      ct.c_double, ct.c_double]
p = para(1.5,2.5,3.5)
c_y0 = (ct.c_double * 2)(4.5,5.5)
ptr = lib.setup(1,2,3,4,p, c_y0,6.5,7.5) 
myobj = ct.c_void_p(ptr)

<强>输出

1 2 3 4 1.500000 2.500000 3.500000 4.500000 5.500000 6.500000 7.500000

使用ctypes在Python中包装c函数

1 个答案: