]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | FUNCTION EPF(A,B,C,D) | |
3 | C CALCULATE TOTALLY ANTISYMMETRIC TENSOR EPSILON CONTRACTED | |
4 | C WITH FOUR 4-VECTORS. | |
5 | #include "isajet/itapes.inc" | |
6 | DIMENSION A(4),B(4),C(4),D(4) | |
7 | #if defined(CERNLIB_DOUBLE) | |
8 | DOUBLE PRECISION EPF | |
9 | DOUBLE PRECISION A,B,C,D,CD,BCD | |
10 | #endif | |
11 | CD(I,J)=C(I)*D(J)-C(J)*D(I) | |
12 | BCD(I,J,K)=B(I)*CD(J,K)-B(J)*CD(I,K)+B(K)*CD(I,J) | |
13 | EPF=A(1)*BCD(2,3,4)-A(2)*BCD(1,3,4)+A(3)*BCD(1,2,4) | |
14 | 1-A(4)*BCD(1,2,3) | |
15 | RETURN | |
16 | END |