]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | DOUBLE PRECISION FUNCTION SSDLAM(A,B,C) | |
3 | C----------------------------------------------------------------------- | |
4 | C Kinematic function | |
5 | C----------------------------------------------------------------------- | |
6 | #if defined(CERNLIB_IMPNONE) | |
7 | IMPLICIT NONE | |
8 | #endif | |
9 | DOUBLE PRECISION A,B,C | |
10 | C Rewrite SSDLAM=A**2+B**2+C**2-2*A*B-2*A*C-2*B*C | |
11 | IF(A.GE.B.AND.A.GE.C) THEN | |
12 | SSDLAM=(A-B-C)**2-4*B*C | |
13 | ELSEIF(B.GE.A.AND.B.GE.C) THEN | |
14 | SSDLAM=(B-A-C)**2-4*A*C | |
15 | ELSE | |
16 | SSDLAM=(C-A-B)**2-4*A*B | |
17 | ENDIF | |
18 | RETURN | |
19 | END |