#! /usr/bin/env python
#
def rnglib_test04 ( ):

#*****************************************************************************80
#
## RNGLIB_TEST04 demonstrates the use of multiple streams.
#
#  Licensing:
#
#    This code is distributed under the GNU LGPL license. 
#
#  Modified:
#
#    27 May 2013
#
#  Author:
#
#    John Burkardt
#
  import platform
  from cgn_set import cgn_set
  from init_generator import init_generator
  from initialize import initialize
  from r4_uni_01 import r4_uni_01

  print ( '' )
  print ( 'RNGLIB_TEST04' )
  print ( '  Python version: %s' % ( platform.python_version ( ) ) )
  print ( '  R4_UNI_01 ( ) returns a random real number' )
  print ( '  in [0,1] using the current generator.' )
#
#  Initialize the package.
#
  print ( '' )
  print ( '  INITIALIZE initializes the random number generator.' )
  print ( '  It only needs to be called once before using the package.' )

  initialize ( )

  print ( '' )
  print ( '  Let us call generators #3, #6 and #9.' )
#
#  Use three separate generators, 3, 6 and 9.
#  Force them to start at their initial seeds.
#
  g = [ 3, 6, 9 ]
  print ( '' )
  for j in range ( 0, 3 ):
    print ( '  Initialize generator %d' % ( g[j] ) )
    cgn_set ( g[j] )
    init_generator ( 0 )
#
#  Call the generators in the order 3, 6, 9.
#
  print ( '' )
  print ( '   I    R4_UNI_01 ( 3 )  R4_UNI_01 ( 6 )  R4_UNI_01 ( 9 )' )
  print ( '' )
  for i in range ( 1, 11 ):
    print ( '  %2d' % ( i ) ),
    for j in range ( 0, 3 ):
      cgn_set ( g[j] )
      u = r4_uni_01 ( )
      print ( '  %14.6g' % ( u ) ),
    print ( '' )
#
#  Restart the generators at their initial seeds.
#
  g = [ 6, 9, 3 ]
  print ( '' )
  for j in range ( 0, 3 ):
    print ( '  Reinitialize generator %d' % ( g[j] ) )
    cgn_set ( g[j] )
    init_generator ( 0 )
#
#  Call them in a different order, same result.
#
  print ( '' )
  print ( '   I    R4_UNI_01 ( 6 )  R4_UNI_01 ( 9 )  R4_UNI_01 ( 3 )' )
  print ( '' )
  for i in range ( 1, 11 ):
    print ( '  %2d' % ( i ) ),
    for j in range ( 0, 3 ):
      cgn_set ( g[j] )
      u = r4_uni_01 ( )
      print ( '  %14.6g' % ( u ) ),
    print ( '' )
#
#  Terminate.
#
  print ( '' )
  print ( 'RNGLIB_TEST04:' )
  print ( '  Normal end of execution.' )
  return

if ( __name__ == '__main__' ):
  rnglib_test04 ( )