scale_random Module Reference

module RANDOM More...

Functions/Subroutines
subroutine, public	random_setup
	Setup. More...
subroutine	random_uniform_1d (var)
	Get uniform random number ( 1D ) More...
subroutine, public	random_finalize
	finalize More...
subroutine, public	random_knuth_shuffle (num, a)

Detailed Description

module RANDOM

Description

random number generation module

Author

Team SCALE

NAMELIST

• PARAM_RANDOM

  name	type	default value	comment
  RANDOM_FIX	logical	.false.

  
  

History Output

No history output

Function/Subroutine Documentation

random_setup()

subroutine, public scale_random::random_setup

Setup.

Definition at line 69 of file scale_random.F90.

    use scale_prc, only: &
       prc_abort
    implicit none
 
    namelist / param_random / &
       random_fix
 
    integer :: nseeds, ierr
    !---------------------------------------------------------------------------
 
    log_newline
    log_info("RANDOM_setup",*) 'Setup'
 
    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_random,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       log_info("RANDOM_setup",*) 'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       log_error("RANDOM_setup",*) 'Not appropriate names in namelist PARAM_RANDOM. Check!'
       call prc_abort
    endif
    log_nml(param_random)
 
    call random_seed
    call random_seed(size=nseeds)
 
    allocate( random_seedvar(nseeds))
 
    log_newline
    log_info("RANDOM_setup",*) 'Array size for random seed:', nseeds
    if ( random_fix ) then
       log_info("RANDOM_setup",*) 'random seed is fixed.'
    endif
 
    call random_reset
 
    return

References scale_io::io_fid_conf, scale_prc::prc_abort(), and scale_prc::prc_universal_myrank.

Referenced by mod_rm_driver::rm_driver(), and mod_rm_prep::rm_prep().

Here is the call graph for this function:

Here is the caller graph for this function:

random_uniform_1d()

subroutine scale_random::random_uniform_1d

(

real(rp), dimension(:), intent(out)

var

)

Get uniform random number ( 1D )

Definition at line 154 of file scale_random.F90.

    implicit none
    real(RP), intent(out) :: var(:)
    !---------------------------------------------------------------------------
 
    call random_number(var)
    !$acc update device(var) if(acc_is_present(var))
 
    return

random_finalize()

subroutine, public scale_random::random_finalize

finalize

Definition at line 229 of file scale_random.F90.

    implicit none
    !---------------------------------------------------------------------------
 
    random_fix = .false.
 
    deallocate( random_seedvar )
 
    return

Referenced by mod_rm_driver::rm_driver(), and mod_rm_prep::rm_prep().

Here is the caller graph for this function:

random_knuth_shuffle()

subroutine, public scale_random::random_knuth_shuffle	(	integer, intent(in)	num,
		integer, dimension(num), intent(out)	a
	)

Definition at line 240 of file scale_random.F90.

    implicit none
    integer, intent(in) :: num
    integer, intent(out) :: a(num)
    integer :: i, randpos, tmp
    real(DP) :: r
 
    do i = 1, num
      a(i) = i
    end do
 
    do i = num, 2, -1
      call random_number(r)
      randpos = int(r * i) + 1
      tmp = a(randpos)
      a(randpos) = a(i)
      a(i) = tmp
    end do
 
    return

References scale_const::const_pi.

Referenced by scale_letkf::letkf_system().

Here is the caller graph for this function: