← Back to branch summary

~maddevelopers/mg5amcnlo/simple_unlops

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/short_ML_SMQCD_default/ddx_ttx/CT_interface.f

  • Committer: olivier-mattelaer
  • Date: 2021-11-12 09:29:31 UTC
  • mfrom: (967.1.15 3.3.0)
  • Revision ID: olivier-mattelaer-20211112092931-4ec9qfzgxkeovqog
version 3.3.0

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/short_ML_SMQCD_default/ddx_ttx/CT_interface.f
152
152
 
153
153
      END
154
154
 
 
155
      SUBROUTINE ML5_0_LOOP_2_2( LID, W1, W2, M1,MP_M1, M2,MP_M2, C1
 
156
     $ ,MP_C1, C2,MP_C2,  RANK, LSYMFACT, LMULTIPLIER, AMPLN, RES,
 
157
     $  STABLE)
 
158
 
 
159
      INTEGER    NEXTERNAL
 
160
      PARAMETER (NEXTERNAL=4)
 
161
      INTEGER    MAXLCOUPLINGS
 
162
      PARAMETER (MAXLCOUPLINGS=4)
 
163
      INTEGER    NLOOPLINE
 
164
      PARAMETER (NLOOPLINE=2)
 
165
      INTEGER    NWAVEFUNCS
 
166
      PARAMETER (NWAVEFUNCS=6)
 
167
      INTEGER    NCOMB
 
168
      PARAMETER (NCOMB=16)
 
169
C     
 
170
C     ARGUMENTS 
 
171
C     
 
172
      INTEGER W1, W2
 
173
      COMPLEX*16 M1, M2
 
174
      COMPLEX*32 MP_M1, MP_M2
 
175
      COMPLEX*16 C1, C2
 
176
      COMPLEX*32 MP_C1, MP_C2
 
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,NCOMB)
 
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
      M2L(1)=M2**2
 
216
      M2L(2)=M1**2
 
217
      ML(1)=M2
 
218
      ML(2)=M2
 
219
      MP_ML(1)=MP_M2
 
220
      MP_ML(2)=MP_M2
 
221
      ML(3)=M1
 
222
      MP_ML(3)=MP_M1
 
223
      ML(4)=M2
 
224
      MP_ML(4)=MP_M2
 
225
      DO I=1,NLOOPLINE
 
226
        PAIRING(I)=1
 
227
      ENDDO
 
228
 
 
229
      LC(1)=C1
 
230
      MP_LC(1)=MP_C1
 
231
      LC(2)=C2
 
232
      MP_LC(2)=MP_C2
 
233
      AMPLNUM=AMPLN
 
234
      ID=LID
 
235
      SYMFACT=LSYMFACT
 
236
      MULTIPLIER=LMULTIPLIER
 
237
      DO I=0,3
 
238
        TEMP=1
 
239
        DO J=1,NLOOPLINE
 
240
          PL(I,J)=0.D0
 
241
          DO K=TEMP,(TEMP+PAIRING(J)-1)
 
242
            PL(I,J)=PL(I,J)-DBLE(W(1+I,WE(K),VALIDH))
 
243
          ENDDO
 
244
          TEMP=TEMP+PAIRING(J)
 
245
        ENDDO
 
246
      ENDDO
 
247
      CALL ML5_0_CTLOOP(NLOOPLINE,PL,M2L,RANK,RES,STABLE)
 
248
 
 
249
      END
 
250
 
 
251
      SUBROUTINE ML5_0_LOOP_3_3( LID, W1, W2, W3, M1,MP_M1, M2,MP_M2,
 
252
     $  M3,MP_M3, C1,MP_C1, C2,MP_C2, C3,MP_C3,  RANK, LSYMFACT,
 
253
     $  LMULTIPLIER, AMPLN, RES, STABLE)
 
254
 
 
255
      INTEGER    NEXTERNAL
 
256
      PARAMETER (NEXTERNAL=4)
 
257
      INTEGER    MAXLCOUPLINGS
 
258
      PARAMETER (MAXLCOUPLINGS=4)
 
259
      INTEGER    NLOOPLINE
 
260
      PARAMETER (NLOOPLINE=3)
 
261
      INTEGER    NWAVEFUNCS
 
262
      PARAMETER (NWAVEFUNCS=6)
 
263
      INTEGER    NCOMB
 
264
      PARAMETER (NCOMB=16)
 
265
C     
 
266
C     ARGUMENTS 
 
267
C     
 
268
      INTEGER W1, W2, W3
 
269
      COMPLEX*16 M1, M2, M3
 
270
      COMPLEX*32 MP_M1, MP_M2, MP_M3
 
271
      COMPLEX*16 C1, C2, C3
 
272
      COMPLEX*32 MP_C1, MP_C2, MP_C3
 
273
 
 
274
      COMPLEX*16 RES(3)
 
275
      INTEGER LID, RANK, LSYMFACT, LMULTIPLIER
 
276
      INTEGER AMPLN
 
277
      LOGICAL STABLE
 
278
C     
 
279
C     LOCAL VARIABLES 
 
280
C     
 
281
      REAL*8 PL(0:3,NLOOPLINE)
 
282
      COMPLEX*16 M2L(NLOOPLINE)
 
283
      INTEGER PAIRING(NLOOPLINE)
 
284
      INTEGER I, J, K, TEMP
 
285
C     
 
286
C     GLOBAL VARIABLES
 
287
C     
 
288
      INTEGER WE(NEXTERNAL)
 
289
      INTEGER ID, SYMFACT, MULTIPLIER, AMPLNUM
 
290
      COMMON/ML5_0_LOOP/WE,ID, SYMFACT, MULTIPLIER,AMPLNUM
 
291
 
 
292
      COMPLEX*16 LC(MAXLCOUPLINGS)
 
293
      COMPLEX*16 ML(NEXTERNAL+2)
 
294
      COMMON/ML5_0_DP_LOOP/LC,ML
 
295
 
 
296
      COMPLEX*32 MP_LC(MAXLCOUPLINGS)
 
297
      COMPLEX*32 MP_ML(NEXTERNAL+2)
 
298
      COMMON/ML5_0_MP_LOOP/MP_LC,MP_ML
 
299
 
 
300
      COMPLEX*16 W(20,NWAVEFUNCS,NCOMB)
 
301
      INTEGER VALIDH
 
302
      COMMON/ML5_0_WFCTS/W
 
303
      COMMON/ML5_0_VALIDH/VALIDH
 
304
 
 
305
C     ----------
 
306
C     BEGIN CODE
 
307
C     ----------
 
308
 
 
309
      WE(1)=W1
 
310
      WE(2)=W2
 
311
      WE(3)=W3
 
312
      M2L(1)=M3**2
 
313
      M2L(2)=M1**2
 
314
      M2L(3)=M2**2
 
315
      ML(1)=M3
 
316
      ML(2)=M3
 
317
      MP_ML(1)=MP_M3
 
318
      MP_ML(2)=MP_M3
 
319
      ML(3)=M1
 
320
      MP_ML(3)=MP_M1
 
321
      ML(4)=M2
 
322
      MP_ML(4)=MP_M2
 
323
      ML(5)=M3
 
324
      MP_ML(5)=MP_M3
 
325
      DO I=1,NLOOPLINE
 
326
        PAIRING(I)=1
 
327
      ENDDO
 
328
 
 
329
      LC(1)=C1
 
330
      MP_LC(1)=MP_C1
 
331
      LC(2)=C2
 
332
      MP_LC(2)=MP_C2
 
333
      LC(3)=C3
 
334
      MP_LC(3)=MP_C3
 
335
      AMPLNUM=AMPLN
 
336
      ID=LID
 
337
      SYMFACT=LSYMFACT
 
338
      MULTIPLIER=LMULTIPLIER
 
339
      DO I=0,3
 
340
        TEMP=1
 
341
        DO J=1,NLOOPLINE
 
342
          PL(I,J)=0.D0
 
343
          DO K=TEMP,(TEMP+PAIRING(J)-1)
 
344
            PL(I,J)=PL(I,J)-DBLE(W(1+I,WE(K),VALIDH))
 
345
          ENDDO
 
346
          TEMP=TEMP+PAIRING(J)
 
347
        ENDDO
 
348
      ENDDO
 
349
      CALL ML5_0_CTLOOP(NLOOPLINE,PL,M2L,RANK,RES,STABLE)
 
350
 
 
351
      END
 
352
 
155
353
      SUBROUTINE ML5_0_LOOP_4_4( LID, W1, W2, W3, W4, M1,MP_M1, M2
156
354
     $ ,MP_M2, M3,MP_M3, M4,MP_M4, C1,MP_C1, C2,MP_C2, C3,MP_C3, C4
157
355
     $ ,MP_C4,  RANK, LSYMFACT, LMULTIPLIER, AMPLN, RES, STABLE)
260
458
 
261
459
      END
262
460
 
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=6)
275
 
      INTEGER    NCOMB
276
 
      PARAMETER (NCOMB=16)
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,NCOMB)
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),VALIDH))
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
 
 
365
 
      SUBROUTINE ML5_0_LOOP_2_2( LID, W1, W2, M1,MP_M1, M2,MP_M2, C1
366
 
     $ ,MP_C1, C2,MP_C2,  RANK, LSYMFACT, LMULTIPLIER, AMPLN, RES,
367
 
     $  STABLE)
368
 
 
369
 
      INTEGER    NEXTERNAL
370
 
      PARAMETER (NEXTERNAL=4)
371
 
      INTEGER    MAXLCOUPLINGS
372
 
      PARAMETER (MAXLCOUPLINGS=4)
373
 
      INTEGER    NLOOPLINE
374
 
      PARAMETER (NLOOPLINE=2)
375
 
      INTEGER    NWAVEFUNCS
376
 
      PARAMETER (NWAVEFUNCS=6)
377
 
      INTEGER    NCOMB
378
 
      PARAMETER (NCOMB=16)
379
 
C     
380
 
C     ARGUMENTS 
381
 
C     
382
 
      INTEGER W1, W2
383
 
      COMPLEX*16 M1, M2
384
 
      COMPLEX*32 MP_M1, MP_M2
385
 
      COMPLEX*16 C1, C2
386
 
      COMPLEX*32 MP_C1, MP_C2
387
 
 
388
 
      COMPLEX*16 RES(3)
389
 
      INTEGER LID, RANK, LSYMFACT, LMULTIPLIER
390
 
      INTEGER AMPLN
391
 
      LOGICAL STABLE
392
 
C     
393
 
C     LOCAL VARIABLES 
394
 
C     
395
 
      REAL*8 PL(0:3,NLOOPLINE)
396
 
      COMPLEX*16 M2L(NLOOPLINE)
397
 
      INTEGER PAIRING(NLOOPLINE)
398
 
      INTEGER I, J, K, TEMP
399
 
C     
400
 
C     GLOBAL VARIABLES
401
 
C     
402
 
      INTEGER WE(NEXTERNAL)
403
 
      INTEGER ID, SYMFACT, MULTIPLIER, AMPLNUM
404
 
      COMMON/ML5_0_LOOP/WE,ID, SYMFACT, MULTIPLIER,AMPLNUM
405
 
 
406
 
      COMPLEX*16 LC(MAXLCOUPLINGS)
407
 
      COMPLEX*16 ML(NEXTERNAL+2)
408
 
      COMMON/ML5_0_DP_LOOP/LC,ML
409
 
 
410
 
      COMPLEX*32 MP_LC(MAXLCOUPLINGS)
411
 
      COMPLEX*32 MP_ML(NEXTERNAL+2)
412
 
      COMMON/ML5_0_MP_LOOP/MP_LC,MP_ML
413
 
 
414
 
      COMPLEX*16 W(20,NWAVEFUNCS,NCOMB)
415
 
      INTEGER VALIDH
416
 
      COMMON/ML5_0_WFCTS/W
417
 
      COMMON/ML5_0_VALIDH/VALIDH
418
 
 
419
 
C     ----------
420
 
C     BEGIN CODE
421
 
C     ----------
422
 
 
423
 
      WE(1)=W1
424
 
      WE(2)=W2
425
 
      M2L(1)=M2**2
426
 
      M2L(2)=M1**2
427
 
      ML(1)=M2
428
 
      ML(2)=M2
429
 
      MP_ML(1)=MP_M2
430
 
      MP_ML(2)=MP_M2
431
 
      ML(3)=M1
432
 
      MP_ML(3)=MP_M1
433
 
      ML(4)=M2
434
 
      MP_ML(4)=MP_M2
435
 
      DO I=1,NLOOPLINE
436
 
        PAIRING(I)=1
437
 
      ENDDO
438
 
 
439
 
      LC(1)=C1
440
 
      MP_LC(1)=MP_C1
441
 
      LC(2)=C2
442
 
      MP_LC(2)=MP_C2
443
 
      AMPLNUM=AMPLN
444
 
      ID=LID
445
 
      SYMFACT=LSYMFACT
446
 
      MULTIPLIER=LMULTIPLIER
447
 
      DO I=0,3
448
 
        TEMP=1
449
 
        DO J=1,NLOOPLINE
450
 
          PL(I,J)=0.D0
451
 
          DO K=TEMP,(TEMP+PAIRING(J)-1)
452
 
            PL(I,J)=PL(I,J)-DBLE(W(1+I,WE(K),VALIDH))
453
 
          ENDDO
454
 
          TEMP=TEMP+PAIRING(J)
455
 
        ENDDO
456
 
      ENDDO
457
 
      CALL ML5_0_CTLOOP(NLOOPLINE,PL,M2L,RANK,RES,STABLE)
458
 
 
459
 
      END
460