我想先声明 Fortran 不是我的技能之一-我接触过它~5次（SO 的其他问题）。
我发现了一堆问题：

• Fortran：
  *参数通过引用传递给函数/子例程（[GNU.GCC]: Argument passing conventions）.这花了我一个小时的时间来弄清楚.在我修复了所有的休息，我得到 * 访问冲突 * 甚至打印一个这样的参数.我在调查期间添加了虚拟函数.无论如何，最简单的版本是指定 value 作为按值传递的值（反之则是从调用代码传递指针）
• 由于我在 Win 上，我必须将 dllexport 属性添加到每个要由 .dll 导出的函数（我也可以通过将 /EXPORT 传递给每个函数的链接器来完成）。我还有另一个挫折，因为当检查函数导出时，我主要使用的工具（检查[SO]: Discover missing module using command-line ("DLL load failed" error) (@CristiFati's answer)）不会显示它：
• • 依赖关系 *：

• • 依赖关系步行者 *：

• 第一个参数是维度，不需要是双精度
• 其他（次要）问题（如重组，添加模块的东西，...）
• Python：
• 参数类型与 Fortran 定义中的参数类型不匹配，从而产生（如问题中的 URL 所述）Undefined Behavior。我说明了传递 NumPy 数组的两种方法（未注解的一种对我来说似乎更好一点）。

请注意，在使用多维数组时，Fortran 和 C 中的轴顺序不同（例如，对于 2D，前者是column major）
设法让示例工作（我不知道 sticksum 的所有操作都应该做什么）。

• dll00.f90：

module dll00

use iso_c_binding
implicit none

contains

#if 0
function dummy(i0, d0) result(res) bind(C, name="dummy")
    !dec$ attributes dllexport :: dummy
    integer(kind=c_int), intent(in), value :: i0
    real(kind=c_double), intent(in), value :: d0
    integer(kind=c_int) :: res

    res = 1618
    print *, "xxx", res
    print *, "yyy", i0
    print *, "zzz", d0
end function dummy
#endif

function sticksum(anz, w, b, a, weight, value, maxTchan, sthresh) result(spek) bind(C, name="sticksum")
    !dec$ attributes dllexport :: sticksum
    integer(kind=c_int), intent(in), value :: anz, maxTchan
    real(kind=c_double), intent(in), value :: b, a, sthresh
    real(kind=c_double), intent(in), dimension(*) :: w, weight, value
    real(kind=c_double) :: spek

    integer :: i, t, j
    real :: stick, maxw

    spek = 123
    i = 1
    t = w(i) * b + a + 1.5
    if (t >= 1) then
        spek = anz
        stick = 0. + value(t)
        maxw = weight(i) * stick
        do i = 2, anz
            t = w(i) * b + a + 1.5
            if (t > maxTchan) exit
            stick = value(t)
            maxw = max(weight(i) * stick, maxw)
            if ((w(i) * w(i) - w(i - 1) * w(i - 1)) > 0.5) then
                spek = spek + maxw
                maxw = 0
            end if
        end do
    end if
end function sticksum

end module dll00

• code00.py：

#!/usr/bin/env python

import ctypes as cts
import sys

import numpy as np

DLL_NAME = "./dll00.{:s}".format("dll" if sys.platform[:3].lower() == "win" else "so")

ELEMENT_TYPE = cts.c_double

DblPtr = cts.POINTER(cts.c_double)

def np_mat_type(dim, element_type=ELEMENT_TYPE):
    return np.ctypeslib.ndpointer(dtype=element_type, shape=(dim,), flags="F_CONTIGUOUS")

def main(*argv):
    dim = 5
    dll = cts.CDLL(DLL_NAME)
    sticksum = dll.sticksum
    sticksum.argtypes = (cts.c_int, np_mat_type(dim), cts.c_double, cts.c_double, np_mat_type(dim), np_mat_type(dim), cts.c_int, cts.c_double)
    #sticksum.argtypes = (cts.c_int, DblPtr, cts.c_double, cts.c_double, DblPtr, DblPtr, cts.c_int, cts.c_double)  # @TODO - cfati: alternative
    sticksum.restype = cts.c_double

    if 0:
        dummy = dll.dummy
        dummy.argtypes = (cts.c_int, cts.c_double)
        dummy.restype = cts.c_int
        print(dummy(cts.c_int(3141593), cts.c_double(2.718282)))

    w = np.arange(dim, dtype=ELEMENT_TYPE)
    a = 1
    b = 2
    weight = np.arange(1, dim + 1, dtype=ELEMENT_TYPE)
    val = np.arange(2, dim + 2, dtype=ELEMENT_TYPE)
    #print(w, weight, val)

    mtc = 3
    stresh = 6

    res = sticksum(dim, w, b, a, weight, val, mtc, stresh)
    #res = sticksum(dim, np.asfortranarray(w, dtype=ELEMENT_TYPE).ctypes.data_as(DblPtr), b, a, np.asfortranarray(weight, dtype=ELEMENT_TYPE).ctypes.data_as(DblPtr), np.asfortranarray(val, dtype=ELEMENT_TYPE).ctypes.data_as(DblPtr), mtc, stresh)  # @TODO - cfati: alternative
    print("\n{:s} returned: {:.3f}".format(sticksum.__name__, res))

if __name__ == "__main__":
    print("Python {:s} {:03d}bit on {:s}\n".format(" ".join(elem.strip() for elem in sys.version.split("\n")),
                                                   64 if sys.maxsize > 0x100000000 else 32, sys.platform))
    rc = main(*sys.argv[1:])
    print("\nDone.\n")
    sys.exit(rc)

输出：

[cfati@CFATI-5510-0:e:\Work\Dev\StackExchange\StackOverflow\q075937942]> sopr.bat
### Set shorter prompt to better fit when pasted in StackOverflow (or other) pages ###

[prompt]> "c:\Install\pc032\Microsoft\VisualStudioCommunity\2019\VC\Auxiliary\Build\vcvarsall.bat" x64 > nul

[prompt]>
[prompt]> dir /b
code00.py
dll00.f90

[prompt]>
[prompt]> "f:\Install\pc032\Intel\OneAPI\Version\compiler\2021.3.0\windows\bin\intel64\ifort.exe" /nologo /fpp /c dll00.f90

[prompt]>
[prompt]> link /NOLOGO /DLL /OUT:dll00.dll /LIBPATH:"f:\Install\pc032\Intel\OneAPI\Version\compiler\2021.3.0\windows\compiler\lib\intel64_win" dll00.obj
   Creating library dll00.lib and object dll00.exp

[prompt]> dir /b
code00.py
dll00.dll
dll00.exp
dll00.f90
dll00.lib
dll00.mod
dll00.obj

[prompt]>
[prompt]> "e:\Work\Dev\VEnvs\py_pc064_03.10_test0\Scripts\python.exe" ./code00.py
Python 3.10.9 (tags/v3.10.9:1dd9be6, Dec  6 2022, 20:01:21) [MSC v.1934 64 bit (AMD64)] 064bit on win32

sticksum returned: 5.000

Done.

读起来可能也很有趣：

• [SO]: What is numpy.ctypeslib.as_ctypes exacty doing (@CristiFati's answer)
• [SO]: How to convert a ctypes C pointer to a Numpy array? (@CristiFati's answer)
• [SO]: How can i cast a double pointer ctype to numpy array? (@CristiFati's answer)