1 ! Check that the real(4) and real(8) random number generators return the same
4 integer, dimension(:), allocatable
:: seed
5 real(kind
=4), dimension(10) :: r4
6 real(kind
=8), dimension(10) :: r8
7 real, parameter :: delta
= 0.0001
10 call random_seed (size
=n
)
12 call random_seed (get
=seed
)
13 ! Test both array valued and scalar routines.
14 call random_number(r4
)
15 call random_number (r4(10))
17 ! Reset the seed and get the real(8) values.
18 call random_seed (put
=seed
)
19 call random_number(r8
)
20 call random_number (r8(10))
22 if (any ((r4
- r8
) .gt
. delta
)) STOP 1