~maddevelopers/mg5amcnlo/nlo_reweighting

« 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-02-27 14:33:01 UTC
  • mfrom: (256.19.84 2.3.4)
  • Revision ID: olivier.mattelaer@uclouvain.be-20160227143301-zjhafyh7cjw3eamq
merge with 2.3.4 (including spin2/ninja branch)

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