rngtest.f90

Go to the documentation of this file.

!#############################################################################
!#                                                                           #
!# fosite - 3D hydrodynamical simulation program                             #
!# module: rngtest.f90                                                       #
!#                                                                           #
!# Copyright (C) 2015-2018                                                   #
!# Manuel Jung <mjung@astrophysik.uni-kiel.de>                               #
!#                                                                           #
!# This program is free software; you can redistribute it and/or modify      #
!# it under the terms of the GNU General Public License as published by      #
!# the Free Software Foundation; either version 2 of the License, or (at     #
!# your option) any later version.                                           #
!#                                                                           #
!# This program is distributed in the hope that it will be useful, but       #
!# WITHOUT ANY WARRANTY; without even the implied warranty of                #
!# MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE, GOOD TITLE or        #
!# NON INFRINGEMENT.  See the GNU General Public License for more            #
!# details.                                                                  #
!#                                                                           #
!# You should have received a copy of the GNU General Public License         #
!# along with this program; if not, write to the Free Software               #
!# Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.                 #
!#                                                                           #
!#############################################################################

!----------------------------------------------------------------------------!
!----------------------------------------------------------------------------!
PROGRAM rngtest
  USE rngs
#ifdef NECSXAURORA
  USE asl_unified
#endif
#include "tap.h"
  IMPLICIT NONE
  !--------------------------------------------------------------------------!
  REAL               :: rave,rmax,rmin,rnew
#ifdef NECSXAURORA
  REAL, DIMENSION(:), POINTER :: r
  INTEGER :: rng, imax
#else
  REAL               :: r
  INTEGER(KIND=I8)   :: i,x,x0,imax
#endif
  !--------------------------------------------------------------------------!

#ifdef NECSXAURORA
  tap_plan(5)

  CALL asl_library_initialize()
  CALL asl_random_create(rng, asl_randommethod_mt19937_64)
  CALL asl_random_distribute_uniform(rng)

  imax = 100000000_i8

  ALLOCATE(r(imax))

  CALL asl_random_generate_d(rng, imax, r)

  rave = sum(r)/imax
  rmin = minval(r)
  rmax = maxval(r)

  CALL asl_random_destroy(rng)
  CALL asl_library_finalize()

#else
  tap_plan(6)

  ! Check Kiss64
  DO i=1, 100000000
    x = kiss64()
  END DO

  x0 = 1666297717051644203_i8
  tap_check(x.EQ.x0,"Kiss64")

  imax = 100000000_i8
  r = 0.
  rmin = 1.
  rmax = 0.
  DO i=1, imax
   rnew = dkiss64()
   r = r + rnew
   rmin = min(rmin,rnew)
   rmax = max(rmax,rnew)
  END DO

  rave = r/imax
#endif


  ! Check if the random numbers are in (0,1)
  ! and if the average is near 0.5
  tap_check_close(rave,0.5,1.e-4,"Average close to 0.5.")
  tap_check_ge(rmin,0.,"All are bigger (or equal) than 0.")
  tap_check_le(rmax,1.,"All are smaller (or equal) than 1.")
  tap_check_close(rmin,0.,1.e-4,"Lower limit is close to 0.")
  tap_check_close(rmax,1.,1.e-4,"Upper limit is close to 1.")

#ifdef NECSXAURORA
  DEALLOCATE(r)
#endif
  ! Check SuperKiss64
  tap_done

END PROGRAM rngtest