Fortran Wiki
random_number

Description

Returns a single pseudorandom number or an array of pseudorandom numbers from the uniform distribution over the range 0x<1.

The runtime-library implements George Marsaglia’s KISS (Keep It Simple Stupid) random number generator (RNG). This RNG combines:

  1. The congruential generator x(n)=69069x(n1)+1327217885 with a period of 2 32,

  2. A 3-shift shift-register generator with a period of 2 321,

  3. Two 16-bit multiply-with-carry generators with a period of 597273182964842497>2 59.

The overall period exceeds 2 123.

Please note, this RNG is thread safe if used within OpenMP directives, i.e., its state will be consistent while called from multiple threads. However, the KISS generator does not create random numbers in parallel from multiple sources, but in sequence from a single source. If an OpenMP-enabled application heavily relies on random numbers, one should consider employing a dedicated parallel random number generator instead.

Standard

Fortran 95 and later

Class

Subroutine

Syntax

random_number(harvest)

Arguments

  • harvest - Shall be a scalar or an array of type real.

Example

program test_random_number
  real :: r(5,5)
  call init_random_seed()         ! see example of random_seed
  call random_number(r)
end program

See also

random_seed

category: intrinsics