Fortran 和 C 混合编程

Question

我是一名理论物理研究生，从事宇宙学研究。在我的研究过程中，我使用了相当庞大的 Fortran 代码库，并使用 C 来满足我的编程需求。

我已经能够在大量测试文件中链接这两个程序，并且它们运行得非常好。但对于他们来说，我一直使用目标文件来链接它们。但是当我尝试通过 C 运行真正的交易时，包括对 Fortran 头文件的引用。它们似乎可以很好地集成和调用，但 Fortran 头文件的格式与 C 编译器不兼容，因此当它跳转到头文件时，它开始抛出无法理解语法的错误。

例如，Fortran 头文件使用 real*8 定义双精度变量，因此当 C 读取它们时会抛出错误。文件中的注释也会发生同样的情况。

那么我想请问有什么办法可以解决这个问题吗？即使 fortran 格式的头文件可以通过 C 读取。

我浏览了互联网并发现了令人困惑的答案，我不知道该遵循哪一个。任何有关此事的帮助将不胜感激:)

Answer 1

抱歉，你太令人困惑了。什么是 Fortran 头文件？例如，您无法使用 C 编译器读取 Fortran 包含文件！两种语言差别太大了。此外，Fortran 包含文件几乎从来都不是与 C 的头文件可比的头文件。

不知道你用的是什么编译器。但如果您选择了最新的 GCC 版本（Gnu 编译器集合），那么其中包含的 Fortran 编译器能够考虑 ISO_C_BINDING 功能，这使得 Fortran-C 的耦合更加容易。

示例：

MODULE my_fortran
  USE iso_c_binding
  IMPLICIT NONE
  CONTAINS
  SUBROUTINE my_subroutine(a,b) BIND(C,name="my_sub")
    INTEGER(c_int),INTENT(in),VALUE :: a
    REAL(C_DOUBLE),INTENT(out) :: b
    ...
  END SUBROUTINE
END MODULE

名为“my_sub.h”的 C 头文件 例如：

void my_sub(int, double *);

C 文件

#include "my_sub.h"

int main(){
  double b;
  int a=3;
  my_sub(a,&b);
  ...
}

Answer 2

Fortran 通常通过引用传递变量（传递指针）。
这意味着您必须在调用的 C 程序中给出地址。

函数结果，可以按值传递，例如以下
代码调用一个名为“gamma”的 FORTRAN 函数：

double   x, y;
   ..................
   x = 2.0; 
   y = gamma_(&x) 

确保
调用程序中变量的大小与大小相同
在 Fortran 例程中：

float  --- REAL     (this is typical, but not always the case)
double --- REAL*8

Fortran 函数必须在 C 调用函数的开头声明：

extern void read_(int *small, float *medium, double *large);

注意我们必须将所有变量作为指针传递给 Fortran。尽管 Fortran 中的函数名称不区分大小写，但它在 C 声明和调用时会获得一个下划线：

 read_(&small, &medium, &large);

Fortran 函数按如下方式接收变量：

SUBROUTINE READ(small,medium,large)

  INTEGER       small
  REAL          medium
  DOUBLE        large

这些变量的精确大小取决于系统的体系结构（32 位与 64 位）位），因此您需要确认您系统上的 C 和 Fortran 中 int、float 和 double 之间的对应关系。