~maddevelopers/mg5amcnlo/2.7.1.3

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/long_ML_SMQCD_optimized/gg_wmtbx/mp_helas_calls_uvct_1.f

  • Committer: olivier Mattelaer
  • Date: 2016-05-12 11:00:18 UTC
  • mfrom: (262.1.150 2.3.4)
  • Revision ID: olivier.mattelaer@uclouvain.be-20160512110018-sevb79f0wm4g8mpp
pass to 2.4.0

Show diffs side-by-side

added added

removed removed

Lines of Context:
1
1
      SUBROUTINE ML5_0_MP_HELAS_CALLS_UVCT_1(P,NHEL,H,IC)
2
2
C     
 
3
      USE ML5_0_POLYNOMIAL_CONSTANTS
3
4
      IMPLICIT NONE
4
5
C     
5
6
C     CONSTANTS
17
18
      PARAMETER (NLOOPAMPS=396)
18
19
      INTEGER    NWAVEFUNCS,NLOOPWAVEFUNCS
19
20
      PARAMETER (NWAVEFUNCS=28,NLOOPWAVEFUNCS=267)
20
 
      INTEGER MAXLWFSIZE
21
 
      PARAMETER (MAXLWFSIZE=4)
22
 
      INTEGER LOOPMAXCOEFS, VERTEXMAXCOEFS
23
 
      PARAMETER (LOOPMAXCOEFS=70, VERTEXMAXCOEFS=5)
24
21
      REAL*16     ZERO
25
22
      PARAMETER (ZERO=0.0E0_16)
26
23
      COMPLEX*32     IZERO
55
52
     $ ,MP_CT_REQ_SO_DONE,LOOP_REQ_SO_DONE,MP_LOOP_REQ_SO_DONE
56
53
     $ ,CTCALL_REQ_SO_DONE,FILTER_SO
57
54
      COMMON/ML5_0_SO_REQS/UVCT_REQ_SO_DONE,MP_UVCT_REQ_SO_DONE
58
 
     $ ,CT_REQ_SO_DONE,MP_CT_REQ_SO_DONE,LOOP_REQ_SO_DONE,MP_LOOP_REQ_S
59
 
     $ O_DONE,CTCALL_REQ_SO_DONE,FILTER_SO
 
55
     $ ,CT_REQ_SO_DONE,MP_CT_REQ_SO_DONE,LOOP_REQ_SO_DONE
 
56
     $ ,MP_LOOP_REQ_SO_DONE,CTCALL_REQ_SO_DONE,FILTER_SO
60
57
 
61
58
      COMPLEX*32 AMP(NBORNAMPS)
62
59
      COMMON/ML5_0_MP_AMPS/AMP
63
60
      COMPLEX*32 W(20,NWAVEFUNCS)
64
61
      COMMON/ML5_0_MP_W/W
65
62
 
66
 
      COMPLEX*32 WL(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE,0:NLOOPWAVEF
67
 
     $ UNCS)
 
63
      COMPLEX*32 WL(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE
 
64
     $ ,0:NLOOPWAVEFUNCS)
68
65
      COMPLEX*32 PL(0:3,0:NLOOPWAVEFUNCS)
69
66
      COMMON/ML5_0_MP_WL/WL,PL
70
67
 
84
81
 
85
82
C     Amplitude(s) for UVCT diagram with ID 135
86
83
      CALL MP_FFV1_0(W(1,5),W(1,7),W(1,6),GC_5,AMPL(1,237))
87
 
      AMPL(1,237)=AMPL(1,237)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
88
 
     $ +1.0D0*UVWFCT_T_0+1.0D0*UVWFCT_B_0)
 
84
      AMPL(1,237)=AMPL(1,237)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1+1.0D0
 
85
     $ *UVWFCT_T_0+1.0D0*UVWFCT_B_0)
89
86
C     Amplitude(s) for UVCT diagram with ID 136
90
87
      CALL MP_FFV1_0(W(1,5),W(1,7),W(1,6),GC_5,AMPL(2,238))
91
 
      AMPL(2,238)=AMPL(2,238)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0*UVWFCT_G_2_1
92
 
     $ EPS)
 
88
      AMPL(2,238)=AMPL(2,238)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0
 
89
     $ *UVWFCT_G_2_1EPS)
93
90
C     Amplitude(s) for UVCT diagram with ID 137
94
91
      CALL MP_FFV1_0(W(1,8),W(1,4),W(1,6),GC_5,AMPL(1,239))
95
 
      AMPL(1,239)=AMPL(1,239)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
96
 
     $ +1.0D0*UVWFCT_T_0+1.0D0*UVWFCT_B_0)
 
92
      AMPL(1,239)=AMPL(1,239)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1+1.0D0
 
93
     $ *UVWFCT_T_0+1.0D0*UVWFCT_B_0)
97
94
C     Amplitude(s) for UVCT diagram with ID 138
98
95
      CALL MP_FFV1_0(W(1,8),W(1,4),W(1,6),GC_5,AMPL(2,240))
99
 
      AMPL(2,240)=AMPL(2,240)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0*UVWFCT_G_2_1
100
 
     $ EPS)
 
96
      AMPL(2,240)=AMPL(2,240)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0
 
97
     $ *UVWFCT_G_2_1EPS)
101
98
C     Amplitude(s) for UVCT diagram with ID 139
102
99
      CALL MP_FFV2_0(W(1,10),W(1,9),W(1,3),GC_47,AMPL(1,241))
103
 
      AMPL(1,241)=AMPL(1,241)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
104
 
     $ +1.0D0*UVWFCT_T_0+1.0D0*UVWFCT_B_0)
 
100
      AMPL(1,241)=AMPL(1,241)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1+1.0D0
 
101
     $ *UVWFCT_T_0+1.0D0*UVWFCT_B_0)
105
102
C     Amplitude(s) for UVCT diagram with ID 140
106
103
      CALL MP_FFV2_0(W(1,10),W(1,9),W(1,3),GC_47,AMPL(2,242))
107
 
      AMPL(2,242)=AMPL(2,242)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0*UVWFCT_G_2_1
108
 
     $ EPS)
 
104
      AMPL(2,242)=AMPL(2,242)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0
 
105
     $ *UVWFCT_G_2_1EPS)
109
106
C     Amplitude(s) for UVCT diagram with ID 141
110
107
      CALL MP_FFV1_0(W(1,8),W(1,9),W(1,2),GC_5,AMPL(1,243))
111
 
      AMPL(1,243)=AMPL(1,243)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
112
 
     $ +1.0D0*UVWFCT_T_0+1.0D0*UVWFCT_B_0)
 
108
      AMPL(1,243)=AMPL(1,243)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1+1.0D0
 
109
     $ *UVWFCT_T_0+1.0D0*UVWFCT_B_0)
113
110
C     Amplitude(s) for UVCT diagram with ID 142
114
111
      CALL MP_FFV1_0(W(1,8),W(1,9),W(1,2),GC_5,AMPL(2,244))
115
 
      AMPL(2,244)=AMPL(2,244)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0*UVWFCT_G_2_1
116
 
     $ EPS)
 
112
      AMPL(2,244)=AMPL(2,244)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0
 
113
     $ *UVWFCT_G_2_1EPS)
117
114
C     Amplitude(s) for UVCT diagram with ID 143
118
115
      CALL MP_FFV2_0(W(1,11),W(1,12),W(1,3),GC_47,AMPL(1,245))
119
 
      AMPL(1,245)=AMPL(1,245)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
120
 
     $ +1.0D0*UVWFCT_T_0+1.0D0*UVWFCT_B_0)
 
116
      AMPL(1,245)=AMPL(1,245)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1+1.0D0
 
117
     $ *UVWFCT_T_0+1.0D0*UVWFCT_B_0)
121
118
C     Amplitude(s) for UVCT diagram with ID 144
122
119
      CALL MP_FFV2_0(W(1,11),W(1,12),W(1,3),GC_47,AMPL(2,246))
123
 
      AMPL(2,246)=AMPL(2,246)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0*UVWFCT_G_2_1
124
 
     $ EPS)
 
120
      AMPL(2,246)=AMPL(2,246)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0
 
121
     $ *UVWFCT_G_2_1EPS)
125
122
C     Amplitude(s) for UVCT diagram with ID 145
126
123
      CALL MP_FFV1_0(W(1,11),W(1,7),W(1,2),GC_5,AMPL(1,247))
127
 
      AMPL(1,247)=AMPL(1,247)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
128
 
     $ +1.0D0*UVWFCT_T_0+1.0D0*UVWFCT_B_0)
 
124
      AMPL(1,247)=AMPL(1,247)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1+1.0D0
 
125
     $ *UVWFCT_T_0+1.0D0*UVWFCT_B_0)
129
126
C     Amplitude(s) for UVCT diagram with ID 146
130
127
      CALL MP_FFV1_0(W(1,11),W(1,7),W(1,2),GC_5,AMPL(2,248))
131
 
      AMPL(2,248)=AMPL(2,248)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0*UVWFCT_G_2_1
132
 
     $ EPS)
 
128
      AMPL(2,248)=AMPL(2,248)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0
 
129
     $ *UVWFCT_G_2_1EPS)
133
130
C     Amplitude(s) for UVCT diagram with ID 147
134
131
      CALL MP_FFV1_0(W(1,8),W(1,12),W(1,1),GC_5,AMPL(1,249))
135
 
      AMPL(1,249)=AMPL(1,249)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
136
 
     $ +1.0D0*UVWFCT_T_0+1.0D0*UVWFCT_B_0)
 
132
      AMPL(1,249)=AMPL(1,249)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1+1.0D0
 
133
     $ *UVWFCT_T_0+1.0D0*UVWFCT_B_0)
137
134
C     Amplitude(s) for UVCT diagram with ID 148
138
135
      CALL MP_FFV1_0(W(1,8),W(1,12),W(1,1),GC_5,AMPL(2,250))
139
 
      AMPL(2,250)=AMPL(2,250)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0*UVWFCT_G_2_1
140
 
     $ EPS)
 
136
      AMPL(2,250)=AMPL(2,250)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0
 
137
     $ *UVWFCT_G_2_1EPS)
141
138
C     Amplitude(s) for UVCT diagram with ID 149
142
139
      CALL MP_FFV1_0(W(1,10),W(1,7),W(1,1),GC_5,AMPL(1,251))
143
 
      AMPL(1,251)=AMPL(1,251)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
144
 
     $ +1.0D0*UVWFCT_T_0+1.0D0*UVWFCT_B_0)
 
140
      AMPL(1,251)=AMPL(1,251)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1+1.0D0
 
141
     $ *UVWFCT_T_0+1.0D0*UVWFCT_B_0)
145
142
C     Amplitude(s) for UVCT diagram with ID 150
146
143
      CALL MP_FFV1_0(W(1,10),W(1,7),W(1,1),GC_5,AMPL(2,252))
147
 
      AMPL(2,252)=AMPL(2,252)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0*UVWFCT_G_2_1
148
 
     $ EPS)
 
144
      AMPL(2,252)=AMPL(2,252)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0
 
145
     $ *UVWFCT_G_2_1EPS)
149
146
 
150
147
      GOTO 1001
151
148
 3000 CONTINUE