小编Gvx*_*ørt的帖子

我正在 Fortran 90/95 的 Monte Carlo 模拟中编写一个检查点函数，我使用的编译器是 ifort 18.0.2，在详细说明我正在使用的伪随机生成器的版本之前：

A C-program for MT19937, with initialization, improved 2002/1/26.
Coded by Takuji Nishimura and Makoto Matsumoto.

Code converted to Fortran 95 by Josi Rui Faustino de Sousa
Date: 2002-02-01

源代码见mt19937。

我的蒙特卡罗模拟代码的一般结构如下：

program montecarlo
 call read_iseed(...)
 call mc_subroutine(...)
end 

内 read_iseed

subroutine read_iseed(...)
  use mt19937

    if (Restart == 'n') then

    call system('od -vAn -N4 -td4 < /dev/urandom > '//trim(IN_ISEED)
    open(unit=7,file=trim(IN_ISEED),status='old')
    read(7,*) i
    close(7)

    !This is only used to initialise the PRNG …