ISAJET/code/epf.F

   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