anchor.diff:  !  Updated: 2019-09-24
   !  Rayleigh distribution with radial part of Box-Muller.
   ! ================================================================================================ !
  -subroutine rnd_rayleigh(seriesID, rndVec, vLen)
  +subroutine rnd_rayleigh_nocut(seriesID, rndVec, vLen)
   
     use crcoall
     use mathlib_bouncer
  
positive.code:     rnd_seriesData(seriesID)%hasNormR   = .true.
    end if
  
  end subroutine rnd_normal
  
  ! ================================================================================================ !
  !  V.K. Berglyd Olsen, BE-ABP-HSS
  
negative.code:     logical          :: hasNormR   = .false. ! A normal distributed remaining number is cached
    end type type_rndRecord
    type(type_rndRecord), private, save :: rnd_seriesData(rnd_nSeries)
  
    interface rnd_normal
      module procedure rnd_normal_cut
      module procedure rnd_normal_nocut
    end interface rnd_normal
  
    interface rnd_rayleigh
      module procedure rnd_rayleigh_nocut
      module procedure rnd_rayleigh_maxcut
      module procedure rnd_rayleigh_maxmincut
    end interface rnd_rayleigh
  
    private :: rnd_normal_cut
    private :: rnd_normal_nocut
    private :: rnd_rayleigh_nocut
    private :: rnd_rayleigh_maxcut
    private :: rnd_rayleigh_maxmincut
  
  contains
  
  