使用r8mat_normal_ab生成大正态分布数组时出现段错误求助

尝试生成尺寸为(150000,35)的正态分布随机数数组，使用r8lib的r8mat_normal_ab函数时，当数组第一维度超过约27000就触发段错误。测试代码如下：

program test2
    integer, parameter :: dp=kind(0.d0)
    integer, parameter :: n_sim = 30000
    integer, parameter :: nt = 35
    real(dp) :: shockARWmat(n_sim,nt+1)
    real(dp) :: sigmaARW=sqrt(.0344705)
    integer :: seed_sim=234567, i

    call r8mat_normal_ab (n_sim, nt+1, 0., sigmaARW, seed_sim, shockARWmat)

    print*, 'success !!'

end program

问题原因 #

段错误的核心原因是静态局部数组默认分配在栈内存中，而栈的默认容量通常很小（多数系统默认栈大小在8MB-16MB之间）。当n_sim=30000时，数组总大小为30000*36*8≈8.2MB，接近栈的上限；当n_sim=150000时，数组大小会超过40MB，直接触发栈溢出。-mcmodel=large参数仅针对全局/静态数据的内存模型，对栈上的局部数组无效，所以无法解决问题。

解决方案 #

1. 使用可分配数组（推荐，跨平台兼容） #

将静态数组改为可分配数组，动态分配到堆内存（堆的容量远大于栈）。修改代码如下：

program test2
    integer, parameter :: dp=kind(0.d0)
    integer, parameter :: n_sim = 150000
    integer, parameter :: nt = 35
    real(dp), allocatable :: shockARWmat(:,:)  ! 声明可分配数组
    real(dp) :: sigmaARW=sqrt(.0344705)
    integer :: seed_sim=234567

    allocate(shockARWmat(n_sim, nt+1))  ! 动态分配堆内存
    call r8mat_normal_ab(n_sim, nt+1, 0.d0, sigmaARW, seed_sim, shockARWmat)

    print*, 'success !!'

    deallocate(shockARWmat)  ! 可选，程序结束时会自动释放内存
end program

2. 调整栈内存大小（不推荐，兼容性差） #

如果必须使用静态数组，可以临时增大栈的容量：

• Linux/macOS：编译运行前执行ulimit -s unlimited解除栈大小限制；或编译时通过链接器参数指定栈大小，比如Intel编译器：ifort -m64 -Wl,--stack,104857600 test2.f90（设置100MB栈）。
• Windows：通过Intel编译器的-heap-arrays参数将大数组分配到堆，比如ifort -heap-arrays 1024 test2.f90（指定大于1024字节的数组分配到堆）。

3. 确认64位编译模式 #

确保编译时启用64位模式（Intel编译器默认启用，也可明确添加-m64参数），避免32位程序的内存地址空间限制。

内容的提问来源于stack exchange，提问作者andrea m.