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

Question

我正在尝试使用 Python 从第 3 方 Fortran DLL 调用函数（冒泡排序）。我的问题是将一个字符传递给函数。我让它与 cffi 一起工作，如下所示，但我想使用 ctypes。

cffi版本：

import numpy as np

import cffi as cf

ffi=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 = 1

m1 = ffi.new("int*", m1)

m2 = 16

m2 = 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 = 0

ifail = 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 np

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 = '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 也是如此。## CTYPESimport 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 = 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)