使用 Char 作为参数从 Fortran DLL 调用函数

Question

我正在尝试使用 Python 从第 3 方 Fortran DLL 调用函数（冒泡排序）。我的问题是将一个字符传递给函数。我让它与 cffi 一起工作，如下所示，但我想使用 ctypes。cffi版本：import numpy as npimport cffi as cfffi=cf.FFI()lib=ffi.dlopen(r"C:\Windows\SysWOW64\DLL20DDS")ffi.cdef("""void M01CAF( double rv[], const int *m1, const int *m2,         const wchar_t *order, int *ifail);""")m1 = 1m1 = ffi.new("int*", m1)m2 = 16m2 = ffi.new("int*", m2)order = ffi.new('wchar_t *', "A")rvx = np.array([1.3, 5.9, 4.1, 2.3, 0.5, 5.8, 1.3, 6.5,                 2.3, 0.5, 6.5, 9.9, 2.1, 1.1, 1.2, 8.6], dtype=float, order='F')rv = ffi.cast("double* ", rvx.__array_interface__['data'][0])ifail = 0ifail = ffi.new('int*', ifail)lib.M01CAF(rv, m1, m2, order, ifail)print(rvx)输出：[0.5 0.5 1.1 1.2 1.3 1.3 2.1 2.3 2.3 4.1 5.8 5.9 6.5 6.5 8.6 9.9]现在我的 ctypes 版本：import numpy as npimport ctypes as ctflib = ct.WinDLL('C:\Windows\SysWOW64\DLL20DDS.dll')func = getattr(flib, "M01CAF")func.restype = Nonem1 = (ct.c_int32)(1)m2 = (ct.c_int32)(16)ifail = ct.c_int32(0)rv = np.array([1.3, 5.9, 4.1, 2.3, 0.5, 5.8, 1.3, 6.5,                 2.3, 0.5, 6.5, 9.9, 2.1, 1.1, 1.2, 8.6], dtype=ct.c_double, order='F')order = 'A'order = ct.c_wchar_p(order)func.argtypes = (np.ctypeslib.ndpointer(dtype=ct.c_double, shape=(m2.value,)),                 ct.POINTER(ct.c_int32), ct.POINTER(ct.c_int32),                 ct.c_wchar_p, ct.POINTER(ct.c_int32))func(rv, m1, m2, order, ifail)错误信息：OSError: exception: access violation writing 0xFFFFFFFC______________________________________________________cffi 遇到的一个问题是：当我再次调用该函数时，它不会对数组进行排序。我需要重新启动内核才能获得正确的结果：rvx = np.array([1.3, 5.9, 4.1, 2.3, 0.5, 5.8, 1.3, 6.5,                 2.3, 0.5, 6.5, 9.9, 2.1, 1.1, 1.2, 8.6], dtype=float, order='F')rv = ffi.cast("double* ", rvx.__array_interface__['data'][0])lib.M01CAF(rv, m1, m2, order, ifail)print(rvx)输出：[1.3 5.9 4.1 2.3 0.5 5.8 1.3 6.5 2.3 0.5 6.5 9.9 2.1 1.1 1.2 8.6]______________________________________________________现在我在 Win10 64 位上使用 Spyder 和 Python 3.8.3 32 位。

Answer 1

通过将 char 的长度作为常规整数（而不是指针）传递到 char 本身的指针之后，解决了所有问题。CFFI 也是如此。

## CTYPES

import ctypes as ct

flib = ct.WinDLL('C:\Windows\SysWOW64\DLL20DDS.dll')

func = getattr(flib, "M01CAF")

func.restype = None

m1 = (ct.c_int32)(1)

m2 = (ct.c_int32)(16)

ifail = ct.c_int32(0)

rv = np.array([1.3, 5.9, 4.1, 2.3, 0.5, 5.8, 1.3, 6.5, 

                2.3, 0.5, 6.5, 9.9, 2.1, 1.1, 1.2, 8.6], dtype=ct.c_double, order='F')

order = ct.c_wchar_p('A')

func.argtypes = (np.ctypeslib.ndpointer(dtype=ct.c_double, shape=(m2.value,)), 

                ct.POINTER(ct.c_int32), ct.POINTER(ct.c_int32), 

                ct.c_wchar_p, ct.c_int32, ct.POINTER(ct.c_int32))

print("BEFORE:")

print(rv)

func(rv, m1, m2, order, 1, ifail)

print("AFTER:")

print(rv)