anchor.diff:  !  unintrupted. This will also hold after a restart from checkpoint.
   !  At 64-bit, tail truncation occurs at 9.42 sigmas, and 6.66 sigmas at 32-bit.
   ! ================================================================================================ !
  -subroutine rnd_normal(seriesID, rndVec, vLen)
  +subroutine rnd_normal_nocut(seriesID, rndVec, vLen)
   
     use crcoall
     use mathlib_bouncer
  
positive.code:     rndVec(i) = sqrt((-two)*log_mb(rndTmp(i)))
    end do
  
  end subroutine rnd_rayleigh
  
  ! ================================================================================================ !
  !  V.K. Berglyd Olsen, BE-ABP-HSS
  !  Created: 2019-09-24
  !  Updated: 2019-09-24
  !  Irwin-Hall approximation of a normal distribution of order N = 1:10
  !  1st order is equivalen to uniform, 2nd order is triangtular, and 3rd order and up approaches a
  !  Gaussian distribution. First few orders should be faster than Box-Muller.
  ! ================================================================================================ !
  subroutine rnd_irwinHall(seriesID, rndVec, vLen, nOrder)
  
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
  
  