From: hristov Date: Mon, 15 Jul 2002 13:53:23 +0000 (+0000) Subject: Use double precision constants X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=commitdiff_plain;h=9e58b465303c25b12412a37135e4fc4cd0fe4bf2;hp=6379c54322f93e5ecfa4b2e0b8ea114a2f9f6f17 Use double precision constants --- 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 )