小编Max*_*sky的帖子

我试图使用Fortran中的内置函数生成在[0,1]中均匀分布的随机数$\xi_i $序列.序列必须是可重现的,所以我想通过索引$ i $(即序列中$\xi_i $的位置)对随机数生成器播种,而不是使用系统时钟作为种子.以下是我的代码:

module rand

contains

  function generate_random(iseed) result(xi1)
  !!
  implicit none
  integer, intent(in) :: iseed

  integer, dimension(:), allocatable :: seed
  integer :: i, j, n
  real :: xi1

  !!-generate a seed first
  call random_seed(size = n)
  allocate(seed(n))
  seed = iseed * (/(i, i=1,n,1)/)
  call random_seed(PUT = seed)
  deallocate(seed)

  call random_number(xi1)

  !!
end function generate_random

end module rand

program test

  use rand
  implicit none
  integer :: i, imax

  imax=100

  do i=1,imax
     print *, generate_random(i)
  enddo

end program test …