Fortran仿R/Python设随机数种子失效原因及set_seed模块实现

异常原因 #

Fortran内置随机数生成逻辑和R、Python的默认行为存在本质差异，相同种子下结果不一致的核心原因有三点：

• Fortran语言标准未强制规定内置随机数生成器的具体算法，不同编译器、同一编译器的不同版本，内置RNG的实现、要求的种子数组长度都可能不同，不存在跨编译器的默认一致性。
• 多数人模仿R/Python直接传入单个整数作为种子的写法是错误的：Fortran的random_seed要求传入的种子是长度与编译器匹配的数组，而非单个整数。如果只传入单元素数组，数组其余位置会保留内存中的随机残留值，根本没有完成完整初始化。
• 不显式完整填充整个种子数组时，每次程序启动的内存初始状态随机，哪怕传入的单个整数参数相同，未被赋值的种子位取值完全随机，自然会出现两次运行结果不一致、传入200的结果不符合预期的问题。

可复用set_seed模块实现 #

下面的模块可以实现和R、Python一致的种子设置逻辑：传入单个整数即可完成全量种子位初始化，同编译环境下相同种子可以稳定生成完全一致的随机数序列。

module set_seed_mod
    implicit none
contains
    subroutine set_seed(input_seed)
        integer, intent(in) :: input_seed
        integer :: seed_size, idx
        integer, allocatable :: seed_array(:)
        ! 先获取当前环境下RNG要求的种子数组长度
        call random_seed(size = seed_size)
        allocate(seed_array(seed_size))

        ! 用固定规则将单输入种子扩展填充到整个种子数组，消除内存残留值影响
        seed_array(1) = input_seed
        do idx = 2, seed_size
            seed_array(idx) = mod(seed_array(idx-1)*1103515245 + 12345, 2147483647)
        end do

        ! 写入完整种子数组完成初始化
        call random_seed(put = seed_array)
        deallocate(seed_array)
    end subroutine set_seed
end module set_seed_mod

使用说明 #

• 程序中调用random_number生成随机数前，先通过use set_seed_mod引入模块，执行call set_seed(你的种子整数值)即可固定随机序列。相同种子下重复调用、重复运行程序，生成的随机数结果完全一致。
• 如果需要跨编译器、跨编程语言复现实验结果，不要依赖Fortran内置RNG，可替换为实现了标准算法（如MT19937、PCG）的第三方Fortran RNG代码。

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