[u/mrichter/AliRoot.git] / GEANT321 / gheisha / gdpdk.F

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1995/10/24 10:21:01  cernlib
* Geant
*
*
#include "geant321/pilot.h"
*CMZ :  3.21/02 29/03/94  15.41.39  by  S.Giani
*-- Author :
      FUNCTION GDPDK(A,B,C)
C
C *** DOUBLE PRECISION VERSION OF "PDK" ***
C *** NVE 29-MAR-1988 CERN GENEVA ***
C
C CALLED BY : PHASP
C ORIGIN : H.FESEFELDT (02-DEC-1986)
C
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
C
C     GDPDK = SQRT(A*A+(B*B-C*C)**2/(A*A) - 2.0*(B*B+C*C))/2.0
C
      A2 = A*A
      B2 = B*B
      C2 = C*C
      IF(A2) 21,21,61
   61 CONTINUE
      ARG=A2+(B2-C2)**2/A2-2.0*(B2+C2)
      IF (ARG) 21,21,31
   21 GDPDK=0.0
      GOTO 41
   31 CONTINUE
      GDPDK = 0.5*SQRT(ABS(A2 + (B2-C2)**2/A2 - 2.0*(B2+C2)))
   41 CONTINUE
      RETURN
      END
Commit	Line	Data
fe4da5cc	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(AA+(BB-CC)2/(AA) - 2.0(BB+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.5SQRT(ABS(A2 + (B2-C2)2/A2 - 2.0(B2+C2)))
	35	41 CONTINUE
	36	RETURN
	37	END