ranl_normal 

Routine to generate normally distributed random numbers using summation of uniformly distributed random numbers

FORMAT

  RANL_NORMAL (s1, s2, nsum vnormal, n)

Arguments

  s1                  integer*4
                      On entry, s1>=1, the first part of the two-integer
                      seed.
                      On exit, s1 is changed according to n*nsum steps of the
                      L'Ecuyer method.

  s2                  integer*4
                      On entry, s2>=1, the second part of the two-integer
                      seed.
                      On exit, s2 is changed according to n*nsum steps of the
                      L'Ecuyer method.

  nsum                integer*4
                      On entry, the number of uniform[0,1] numbers to sum to
                      get each N(0,1) output.

  vnormal             integer*4
                      On exit, a vector vnormal(1),...,vnormal(n) of N(0,1)
                      normally distributed numbers.

  n                   integer*4
                      On entry, n>=1 specifies the number of N(0,1) results
                      to return in vnormal(1),...,vnormal(n).

Description

  The RANL_NORMAL routine returns a vector of N(0,1) normally distributed
  random numbers.

  For parallel applications using the RANL_NORMAL routine, CXML provides two
  auxiliary, input programs.  Refer to the descriptions of RANL_SKIP64 and
  RANL_SKIP2.

Example

  See the example for RANL_SKIP64.