Use double precision constants
authorhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 15 Jul 2002 13:53:23 +0000 (13:53 +0000)
committerhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 15 Jul 2002 13:53:23 +0000 (13:53 +0000)
GEANT321/peanut/bimsel.F

index 805d26a2843ec1ac4cf3aab28f18e4daab33acc6..78dd23293ecbc016ef700b1c7637cfde556d9a07 100644 (file)
@@ -2,6 +2,9 @@
 * $Id$
 *
 * $Log$
+* Revision 1.5  2002/06/25 08:17:33  hristov
+* Additional protection
+*
 * Revision 1.4  2002/06/22 10:50:18  hristov
 * Better protection
 *
             IF ( BIMPTR .GE. RADTOT ) THEN
                X1 = BIMPTR - RADTOT
                DUMMY = 2.D+00 * X1 / ( RADPRO + X1 )
-               IF (ABS(DUMMY).GT.1.0) THEN
+               IF (ABS(DUMMY).GT.1.D0) THEN
                   PRINT *,"Warning in GEANT321/peanut/bimsel.F "
                   PRINT *,"Illegal ACOS argument ",DUMMY
-                  DUMMY = SIGN(1.0, DUMMY )
+                  DUMMY = SIGN(1.D0, DUMMY )
                ENDIF
                ANGRED = ACOS ( DUMMY ) / PI
                X1 = X1 / ( R0PROT * RPRONU * RADCO2 )
                DSKRED = ( 0.5D+00 * X1 * X1 + X1 + 1.D+00 ) * EXP (-X1)
      &                * ANGRED
-               IF (DSKRED.EQ.0.0) DSKRED=1.D+00
+               IF (DSKRED.EQ.0.D0) DSKRED=1.D+00
             ELSE
                X1 = RADPRO + BIMPTR - RADTOT
                DUMMY = 2.D+00 * X1 / ( RADPRO + X1 )
-               IF (ABS(DUMMY).GT.1.0) THEN
+               IF (ABS(DUMMY).GT.1.D0) THEN
                   PRINT *,"Warning in GEANT321/peanut/bimsel.F "
                   PRINT *,"Illegal ACOS argument ",DUMMY
-                  DUMMY = SIGN(1.0, DUMMY )
+                  DUMMY = SIGN(1.D0, DUMMY )
                ENDIF
                ANGRED = ACOS ( DUMMY ) / PI
                X1 = X1 / ( R0PROT * RPRONU * RADCO2 )