GEANT321/gtrak/grndm.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.2  1998/03/02 15:33:35  japost
   6 *    A significant error on Linux is corrected:
   7 *
   8 *     Two lines of grndm were not defined for any system that was not listed.
   9 *   In such a system (and Linux was one) the behavior of grndm was incorrect.
  10 *
  11 *     A default behavior is now defined for all systems. Any system that has
  12 *   a problem with this behavior (using integer division) must be added to
  13 *   the #ifdef together with IBMALL etc.
  14 *
  15 * Revision 1.1.1.1  1995/10/24 10:21:42  cernlib
  16 * Geant
  17 *
  18 *
  19 #include "geant321/pilot.h"
  20 *CMZ :  3.21/02 29/03/94  15.41.23  by  S.Giani
  21 *-- Author :
  22       SUBROUTINE GRNDM(RVEC,LEN)
  23 C.
  24 C.    ******************************************************************
  25 C.    *                                                                *
  26 C.    *       To generate a vector RVECV of LEN random numbers         *
  27 C.    *         Copy of the CERN Library routine RANECU                *
  28 C.    *                                                                *
  29 C.    *    ==>Called by : <USER>, many GEANT routines                  *
  30 C.    *       Author    F.Carminati *********                          *
  31 C.    *                                                                *
  32 C.    ******************************************************************
  33 C.
  34 #if defined(CERNLIB_CRAY)
  35 CDIR$ STACK
  36 #endif
  37       DIMENSION RVEC(*)
  38 #if defined(CERNLIB_CRAY)
  39       REAL ISEED, ISEED1, ISEED2, K, C1, C2, IZ
  40 #endif
  41 #if defined(CERNLIB_IBMALL)||defined(CERNLIB_IBMRT)||defined(CERNLIB_APO10K)||defined(CERNLIB_HPUX)||defined(CERNLIB_MSDOS)||defined(CERNLIB_WINNT)
  42       DOUBLE PRECISION ISEED, ISEED1, ISEED2, K, C1, C2, IZ
  43 #endif
  44 #if defined(CERNLIB_IBMALL)||defined(CERNLIB_IBMRT)||defined(CERNLIB_CRAY)||defined(CERNLIB_APO10K)||defined(CERNLIB_HPUX)||defined(CERNLIB_MSDOS)||defined(CERNLIB_WINNT)
  45       PARAMETER (C1=1./53668.,C2=1./52774.)
  46 #endif
  47       PARAMETER  (MAXSEQ=215)
  48       COMMON / GCRNDM / JSEQ(2), ISEED(2,MAXSEQ)
  49 C
  50       ISEQ   = JSEQ(1)
  51       ISEED1 = ISEED(1,ISEQ)
  52       ISEED2 = ISEED(2,ISEQ)
  53       DO 100 I= 1, LEN
  54 #if defined(CERNLIB_IBMALL)||defined(CERNLIB_IBMRT)||defined(CERNLIB_CRAY)||defined(CERNLIB_APO10K)||defined(CERNLIB_HPUX)||defined(CERNLIB_MSDOS)||defined(CERNLIB_WINNT)
  55         K = INT(ISEED1*C1)
  56 #else
  57         K = ISEED1/53668
  58 #endif
  59         ISEED1 = 40014*(ISEED1 - K*53668) - K*12211
  60         IF (ISEED1 .LT. 0) ISEED1=ISEED1+2147483563
  61 C
  62 #if defined(CERNLIB_IBMALL)||defined(CERNLIB_IBMRT)||defined(CERNLIB_CRAY)||defined(CERNLIB_APO10K)||defined(CERNLIB_HPUX)||defined(CERNLIB_MSDOS)||defined(CERNLIB_WINNT)
  63         K = INT(ISEED2*C2)
  64 #else
  65         K = ISEED2/52774
  66 #endif
  67         ISEED2 = 40692*(ISEED2 - K*52774) - K* 3791
  68         IF (ISEED2 .LT. 0) ISEED2=ISEED2+2147483399
  69 C
  70         IZ = ISEED1 - ISEED2
  71         IF (IZ .LE. 0)  IZ = IZ + 2147483562
  72 C
  73         RVEC(I) = IZ * 4.6566128E-10
  74   100 CONTINUE
  75       ISEED(1,ISEQ) = ISEED1
  76       ISEED(2,ISEQ) = ISEED2
  77       END