传递可分配数组至Fortran子程序触发SIGSEGV段错误求助

问题现象 #

• 初始未分配的二维可分配数组sort_cent，调用语句call rbind(sort_cent,tmparray([arg],:))修改数组值与大小后，程序触发signal SIGSEGV, Segmentation fault错误中断
• 子程序rbind能执行到末尾，输出End of the function，推测问题出在数组返回环节

相关代码 #

主程序调用片段 #

real, allocatable :: sort_cent(:,:), tmparray(:,:)
! 此处省略tmparray的初始化与赋值逻辑
call rbind(sort_cent, tmparray([arg],:))

rbind子程序实现 #

subroutine rbind(arr, new_row)
    real, allocatable, intent(inout) :: arr(:,:)
    real, intent(in) :: new_row(:)
    real, allocatable :: temp(:,:)
    integer :: n_rows, n_cols

    if (allocated(arr)) then
        n_rows = size(arr, 1)
        n_cols = size(arr, 2)
        if (n_cols /= size(new_row)) then
            error stop "rbind: 新行维度与原数组列数不匹配"
        end if
        allocate(temp(n_rows+1, n_cols))
        temp(1:n_rows, :) = arr
        temp(n_rows+1, :) = new_row
        deallocate(arr)
        call move_alloc(temp, arr)
    else
        n_cols = size(new_row)
        allocate(arr(1, n_cols))
        arr(1, :) = new_row
    end if
    print *, "End of the function"
end subroutine rbind

程序输出 #

End of the function
signal SIGSEGV, Segmentation fault

解决思路 #

1. 确认实参维度合法性

   • 检查tmparray([arg],:)的返回维度：确保arg是合法的索引，不会导致new_row的维度与后续数组列数不匹配。第一次调用时new_row的列数会决定sort_cent的列数，后续调用必须严格匹配该列数。
   • 可以在子程序开头添加维度检查，比如if (rank(new_row)/=1) error stop "rbind: 输入必须是一维数组"。

2. 修正内存分配与移动逻辑

   • 确保move_alloc的使用完全合规：move_alloc(temp, arr)会将temp的内存所有权转移给arr，同时自动释放temp。如果arr在deallocate后仍处于异常状态，会导致后续操作出错。可以在deallocate(arr)后添加检查：

     deallocate(arr)
     if (allocated(arr)) error stop "rbind: 原数组释放失败"
     

3. 添加显式接口

   • Fortran中传递可分配数组时，必须使用显式接口，否则会因隐式类型匹配错误导致内存访问异常。解决方式：
     • 将rbind子程序放入模块中，主程序通过use语句引用该模块；
     • 或者在主程序中声明显式接口块：

       interface
           subroutine rbind(arr, new_row)
               real, allocatable, intent(inout) :: arr(:,:)
               real, intent(in) :: new_row(:)
           end subroutine rbind
       end interface
       

4. 使用内存调试工具定位

   • 用valgrind工具运行程序，获取精确的内存错误位置：

     valgrind --leak-check=full --track-origins=yes ./your_program
     

     工具会输出具体的非法内存访问点，帮助确认是move_alloc后的数组访问错误，还是其他隐含问题。

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