如何让含递归子程序的Fortran快排程序在Windows gfortran下支持大数组？

问题根源 #

Rosetta Code上的示例快排代码默认使用栈内存存储数组，且采用递归实现。Windows系统默认栈空间仅约1MB：

• 大数组直接占用栈内存，容易触发栈溢出；
• 快排最坏递归深度为O(n)，500层递归的栈帧会快速耗尽系统栈，导致程序挂起。

具体解决方案 #

1. 改用堆分配数组（最直接有效） #

将固定大小的栈数组替换为可分配数组，堆内存无栈空间限制：
原代码中类似integer :: arr(nl*10)的固定数组定义，修改为：

program quicksort
  integer, parameter :: nl = 100
  integer, allocatable :: arr(:)
  ! 分配500个元素的堆数组
  allocate(arr(nl*5))
  ! 数组初始化、调用快排逻辑...
  call quicksort(arr, 1, size(arr))
  ! 使用完释放内存
  deallocate(arr)
end program

2. 替换为迭代版快排（避免递归栈耗尽） #

递归版快排的系统栈占用随递归深度增长，换成迭代版用自定义栈模拟递归，完全避开系统栈限制：

subroutine quicksort_iterative(arr, left, right)
  integer, intent(inout) :: arr(:)
  integer, intent(in) :: left, right
  integer, allocatable :: stack(:)
  integer :: top, l, r, pivot, temp, i, j

  ! 自定义栈，大小足够存储所有待处理区间的左右边界
  allocate(stack(2*(right-left+1)))
  top = 1
  stack(top*2-1) = left
  stack(top*2) = right

  do while (top > 0)
    r = stack(top*2)
    l = stack(top*2-1)
    top = top - 1

    pivot = arr((l+r)/2)
    i = l
    j = r

    ! 分区操作
    do while (i <= j)
      do while (arr(i) < pivot)
        i = i + 1
      end do
      do while (arr(j) > pivot)
        j = j - 1
      end do
      if (i <= j) then
        temp = arr(i)
        arr(i) = arr(j)
        arr(j) = temp
        i = i + 1
        j = j - 1
      end if
    end do

    ! 将未排序的区间压入自定义栈
    if (l < j) then
      top = top + 1
      stack(top*2-1) = l
      stack(top*2) = j
    end if
    if (i < r) then
      top = top + 1
      stack(top*2-1) = i
      stack(top*2) = r
    end if
  end do
  deallocate(stack)
end subroutine quicksort_iterative

3. 增大系统栈空间（临时应急方案） #

如果坚持使用递归版和栈数组，可通过编译链接选项直接增大程序栈空间。例如设置栈大小为8MB：

gfortran -frecursive -Wl,--stack,8388608 your_code.f90 -o quicksort.exe

注：该方法仅缓解栈溢出问题，长期来看堆分配+迭代版是更稳健的方案。

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