Fortran随机数生成的种子精确定义及使用方法问询

(a) 如何调用random_seed() #

random_seed是Fortran内置的随机数种子控制子程序，有两种调用方式：

• 无参数调用：call random_seed()，自动用系统随机值（如系统时间）初始化种子，每次运行得到的随机序列不同，适合非重复模拟场景。
• 带参数调用：通过SIZE、GET、PUT三个可选参数精确控制种子状态，这是固定种子、复现模拟结果的核心方式。

(b) 如何初始化精确整数种子（如123456789） #

Fortran随机数生成器需要整数数组作为种子状态，单个整数需先映射到数组中。示例代码：

program fixed_seed_demo
    implicit none
    integer :: seed_len, i
    integer, allocatable :: seed(:)
    integer, parameter :: target_seed = 123456789

    ! 获取种子数组长度
    call random_seed(size=seed_len)
    allocate(seed(seed_len))
    ! 将目标整数映射到数组（避免全同元素，部分编译器不支持）
    do i = 1, seed_len
        seed(i) = target_seed + i*1000
    end do
    ! 设置种子
    call random_seed(put=seed)

    ! 测试生成随机数
    real :: r
    call random_number(r)
    print *, "第一个随机数：", r
end program fixed_seed_demo

该初始化步骤必须放在所有random_number调用之前，且仅执行一次。

(c) random_number()的调用时机 #

• 必须在random_seed初始化完成后调用。
• 针对布朗动力学模拟：
  1. 程序启动时完成种子初始化（固定或自动），仅执行一次。
  2. 每个时间步调用Box-Muller子程序时，由该子程序内部调用random_number()生成均匀随机数，再转换为高斯分布。
  3. 禁止在每个时间步重新初始化random_seed，否则会导致随机序列重复，破坏模拟随机性。

(d) random_seed()需传入哪些参数 #

random_seed包含三个可选关键字参数（均为整数类型）：

• SIZE：获取当前随机数生成器所需的种子数组长度。
• GET：获取当前种子状态，保存到整数数组中，用于后续恢复模拟。
• PUT：传入整数数组设置种子状态，实现固定种子或恢复模拟进度。
  三个参数可单独使用，通常不同时用GET和PUT。

(e) SIZE、GET、PUT参数的用法 #

SIZE参数

用于查询种子数组长度（不同编译器要求的长度可能不同）：

integer :: seed_len
call random_seed(size=seed_len)
! seed_len即为当前编译器所需的种子数组元素个数

GET参数

用于保存当前随机数生成器状态，方便后续恢复模拟：

integer, allocatable :: saved_seed(:)
integer :: seed_len

call random_seed(size=seed_len)
allocate(saved_seed(seed_len))
call random_seed(get=saved_seed)
! 可将saved_seed保存到文件或变量，后续用PUT恢复

PUT参数

用于设置种子状态，分两种场景：

1. 固定种子：如问题(b)示例，传入构造好的整数数组，确保每次运行随机序列完全一致。
2. 恢复状态：用之前GET保存的种子数组，回到之前的模拟进度：

call random_seed(put=saved_seed)
! 后续调用random_number会从保存的状态继续生成随机数

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