← Back to branch summary

~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/LoopSquaredOrder_IOTest/Loop_sqso_uux_ddx/loop_matrix_QCDQEDpert_default.f

  • Committer: Paolo Torrielli
  • Date: 2016-02-18 17:03:15 UTC
  • mfrom: (262.1.71 2.3.4)
  • mto: (262.1.75 2.3.4)
  • mto: This revision was merged to the branch mainline in revision 263.
  • Revision ID: torriell@physik.uzh.ch-20160218170315-gjajinyn9929pef3
merged with newest 2.3.4

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/LoopSquaredOrder_IOTest/Loop_sqso_uux_ddx/loop_matrix_QCDQEDpert_default.f
61
61
      INTEGER MAXLWFSIZE
62
62
      PARAMETER (MAXLWFSIZE=4)
63
63
      INTEGER LOOPMAXCOEFS, VERTEXMAXCOEFS
64
 
      PARAMETER (LOOPMAXCOEFS=70, VERTEXMAXCOEFS=15)
 
64
      PARAMETER (LOOPMAXCOEFS=15, VERTEXMAXCOEFS=5)
65
65
      INTEGER    NCOMB
66
66
      PARAMETER (NCOMB=16)
67
67
      REAL*8     ZERO
832
832
          CALL FFV1P0_3(W(1,1),W(1,2),GC_2,ZERO,ZERO,W(1,6))
833
833
C         Amplitude(s) for born diagram with ID 2
834
834
          CALL FFV1_0(W(1,4),W(1,3),W(1,6),GC_1,AMP(2))
835
 
          CALL FFV2_5_3(W(1,1),W(1,2),GC_69,GC_76,MDL_MZ,MDL_WZ,W(1,7))
 
835
          CALL FFV2_5P0_3(W(1,1),W(1,2),GC_69,GC_76,MDL_MZ,MDL_WZ,W(1
 
836
     $     ,7))
836
837
C         Amplitude(s) for born diagram with ID 3
837
838
          CALL FFV2_3_0(W(1,4),W(1,3),W(1,7),GC_68,GC_76,AMP(3))
838
839
          CALL FFS2_3(W(1,1),W(1,2),GC_123,MDL_MZ,MDL_WZ,W(1,8))
841
842
          CALL FFS4_3(W(1,1),W(1,2),GC_122,MDL_MH,MDL_WH,W(1,9))
842
843
C         Amplitude(s) for born diagram with ID 5
843
844
          CALL FFS4_0(W(1,4),W(1,3),W(1,9),GC_105,AMP(5))
844
 
          CALL FFV2_3(W(1,1),W(1,3),GC_124,MDL_MW,MDL_WW,W(1,10))
 
845
          CALL FFV2P0_3(W(1,1),W(1,3),GC_124,MDL_MW,MDL_WW,W(1,10))
845
846
C         Amplitude(s) for born diagram with ID 6
846
847
          CALL FFV2_0(W(1,4),W(1,2),W(1,10),GC_59,AMP(6))
847
848
          CALL FFS1_3_3(W(1,1),W(1,3),GC_13,GC_22,MDL_MW,MDL_WW,W(1
905
906
C         Counter-term amplitude(s) for loop diagram number 27
906
907
          CALL FFS4_2_0(W(1,4),W(1,2),W(1,11),R2_DXUGM,R2_DXUGMA
907
908
     $     ,AMPL(1,33))
908
 
          CALL FFV2_3(W(1,4),W(1,2),GC_59,MDL_MW,MDL_WW,W(1,13))
 
909
          CALL FFV2P0_3(W(1,4),W(1,2),GC_59,MDL_MW,MDL_WW,W(1,13))
909
910
C         Counter-term amplitude(s) for loop diagram number 28
910
911
          CALL FFV2_0(W(1,1),W(1,3),W(1,13),R2_UXDW,AMPL(1,34))
911
912
          CALL FFS1_3_3(W(1,4),W(1,2),GC_31,GC_40,MDL_MW,MDL_WW,W(1
916
917
          CALL FFV1P0_3(W(1,4),W(1,3),GC_1,ZERO,ZERO,W(1,15))
917
918
C         Counter-term amplitude(s) for loop diagram number 32
918
919
          CALL FFV1_0(W(1,1),W(1,2),W(1,15),R2_UUA,AMPL(1,36))
919
 
          CALL FFV2_3_3(W(1,4),W(1,3),GC_68,GC_76,MDL_MZ,MDL_WZ,W(1
 
920
          CALL FFV2_3P0_3(W(1,4),W(1,3),GC_68,GC_76,MDL_MZ,MDL_WZ,W(1
920
921
     $     ,16))
921
922
C         Counter-term amplitude(s) for loop diagram number 33
922
923
          CALL FFV2_5_0(W(1,1),W(1,2),W(1,16),R2_UUZ_V2,R2_UUZ_V5
1197
1198
     $     ,43))
1198
1199
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,43),2,4,15,1,1,73,H)
1199
1200
C         Coefficient construction for loop diagram with ID 30
1200
 
          CALL FFV2L2_3(PL(0,10),W(1,4),GC_59,MDL_MW,MDL_WW,PL(0,44)
 
1201
          CALL FFV2L2P0_3(PL(0,10),W(1,4),GC_59,MDL_MW,MDL_WW,PL(0,44)
1201
1202
     $     ,COEFS)
1202
 
          CALL ML5_0_UPDATE_WL_1_2(WL(1,0,1,10),4,COEFS,4,4,WL(1,0,1
 
1203
          CALL ML5_0_UPDATE_WL_1_0(WL(1,0,1,10),4,COEFS,4,4,WL(1,0,1
1203
1204
     $     ,44))
1204
1205
          CALL FFV2L3_1(PL(0,44),W(1,3),GC_124,MDL_MU,ZERO,PL(0,45)
1205
1206
     $     ,COEFS)
1206
 
          CALL ML5_0_UPDATE_WL_3_1(WL(1,0,1,44),4,COEFS,4,4,WL(1,0,1
 
1207
          CALL ML5_0_UPDATE_WL_1_1(WL(1,0,1,44),4,COEFS,4,4,WL(1,0,1
1207
1208
     $     ,45))
1208
 
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,45),4,4,16,1,1,74,H)
 
1209
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,45),2,4,16,1,1,74,H)
1209
1210
C         Coefficient construction for loop diagram with ID 31
1210
1211
          CALL FFS1_3L2_3(PL(0,10),W(1,4),GC_31,GC_40,MDL_MW,MDL_WW
1211
1212
     $     ,PL(0,46),COEFS)
1255
1256
     $     ,54))
1256
1257
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,54),2,4,6,1,1,80,H)
1257
1258
C         Coefficient construction for loop diagram with ID 37
1258
 
          CALL FFV2_5L1_3(PL(0,12),W(1,2),GC_69,GC_76,MDL_MZ,MDL_WZ
 
1259
          CALL FFV2_5L1P0_3(PL(0,12),W(1,2),GC_69,GC_76,MDL_MZ,MDL_WZ
1259
1260
     $     ,PL(0,55),COEFS)
1260
 
          CALL ML5_0_UPDATE_WL_1_2(WL(1,0,1,12),4,COEFS,4,4,WL(1,0,1
 
1261
          CALL ML5_0_UPDATE_WL_1_0(WL(1,0,1,12),4,COEFS,4,4,WL(1,0,1
1261
1262
     $     ,55))
1262
1263
          CALL FFV2_3L3_2(PL(0,55),W(1,4),GC_68,GC_76,MDL_MD,ZERO,PL(0
1263
1264
     $     ,56),COEFS)
1264
 
          CALL ML5_0_UPDATE_WL_3_1(WL(1,0,1,55),4,COEFS,4,4,WL(1,0,1
 
1265
          CALL ML5_0_UPDATE_WL_1_1(WL(1,0,1,55),4,COEFS,4,4,WL(1,0,1
1265
1266
     $     ,56))
1266
1267
          CALL FFV1L1P0_3(PL(0,56),W(1,3),GC_11,ZERO,ZERO,PL(0,57)
1267
1268
     $     ,COEFS)
1268
 
          CALL ML5_0_UPDATE_WL_4_0(WL(1,0,1,56),4,COEFS,4,4,WL(1,0,1
 
1269
          CALL ML5_0_UPDATE_WL_2_0(WL(1,0,1,56),4,COEFS,4,4,WL(1,0,1
1269
1270
     $     ,57))
1270
 
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,57),4,4,17,1,1,81,H)
 
1271
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,57),2,4,17,1,1,81,H)
1271
1272
C         Coefficient construction for loop diagram with ID 38
1272
1273
          CALL FFS2L1_3(PL(0,12),W(1,2),GC_123,MDL_MZ,MDL_WZ,PL(0,58)
1273
1274
     $     ,COEFS)
1309
1310
C         Coefficient construction for loop diagram with ID 41
1310
1311
          CALL FFV2_3L3_1(PL(0,55),W(1,3),GC_68,GC_76,MDL_MD,ZERO,PL(0
1311
1312
     $     ,66),COEFS)
1312
 
          CALL ML5_0_UPDATE_WL_3_1(WL(1,0,1,55),4,COEFS,4,4,WL(1,0,1
 
1313
          CALL ML5_0_UPDATE_WL_1_1(WL(1,0,1,55),4,COEFS,4,4,WL(1,0,1
1313
1314
     $     ,66))
1314
1315
          CALL FFV1L2P0_3(PL(0,66),W(1,4),GC_11,ZERO,ZERO,PL(0,67)
1315
1316
     $     ,COEFS)
1316
 
          CALL ML5_0_UPDATE_WL_4_0(WL(1,0,1,66),4,COEFS,4,4,WL(1,0,1
 
1317
          CALL ML5_0_UPDATE_WL_2_0(WL(1,0,1,66),4,COEFS,4,4,WL(1,0,1
1317
1318
     $     ,67))
1318
 
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,67),4,4,19,1,1,85,H)
 
1319
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,67),2,4,19,1,1,85,H)
1319
1320
C         Coefficient construction for loop diagram with ID 42
1320
1321
          CALL FFS2L3_1(PL(0,58),W(1,3),GC_104,MDL_MD,ZERO,PL(0,68)
1321
1322
     $     ,COEFS)
1337
1338
     $     ,71))
1338
1339
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,71),2,4,20,1,1,87,H)
1339
1340
C         Coefficient construction for loop diagram with ID 44
1340
 
          CALL FFV2L1_3(PL(0,12),W(1,3),GC_124,MDL_MW,MDL_WW,PL(0,72)
1341
 
     $     ,COEFS)
1342
 
          CALL ML5_0_UPDATE_WL_1_2(WL(1,0,1,12),4,COEFS,4,4,WL(1,0,1
 
1341
          CALL FFV2L1P0_3(PL(0,12),W(1,3),GC_124,MDL_MW,MDL_WW,PL(0
 
1342
     $     ,72),COEFS)
 
1343
          CALL ML5_0_UPDATE_WL_1_0(WL(1,0,1,12),4,COEFS,4,4,WL(1,0,1
1343
1344
     $     ,72))
1344
1345
          CALL FFV2L3_1(PL(0,72),W(1,2),GC_59,MDL_MD,ZERO,PL(0,73)
1345
1346
     $     ,COEFS)
1346
 
          CALL ML5_0_UPDATE_WL_3_1(WL(1,0,1,72),4,COEFS,4,4,WL(1,0,1
 
1347
          CALL ML5_0_UPDATE_WL_1_1(WL(1,0,1,72),4,COEFS,4,4,WL(1,0,1
1347
1348
     $     ,73))
1348
1349
          CALL FFV1L2P0_3(PL(0,73),W(1,4),GC_11,ZERO,ZERO,PL(0,74)
1349
1350
     $     ,COEFS)
1350
 
          CALL ML5_0_UPDATE_WL_4_0(WL(1,0,1,73),4,COEFS,4,4,WL(1,0,1
 
1351
          CALL ML5_0_UPDATE_WL_2_0(WL(1,0,1,73),4,COEFS,4,4,WL(1,0,1
1351
1352
     $     ,74))
1352
 
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,74),4,4,21,1,1,88,H)
 
1353
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,74),2,4,21,1,1,88,H)
1353
1354
C         Coefficient construction for loop diagram with ID 45
1354
1355
          CALL FFS1_3L1_3(PL(0,12),W(1,3),GC_13,GC_22,MDL_MW,MDL_WW
1355
1356
     $     ,PL(0,75),COEFS)
1394
1395
     $     ,COEFS)
1395
1396
          CALL ML5_0_UPDATE_WL_1_1(WL(1,0,1,83),4,COEFS,4,4,WL(1,0,1
1396
1397
     $     ,84))
1397
 
          CALL FFV2_3L2_3(PL(0,84),W(1,4),GC_68,GC_76,MDL_MZ,MDL_WZ
 
1398
          CALL FFV2_3L2P0_3(PL(0,84),W(1,4),GC_68,GC_76,MDL_MZ,MDL_WZ
1398
1399
     $     ,PL(0,85),COEFS)
1399
 
          CALL ML5_0_UPDATE_WL_2_2(WL(1,0,1,84),4,COEFS,4,4,WL(1,0,1
 
1400
          CALL ML5_0_UPDATE_WL_2_0(WL(1,0,1,84),4,COEFS,4,4,WL(1,0,1
1400
1401
     $     ,85))
1401
 
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,85),4,4,22,1,1,91,H)
 
1402
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,85),2,4,22,1,1,91,H)
1402
1403
C         Coefficient construction for loop diagram with ID 48
1403
1404
          CALL FFS2L3_2(PL(0,0),W(1,1),GC_123,MDL_MU,ZERO,PL(0,86)
1404
1405
     $     ,COEFS)
1448
1449
     $     ,COEFS)
1449
1450
          CALL ML5_0_UPDATE_WL_1_1(WL(1,0,1,95),4,COEFS,4,4,WL(1,0,1
1450
1451
     $     ,96))
1451
 
          CALL FFV2L2_3(PL(0,96),W(1,4),GC_59,MDL_MW,MDL_WW,PL(0,97)
 
1452
          CALL FFV2L2P0_3(PL(0,96),W(1,4),GC_59,MDL_MW,MDL_WW,PL(0,97)
1452
1453
     $     ,COEFS)
1453
 
          CALL ML5_0_UPDATE_WL_2_2(WL(1,0,1,96),4,COEFS,4,4,WL(1,0,1
 
1454
          CALL ML5_0_UPDATE_WL_2_0(WL(1,0,1,96),4,COEFS,4,4,WL(1,0,1
1454
1455
     $     ,97))
1455
 
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,97),4,4,24,1,1,94,H)
 
1456
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,97),2,4,24,1,1,94,H)
1456
1457
C         Coefficient construction for loop diagram with ID 51
1457
1458
          CALL FFS1_3L3_2(PL(0,0),W(1,1),GC_13,GC_22,MDL_MD,ZERO,PL(0
1458
1459
     $     ,98),COEFS)
1486
1487
     $     ,COEFS)
1487
1488
          CALL ML5_0_UPDATE_WL_1_1(WL(1,0,1,83),4,COEFS,4,4,WL(1,0,1
1488
1489
     $     ,104))
1489
 
          CALL FFV2_3L1_3(PL(0,104),W(1,3),GC_68,GC_76,MDL_MZ,MDL_WZ
 
1490
          CALL FFV2_3L1P0_3(PL(0,104),W(1,3),GC_68,GC_76,MDL_MZ,MDL_WZ
1490
1491
     $     ,PL(0,105),COEFS)
1491
 
          CALL ML5_0_UPDATE_WL_2_2(WL(1,0,1,104),4,COEFS,4,4,WL(1,0,1
 
1492
          CALL ML5_0_UPDATE_WL_2_0(WL(1,0,1,104),4,COEFS,4,4,WL(1,0,1
1492
1493
     $     ,105))
1493
 
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,105),4,4,25,1,1,97,H)
 
1494
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,105),2,4,25,1,1,97,H)
1494
1495
C         Coefficient construction for loop diagram with ID 54
1495
1496
          CALL FFV1L3_2(PL(0,87),W(1,4),GC_11,MDL_MD,ZERO,PL(0,106)
1496
1497
     $     ,COEFS)
1512
1513
     $     ,109))
1513
1514
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,109),2,1,26,1,1,99,H)
1514
1515
C         Coefficient construction for loop diagram with ID 56
1515
 
          CALL FFV2L2_3(PL(0,0),W(1,1),GC_124,MDL_MW,MDL_WW,PL(0,110)
1516
 
     $     ,COEFS)
1517
 
          CALL ML5_0_UPDATE_WL_0_2(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
 
1516
          CALL FFV2L2P0_3(PL(0,0),W(1,1),GC_124,MDL_MW,MDL_WW,PL(0
 
1517
     $     ,110),COEFS)
 
1518
          CALL ML5_0_UPDATE_WL_0_0(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
1518
1519
     $     ,110))
1519
1520
          CALL FFV2L3_1(PL(0,110),W(1,2),GC_59,MDL_MD,ZERO,PL(0,111)
1520
1521
     $     ,COEFS)
1521
 
          CALL ML5_0_UPDATE_WL_2_1(WL(1,0,1,110),4,COEFS,4,4,WL(1,0,1
 
1522
          CALL ML5_0_UPDATE_WL_0_1(WL(1,0,1,110),4,COEFS,4,4,WL(1,0,1
1522
1523
     $     ,111))
1523
1524
          CALL FFV1L2P0_3(PL(0,111),W(1,4),GC_11,ZERO,ZERO,PL(0,112)
1524
1525
     $     ,COEFS)
1525
 
          CALL ML5_0_UPDATE_WL_3_0(WL(1,0,1,111),4,COEFS,4,4,WL(1,0,1
 
1526
          CALL ML5_0_UPDATE_WL_1_0(WL(1,0,1,111),4,COEFS,4,4,WL(1,0,1
1526
1527
     $     ,112))
1527
1528
          CALL FFV1L3_1(PL(0,112),W(1,3),GC_11,MDL_MD,ZERO,PL(0,113)
1528
1529
     $     ,COEFS)
1529
 
          CALL ML5_0_UPDATE_WL_3_1(WL(1,0,1,112),4,COEFS,4,4,WL(1,0,1
 
1530
          CALL ML5_0_UPDATE_WL_1_1(WL(1,0,1,112),4,COEFS,4,4,WL(1,0,1
1530
1531
     $     ,113))
1531
 
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,113),4,4,27,1,1,100,H)
 
1532
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,113),2,4,27,1,1,100,H)
1532
1533
C         Coefficient construction for loop diagram with ID 57
1533
1534
          CALL FFS1_3L2_3(PL(0,0),W(1,1),GC_13,GC_22,MDL_MW,MDL_WW
1534
1535
     $     ,PL(0,114),COEFS)
1548
1549
     $     ,117))
1549
1550
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,117),2,4,27,1,1,101,H)
1550
1551
C         Coefficient construction for loop diagram with ID 58
1551
 
          CALL FFV2L1_3(PL(0,0),W(1,3),GC_124,MDL_MW,MDL_WW,PL(0,118)
1552
 
     $     ,COEFS)
1553
 
          CALL ML5_0_UPDATE_WL_0_2(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
 
1552
          CALL FFV2L1P0_3(PL(0,0),W(1,3),GC_124,MDL_MW,MDL_WW,PL(0
 
1553
     $     ,118),COEFS)
 
1554
          CALL ML5_0_UPDATE_WL_0_0(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
1554
1555
     $     ,118))
1555
1556
          CALL FFV2L3_2(PL(0,118),W(1,4),GC_59,MDL_MU,ZERO,PL(0,119)
1556
1557
     $     ,COEFS)
1557
 
          CALL ML5_0_UPDATE_WL_2_1(WL(1,0,1,118),4,COEFS,4,4,WL(1,0,1
 
1558
          CALL ML5_0_UPDATE_WL_0_1(WL(1,0,1,118),4,COEFS,4,4,WL(1,0,1
1558
1559
     $     ,119))
1559
1560
          CALL FFV1L1_2(PL(0,119),W(1,5),GC_11,MDL_MU,ZERO,PL(0,120)
1560
1561
     $     ,COEFS)
1561
 
          CALL ML5_0_UPDATE_WL_3_1(WL(1,0,1,119),4,COEFS,4,4,WL(1,0,1
 
1562
          CALL ML5_0_UPDATE_WL_1_1(WL(1,0,1,119),4,COEFS,4,4,WL(1,0,1
1562
1563
     $     ,120))
1563
 
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,120),4,4,28,1,1,102,H)
 
1564
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,120),2,4,28,1,1,102,H)
1564
1565
C         Coefficient construction for loop diagram with ID 59
1565
1566
          CALL FFS1_3L1_3(PL(0,0),W(1,3),GC_13,GC_22,MDL_MW,MDL_WW
1566
1567
     $     ,PL(0,121),COEFS)
1590
1591
     $     ,126))
1591
1592
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,126),2,4,4,1,1,104,H)
1592
1593
C         Coefficient construction for loop diagram with ID 61
1593
 
          CALL FFV2_5L2_3(PL(0,0),W(1,1),GC_69,GC_76,MDL_MZ,MDL_WZ
 
1594
          CALL FFV2_5L2P0_3(PL(0,0),W(1,1),GC_69,GC_76,MDL_MZ,MDL_WZ
1594
1595
     $     ,PL(0,127),COEFS)
1595
 
          CALL ML5_0_UPDATE_WL_0_2(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
 
1596
          CALL ML5_0_UPDATE_WL_0_0(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
1596
1597
     $     ,127))
1597
1598
          CALL FFV2_5L3_1(PL(0,127),W(1,2),GC_69,GC_76,MDL_MU,ZERO
1598
1599
     $     ,PL(0,128),COEFS)
1599
 
          CALL ML5_0_UPDATE_WL_2_1(WL(1,0,1,127),4,COEFS,4,4,WL(1,0,1
 
1600
          CALL ML5_0_UPDATE_WL_0_1(WL(1,0,1,127),4,COEFS,4,4,WL(1,0,1
1600
1601
     $     ,128))
1601
1602
          CALL FFV1L2_1(PL(0,128),W(1,12),GC_11,MDL_MU,ZERO,PL(0,129)
1602
1603
     $     ,COEFS)
1603
 
          CALL ML5_0_UPDATE_WL_3_1(WL(1,0,1,128),4,COEFS,4,4,WL(1,0,1
 
1604
          CALL ML5_0_UPDATE_WL_1_1(WL(1,0,1,128),4,COEFS,4,4,WL(1,0,1
1604
1605
     $     ,129))
1605
 
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,129),4,4,29,1,1,105,H)
 
1606
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,129),2,4,29,1,1,105,H)
1606
1607
C         Coefficient construction for loop diagram with ID 62
1607
1608
          CALL FFS2L2_3(PL(0,0),W(1,1),GC_123,MDL_MZ,MDL_WZ,PL(0,130)
1608
1609
     $     ,COEFS)
1632
1633
     $     ,135))
1633
1634
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,135),2,4,30,1,1,107,H)
1634
1635
C         Coefficient construction for loop diagram with ID 64
1635
 
          CALL FFV2L1_3(PL(0,0),W(1,3),GC_127,MDL_MW,MDL_WW,PL(0,136)
1636
 
     $     ,COEFS)
1637
 
          CALL ML5_0_UPDATE_WL_0_2(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
 
1636
          CALL FFV2L1P0_3(PL(0,0),W(1,3),GC_127,MDL_MW,MDL_WW,PL(0
 
1637
     $     ,136),COEFS)
 
1638
          CALL ML5_0_UPDATE_WL_0_0(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
1638
1639
     $     ,136))
1639
1640
          CALL FFV2L3_2(PL(0,136),W(1,4),GC_62,MDL_MC,ZERO,PL(0,137)
1640
1641
     $     ,COEFS)
1641
 
          CALL ML5_0_UPDATE_WL_2_1(WL(1,0,1,136),4,COEFS,4,4,WL(1,0,1
 
1642
          CALL ML5_0_UPDATE_WL_0_1(WL(1,0,1,136),4,COEFS,4,4,WL(1,0,1
1642
1643
     $     ,137))
1643
1644
          CALL FFV1L1_2(PL(0,137),W(1,5),GC_11,MDL_MC,ZERO,PL(0,138)
1644
1645
     $     ,COEFS)
1645
 
          CALL ML5_0_UPDATE_WL_3_1(WL(1,0,1,137),4,COEFS,4,4,WL(1,0,1
 
1646
          CALL ML5_0_UPDATE_WL_1_1(WL(1,0,1,137),4,COEFS,4,4,WL(1,0,1
1646
1647
     $     ,138))
1647
 
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,138),4,4,31,1,1,108,H)
 
1648
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,138),2,4,31,1,1,108,H)
1648
1649
C         Coefficient construction for loop diagram with ID 65
1649
1650
          CALL FFS1_3L1_3(PL(0,0),W(1,3),GC_14,GC_23,MDL_MW,MDL_WW
1650
1651
     $     ,PL(0,139),COEFS)
1660
1661
     $     ,141))
1661
1662
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,141),2,4,31,1,1,109,H)
1662
1663
C         Coefficient construction for loop diagram with ID 66
1663
 
          CALL FFV2L1_3(PL(0,0),W(1,3),GC_130,MDL_MW,MDL_WW,PL(0,142)
1664
 
     $     ,COEFS)
1665
 
          CALL ML5_0_UPDATE_WL_0_2(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
 
1664
          CALL FFV2L1P0_3(PL(0,0),W(1,3),GC_130,MDL_MW,MDL_WW,PL(0
 
1665
     $     ,142),COEFS)
 
1666
          CALL ML5_0_UPDATE_WL_0_0(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
1666
1667
     $     ,142))
1667
1668
          CALL FFV2L3_2(PL(0,142),W(1,4),GC_65,MDL_MT,MDL_WT,PL(0,143)
1668
1669
     $     ,COEFS)
1669
 
          CALL ML5_0_UPDATE_WL_2_1(WL(1,0,1,142),4,COEFS,4,4,WL(1,0,1
 
1670
          CALL ML5_0_UPDATE_WL_0_1(WL(1,0,1,142),4,COEFS,4,4,WL(1,0,1
1670
1671
     $     ,143))
1671
1672
          CALL FFV1L1_2(PL(0,143),W(1,5),GC_11,MDL_MT,MDL_WT,PL(0,144)
1672
1673
     $     ,COEFS)
1673
 
          CALL ML5_0_UPDATE_WL_3_1(WL(1,0,1,143),4,COEFS,4,4,WL(1,0,1
 
1674
          CALL ML5_0_UPDATE_WL_1_1(WL(1,0,1,143),4,COEFS,4,4,WL(1,0,1
1674
1675
     $     ,144))
1675
 
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,144),4,4,32,1,1,110,H)
 
1676
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,144),2,4,32,1,1,110,H)
1676
1677
C         Coefficient construction for loop diagram with ID 67
1677
1678
          CALL FFS1_3L1_3(PL(0,0),W(1,3),GC_15,GC_24,MDL_MW,MDL_WW
1678
1679
     $     ,PL(0,145),COEFS)
1702
1703
     $     ,150))
1703
1704
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,150),2,4,3,1,1,112,H)
1704
1705
C         Coefficient construction for loop diagram with ID 69
1705
 
          CALL FFV2_3L1_3(PL(0,0),W(1,3),GC_68,GC_76,MDL_MZ,MDL_WZ
 
1706
          CALL FFV2_3L1P0_3(PL(0,0),W(1,3),GC_68,GC_76,MDL_MZ,MDL_WZ
1706
1707
     $     ,PL(0,151),COEFS)
1707
 
          CALL ML5_0_UPDATE_WL_0_2(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
 
1708
          CALL ML5_0_UPDATE_WL_0_0(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
1708
1709
     $     ,151))
1709
1710
          CALL FFV2_3L3_2(PL(0,151),W(1,4),GC_68,GC_76,MDL_MD,ZERO
1710
1711
     $     ,PL(0,152),COEFS)
1711
 
          CALL ML5_0_UPDATE_WL_2_1(WL(1,0,1,151),4,COEFS,4,4,WL(1,0,1
 
1712
          CALL ML5_0_UPDATE_WL_0_1(WL(1,0,1,151),4,COEFS,4,4,WL(1,0,1
1712
1713
     $     ,152))
1713
1714
          CALL FFV1L1_2(PL(0,152),W(1,5),GC_11,MDL_MD,ZERO,PL(0,153)
1714
1715
     $     ,COEFS)
1715
 
          CALL ML5_0_UPDATE_WL_3_1(WL(1,0,1,152),4,COEFS,4,4,WL(1,0,1
 
1716
          CALL ML5_0_UPDATE_WL_1_1(WL(1,0,1,152),4,COEFS,4,4,WL(1,0,1
1716
1717
     $     ,153))
1717
 
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,153),4,4,33,1,1,113,H)
 
1718
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,153),2,4,33,1,1,113,H)
1718
1719
C         Coefficient construction for loop diagram with ID 70
1719
1720
          CALL FFS2L1_3(PL(0,0),W(1,3),GC_104,MDL_MZ,MDL_WZ,PL(0,154)
1720
1721
     $     ,COEFS)
1746
1747
C         Coefficient construction for loop diagram with ID 72
1747
1748
          CALL FFV1L2_1(PL(0,111),W(1,12),GC_11,MDL_MD,ZERO,PL(0,160)
1748
1749
     $     ,COEFS)
1749
 
          CALL ML5_0_UPDATE_WL_3_1(WL(1,0,1,111),4,COEFS,4,4,WL(1,0,1
 
1750
          CALL ML5_0_UPDATE_WL_1_1(WL(1,0,1,111),4,COEFS,4,4,WL(1,0,1
1750
1751
     $     ,160))
1751
 
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,160),4,4,35,1,1,116,H)
 
1752
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,160),2,4,35,1,1,116,H)
1752
1753
C         Coefficient construction for loop diagram with ID 73
1753
1754
          CALL FFV1L2_1(PL(0,115),W(1,12),GC_11,MDL_MD,ZERO,PL(0,161)
1754
1755
     $     ,COEFS)
1756
1757
     $     ,161))
1757
1758
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,161),2,4,35,1,1,117,H)
1758
1759
C         Coefficient construction for loop diagram with ID 74
1759
 
          CALL FFV2L2_3(PL(0,0),W(1,1),GC_125,MDL_MW,MDL_WW,PL(0,162)
1760
 
     $     ,COEFS)
1761
 
          CALL ML5_0_UPDATE_WL_0_2(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
 
1760
          CALL FFV2L2P0_3(PL(0,0),W(1,1),GC_125,MDL_MW,MDL_WW,PL(0
 
1761
     $     ,162),COEFS)
 
1762
          CALL ML5_0_UPDATE_WL_0_0(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
1762
1763
     $     ,162))
1763
1764
          CALL FFV2L3_1(PL(0,162),W(1,2),GC_60,MDL_MS,ZERO,PL(0,163)
1764
1765
     $     ,COEFS)
1765
 
          CALL ML5_0_UPDATE_WL_2_1(WL(1,0,1,162),4,COEFS,4,4,WL(1,0,1
 
1766
          CALL ML5_0_UPDATE_WL_0_1(WL(1,0,1,162),4,COEFS,4,4,WL(1,0,1
1766
1767
     $     ,163))
1767
1768
          CALL FFV1L2_1(PL(0,163),W(1,12),GC_11,MDL_MS,ZERO,PL(0,164)
1768
1769
     $     ,COEFS)
1769
 
          CALL ML5_0_UPDATE_WL_3_1(WL(1,0,1,163),4,COEFS,4,4,WL(1,0,1
 
1770
          CALL ML5_0_UPDATE_WL_1_1(WL(1,0,1,163),4,COEFS,4,4,WL(1,0,1
1770
1771
     $     ,164))
1771
 
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,164),4,4,36,1,1,118,H)
 
1772
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,164),2,4,36,1,1,118,H)
1772
1773
C         Coefficient construction for loop diagram with ID 75
1773
1774
          CALL FFS1_3L2_3(PL(0,0),W(1,1),GC_16,GC_25,MDL_MW,MDL_WW
1774
1775
     $     ,PL(0,165),COEFS)
1784
1785
     $     ,167))
1785
1786
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,167),2,4,36,1,1,119,H)
1786
1787
C         Coefficient construction for loop diagram with ID 76
1787
 
          CALL FFV2L2_3(PL(0,0),W(1,1),GC_126,MDL_MW,MDL_WW,PL(0,168)
1788
 
     $     ,COEFS)
1789
 
          CALL ML5_0_UPDATE_WL_0_2(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
 
1788
          CALL FFV2L2P0_3(PL(0,0),W(1,1),GC_126,MDL_MW,MDL_WW,PL(0
 
1789
     $     ,168),COEFS)
 
1790
          CALL ML5_0_UPDATE_WL_0_0(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
1790
1791
     $     ,168))
1791
1792
          CALL FFV2L3_1(PL(0,168),W(1,2),GC_61,MDL_MB,ZERO,PL(0,169)
1792
1793
     $     ,COEFS)
1793
 
          CALL ML5_0_UPDATE_WL_2_1(WL(1,0,1,168),4,COEFS,4,4,WL(1,0,1
 
1794
          CALL ML5_0_UPDATE_WL_0_1(WL(1,0,1,168),4,COEFS,4,4,WL(1,0,1
1794
1795
     $     ,169))
1795
1796
          CALL FFV1L2_1(PL(0,169),W(1,12),GC_11,MDL_MB,ZERO,PL(0,170)
1796
1797
     $     ,COEFS)
1797
 
          CALL ML5_0_UPDATE_WL_3_1(WL(1,0,1,169),4,COEFS,4,4,WL(1,0,1
 
1798
          CALL ML5_0_UPDATE_WL_1_1(WL(1,0,1,169),4,COEFS,4,4,WL(1,0,1
1798
1799
     $     ,170))
1799
 
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,170),4,4,37,1,1,120,H)
 
1800
          CALL ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,170),2,4,37,1,1,120,H)
1800
1801
C         Coefficient construction for loop diagram with ID 77
1801
1802
          CALL FFS1_3L2_3(PL(0,0),W(1,1),GC_19,GC_28,MDL_MW,MDL_WW
1802
1803
     $     ,PL(0,171),COEFS)
1839
1840
C      recomputes the same PS point with a different CTMode
1840
1841
 300  CONTINUE
1841
1842
 
 
1843
C     Make sure that the loop calls are performed since this is new
 
1844
C      evaluation.
 
1845
      CTCALL_REQ_SO_DONE=.FALSE.
 
1846
 
1842
1847
C     Free cache when using IREGI
1843
1848
      IF(IREGIRECY.AND.MLREDUCTIONLIB(I_LIB).EQ.3) THEN
1844
1849
        CALL IREGI_FREE_PS()
1926
1931
     $   _MU),2,I_SO,15)
1927
1932
C       CutTools call for loop numbers 23,24
1928
1933
        CALL ML5_0_LOOP_4(1,2,4,3,DCMPLX(ZERO),DCMPLX(MDL_MU)
1929
 
     $   ,DCMPLX(MDL_MW),DCMPLX(MDL_MU),4,I_SO,16)
 
1934
     $   ,DCMPLX(MDL_MW),DCMPLX(MDL_MU),2,I_SO,16)
1930
1935
C       CutTools call for loop numbers 30,31
1931
1936
        CALL ML5_0_LOOP_4(1,2,4,3,DCMPLX(MDL_MU),DCMPLX(MDL_MZ)
1932
 
     $   ,DCMPLX(MDL_MD),DCMPLX(ZERO),4,I_SO,17)
 
1937
     $   ,DCMPLX(MDL_MD),DCMPLX(ZERO),2,I_SO,17)
1933
1938
C       CutTools call for loop numbers 32
1934
1939
        CALL ML5_0_LOOP_4(1,2,4,3,DCMPLX(MDL_MU),DCMPLX(MDL_MH)
1935
1940
     $   ,DCMPLX(MDL_MD),DCMPLX(ZERO),2,I_SO,18)
1936
1941
C       CutTools call for loop numbers 34,35
1937
1942
        CALL ML5_0_LOOP_4(1,2,3,4,DCMPLX(MDL_MU),DCMPLX(MDL_MZ)
1938
 
     $   ,DCMPLX(MDL_MD),DCMPLX(ZERO),4,I_SO,19)
 
1943
     $   ,DCMPLX(MDL_MD),DCMPLX(ZERO),2,I_SO,19)
1939
1944
C       CutTools call for loop numbers 36
1940
1945
        CALL ML5_0_LOOP_4(1,2,3,4,DCMPLX(MDL_MU),DCMPLX(MDL_MH)
1941
1946
     $   ,DCMPLX(MDL_MD),DCMPLX(ZERO),2,I_SO,20)
1942
1947
C       CutTools call for loop numbers 37,38
1943
1948
        CALL ML5_0_LOOP_4(1,3,2,4,DCMPLX(MDL_MU),DCMPLX(MDL_MW)
1944
 
     $   ,DCMPLX(MDL_MD),DCMPLX(ZERO),4,I_SO,21)
 
1949
     $   ,DCMPLX(MDL_MD),DCMPLX(ZERO),2,I_SO,21)
1945
1950
C       CutTools call for loop numbers 40,41
1946
1951
        CALL ML5_0_LOOP_4(1,2,3,4,DCMPLX(MDL_MU),DCMPLX(ZERO)
1947
 
     $   ,DCMPLX(MDL_MD),DCMPLX(MDL_MZ),4,I_SO,22)
 
1952
     $   ,DCMPLX(MDL_MD),DCMPLX(MDL_MZ),2,I_SO,22)
1948
1953
C       CutTools call for loop numbers 42
1949
1954
        CALL ML5_0_LOOP_4(1,2,3,4,DCMPLX(MDL_MU),DCMPLX(ZERO)
1950
1955
     $   ,DCMPLX(MDL_MD),DCMPLX(MDL_MH),2,I_SO,23)
1951
1956
C       CutTools call for loop numbers 43,44
1952
1957
        CALL ML5_0_LOOP_4(1,3,2,4,DCMPLX(MDL_MD),DCMPLX(ZERO)
1953
 
     $   ,DCMPLX(MDL_MU),DCMPLX(MDL_MW),4,I_SO,24)
 
1958
     $   ,DCMPLX(MDL_MU),DCMPLX(MDL_MW),2,I_SO,24)
1954
1959
C       CutTools call for loop numbers 46,47
1955
1960
        CALL ML5_0_LOOP_4(1,2,4,3,DCMPLX(MDL_MU),DCMPLX(ZERO)
1956
 
     $   ,DCMPLX(MDL_MD),DCMPLX(MDL_MZ),4,I_SO,25)
 
1961
     $   ,DCMPLX(MDL_MD),DCMPLX(MDL_MZ),2,I_SO,25)
1957
1962
C       CutTools call for loop numbers 48
1958
1963
        CALL ML5_0_LOOP_4(1,2,4,3,DCMPLX(MDL_MU),DCMPLX(ZERO)
1959
1964
     $   ,DCMPLX(MDL_MD),DCMPLX(MDL_MH),2,I_SO,26)
1960
1965
C       CutTools call for loop numbers 49,50
1961
1966
        CALL ML5_0_LOOP_4(1,2,4,3,DCMPLX(MDL_MW),DCMPLX(MDL_MD)
1962
 
     $   ,DCMPLX(ZERO),DCMPLX(MDL_MD),4,I_SO,27)
 
1967
     $   ,DCMPLX(ZERO),DCMPLX(MDL_MD),2,I_SO,27)
1963
1968
C       CutTools call for loop numbers 51,52
1964
1969
        CALL ML5_0_LOOP_3(3,4,5,DCMPLX(MDL_MW),DCMPLX(MDL_MU)
1965
 
     $   ,DCMPLX(MDL_MU),4,I_SO,28)
 
1970
     $   ,DCMPLX(MDL_MU),2,I_SO,28)
1966
1971
C       CutTools call for loop numbers 54,55
1967
1972
        CALL ML5_0_LOOP_3(1,2,12,DCMPLX(MDL_MZ),DCMPLX(MDL_MU)
1968
 
     $   ,DCMPLX(MDL_MU),4,I_SO,29)
 
1973
     $   ,DCMPLX(MDL_MU),2,I_SO,29)
1969
1974
C       CutTools call for loop numbers 56
1970
1975
        CALL ML5_0_LOOP_3(1,2,12,DCMPLX(MDL_MH),DCMPLX(MDL_MU)
1971
1976
     $   ,DCMPLX(MDL_MU),2,I_SO,30)
1972
1977
C       CutTools call for loop numbers 57,58
1973
1978
        CALL ML5_0_LOOP_3(3,4,5,DCMPLX(MDL_MW),DCMPLX(MDL_MC)
1974
 
     $   ,DCMPLX(MDL_MC),4,I_SO,31)
 
1979
     $   ,DCMPLX(MDL_MC),2,I_SO,31)
1975
1980
C       CutTools call for loop numbers 59,60
1976
1981
        CALL ML5_0_LOOP_3(3,4,5,DCMPLX(MDL_MW),DCMPLX(MDL_MT)
1977
 
     $   ,DCMPLX(MDL_MT),4,I_SO,32)
 
1982
     $   ,DCMPLX(MDL_MT),2,I_SO,32)
1978
1983
C       CutTools call for loop numbers 62,63
1979
1984
        CALL ML5_0_LOOP_3(3,4,5,DCMPLX(MDL_MZ),DCMPLX(MDL_MD)
1980
 
     $   ,DCMPLX(MDL_MD),4,I_SO,33)
 
1985
     $   ,DCMPLX(MDL_MD),2,I_SO,33)
1981
1986
C       CutTools call for loop numbers 64
1982
1987
        CALL ML5_0_LOOP_3(3,4,5,DCMPLX(MDL_MH),DCMPLX(MDL_MD)
1983
1988
     $   ,DCMPLX(MDL_MD),2,I_SO,34)
1984
1989
C       CutTools call for loop numbers 65,66
1985
1990
        CALL ML5_0_LOOP_3(1,2,12,DCMPLX(MDL_MW),DCMPLX(MDL_MD)
1986
 
     $   ,DCMPLX(MDL_MD),4,I_SO,35)
 
1991
     $   ,DCMPLX(MDL_MD),2,I_SO,35)
1987
1992
C       CutTools call for loop numbers 67,68
1988
1993
        CALL ML5_0_LOOP_3(1,2,12,DCMPLX(MDL_MW),DCMPLX(MDL_MS)
1989
 
     $   ,DCMPLX(MDL_MS),4,I_SO,36)
 
1994
     $   ,DCMPLX(MDL_MS),2,I_SO,36)
1990
1995
C       CutTools call for loop numbers 69,70
1991
1996
        CALL ML5_0_LOOP_3(1,2,12,DCMPLX(MDL_MW),DCMPLX(MDL_MB)
1992
 
     $   ,DCMPLX(MDL_MB),4,I_SO,37)
 
1997
     $   ,DCMPLX(MDL_MB),2,I_SO,37)
1993
1998
C       At this point, all reductions needed for (QCD=4 QED=2), i.e.
1994
1999
C        of split order ID=2, are computed.
1995
2000
        IF(FILTER_SO.AND.SQSO_TARGET.EQ.2) GOTO 5000