1 #include "isajet/pilot.h"
3 C Calculate com momentum for A-->B+C with double precision.
4 C Needed to fix bug on 32-bit machines at high energy.
5 C Ver. 7.27: Rewrite order and then take abs value to be sure.
6 #include "isajet/itapes.inc"
7 #if defined(CERNLIB_DOUBLE)
8 DOUBLE PRECISION DA,DB,DC,DVAL
10 C Convert to double precision
14 DVAL=(DA-(DB+DC))*(DA+(DB+DC))*(DA-(DB-DC))*(DA+(DB-DC))
15 C Convert back to single precision
17 DBLPCM=SQRT(ABS(VAL))/(2.*A)