GEANT321/gheisha/gdpdk.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:21:01  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.39  by  S.Giani
  11 *-- Author :
  12       FUNCTION GDPDK(A,B,C)
  13 C
  14 C *** DOUBLE PRECISION VERSION OF "PDK" ***
  15 C *** NVE 29-MAR-1988 CERN GENEVA ***
  16 C
  17 C CALLED BY : PHASP
  18 C ORIGIN : H.FESEFELDT (02-DEC-1986)
  19 C
  20       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
  21 C
  22 C     GDPDK = SQRT(A*A+(B*B-C*C)**2/(A*A) - 2.0*(B*B+C*C))/2.0
  23 C
  24       A2 = A*A
  25       B2 = B*B
  26       C2 = C*C
  27       IF(A2) 21,21,61
  28    61 CONTINUE
  29       ARG=A2+(B2-C2)**2/A2-2.0*(B2+C2)
  30       IF (ARG) 21,21,31
  31    21 GDPDK=0.0
  32       GOTO 41
  33    31 CONTINUE
  34       GDPDK = 0.5*SQRT(ABS(A2 + (B2-C2)**2/A2 - 2.0*(B2+C2)))
  35    41 CONTINUE
  36       RETURN
  37       END