如何让Fortran生成不同的随机数

标签：io   for   sp   on   时间   c   size   table   不同的   

用Fortran生成随机数的方法很简单，就是：

call random_seed ()
call random_number (rd)

生成随机数组可以这样：

do k = 1,10

    call random_seed ()

    call random_number (rd)

    x(k) = rd

    !do something

    ....

end do

但如果循环过程中do something这里所消耗的时间比较短，那么

所生成的随机数很可能是一样的。在小木虫上一位网友给出了相应的解决方法，请见：

program test
   real :: r
   
!  call random_seed()
   call init_random_seed()
   
   do i=1,10
      call random_number(r)
      write(*,*) r
   end do
   
end program test
subroutine init_random_seed()
   integer :: i, n, clock
   integer, dimension(:), allocatable :: seed
   
   call random_seed(size = n)
   allocate(seed(n))
   
   call system_clock(count=clock)
   
   seed = clock + 37 * (/ (i - 1, i = 1, n) /)
   call random_seed(put = seed)
   
   deallocate(seed)
end subroutine init_random_seed

如何让Fortran生成不同的随机数

标签：io   for   sp   on   时间   c   size   table   不同的   

原文地址：http://www.cnblogs.com/scicalweb/p/3949196.html