Question

我正在尝试将单个维度数组从FORTRAN程序传递给C.

调用C函数，但它保存的值是垃圾。但是，如果我尝试使用整数变量调用相同的函数，我可以传递所需的值。任何人都可以帮我解决这个问题吗？

我使用的代码与此类似

档案：fortran_prog.f

program test
    real*4 :: a(4)
    data a / 1,2,3,4 /
    call test_func(a)
end program test

文件：c_prog.c

int test_func(double a[]) {
    int i;

    for(i=0;i<4;i++) {
        printf("%f\n",a[i]);
    }

    return 0;
}

Answer 1

program test_Cfunc

   use iso_c_binding
   implicit none

   interface
      function test_func (a) bind (C, name="test_func")
         import
         integer (c_int) :: test_func
         real (c_double), dimension (1:4), intent (in) :: a
      end function test_func
   end interface

   real (c_double), dimension (1:4) :: a = [ 2.3, 3.4, 4.5, 5.6 ]
   integer (c_int) :: result

   result = test_func (a)
   write (*, *) result

end program test_Cfunc

使用Fortran的ISO C绑定，该解决方案可以移植到同一供应商的成对编译器，或Fortran编译器供应商支持的组合。您不必了解特定编译器的传递约定，也不必处理Fortran编译器的名称修改（由name的{{1}}子句覆盖）。使用bind块描述Fortran的C例程，指定在ISO C Binding中提供的带有Fortran类值的C类型。 “内部模块”一章中的gfortran手册中列出了各种类型。另请参阅“混合语言编程”一章。由于ISO C Binding是语言标准的一部分，因此本文档比gfortran更为通用。

Answer 2

在Fortran和C之间传递数组是一个非常重要的问题。特定的C和Fortran编译器很重要。

我看到的第一个问题是您指定double以匹配real*4。这几乎在所有平台上都无效。将C函数声明为：

int test_func (float *a)

这可以在某些平台上运行，尽管许多Fortran编译器传递“数组描述符”的地址而不是数组本身。查看Fortran编译器的文档。

无法将数组从FORTRAN传递到C.

2 个答案: