← Back to branch summary

~maddevelopers/mg5amcnlo/2.9.4

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/short_ML_SMQCD_LoopInduced/gg_hh/CT_interface.f

  • Committer: olivier-mattelaer
  • Date: 2020-10-31 07:50:46 UTC
  • mfrom: (290.1.15 2.8.2.py3)
  • Revision ID: olivier-mattelaer-20201031075046-e4y4pft8ugi89svi
pass to 2.8.2

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/short_ML_SMQCD_LoopInduced/gg_hh/CT_interface.f
152
152
 
153
153
      END
154
154
 
155
 
      SUBROUTINE ML5_0_LOOP_3_3( LID, W1, W2, W3, M1,MP_M1, M2,MP_M2,
156
 
     $  M3,MP_M3, C1,MP_C1, C2,MP_C2, C3,MP_C3,  RANK, LSYMFACT,
157
 
     $  LMULTIPLIER, AMPLN, RES, STABLE)
158
 
 
159
 
      INTEGER    NEXTERNAL
160
 
      PARAMETER (NEXTERNAL=4)
161
 
      INTEGER    MAXLCOUPLINGS
162
 
      PARAMETER (MAXLCOUPLINGS=4)
163
 
      INTEGER    NLOOPLINE
164
 
      PARAMETER (NLOOPLINE=3)
165
 
      INTEGER    NWAVEFUNCS
166
 
      PARAMETER (NWAVEFUNCS=5)
167
 
      INTEGER    NCOMB
168
 
      PARAMETER (NCOMB=4)
169
 
C     
170
 
C     ARGUMENTS 
171
 
C     
172
 
      INTEGER W1, W2, W3
173
 
      COMPLEX*16 M1, M2, M3
174
 
      COMPLEX*32 MP_M1, MP_M2, MP_M3
175
 
      COMPLEX*16 C1, C2, C3
176
 
      COMPLEX*32 MP_C1, MP_C2, MP_C3
177
 
 
178
 
      COMPLEX*16 RES(3)
179
 
      INTEGER LID, RANK, LSYMFACT, LMULTIPLIER
180
 
      INTEGER AMPLN
181
 
      LOGICAL STABLE
182
 
C     
183
 
C     LOCAL VARIABLES 
184
 
C     
185
 
      REAL*8 PL(0:3,NLOOPLINE)
186
 
      COMPLEX*16 M2L(NLOOPLINE)
187
 
      INTEGER PAIRING(NLOOPLINE)
188
 
      INTEGER I, J, K, TEMP
189
 
C     
190
 
C     GLOBAL VARIABLES
191
 
C     
192
 
      INTEGER WE(NEXTERNAL)
193
 
      INTEGER ID, SYMFACT, MULTIPLIER, AMPLNUM
194
 
      COMMON/ML5_0_LOOP/WE,ID, SYMFACT, MULTIPLIER,AMPLNUM
195
 
 
196
 
      COMPLEX*16 LC(MAXLCOUPLINGS)
197
 
      COMPLEX*16 ML(NEXTERNAL+2)
198
 
      COMMON/ML5_0_DP_LOOP/LC,ML
199
 
 
200
 
      COMPLEX*32 MP_LC(MAXLCOUPLINGS)
201
 
      COMPLEX*32 MP_ML(NEXTERNAL+2)
202
 
      COMMON/ML5_0_MP_LOOP/MP_LC,MP_ML
203
 
 
204
 
      COMPLEX*16 W(20,NWAVEFUNCS)
205
 
      INTEGER VALIDH
206
 
      COMMON/ML5_0_WFCTS/W
207
 
      COMMON/ML5_0_VALIDH/VALIDH
208
 
 
209
 
C     ----------
210
 
C     BEGIN CODE
211
 
C     ----------
212
 
 
213
 
      WE(1)=W1
214
 
      WE(2)=W2
215
 
      WE(3)=W3
216
 
      M2L(1)=M3**2
217
 
      M2L(2)=M1**2
218
 
      M2L(3)=M2**2
219
 
      ML(1)=M3
220
 
      ML(2)=M3
221
 
      MP_ML(1)=MP_M3
222
 
      MP_ML(2)=MP_M3
223
 
      ML(3)=M1
224
 
      MP_ML(3)=MP_M1
225
 
      ML(4)=M2
226
 
      MP_ML(4)=MP_M2
227
 
      ML(5)=M3
228
 
      MP_ML(5)=MP_M3
229
 
      DO I=1,NLOOPLINE
230
 
        PAIRING(I)=1
231
 
      ENDDO
232
 
 
233
 
      LC(1)=C1
234
 
      MP_LC(1)=MP_C1
235
 
      LC(2)=C2
236
 
      MP_LC(2)=MP_C2
237
 
      LC(3)=C3
238
 
      MP_LC(3)=MP_C3
239
 
      AMPLNUM=AMPLN
240
 
      ID=LID
241
 
      SYMFACT=LSYMFACT
242
 
      MULTIPLIER=LMULTIPLIER
243
 
      DO I=0,3
244
 
        TEMP=1
245
 
        DO J=1,NLOOPLINE
246
 
          PL(I,J)=0.D0
247
 
          DO K=TEMP,(TEMP+PAIRING(J)-1)
248
 
            PL(I,J)=PL(I,J)-DBLE(W(1+I,WE(K)))
249
 
          ENDDO
250
 
          TEMP=TEMP+PAIRING(J)
251
 
        ENDDO
252
 
      ENDDO
253
 
      CALL ML5_0_CTLOOP(NLOOPLINE,PL,M2L,RANK,RES,STABLE)
254
 
 
255
 
      END
256
 
 
257
155
      SUBROUTINE ML5_0_LOOP_4_4( LID, W1, W2, W3, W4, M1,MP_M1, M2
258
156
     $ ,MP_M2, M3,MP_M3, M4,MP_M4, C1,MP_C1, C2,MP_C2, C3,MP_C3, C4
259
157
     $ ,MP_C4,  RANK, LSYMFACT, LMULTIPLIER, AMPLN, RES, STABLE)
362
260
 
363
261
      END
364
262
 
 
263
      SUBROUTINE ML5_0_LOOP_3_3( LID, W1, W2, W3, M1,MP_M1, M2,MP_M2,
 
264
     $  M3,MP_M3, C1,MP_C1, C2,MP_C2, C3,MP_C3,  RANK, LSYMFACT,
 
265
     $  LMULTIPLIER, AMPLN, RES, STABLE)
 
266
 
 
267
      INTEGER    NEXTERNAL
 
268
      PARAMETER (NEXTERNAL=4)
 
269
      INTEGER    MAXLCOUPLINGS
 
270
      PARAMETER (MAXLCOUPLINGS=4)
 
271
      INTEGER    NLOOPLINE
 
272
      PARAMETER (NLOOPLINE=3)
 
273
      INTEGER    NWAVEFUNCS
 
274
      PARAMETER (NWAVEFUNCS=5)
 
275
      INTEGER    NCOMB
 
276
      PARAMETER (NCOMB=4)
 
277
C     
 
278
C     ARGUMENTS 
 
279
C     
 
280
      INTEGER W1, W2, W3
 
281
      COMPLEX*16 M1, M2, M3
 
282
      COMPLEX*32 MP_M1, MP_M2, MP_M3
 
283
      COMPLEX*16 C1, C2, C3
 
284
      COMPLEX*32 MP_C1, MP_C2, MP_C3
 
285
 
 
286
      COMPLEX*16 RES(3)
 
287
      INTEGER LID, RANK, LSYMFACT, LMULTIPLIER
 
288
      INTEGER AMPLN
 
289
      LOGICAL STABLE
 
290
C     
 
291
C     LOCAL VARIABLES 
 
292
C     
 
293
      REAL*8 PL(0:3,NLOOPLINE)
 
294
      COMPLEX*16 M2L(NLOOPLINE)
 
295
      INTEGER PAIRING(NLOOPLINE)
 
296
      INTEGER I, J, K, TEMP
 
297
C     
 
298
C     GLOBAL VARIABLES
 
299
C     
 
300
      INTEGER WE(NEXTERNAL)
 
301
      INTEGER ID, SYMFACT, MULTIPLIER, AMPLNUM
 
302
      COMMON/ML5_0_LOOP/WE,ID, SYMFACT, MULTIPLIER,AMPLNUM
 
303
 
 
304
      COMPLEX*16 LC(MAXLCOUPLINGS)
 
305
      COMPLEX*16 ML(NEXTERNAL+2)
 
306
      COMMON/ML5_0_DP_LOOP/LC,ML
 
307
 
 
308
      COMPLEX*32 MP_LC(MAXLCOUPLINGS)
 
309
      COMPLEX*32 MP_ML(NEXTERNAL+2)
 
310
      COMMON/ML5_0_MP_LOOP/MP_LC,MP_ML
 
311
 
 
312
      COMPLEX*16 W(20,NWAVEFUNCS)
 
313
      INTEGER VALIDH
 
314
      COMMON/ML5_0_WFCTS/W
 
315
      COMMON/ML5_0_VALIDH/VALIDH
 
316
 
 
317
C     ----------
 
318
C     BEGIN CODE
 
319
C     ----------
 
320
 
 
321
      WE(1)=W1
 
322
      WE(2)=W2
 
323
      WE(3)=W3
 
324
      M2L(1)=M3**2
 
325
      M2L(2)=M1**2
 
326
      M2L(3)=M2**2
 
327
      ML(1)=M3
 
328
      ML(2)=M3
 
329
      MP_ML(1)=MP_M3
 
330
      MP_ML(2)=MP_M3
 
331
      ML(3)=M1
 
332
      MP_ML(3)=MP_M1
 
333
      ML(4)=M2
 
334
      MP_ML(4)=MP_M2
 
335
      ML(5)=M3
 
336
      MP_ML(5)=MP_M3
 
337
      DO I=1,NLOOPLINE
 
338
        PAIRING(I)=1
 
339
      ENDDO
 
340
 
 
341
      LC(1)=C1
 
342
      MP_LC(1)=MP_C1
 
343
      LC(2)=C2
 
344
      MP_LC(2)=MP_C2
 
345
      LC(3)=C3
 
346
      MP_LC(3)=MP_C3
 
347
      AMPLNUM=AMPLN
 
348
      ID=LID
 
349
      SYMFACT=LSYMFACT
 
350
      MULTIPLIER=LMULTIPLIER
 
351
      DO I=0,3
 
352
        TEMP=1
 
353
        DO J=1,NLOOPLINE
 
354
          PL(I,J)=0.D0
 
355
          DO K=TEMP,(TEMP+PAIRING(J)-1)
 
356
            PL(I,J)=PL(I,J)-DBLE(W(1+I,WE(K)))
 
357
          ENDDO
 
358
          TEMP=TEMP+PAIRING(J)
 
359
        ENDDO
 
360
      ENDDO
 
361
      CALL ML5_0_CTLOOP(NLOOPLINE,PL,M2L,RANK,RES,STABLE)
 
362
 
 
363
      END
 
364