如何实现接收Fortran函数的C积分函数的Fortran-C互操作接口

需求说明 #

现有一个C库函数Calculate_Integral，功能是对传入的C浮点函数在指定区间执行定积分，函数签名已给定。需要从Fortran中调用该函数，并且传入Fortran编写的被积函数，让C函数能回调这个Fortran函数完成积分计算。

一、C端代码（保持原有实现） #

C函数的签名和实现无需修改，确保函数指针参数的类型正确，以下是梯形法实现的示例：

#include <math.h>

float Calculate_Integral(float (*func)(float), float a, float b, int n) {
    if (n <= 0 || a >= b) return 0.0f;
    
    float h = (b - a) / n;
    float sum = 0.5f * (func(a) + func(b));
    
    for (int i = 1; i < n; i++) {
        float x = a + i * h;
        sum += func(x);
    }
    
    return sum * h;
}

二、Fortran端接口与调用实现 #

1. 正确的Fortran接口定义 #

使用iso_c_binding模块实现类型和接口的匹配，修正原代码中的返回值命名错误后，完整接口如下：

module integral_interface
    use, intrinsic :: iso_c_binding
    implicit none
    
    interface
        function FCalculate_Integral(func, a, b, n) bind(C, name="Calculate_Integral")
            import :: C_FLOAT, C_INT
            implicit none
            real(C_FLOAT) :: FCalculate_Integral  ! 返回值变量名需与函数名一致
            real(C_FLOAT), value :: a, b
            integer(C_INT), value :: n
            ! 嵌套接口定义C函数指针对应的Fortran函数签名
            interface
                function func(x) bind(C)
                    import :: C_FLOAT
                    implicit none
                    real(C_FLOAT) :: func
                    real(C_FLOAT), value :: x
                end function func
            end interface
        end function FCalculate_Integral
    end interface
end module integral_interface

2. Fortran被积函数定义 #

要被C回调的Fortran函数必须添加bind(C)属性，确保签名和C的函数指针完全匹配：

! 示例被积函数：f(x) = x²
function my_func(x) bind(C)
    use, intrinsic :: iso_c_binding
    implicit none
    real(C_FLOAT) :: my_func
    real(C_FLOAT), value :: x
    my_func = x * x
end function my_func

3. Fortran主程序调用 #

在主程序中使用上述接口和被积函数完成调用：

program test_integral
    use integral_interface
    use, intrinsic :: iso_c_binding
    implicit none
    real(C_FLOAT) :: result
    real(C_FLOAT), parameter :: a = 0.0_C_FLOAT, b = 2.0_C_FLOAT
    integer(C_INT), parameter :: n = 1000_C_INT
    
    result = FCalculate_Integral(my_func, a, b, n)
    print *, "积分结果：", result
    ! 理论值为8/3 ≈ 2.6666667
end program test_integral

关键注意事项 #

• 类型严格匹配：必须使用iso_c_binding模块提供的C_FLOAT、C_INT，避免Fortran和C的基础类型长度不匹配导致错误。
• 传值属性：C默认传值，Fortran默认传引用，因此所有C函数的参数（包括回调函数的参数）都要在Fortran接口中添加value属性。
• 回调函数的bind(C)：Fortran被积函数必须加bind(C)，这样编译器会生成符合C调用约定的函数，保证C能正确调用。
• 接口命名一致性：Fortran接口的返回值变量名必须和函数名一致（原代码中返回值写的Calculate_Integral是错误的，应改为FCalculate_Integral）。
• 避免名字混淆：bind(C, name="Calculate_Integral")明确指定C函数的名字，防止Fortran编译器对函数名进行命名修饰（比如添加下划线）。

内容的提问来源于stack exchange，提问作者Akhil Akkapelli