~maddevelopers/mg5amcnlo/WWW5_caching

« back to all changes in this revision

Viewing changes to users/mardelcourt/PROC_242195/PROC_242195/Source/DHELAS/VVV1_0.f

  • Committer: John Doe
  • Date: 2013-03-25 20:27:02 UTC
  • Revision ID: john.doe@gmail.com-20130325202702-5sk3t1r8h33ca4p4
first clean version

Show diffs side-by-side

added added

removed removed

Lines of Context:
 
1
C     This File is Automatically generated by ALOHA 
 
2
C     The process calculated in this file is: 
 
3
C     P(3,1)*Metric(1,2) - P(3,2)*Metric(1,2) - P(2,1)*Metric(1,3) +
 
4
C      P(2,3)*Metric(1,3) + P(1,2)*Metric(2,3) - P(1,3)*Metric(2,3)
 
5
C     
 
6
      SUBROUTINE VVV1_0(V1, V2, V3, COUP,VERTEX)
 
7
      IMPLICIT NONE
 
8
      COMPLEX*16 CI
 
9
      PARAMETER (CI=(0D0,1D0))
 
10
      COMPLEX*16 V2(*)
 
11
      COMPLEX*16 TMP12
 
12
      COMPLEX*16 V3(*)
 
13
      COMPLEX*16 TMP1
 
14
      REAL*8 P1(0:3)
 
15
      COMPLEX*16 TMP10
 
16
      REAL*8 P2(0:3)
 
17
      COMPLEX*16 TMP7
 
18
      REAL*8 P3(0:3)
 
19
      COMPLEX*16 VERTEX
 
20
      COMPLEX*16 TMP13
 
21
      COMPLEX*16 TMP11
 
22
      COMPLEX*16 COUP
 
23
      COMPLEX*16 TMP9
 
24
      COMPLEX*16 V1(*)
 
25
      COMPLEX*16 TMP3
 
26
      COMPLEX*16 TMP8
 
27
      P1(0) = DBLE(V1(1))
 
28
      P1(1) = DBLE(V1(2))
 
29
      P1(2) = DIMAG(V1(2))
 
30
      P1(3) = DIMAG(V1(1))
 
31
      P2(0) = DBLE(V2(1))
 
32
      P2(1) = DBLE(V2(2))
 
33
      P2(2) = DIMAG(V2(2))
 
34
      P2(3) = DIMAG(V2(1))
 
35
      P3(0) = DBLE(V3(1))
 
36
      P3(1) = DBLE(V3(2))
 
37
      P3(2) = DIMAG(V3(2))
 
38
      P3(3) = DIMAG(V3(1))
 
39
      TMP1 = (V3(3)*V1(3)-V3(4)*V1(4)-V3(5)*V1(5)-V3(6)*V1(6))
 
40
      TMP9 = (V3(3)*P2(0)-V3(4)*P2(1)-V3(5)*P2(2)-V3(6)*P2(3))
 
41
      TMP8 = (V3(3)*P1(0)-V3(4)*P1(1)-V3(5)*P1(2)-V3(6)*P1(3))
 
42
      TMP3 = (V2(3)*V1(3)-V2(4)*V1(4)-V2(5)*V1(5)-V2(6)*V1(6))
 
43
      TMP7 = (V2(3)*V3(3)-V2(4)*V3(4)-V2(5)*V3(5)-V2(6)*V3(6))
 
44
      TMP11 = (V2(3)*P3(0)-V2(4)*P3(1)-V2(5)*P3(2)-V2(6)*P3(3))
 
45
      TMP10 = (V2(3)*P1(0)-V2(4)*P1(1)-V2(5)*P1(2)-V2(6)*P1(3))
 
46
      TMP13 = (V1(3)*P3(0)-V1(4)*P3(1)-V1(5)*P3(2)-V1(6)*P3(3))
 
47
      TMP12 = (V1(3)*P2(0)-V1(4)*P2(1)-V1(5)*P2(2)-V1(6)*P2(3))
 
48
      VERTEX = COUP*(TMP1*(-CI*(TMP11)+CI*(TMP10))+(TMP3*(-CI*(TMP8)
 
49
     $ +CI*(TMP9))+TMP7*(-CI*(TMP12)+CI*(TMP13))))
 
50
      END
 
51
 
 
52