Fortran与C接口适配任意维度数组及标量的问题咨询

一、为啥调用二维数组会报错？Rank不匹配的锅！ #

你碰到的这个报错，核心原因是Fortran通用接口的严格匹配规则：当你写call cset(array2d, size(array2d))时，array2d是个二维数组（Rank=2），但你通用接口里的cset_array子程序，它的哑元array(*)是个一维的假定大小数组（Rank=1）。Fortran在找通用接口对应的特定子程序时，会比对实参和哑元的维度数、类型这些关键信息，维度数对不上，自然就找不到合适的子程序，直接报错了。

那为啥直接调用cset_array(array2d, ...)或者传array2d(1,1)能跑通？

• 直接调用cset_array时，Fortran会触发序列关联：多维数组在内存里是连续存储的（Fortran默认列优先），编译器会自动把它当成一维的数据块传给array(*)，这是Fortran允许的显式兼容操作，所以能正常调用C函数。
• 传array2d(1,1)的时候，你其实传的是数组第一个元素的标量引用，这时候会匹配cset_scalar子程序；而C函数里拿着这个指针遍历整个数组的元素数，刚好覆盖了多维数组的连续内存区域，所以能“凑巧”工作，但这种写法依赖内存布局，不够规范，不推荐这么用。

二、更规范的通用接口写法：一招兼容所有维度和标量 #

要让通用接口支持任意维度的数组和标量，我们可以用Fortran的**延迟形状数组（array(..)）**写一个统一的包装子程序，不用为每个维度单独写特定接口。具体实现如下：

1. 原C代码不用改 #

你的C函数逻辑没问题，只需要连续内存块的首地址和总元素数，完全可以保留：

#include <stdlib.h>
void cset( int *array, int size ) {
    for (size_t i = 0; i < size; i++) {
        array[i] = 1;
    }
    return;
}

2. 优化后的Fortran代码 #

我们写一个通用接口，用一个子程序统一处理标量和任意维度的数组，底层再调用绑定到C的函数：

program main
    use iso_c_binding
    implicit none

    ! 通用接口：一个子程序搞定所有情况
    interface cset
        module procedure cset_handle_all
    end interface cset

    ! 测试变量
    integer :: scalar = 0
    integer :: array(3) = 0
    integer :: array2d(3,3) = 0

    ! 现在可以直接调用所有类型的参数
    call cset(scalar)
    print *, scalar

    call cset(array)
    print *, array

    call cset(array2d)
    print *, array2d

contains

    ! 统一处理标量和任意维度数组的包装子程序
    subroutine cset_handle_all(x)
        use iso_c_binding
        integer(c_int), intent(inout) :: x(..)  ! 延迟形状数组，支持任意维度（包括标量）
        type(c_ptr) :: x_ptr
        integer(c_int) :: total_size

        x_ptr = c_loc(x)  ! 获取变量/数组的首地址
        total_size = int(size(x), c_int)  ! 获取总元素数（标量的size是1）

        ! 调用底层的C绑定子程序
        call c_cset(x_ptr, total_size)
    end subroutine cset_handle_all

    ! 直接绑定到C函数的底层子程序
    subroutine c_cset(array_ptr, size) bind(C, NAME="cset")
        use iso_c_binding
        type(c_ptr), value, intent(in) :: array_ptr
        integer(c_int), value, intent(in) :: size
    end subroutine c_cset

end program main

这个写法的好处：

• 延迟形状数组x(..)自动适配任意维度的数组，甚至标量（Fortran 2008及以后把标量视为0维数组）。
• 用c_loc拿首地址、size(x)拿总元素数，完全不关心变量是标量还是几维数组，逻辑统一。
• 通用接口只需要绑定一个子程序，代码更简洁，后期维护更方便。

三、你当前的标量传递方式合理吗？没问题，但可以更简洁 #

你现在的标量传递方式是合理且可行的：

• Fortran里标量默认按引用传递，所以cset_scalar的哑元会被传成指针，和C函数的int *array参数完全兼容。
• 当你传入size=1时，C函数只修改第一个元素（也就是标量本身），逻辑完全正确。

不过上面的统一包装方式更优，不需要单独为标量写特定子程序，直接把标量当成特殊的“0维数组”处理，代码更简洁，也避免了重复逻辑。

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