From 9e58b465303c25b12412a37135e4fc4cd0fe4bf2 Mon Sep 17 00:00:00 2001 From: hristov Date: Mon, 15 Jul 2002 13:53:23 +0000 Subject: [PATCH] Use double precision constants --- GEANT321/peanut/bimsel.F | 13 ++++++++----- 1 file changed, 8 insertions(+), 5 deletions(-) diff --git a/GEANT321/peanut/bimsel.F b/GEANT321/peanut/bimsel.F index 805d26a2843..78dd23293ec 100644 --- a/GEANT321/peanut/bimsel.F +++ b/GEANT321/peanut/bimsel.F @@ -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 * @@ -296,23 +299,23 @@ 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 ) -- 2.31.1