5 * Revision 1.1.1.1 1995/10/24 10:21:01 cernlib
9 #include "geant321/pilot.h"
10 *CMZ : 3.21/02 29/03/94 15.41.39 by S.Giani
14 C *** DOUBLE PRECISION VERSION OF "PDK" ***
15 C *** NVE 29-MAR-1988 CERN GENEVA ***
18 C ORIGIN : H.FESEFELDT (02-DEC-1986)
20 IMPLICIT DOUBLE PRECISION (A-H,O-Z)
22 C GDPDK = SQRT(A*A+(B*B-C*C)**2/(A*A) - 2.0*(B*B+C*C))/2.0
29 ARG=A2+(B2-C2)**2/A2-2.0*(B2+C2)
34 GDPDK = 0.5*SQRT(ABS(A2 + (B2-C2)**2/A2 - 2.0*(B2+C2)))