从 C 调用 Fortran 77

Question

我在 fortran 77 中有一个 ema 函数。我如何从 C 调用 fortran 函数。我知道在 fortran 2003 中有一个 iso_c_bindings 但我需要一些针对旧 77 标准的解决方案。

subroutine ema(inarray,leninarray,n,outarray,lenoutarray)

integer leninarray, lenoutarray, n, i
double precision inarray(leninarray), outarray(lenoutarray)

do 10 i = n+1, leninarray
    outarray = inarray(i) * (1/n) + outarray(i-1) * (1-(1/n))

continue
end

感谢您的帮助 问候

Answer 1

Fortran 的调用约定与 C 几乎相同，但有以下差异：

• 一切都通过引用传递：您需要将地址传递给 fortran 例程
• 参数的顺序是相同的（从从右到左）
• 调用 fortran 函数（不是 fortran 子例程）时，需要将返回值的地址压入堆栈中的参数后面。在相应的 C 原型中，这相当于将函数声明为 void，并具有一个额外的 first 参数，该参数是指向返回值的指针。
• 有时，您需要在 Fortran 名称后附加一个额外的下划线。这需要一些实验。

例如，如果您有一个 Fortran 例程，

subroutine ema(inarray,leninarray,n,outarray,lenoutarray)

它会变成

void ema(double*, int*, int*, double*, int*)

或

void ema_(double*, int*, int*, double*, int*)

取决于您的机器。在 C++ 中，

extern "C" void ema(double*, int*, int*, double*, int*)

您希望或不要忘记 fortran 数组从 1 开始。

extern "C" void ema_(double*, int*, int*, double*, int*)

当您将索引传递给例程时，