关于在 iso_c_binding 中使用 C 中的 Fortran 函数

Question

学习了几天Fortran的C互操作性，从C调用一个Fortran函数DLL。我在这里找到了这个 link：

我尝试像这样创建一个 Fortran DLL，我的编译器是 Intel Fortran 编译器：

module integration
  implicit none

contains

  function Integrate(func, a,b, intsteps) result(integral)
!DEC$ ATTRIBUTES DLLEXPORT :: Integrate
    interface
      real function func(x)
        real, intent(in) :: x
      end function func
    end interface

    real :: integral, a, b
    integer :: intsteps
    intent(in) :: a, b, intsteps
    optional :: intsteps

    real :: x, dx
    integer :: i,n
    integer, parameter :: rk = kind(x)

    n = 1000
    if (present(intsteps)) n = intsteps

    dx = (b-a)/n

    integral = 0.0_rk
    do i = 1,n
      x = a + (1.0_rk * i - 0.5_rk) * dx
      integral = integral + func(x)
    end do

    integral = integral * dx
  end function

end module integration


    real(c_float) function wrapper_integrate(func, a, b, intsteps) result(integral) bind(C, name='integrate')
  use iso_c_binding
  use integration

 interface
    function iFunc(x) bind(C)
      use, intrinsic :: iso_c_binding
      real(c_float) :: iFunc
      real(c_float), intent(in) :: x
    end function iFunc
  end interface

  type(C_FUNPTR), INTENT(IN), VALUE :: func
  real(c_float) :: a,b
  integer(c_int),intent(in) :: intsteps
  optional :: intsteps

  procedure(iFunc),pointer :: myfunc
  call c_f_procpointer(func, myfunc)

  if (present(intsteps)) then
    integral = Integrate(myfunc,a,b,intsteps)   <==error #8128
  else
    integral = Integrate(myfunc,a,b)          <==error #8128
  endif

end function wrapper_integrate

当我构建 DLL 时，错误显示如下

error #8128: The BIND attribute of the associated actual procedure differs from the BIND attribute of the dummy procedure.   [MYFUNC]

代码中指出的错误行，看起来myfunc不等于模块集成中定义的func。也许我可以使用 iso_c_binding 修改模块 integration.

有时一些 Fortrans 函数有很多参数要在函数中分配和使用，在 Fortran 中使用 iso_c_binding 并不容易。那么如何解决MYFUNC中的错误？

Answer 1

如果虚拟过程有一个显式接口，那么它的特征必须与实际过程的特征相同，除了一些在这里不相关的东西。

过程的特征包括是否具有BIND属性。

为了使特性一致，您可以使用一个没有 BIND(C) 的小中间过程来调用您提供的 BIND(C) 过程。

例如，进行一些其他更改以避免围绕 c_float 和 c_int 是否指定默认实数和整数类型的假设：

function wrapper_integrate(func, a, b, intsteps)  &
    result(integral) bind(C, name='integrate')
  use iso_c_binding
  use integration
  implicit none

  interface
    function func(x) bind(C)
      use, intrinsic :: iso_c_binding
      implicit none
      real(c_float), intent(in) :: x
      real(c_float) :: func
    end function func
  end interface

  real(c_float), intent(in) :: a,b
  integer(c_int), intent(in), optional :: intsteps
  real(c_float) :: integral

  real :: local_a, local_b
  integer :: local_intsteps

  local_a = a
  local_b = b
  if (present(intsteps)) then
    local_intsteps = intsteps
    integral = Integrate(local_func, a, b, local_intsteps)
  else
    integral = Integrate(local_func, a, b)
  end if
contains
  function local_func(x)
    real, intent(in) :: x
    real :: local_func

    real(c_float) :: local_x

    local_x = x
    local_func = func(local_x)
  end function local_func
end function wrapper_integrate

请注意，可互操作过程中存在可选参数是 Fortran 2015 的一项功能。

关于在 iso_c_binding 中使用 C 中的 Fortran 函数

About using Fortran function in C with iso_c_binding

fortran

visual-c++

fortran-iso-c-binding

intel-fortran