← Back to branch summary

~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/long_ML_SMQCD_default/gg_wmtbx/loop_matrix.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:
tests/input_files/IOTestsComparison/long_ML_SMQCD_default/gg_wmtbx/loop_matrix.f
47
47
C     ----------
48
48
      IF(ABS(REFERENCE_VALUE).EQ.0.0D0) THEN
49
49
        ML5_0_ISZERO=.FALSE.
50
 
        WRITE(*,*) '##E02 ERRROR Reference value for comparison i'
51
 
     $   //'s zero.'
 
50
        WRITE(*,*) '##E02 ERRROR Reference value for comparison is'
 
51
     $   //' zero.'
52
52
        STOP
53
53
      ELSE
54
54
        ML5_0_ISZERO=((ABS(TOTEST)/ABS(REFERENCE_VALUE)).LT.ZEROTHRES)
55
55
      ENDIF
56
56
      IF(AMPLN.NE.-1) THEN
57
57
        IF((.NOT.ML5_0_ISZERO).AND.(.NOT.S(AMPLN))) THEN
58
 
          WRITE(*,*) '##W01 WARNING Contribution ',AMPLN,' is detecte'
59
 
     $     //'d as contributing with CR=',(ABS(TOTEST)/ABS(REFERENCE_VA
60
 
     $     LUE)),' but is unstable.'
 
58
          WRITE(*,*) '##W01 WARNING Contribution ',AMPLN,' is detected'
 
59
     $     //' as contributing with CR=',(ABS(TOTEST)
 
60
     $     /ABS(REFERENCE_VALUE)),' but is unstable.'
61
61
        ENDIF
62
62
      ENDIF
63
63
 
73
73
C     and helicities for the point in phase space P(0:3,NEXTERNAL)
74
74
C     and external lines W(0:6,NEXTERNAL)
75
75
C     
76
 
C     Process: g g > w- t b~ QED=1 QCD=2 [ virt = QCD ]
 
76
C     Process: g g > w- t b~ QED<=1 QCD<=2 [ virt = QCD ]
77
77
C     
78
78
      IMPLICIT NONE
79
79
C     
167
167
      DATA NEXTREF/ZERO/
168
168
      INTEGER NPSPOINTS
169
169
      DATA NPSPOINTS/0/
 
170
      LOGICAL FOUND_VALID_REDUCTION_METHOD
 
171
      DATA FOUND_VALID_REDUCTION_METHOD/.FALSE./
170
172
 
171
173
      REAL*8 ACC
172
174
      REAL*8 DP_RES(3,MAXSTABILITYLENGTH)
205
207
C     
206
208
C     GLOBAL VARIABLES
207
209
C     
 
210
      INCLUDE 'process_info.inc'
208
211
      INCLUDE 'coupl.inc'
209
212
      INCLUDE 'mp_coupl.inc'
210
213
      INCLUDE 'MadLoopParams.inc'
307
310
      LOGICAL ML_INIT
308
311
      COMMON/ML_INIT/ML_INIT
309
312
 
 
313
C     This variable controls the *local* initialization of this
 
314
C      particular SubProcess.
 
315
C     For example, the reading of the filters must be done
 
316
C      independently by each SubProcess.
 
317
      LOGICAL LOCAL_ML_INIT
 
318
      DATA LOCAL_ML_INIT/.TRUE./
 
319
 
 
320
C     Variables related to turning off the Lorentz rotation test when
 
321
C      spin-2 particles are external
 
322
      LOGICAL WARNED_LORENTZ_STAB_TEST_OFF
 
323
      DATA WARNED_LORENTZ_STAB_TEST_OFF/.FALSE./
 
324
      INTEGER NROTATIONS_DP_BU,NROTATIONS_QP_BU
 
325
 
310
326
C     ----------
311
327
C     BEGIN CODE
312
328
C     ----------
321
337
          DOUBLECHECKHELICITYFILTER = .FALSE.
322
338
        ENDIF
323
339
        ML_INIT = .FALSE.
324
 
      ENDIF
325
 
 
326
 
C     Setup the file paths
327
 
      CALL JOINPATH(MLPATH,PARAMFNAME,PARAMFN)
328
 
      CALL JOINPATH(MLPATH,PROC_PREFIX,TMP)
329
 
      CALL JOINPATH(TMP,HELCONFIGFNAME,HELCONFIGFN)
330
 
      CALL JOINPATH(TMP,LOOPFILTERFNAME,LOOPFILTERFN)
331
 
      CALL JOINPATH(TMP,COLORNUMFNAME,COLORNUMFN)
332
 
      CALL JOINPATH(TMP,COLORDENOMFNAME,COLORDENOMFN)
333
 
      CALL JOINPATH(TMP,HELFILTERFNAME,HELFILTERFN)
334
 
 
335
 
      OPEN(1, FILE=COLORNUMFN, ERR=104, STATUS='OLD',          
336
 
     $  ACTION='READ')
337
 
      DO I=1,NLOOPAMPS
338
 
        READ(1,*,END=105) (CF_N(I,J),J=1,NBORNAMPS)
339
 
      ENDDO
340
 
      GOTO 105
341
 
 104  CONTINUE
342
 
      STOP 'Color factors could not be initialized from file ML5_0_Col'
343
 
     $ //'orNumFactors.dat. File not found'
344
 
 105  CONTINUE
345
 
      CLOSE(1)
346
 
      OPEN(1, FILE=COLORDENOMFN, ERR=106, STATUS='OLD',          
347
 
     $  ACTION='READ')
348
 
      DO I=1,NLOOPAMPS
349
 
        READ(1,*,END=107) (CF_D(I,J),J=1,NBORNAMPS)
350
 
      ENDDO
351
 
      GOTO 107
352
 
 106  CONTINUE
353
 
      STOP 'Color factors could not be initialized from file ML5_0_Col'
354
 
     $ //'orDenomFactors.dat. File not found'
355
 
 107  CONTINUE
356
 
      CLOSE(1)
357
 
      OPEN(1, FILE=HELCONFIGFN, ERR=108, STATUS='OLD',                
358
 
     $   ACTION='READ')
359
 
      DO H=1,NCOMB
360
 
        READ(1,*,END=109) (HELC(I,H),I=1,NEXTERNAL)
361
 
      ENDDO
362
 
      GOTO 109
363
 
 108  CONTINUE
364
 
      STOP 'Color helictiy configurations could not be initialize'
365
 
     $ //'d from file ML5_0_HelConfigs.dat. File not found'
366
 
 109  CONTINUE
367
 
      CLOSE(1)
368
 
      IF(BOOTANDSTOP) THEN
369
 
        WRITE(*,*) '##Stopped by user request.'
370
 
        STOP
 
340
C       For now only CutTools is interfaced in the default mode.
 
341
C        Samurai could follow.
 
342
        DO I=1,SIZE(MLREDUCTIONLIB)
 
343
          IF (MLREDUCTIONLIB(I).EQ.1) THEN
 
344
            FOUND_VALID_REDUCTION_METHOD = .TRUE.
 
345
          ENDIF
 
346
        ENDDO
 
347
        IF (.NOT.FOUND_VALID_REDUCTION_METHOD) THEN
 
348
          WRITE(*,*) 'ERROR:: For now, only CutTools is interfaced to'
 
349
     $     //' MadLoop in the non-optimized output.'
 
350
          WRITE(*,*) 'ERROR:: Make sure to include 1 in the parameter'
 
351
     $     //' MLReductionLib of the card MadLoopParams.dat'
 
352
          STOP 1
 
353
        ENDIF
 
354
      ENDIF
 
355
      IF (LOCAL_ML_INIT) THEN
 
356
C       Setup the file paths
 
357
        CALL JOINPATH(MLPATH,PARAMFNAME,PARAMFN)
 
358
        CALL JOINPATH(MLPATH,PROC_PREFIX,TMP)
 
359
        CALL JOINPATH(TMP,HELCONFIGFNAME,HELCONFIGFN)
 
360
        CALL JOINPATH(TMP,LOOPFILTERFNAME,LOOPFILTERFN)
 
361
        CALL JOINPATH(TMP,COLORNUMFNAME,COLORNUMFN)
 
362
        CALL JOINPATH(TMP,COLORDENOMFNAME,COLORDENOMFN)
 
363
        CALL JOINPATH(TMP,HELFILTERFNAME,HELFILTERFN)
 
364
 
 
365
C       Make sure that the loop filter is disabled when there is
 
366
C        spin-2 particles for 2>1 or 1>2 processes
 
367
        IF(MAX_SPIN_EXTERNAL_PARTICLE.GT.3.AND.(NEXTERNAL.LE.3.AND.HELI
 
368
     $CITYFILTERLEVEL.NE.0)) THEN
 
369
          WRITE(*,*) '##INFO: Helicity filter deactivated for 2>1'
 
370
     $     //' processes involving spin 2 particles.'
 
371
          HELICITYFILTERLEVEL = 0
 
372
C         We write a dummy filter for structural reasons here
 
373
          OPEN(1, FILE=HELFILTERFN, ERR=6116, STATUS='NEW'
 
374
     $     ,ACTION='WRITE')
 
375
          DO I=1,NCOMB
 
376
            WRITE(1,*) 'T'
 
377
          ENDDO
 
378
 6116     CONTINUE
 
379
          CLOSE(1)
 
380
        ENDIF
 
381
 
 
382
        OPEN(1, FILE=COLORNUMFN, ERR=104, STATUS='OLD',          
 
383
     $    ACTION='READ')
 
384
        DO I=1,NLOOPAMPS
 
385
          READ(1,*,END=105) (CF_N(I,J),J=1,NBORNAMPS)
 
386
        ENDDO
 
387
        GOTO 105
 
388
 104    CONTINUE
 
389
        STOP 'Color factors could not be initialized from file'
 
390
     $   //' ML5_0_ColorNumFactors.dat. File not found'
 
391
 105    CONTINUE
 
392
        CLOSE(1)
 
393
        OPEN(1, FILE=COLORDENOMFN, ERR=106, STATUS='OLD',          
 
394
     $    ACTION='READ')
 
395
        DO I=1,NLOOPAMPS
 
396
          READ(1,*,END=107) (CF_D(I,J),J=1,NBORNAMPS)
 
397
        ENDDO
 
398
        GOTO 107
 
399
 106    CONTINUE
 
400
        STOP 'Color factors could not be initialized from file'
 
401
     $   //' ML5_0_ColorDenomFactors.dat. File not found'
 
402
 107    CONTINUE
 
403
        CLOSE(1)
 
404
        OPEN(1, FILE=HELCONFIGFN, ERR=108, STATUS='OLD',              
 
405
     $       ACTION='READ')
 
406
        DO H=1,NCOMB
 
407
          READ(1,*,END=109) (HELC(I,H),I=1,NEXTERNAL)
 
408
        ENDDO
 
409
        GOTO 109
 
410
 108    CONTINUE
 
411
        STOP 'Color helictiy configurations could not be initialized'
 
412
     $   //' from file ML5_0_HelConfigs.dat. File not found'
 
413
 109    CONTINUE
 
414
        CLOSE(1)
 
415
        IF(BOOTANDSTOP) THEN
 
416
          WRITE(*,*) '##Stopped by user request.'
 
417
          STOP
 
418
        ENDIF
 
419
        LOCAL_ML_INIT = .FALSE.
 
420
      ENDIF
 
421
 
 
422
C     Make sure that lorentz rotation tests are not used if there is
 
423
C      external loop wavefunction of spin 2 and that one specific
 
424
C      helicity is asked
 
425
      NROTATIONS_DP_BU = NROTATIONS_DP
 
426
      NROTATIONS_QP_BU = NROTATIONS_QP
 
427
      IF(MAX_SPIN_EXTERNAL_PARTICLE.GT.3.AND.USERHEL.NE.-1) THEN
 
428
        IF(.NOT.WARNED_LORENTZ_STAB_TEST_OFF) THEN
 
429
          WRITE(*,*) '##WARNING: Evaluation of a specific helicity was'
 
430
     $     //' asked for this PS point, and there is a spin-2 (or'
 
431
     $     //' higher) particle in the external states.'
 
432
          WRITE(*,*) '##WARNING: As a result, MadLoop disabled the'
 
433
     $     //' Lorentz rotation test for this phase-space point only.'
 
434
          WRITE(*,*) '##WARNING: Further warning of that type'
 
435
     $     //' suppressed.'
 
436
          WARNED_LORENTZ_STAB_TEST_OFF = .FALSE.
 
437
        ENDIF
 
438
        NROTATIONS_QP=0
 
439
        NROTATIONS_DP=0
371
440
      ENDIF
372
441
 
373
442
      IF(NTRY.EQ.0) THEN
374
443
        CALL ML5_0_SET_N_EVALS(N_DP_EVAL,N_QP_EVAL)
375
 
        HELDOUBLECHECKED=(.NOT.DOUBLECHECKHELICITYFILTER).OR.(HELICITYF
376
 
     $   ILTERLEVEL.EQ.0)
 
444
        HELDOUBLECHECKED=(.NOT.DOUBLECHECKHELICITYFILTER)
 
445
     $   .OR.(HELICITYFILTERLEVEL.EQ.0)
377
446
        DO J=1,NCOMB
378
447
          DO I=1,NCTAMPS
379
448
            GOODAMP(I,J)=.TRUE.
452
521
        DONEHELDOUBLECHECK=.FALSE.
453
522
      ENDIF
454
523
 
455
 
      CHECKPHASE=(NTRY.LE.CHECKCYCLE).AND.(((.NOT.FOUNDLOOPFILTER
456
 
     $ ).AND.USELOOPFILTER).OR.(.NOT.FOUNDHELFILTER))
 
524
      CHECKPHASE=(NTRY.LE.CHECKCYCLE).AND.(((.NOT.FOUNDLOOPFILTER)
 
525
     $ .AND.USELOOPFILTER).OR.(.NOT.FOUNDHELFILTER))
457
526
 
458
527
      IF (WRITEOUTFILTERS) THEN
459
528
        IF ((.NOT. CHECKPHASE).AND.(.NOT.FOUNDHELFILTER)) THEN
460
 
          OPEN(1, FILE=HELFILTERFN, ERR=110, STATUS='NEW',ACTION='WRIT'
461
 
     $     //'E')
 
529
          OPEN(1, FILE=HELFILTERFN, ERR=110, STATUS='NEW'
 
530
     $     ,ACTION='WRITE')
462
531
          WRITE(1,*) (GOODHEL(I),I=1,NCOMB)
463
532
 110      CONTINUE
464
533
          CLOSE(1)
465
534
          FOUNDHELFILTER=.TRUE.
466
535
        ENDIF
467
536
 
468
 
        IF ((.NOT. CHECKPHASE).AND.(.NOT.FOUNDLOOPFILTER).AND.USELOOPFI
469
 
     $   LTER) THEN
470
 
          OPEN(1, FILE=LOOPFILTERFN, ERR=111, STATUS='NEW',ACTION='WRI'
471
 
     $     //'TE')
 
537
        IF ((.NOT. CHECKPHASE).AND.(.NOT.FOUNDLOOPFILTER)
 
538
     $   .AND.USELOOPFILTER) THEN
 
539
          OPEN(1, FILE=LOOPFILTERFN, ERR=111, STATUS='NEW'
 
540
     $     ,ACTION='WRITE')
472
541
          DO J=1,NCOMB
473
542
            WRITE(1,*) (GOODAMP(I,J),I=NCTAMPS+1,NLOOPAMPS)
474
543
          ENDDO
522
591
      ENDDO
523
592
 
524
593
      IF (IMPROVEPSPOINT.GE.0) THEN
525
 
C       Make the input PS more precise (exact onshell and energy-moment
526
 
C       um conservation)
 
594
C       Make the input PS more precise (exact onshell and
 
595
C        energy-momentum conservation)
527
596
        CALL ML5_0_IMPROVE_PS_POINT_PRECISION(PS)
528
597
      ENDIF
529
598
 
564
633
      VALIDH=-1
565
634
      DO H=1,NCOMB
566
635
        IF ((HELPICKED.EQ.H).OR.((HELPICKED.EQ.-1).AND.(CHECKPHASE.OR.(
567
 
     $   .NOT.HELDOUBLECHECKED).OR.GOODHEL(H)))) THEN
 
636
     $.NOT.HELDOUBLECHECKED).OR.GOODHEL(H)))) THEN
568
637
          IF (VALIDH.EQ.-1) VALIDH=H
569
638
          DO I=1,NEXTERNAL
570
639
            NHEL(I)=HELC(I,H)
627
696
          CALL R2_GG_1_0(W(1,6,H),W(1,14,H),R2_GGQ,AMPL(1,7))
628
697
          CALL R2_GG_1_0(W(1,6,H),W(1,14,H),R2_GGQ,AMPL(1,8))
629
698
C         Counter-term amplitude(s) for loop diagram number 11
630
 
          CALL VVV1_0(W(1,1,H),W(1,2,H),W(1,13,H),UV_3GB_1EPS,AMPL(2
631
 
     $     ,9))
 
699
          CALL VVV1_0(W(1,1,H),W(1,2,H),W(1,13,H),UV_3GB_1EPS,AMPL(2,9)
 
700
     $     )
632
701
          CALL VVV1_0(W(1,1,H),W(1,2,H),W(1,13,H),UV_3GB_1EPS,AMPL(2
633
702
     $     ,10))
634
703
          CALL VVV1_0(W(1,1,H),W(1,2,H),W(1,13,H),UV_3GB_1EPS,AMPL(2
1150
1219
     $ ,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),GC_5,MP__GC_5,GC_5
1151
1220
     $ ,MP__GC_5,2,1,4,254,AMPL(1,254),S(254))
1152
1221
C     Loop amplitude for loop diagram with ID 11
1153
 
      CALL ML5_0_LOOP_3_3(2,1,2,13,DCMPLX(ZERO),CMPLX(MP__ZERO
1154
 
     $ ,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(ZERO)
 
1222
      CALL ML5_0_LOOP_3_3(2,1,2,13,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
 
1223
     $ ,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(ZERO)
1155
1224
     $ ,CMPLX(MP__ZERO,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5
1156
1225
     $ ,MP__GC_5,3,1,4,255,AMPL(1,255),S(255))
1157
1226
C     Loop amplitude for loop diagram with ID 12
1158
 
      CALL ML5_0_LOOP_3_3(2,1,2,14,DCMPLX(ZERO),CMPLX(MP__ZERO
1159
 
     $ ,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(ZERO)
 
1227
      CALL ML5_0_LOOP_3_3(2,1,2,14,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
 
1228
     $ ,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(ZERO)
1160
1229
     $ ,CMPLX(MP__ZERO,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5
1161
1230
     $ ,MP__GC_5,3,1,4,256,AMPL(1,256),S(256))
1162
1231
C     Loop amplitude for loop diagram with ID 13
1163
 
      CALL ML5_0_LOOP_3_3(3,1,2,13,DCMPLX(ZERO),CMPLX(MP__ZERO
1164
 
     $ ,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(ZERO)
 
1232
      CALL ML5_0_LOOP_3_3(3,1,2,13,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
 
1233
     $ ,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(ZERO)
1165
1234
     $ ,CMPLX(MP__ZERO,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5
1166
1235
     $ ,MP__GC_5,3,1,4,257,AMPL(1,257),S(257))
1167
1236
C     Loop amplitude for loop diagram with ID 14
1168
 
      CALL ML5_0_LOOP_3_3(3,1,2,14,DCMPLX(ZERO),CMPLX(MP__ZERO
1169
 
     $ ,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(ZERO)
 
1237
      CALL ML5_0_LOOP_3_3(3,1,2,14,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
 
1238
     $ ,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(ZERO)
1170
1239
     $ ,CMPLX(MP__ZERO,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5
1171
1240
     $ ,MP__GC_5,3,1,4,258,AMPL(1,258),S(258))
1172
1241
C     Loop amplitude for loop diagram with ID 15
1184
1253
     $ ,MP__GC_5,2,1,1,261,AMPL(1,261),S(261))
1185
1254
C     Loop amplitude for loop diagram with ID 18
1186
1255
      CALL ML5_0_LOOP_3_3(6,5,6,7,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
1187
 
     $ ,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__M
1188
 
     $ DL_MB,KIND=16),GC_5,MP__GC_5,GC_4,MP__GC_4,GC_5,MP__GC_5,2,1,1
1189
 
     $ ,262,AMPL(1,262),S(262))
 
1256
     $ ,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MB)
 
1257
     $ ,CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_4,MP__GC_4,GC_5
 
1258
     $ ,MP__GC_5,2,1,1,262,AMPL(1,262),S(262))
1190
1259
C     Loop amplitude for loop diagram with ID 19
1191
1260
      CALL ML5_0_LOOP_2_2(1,6,14,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB
1192
1261
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),GC_5
1220
1289
     $ ,KIND=16),GC_47,MP__GC_47,GC_5,MP__GC_5,GC_4,MP__GC_4,GC_5
1221
1290
     $ ,MP__GC_5,3,1,1,268,AMPL(1,268),S(268))
1222
1291
C     Loop amplitude for loop diagram with ID 25
1223
 
      CALL ML5_0_LOOP_2_2(12,10,17,DCMPLX(ZERO),CMPLX(MP__ZERO
1224
 
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),GC_5
1225
 
     $ ,MP__GC_5,GC_5,MP__GC_5,1,1,1,269,AMPL(1,269),S(269))
 
1292
      CALL ML5_0_LOOP_2_2(12,10,17,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
 
1293
     $ ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_5
 
1294
     $ ,MP__GC_5,1,1,1,269,AMPL(1,269),S(269))
1226
1295
C     Loop amplitude for loop diagram with ID 26
1227
1296
      CALL ML5_0_LOOP_3_3(7,3,9,10,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT
1228
1297
     $ ,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MB)
1246
1315
     $ ,MP__GC_5,2,1,1,273,AMPL(1,273),S(273))
1247
1316
C     Loop amplitude for loop diagram with ID 30
1248
1317
      CALL ML5_0_LOOP_4_4(15,2,5,3,9,DCMPLX(ZERO),CMPLX(MP__ZERO
1249
 
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(MDL_MT
1250
 
     $ ),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO
1251
 
     $ ,KIND=16),GC_4,MP__GC_4,GC_5,MP__GC_5,GC_47,MP__GC_47,GC_5
1252
 
     $ ,MP__GC_5,3,1,1,274,AMPL(1,274),S(274))
 
1318
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16)
 
1319
     $ ,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(ZERO)
 
1320
     $ ,CMPLX(MP__ZERO,KIND=16),GC_4,MP__GC_4,GC_5,MP__GC_5,GC_47
 
1321
     $ ,MP__GC_47,GC_5,MP__GC_5,3,1,1,274,AMPL(1,274),S(274))
1253
1322
C     Loop amplitude for loop diagram with ID 31
1254
1323
      CALL ML5_0_LOOP_4_4(16,2,3,5,9,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT
1255
1324
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(ZERO)
1262
1331
     $ ,CMPLX(MP__ZERO,KIND=16),GC_4,MP__GC_4,GC_5,MP__GC_5,GC_5
1263
1332
     $ ,MP__GC_5,2,1,1,276,AMPL(1,276),S(276))
1264
1333
C     Loop amplitude for loop diagram with ID 33
1265
 
      CALL ML5_0_LOOP_2_2(12,11,19,DCMPLX(ZERO),CMPLX(MP__ZERO
1266
 
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),GC_5
1267
 
     $ ,MP__GC_5,GC_5,MP__GC_5,1,1,1,277,AMPL(1,277),S(277))
 
1334
      CALL ML5_0_LOOP_2_2(12,11,19,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
 
1335
     $ ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_5
 
1336
     $ ,MP__GC_5,1,1,1,277,AMPL(1,277),S(277))
1268
1337
C     Loop amplitude for loop diagram with ID 34
1269
1338
      CALL ML5_0_LOOP_3_3(7,3,12,11,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT
1270
1339
     $ ,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MB)
1286
1355
     $ ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_5
1287
1356
     $ ,MP__GC_5,1,1,1,281,AMPL(1,281),S(281))
1288
1357
C     Loop amplitude for loop diagram with ID 38
1289
 
      CALL ML5_0_LOOP_2_2(12,11,21,DCMPLX(ZERO),CMPLX(MP__ZERO
1290
 
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),GC_5
1291
 
     $ ,MP__GC_5,GC_5,MP__GC_5,1,1,1,282,AMPL(1,282),S(282))
 
1358
      CALL ML5_0_LOOP_2_2(12,11,21,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
 
1359
     $ ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_5
 
1360
     $ ,MP__GC_5,1,1,1,282,AMPL(1,282),S(282))
1292
1361
C     Loop amplitude for loop diagram with ID 39
1293
1362
      CALL ML5_0_LOOP_3_3(19,2,7,11,DCMPLX(ZERO),CMPLX(MP__ZERO
1294
1363
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(ZERO)
1307
1376
     $ ,MP__GC_5,3,1,1,285,AMPL(1,285),S(285))
1308
1377
C     Loop amplitude for loop diagram with ID 42
1309
1378
      CALL ML5_0_LOOP_4_4(20,2,4,3,11,DCMPLX(ZERO),CMPLX(MP__ZERO
1310
 
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(MDL_MB
1311
 
     $ ),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO
1312
 
     $ ,KIND=16),GC_4,MP__GC_4,GC_5,MP__GC_5,GC_47,MP__GC_47,GC_5
1313
 
     $ ,MP__GC_5,3,1,1,286,AMPL(1,286),S(286))
 
1379
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16)
 
1380
     $ ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(ZERO)
 
1381
     $ ,CMPLX(MP__ZERO,KIND=16),GC_4,MP__GC_4,GC_5,MP__GC_5,GC_47
 
1382
     $ ,MP__GC_47,GC_5,MP__GC_5,3,1,1,286,AMPL(1,286),S(286))
1314
1383
C     Loop amplitude for loop diagram with ID 43
1315
1384
      CALL ML5_0_LOOP_3_3(13,1,5,19,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB
1316
1385
     $ ,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MB)
1335
1404
     $ ,MP__GC_5,3,1,1,290,AMPL(1,290),S(290))
1336
1405
C     Loop amplitude for loop diagram with ID 47
1337
1406
      CALL ML5_0_LOOP_5_5(21,1,2,3,4,5,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB
1338
 
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(MDL_MT
1339
 
     $ ),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO
1340
 
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),GC_5
1341
 
     $ ,MP__GC_5,GC_5,MP__GC_5,GC_47,MP__GC_47,GC_5,MP__GC_5,GC_5
1342
 
     $ ,MP__GC_5,4,1,1,291,AMPL(1,291),S(291))
 
1407
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16)
 
1408
     $ ,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(ZERO)
 
1409
     $ ,CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB
 
1410
     $ ,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5,GC_47,MP__GC_47,GC_5
 
1411
     $ ,MP__GC_5,GC_5,MP__GC_5,4,1,1,291,AMPL(1,291),S(291))
1343
1412
C     Loop amplitude for loop diagram with ID 48
1344
1413
      CALL ML5_0_LOOP_3_3(2,1,2,13,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB
1345
 
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(MDL_MB
1346
 
     $ ),CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5
1347
 
     $ ,MP__GC_5,3,1,1,292,AMPL(1,292),S(292))
 
1414
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16)
 
1415
     $ ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_5
 
1416
     $ ,MP__GC_5,GC_5,MP__GC_5,3,1,1,292,AMPL(1,292),S(292))
1348
1417
C     Loop amplitude for loop diagram with ID 49
1349
1418
      CALL ML5_0_LOOP_4_4(22,1,2,7,5,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB
1350
1419
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(ZERO)
1353
1422
     $ ,MP__GC_5,3,1,1,293,AMPL(1,293),S(293))
1354
1423
C     Loop amplitude for loop diagram with ID 50
1355
1424
      CALL ML5_0_LOOP_3_3(2,1,2,14,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB
1356
 
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(MDL_MB
1357
 
     $ ),CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5
1358
 
     $ ,MP__GC_5,3,1,1,294,AMPL(1,294),S(294))
 
1425
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16)
 
1426
     $ ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_5
 
1427
     $ ,MP__GC_5,GC_5,MP__GC_5,3,1,1,294,AMPL(1,294),S(294))
1359
1428
C     Loop amplitude for loop diagram with ID 51
1360
1429
      CALL ML5_0_LOOP_3_3(3,1,2,13,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB
1361
 
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(MDL_MB
1362
 
     $ ),CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5
1363
 
     $ ,MP__GC_5,3,1,1,295,AMPL(1,295),S(295))
 
1430
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16)
 
1431
     $ ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_5
 
1432
     $ ,MP__GC_5,GC_5,MP__GC_5,3,1,1,295,AMPL(1,295),S(295))
1364
1433
C     Loop amplitude for loop diagram with ID 52
1365
1434
      CALL ML5_0_LOOP_4_4(23,1,2,5,7,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB
1366
1435
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(ZERO)
1380
1449
     $ ,MP__GC_5,3,1,1,298,AMPL(1,298),S(298))
1381
1450
C     Loop amplitude for loop diagram with ID 55
1382
1451
      CALL ML5_0_LOOP_3_3(3,1,2,14,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB
1383
 
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(MDL_MB
1384
 
     $ ),CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5
1385
 
     $ ,MP__GC_5,3,1,1,299,AMPL(1,299),S(299))
 
1452
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16)
 
1453
     $ ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_5
 
1454
     $ ,MP__GC_5,GC_5,MP__GC_5,3,1,1,299,AMPL(1,299),S(299))
1386
1455
C     Loop amplitude for loop diagram with ID 56
1387
1456
      CALL ML5_0_LOOP_5_5(25,1,2,5,4,3,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB
1388
1457
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(ZERO)
1392
1461
     $ ,MP__GC_47,4,1,1,300,AMPL(1,300),S(300))
1393
1462
C     Loop amplitude for loop diagram with ID 57
1394
1463
      CALL ML5_0_LOOP_5_5(26,1,3,2,4,5,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB
1395
 
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(MDL_MT
1396
 
     $ ),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO
1397
 
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),GC_5
1398
 
     $ ,MP__GC_5,GC_47,MP__GC_47,GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5
1399
 
     $ ,MP__GC_5,4,1,1,301,AMPL(1,301),S(301))
 
1464
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16)
 
1465
     $ ,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(ZERO)
 
1466
     $ ,CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB
 
1467
     $ ,KIND=16),GC_5,MP__GC_5,GC_47,MP__GC_47,GC_5,MP__GC_5,GC_5
 
1468
     $ ,MP__GC_5,GC_5,MP__GC_5,4,1,1,301,AMPL(1,301),S(301))
1400
1469
C     Loop amplitude for loop diagram with ID 58
1401
1470
      CALL ML5_0_LOOP_5_5(27,1,3,4,2,5,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB
1402
1471
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(ZERO)
1417
1486
     $ ,MP__GC_5,3,1,1,304,AMPL(1,304),S(304))
1418
1487
C     Loop amplitude for loop diagram with ID 61
1419
1488
      CALL ML5_0_LOOP_4_4(15,1,5,3,12,DCMPLX(ZERO),CMPLX(MP__ZERO
1420
 
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(MDL_MT
1421
 
     $ ),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO
1422
 
     $ ,KIND=16),GC_4,MP__GC_4,GC_5,MP__GC_5,GC_47,MP__GC_47,GC_5
1423
 
     $ ,MP__GC_5,3,1,1,305,AMPL(1,305),S(305))
 
1489
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16)
 
1490
     $ ,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(ZERO)
 
1491
     $ ,CMPLX(MP__ZERO,KIND=16),GC_4,MP__GC_4,GC_5,MP__GC_5,GC_47
 
1492
     $ ,MP__GC_47,GC_5,MP__GC_5,3,1,1,305,AMPL(1,305),S(305))
1424
1493
C     Loop amplitude for loop diagram with ID 62
1425
1494
      CALL ML5_0_LOOP_3_3(17,1,5,19,DCMPLX(ZERO),CMPLX(MP__ZERO
1426
1495
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(ZERO)
1431
1500
     $ ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_5
1432
1501
     $ ,MP__GC_5,1,1,1,307,AMPL(1,307),S(307))
1433
1502
C     Loop amplitude for loop diagram with ID 64
1434
 
      CALL ML5_0_LOOP_2_2(12,10,24,DCMPLX(ZERO),CMPLX(MP__ZERO
1435
 
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),GC_5
1436
 
     $ ,MP__GC_5,GC_5,MP__GC_5,1,1,1,308,AMPL(1,308),S(308))
 
1503
      CALL ML5_0_LOOP_2_2(12,10,24,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
 
1504
     $ ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_5
 
1505
     $ ,MP__GC_5,1,1,1,308,AMPL(1,308),S(308))
1437
1506
C     Loop amplitude for loop diagram with ID 65
1438
1507
      CALL ML5_0_LOOP_3_3(19,1,7,10,DCMPLX(ZERO),CMPLX(MP__ZERO
1439
1508
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(ZERO)
1452
1521
     $ ,MP__GC_5,3,1,1,311,AMPL(1,311),S(311))
1453
1522
C     Loop amplitude for loop diagram with ID 68
1454
1523
      CALL ML5_0_LOOP_4_4(20,1,4,3,10,DCMPLX(ZERO),CMPLX(MP__ZERO
1455
 
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(MDL_MB
1456
 
     $ ),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO
1457
 
     $ ,KIND=16),GC_4,MP__GC_4,GC_5,MP__GC_5,GC_47,MP__GC_47,GC_5
1458
 
     $ ,MP__GC_5,3,1,1,312,AMPL(1,312),S(312))
 
1524
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16)
 
1525
     $ ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(ZERO)
 
1526
     $ ,CMPLX(MP__ZERO,KIND=16),GC_4,MP__GC_4,GC_5,MP__GC_5,GC_47
 
1527
     $ ,MP__GC_47,GC_5,MP__GC_5,3,1,1,312,AMPL(1,312),S(312))
1459
1528
C     Loop amplitude for loop diagram with ID 69
1460
1529
      CALL ML5_0_LOOP_3_3(13,2,5,24,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB
1461
1530
     $ ,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MB)
1463
1532
     $ ,MP__GC_5,2,1,1,313,AMPL(1,313),S(313))
1464
1533
C     Loop amplitude for loop diagram with ID 70
1465
1534
      CALL ML5_0_LOOP_4_4(28,1,5,2,7,DCMPLX(ZERO),CMPLX(MP__ZERO
1466
 
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(MDL_MB
1467
 
     $ ),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO
1468
 
     $ ,KIND=16),GC_4,MP__GC_4,GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5
1469
 
     $ ,MP__GC_5,3,1,1,314,AMPL(1,314),S(314))
 
1535
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16)
 
1536
     $ ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(ZERO)
 
1537
     $ ,CMPLX(MP__ZERO,KIND=16),GC_4,MP__GC_4,GC_5,MP__GC_5,GC_5
 
1538
     $ ,MP__GC_5,GC_5,MP__GC_5,3,1,1,314,AMPL(1,314),S(314))
1470
1539
C     Loop amplitude for loop diagram with ID 71
1471
1540
      CALL ML5_0_LOOP_5_5(29,1,3,2,5,4,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT
1472
 
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(MDL_MB
1473
 
     $ ),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO
1474
 
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),GC_5
1475
 
     $ ,MP__GC_5,GC_47,MP__GC_47,GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5
1476
 
     $ ,MP__GC_5,4,1,1,315,AMPL(1,315),S(315))
 
1541
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16)
 
1542
     $ ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(ZERO)
 
1543
     $ ,CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT
 
1544
     $ ,KIND=16),GC_5,MP__GC_5,GC_47,MP__GC_47,GC_5,MP__GC_5,GC_5
 
1545
     $ ,MP__GC_5,GC_5,MP__GC_5,4,1,1,315,AMPL(1,315),S(315))
1477
1546
C     Loop amplitude for loop diagram with ID 72
1478
1547
      CALL ML5_0_LOOP_5_5(30,1,4,3,2,5,DCMPLX(ZERO),CMPLX(MP__ZERO
1479
 
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(MDL_MB
1480
 
     $ ),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB
1481
 
     $ ,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),GC_4,MP__GC_4
1482
 
     $ ,GC_5,MP__GC_5,GC_47,MP__GC_47,GC_5,MP__GC_5,GC_5,MP__GC_5,4,1
1483
 
     $ ,1,316,AMPL(1,316),S(316))
 
1548
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16)
 
1549
     $ ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(MDL_MB)
 
1550
     $ ,CMPLX(MP__MDL_MB,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
 
1551
     $ ,GC_4,MP__GC_4,GC_5,MP__GC_5,GC_47,MP__GC_47,GC_5,MP__GC_5,GC_5
 
1552
     $ ,MP__GC_5,4,1,1,316,AMPL(1,316),S(316))
1484
1553
C     Loop amplitude for loop diagram with ID 73
1485
1554
      CALL ML5_0_LOOP_3_3(17,2,5,24,DCMPLX(ZERO),CMPLX(MP__ZERO
1486
1555
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(ZERO)
1504
1573
     $ ,GC_4,MP__GC_4,GC_4,MP__GC_4,GC_5,MP__GC_5,GC_5,MP__GC_5,3,1,1
1505
1574
     $ ,320,AMPL(1,320),S(320))
1506
1575
C     Loop amplitude for loop diagram with ID 77
1507
 
      CALL ML5_0_LOOP_3_4_3(33,1,2,1,5,1,2,7,DCMPLX(ZERO),CMPLX(MP__ZER
1508
 
     $ O,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MB)
1509
 
     $ ,CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_6,MP__GC_6,GC_5
1510
 
     $ ,MP__GC_5,1,1,1,321,AMPL(1,321),S(321))
1511
 
      CALL ML5_0_LOOP_3_4_3(34,1,2,1,5,1,2,7,DCMPLX(ZERO),CMPLX(MP__ZER
1512
 
     $ O,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MB)
1513
 
     $ ,CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_6,MP__GC_6,GC_5
1514
 
     $ ,MP__GC_5,1,1,1,322,AMPL(1,322),S(322))
1515
 
      CALL ML5_0_LOOP_3_4_3(35,1,2,1,5,1,2,7,DCMPLX(ZERO),CMPLX(MP__ZER
1516
 
     $ O,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MB)
1517
 
     $ ,CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_6,MP__GC_6,GC_5
1518
 
     $ ,MP__GC_5,1,1,1,323,AMPL(1,323),S(323))
 
1576
      CALL ML5_0_LOOP_3_4_3(33,1,2,1,5,1,2,7,DCMPLX(ZERO)
 
1577
     $ ,CMPLX(MP__ZERO,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
 
1578
     $ ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_6
 
1579
     $ ,MP__GC_6,GC_5,MP__GC_5,1,1,1,321,AMPL(1,321),S(321))
 
1580
      CALL ML5_0_LOOP_3_4_3(34,1,2,1,5,1,2,7,DCMPLX(ZERO)
 
1581
     $ ,CMPLX(MP__ZERO,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
 
1582
     $ ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_6
 
1583
     $ ,MP__GC_6,GC_5,MP__GC_5,1,1,1,322,AMPL(1,322),S(322))
 
1584
      CALL ML5_0_LOOP_3_4_3(35,1,2,1,5,1,2,7,DCMPLX(ZERO)
 
1585
     $ ,CMPLX(MP__ZERO,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
 
1586
     $ ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_6
 
1587
     $ ,MP__GC_6,GC_5,MP__GC_5,1,1,1,323,AMPL(1,323),S(323))
1519
1588
C     Loop amplitude for loop diagram with ID 78
1520
1589
      CALL ML5_0_LOOP_5_5(36,1,3,5,2,4,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT
1521
1590
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(ZERO)
1532
1601
     $ ,MP__GC_47,4,1,1,325,AMPL(1,325),S(325))
1533
1602
C     Loop amplitude for loop diagram with ID 80
1534
1603
      CALL ML5_0_LOOP_5_5(38,1,4,2,3,5,DCMPLX(ZERO),CMPLX(MP__ZERO
1535
 
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(MDL_MT
1536
 
     $ ),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB
1537
 
     $ ,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),GC_4,MP__GC_4
1538
 
     $ ,GC_5,MP__GC_5,GC_5,MP__GC_5,GC_47,MP__GC_47,GC_5,MP__GC_5,4,1
1539
 
     $ ,1,326,AMPL(1,326),S(326))
 
1604
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16)
 
1605
     $ ,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(MDL_MB)
 
1606
     $ ,CMPLX(MP__MDL_MB,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
 
1607
     $ ,GC_4,MP__GC_4,GC_5,MP__GC_5,GC_5,MP__GC_5,GC_47,MP__GC_47,GC_5
 
1608
     $ ,MP__GC_5,4,1,1,326,AMPL(1,326),S(326))
1540
1609
C     Loop amplitude for loop diagram with ID 81
1541
1610
      CALL ML5_0_LOOP_5_5(39,1,2,4,3,5,DCMPLX(ZERO),CMPLX(MP__ZERO
1542
1611
     $ ,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MT)
1546
1615
     $ ,1,327,AMPL(1,327),S(327))
1547
1616
C     Loop amplitude for loop diagram with ID 82
1548
1617
      CALL ML5_0_LOOP_5_5(40,1,2,3,5,4,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT
1549
 
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(MDL_MB
1550
 
     $ ),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO
1551
 
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),GC_5
1552
 
     $ ,MP__GC_5,GC_5,MP__GC_5,GC_47,MP__GC_47,GC_5,MP__GC_5,GC_5
1553
 
     $ ,MP__GC_5,4,1,1,328,AMPL(1,328),S(328))
 
1618
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16)
 
1619
     $ ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(ZERO)
 
1620
     $ ,CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT
 
1621
     $ ,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5,GC_47,MP__GC_47,GC_5
 
1622
     $ ,MP__GC_5,GC_5,MP__GC_5,4,1,1,328,AMPL(1,328),S(328))
1554
1623
C     Loop amplitude for loop diagram with ID 83
1555
1624
      CALL ML5_0_LOOP_5_5(41,1,2,5,3,4,DCMPLX(ZERO),CMPLX(MP__ZERO
1556
1625
     $ ,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MB)
1561
1630
C     Loop amplitude for loop diagram with ID 84
1562
1631
      CALL ML5_0_LOOP_4_5_4(42,1,1,2,1,3,4,1,2,5,DCMPLX(MDL_MT)
1563
1632
     $ ,CMPLX(MP__MDL_MT,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
1564
 
     $ ,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__M
1565
 
     $ DL_MB,KIND=16),GC_47,MP__GC_47,GC_5,MP__GC_5,GC_6,MP__GC_6,GC_5
1566
 
     $ ,MP__GC_5,2,1,1,330,AMPL(1,330),S(330))
 
1633
     $ ,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MB)
 
1634
     $ ,CMPLX(MP__MDL_MB,KIND=16),GC_47,MP__GC_47,GC_5,MP__GC_5,GC_6
 
1635
     $ ,MP__GC_6,GC_5,MP__GC_5,2,1,1,330,AMPL(1,330),S(330))
1567
1636
      CALL ML5_0_LOOP_4_5_4(43,1,1,2,1,3,4,1,2,5,DCMPLX(MDL_MT)
1568
1637
     $ ,CMPLX(MP__MDL_MT,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
1569
 
     $ ,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__M
1570
 
     $ DL_MB,KIND=16),GC_47,MP__GC_47,GC_5,MP__GC_5,GC_6,MP__GC_6,GC_5
1571
 
     $ ,MP__GC_5,2,1,1,331,AMPL(1,331),S(331))
 
1638
     $ ,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MB)
 
1639
     $ ,CMPLX(MP__MDL_MB,KIND=16),GC_47,MP__GC_47,GC_5,MP__GC_5,GC_6
 
1640
     $ ,MP__GC_6,GC_5,MP__GC_5,2,1,1,331,AMPL(1,331),S(331))
1572
1641
      CALL ML5_0_LOOP_4_5_4(44,1,1,2,1,3,4,1,2,5,DCMPLX(MDL_MT)
1573
1642
     $ ,CMPLX(MP__MDL_MT,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
1574
 
     $ ,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__M
1575
 
     $ DL_MB,KIND=16),GC_47,MP__GC_47,GC_5,MP__GC_5,GC_6,MP__GC_6,GC_5
1576
 
     $ ,MP__GC_5,2,1,1,332,AMPL(1,332),S(332))
 
1643
     $ ,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MB)
 
1644
     $ ,CMPLX(MP__MDL_MB,KIND=16),GC_47,MP__GC_47,GC_5,MP__GC_5,GC_6
 
1645
     $ ,MP__GC_6,GC_5,MP__GC_5,2,1,1,332,AMPL(1,332),S(332))
1577
1646
C     Loop amplitude for loop diagram with ID 85
1578
1647
      CALL ML5_0_LOOP_2_2(1,6,13,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT
1579
1648
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),GC_5
1592
1661
     $ ,CMPLX(MP__ZERO,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5
1593
1662
     $ ,MP__GC_5,2,1,1,336,AMPL(1,336),S(336))
1594
1663
C     Loop amplitude for loop diagram with ID 89
1595
 
      CALL ML5_0_LOOP_3_3(46,4,6,8,DCMPLX(ZERO),CMPLX(MP__ZERO
1596
 
     $ ,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MT)
 
1664
      CALL ML5_0_LOOP_3_3(46,4,6,8,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
 
1665
     $ ,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MT)
1597
1666
     $ ,CMPLX(MP__MDL_MT,KIND=16),GC_5,MP__GC_5,GC_4,MP__GC_4,GC_5
1598
1667
     $ ,MP__GC_5,2,1,1,337,AMPL(1,337),S(337))
1599
1668
C     Loop amplitude for loop diagram with ID 90
1614
1683
     $ ,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),GC_5,MP__GC_5,GC_5
1615
1684
     $ ,MP__GC_5,1,1,1,341,AMPL(1,341),S(341))
1616
1685
C     Loop amplitude for loop diagram with ID 94
1617
 
      CALL ML5_0_LOOP_3_3(19,2,9,8,DCMPLX(ZERO),CMPLX(MP__ZERO
1618
 
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(ZERO)
 
1686
      CALL ML5_0_LOOP_3_3(19,2,9,8,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
 
1687
     $ ,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(ZERO)
1619
1688
     $ ,CMPLX(MP__ZERO,KIND=16),GC_4,MP__GC_4,GC_5,MP__GC_5,GC_5
1620
1689
     $ ,MP__GC_5,2,1,1,342,AMPL(1,342),S(342))
1621
1690
C     Loop amplitude for loop diagram with ID 95
1644
1713
     $ ,MP__GC_5,2,1,1,347,AMPL(1,347),S(347))
1645
1714
C     Loop amplitude for loop diagram with ID 100
1646
1715
      CALL ML5_0_LOOP_3_3(2,1,2,13,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT
1647
 
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(MDL_MT
1648
 
     $ ),CMPLX(MP__MDL_MT,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5
1649
 
     $ ,MP__GC_5,3,1,1,348,AMPL(1,348),S(348))
 
1716
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16)
 
1717
     $ ,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),GC_5,MP__GC_5,GC_5
 
1718
     $ ,MP__GC_5,GC_5,MP__GC_5,3,1,1,348,AMPL(1,348),S(348))
1650
1719
C     Loop amplitude for loop diagram with ID 101
1651
1720
      CALL ML5_0_LOOP_3_3(2,1,2,14,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT
1652
 
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(MDL_MT
1653
 
     $ ),CMPLX(MP__MDL_MT,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5
1654
 
     $ ,MP__GC_5,3,1,1,349,AMPL(1,349),S(349))
 
1721
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16)
 
1722
     $ ,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),GC_5,MP__GC_5,GC_5
 
1723
     $ ,MP__GC_5,GC_5,MP__GC_5,3,1,1,349,AMPL(1,349),S(349))
1655
1724
C     Loop amplitude for loop diagram with ID 102
1656
1725
      CALL ML5_0_LOOP_4_4(22,1,2,4,8,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT
1657
1726
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(ZERO)
1660
1729
     $ ,MP__GC_5,3,1,1,350,AMPL(1,350),S(350))
1661
1730
C     Loop amplitude for loop diagram with ID 103
1662
1731
      CALL ML5_0_LOOP_3_3(3,1,2,13,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT
1663
 
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(MDL_MT
1664
 
     $ ),CMPLX(MP__MDL_MT,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5
1665
 
     $ ,MP__GC_5,3,1,1,351,AMPL(1,351),S(351))
 
1732
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16)
 
1733
     $ ,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),GC_5,MP__GC_5,GC_5
 
1734
     $ ,MP__GC_5,GC_5,MP__GC_5,3,1,1,351,AMPL(1,351),S(351))
1666
1735
C     Loop amplitude for loop diagram with ID 104
1667
1736
      CALL ML5_0_LOOP_3_3(3,1,2,14,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT
1668
 
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(MDL_MT
1669
 
     $ ),CMPLX(MP__MDL_MT,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5
1670
 
     $ ,MP__GC_5,3,1,1,352,AMPL(1,352),S(352))
 
1737
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16)
 
1738
     $ ,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),GC_5,MP__GC_5,GC_5
 
1739
     $ ,MP__GC_5,GC_5,MP__GC_5,3,1,1,352,AMPL(1,352),S(352))
1671
1740
C     Loop amplitude for loop diagram with ID 105
1672
1741
      CALL ML5_0_LOOP_4_4(23,1,2,8,4,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT
1673
1742
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(ZERO)
1710
1779
     $ ,MP__GC_5,2,1,1,360,AMPL(1,360),S(360))
1711
1780
C     Loop amplitude for loop diagram with ID 113
1712
1781
      CALL ML5_0_LOOP_4_4(48,1,4,2,8,DCMPLX(ZERO),CMPLX(MP__ZERO
1713
 
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(MDL_MT
1714
 
     $ ),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO
1715
 
     $ ,KIND=16),GC_4,MP__GC_4,GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5
1716
 
     $ ,MP__GC_5,3,1,1,361,AMPL(1,361),S(361))
 
1782
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16)
 
1783
     $ ,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(ZERO)
 
1784
     $ ,CMPLX(MP__ZERO,KIND=16),GC_4,MP__GC_4,GC_5,MP__GC_5,GC_5
 
1785
     $ ,MP__GC_5,GC_5,MP__GC_5,3,1,1,361,AMPL(1,361),S(361))
1717
1786
C     Loop amplitude for loop diagram with ID 114
1718
1787
      CALL ML5_0_LOOP_3_3(19,2,4,28,DCMPLX(ZERO),CMPLX(MP__ZERO
1719
1788
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(ZERO)
1737
1806
     $ ,GC_4,MP__GC_4,GC_4,MP__GC_4,GC_5,MP__GC_5,GC_5,MP__GC_5,3,1,1
1738
1807
     $ ,365,AMPL(1,365),S(365))
1739
1808
C     Loop amplitude for loop diagram with ID 118
1740
 
      CALL ML5_0_LOOP_3_4_3(49,1,2,1,4,1,2,8,DCMPLX(ZERO),CMPLX(MP__ZER
1741
 
     $ O,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MT)
1742
 
     $ ,CMPLX(MP__MDL_MT,KIND=16),GC_5,MP__GC_5,GC_6,MP__GC_6,GC_5
1743
 
     $ ,MP__GC_5,1,1,1,366,AMPL(1,366),S(366))
1744
 
      CALL ML5_0_LOOP_3_4_3(50,1,2,1,4,1,2,8,DCMPLX(ZERO),CMPLX(MP__ZER
1745
 
     $ O,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MT)
1746
 
     $ ,CMPLX(MP__MDL_MT,KIND=16),GC_5,MP__GC_5,GC_6,MP__GC_6,GC_5
1747
 
     $ ,MP__GC_5,1,1,1,367,AMPL(1,367),S(367))
1748
 
      CALL ML5_0_LOOP_3_4_3(51,1,2,1,4,1,2,8,DCMPLX(ZERO),CMPLX(MP__ZER
1749
 
     $ O,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MT)
1750
 
     $ ,CMPLX(MP__MDL_MT,KIND=16),GC_5,MP__GC_5,GC_6,MP__GC_6,GC_5
1751
 
     $ ,MP__GC_5,1,1,1,368,AMPL(1,368),S(368))
 
1809
      CALL ML5_0_LOOP_3_4_3(49,1,2,1,4,1,2,8,DCMPLX(ZERO)
 
1810
     $ ,CMPLX(MP__ZERO,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
 
1811
     $ ,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),GC_5,MP__GC_5,GC_6
 
1812
     $ ,MP__GC_6,GC_5,MP__GC_5,1,1,1,366,AMPL(1,366),S(366))
 
1813
      CALL ML5_0_LOOP_3_4_3(50,1,2,1,4,1,2,8,DCMPLX(ZERO)
 
1814
     $ ,CMPLX(MP__ZERO,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
 
1815
     $ ,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),GC_5,MP__GC_5,GC_6
 
1816
     $ ,MP__GC_6,GC_5,MP__GC_5,1,1,1,367,AMPL(1,367),S(367))
 
1817
      CALL ML5_0_LOOP_3_4_3(51,1,2,1,4,1,2,8,DCMPLX(ZERO)
 
1818
     $ ,CMPLX(MP__ZERO,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
 
1819
     $ ,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),GC_5,MP__GC_5,GC_6
 
1820
     $ ,MP__GC_6,GC_5,MP__GC_5,1,1,1,368,AMPL(1,368),S(368))
1752
1821
C     Loop amplitude for loop diagram with ID 119
1753
1822
      CALL ML5_0_LOOP_2_2(52,6,13,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
1754
1823
     $ ,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),GC_4,MP__GC_4,GC_4
1887
1956
                WRITE(*,*) '##W02A WARNING Inconsistent helicity '
1888
1957
     $           ,HELPICKED
1889
1958
                IF(HELINITSTARTOVER) THEN
1890
 
                  WRITE(*,*) '##I01 INFO Initialization starting ove'
1891
 
     $             //'r because of inconsistency in the helicit'
1892
 
     $             //'y filter setup.'
 
1959
                  WRITE(*,*) '##I01 INFO Initialization starting over'
 
1960
     $             //' because of inconsistency in the helicity filter'
 
1961
     $             //' setup.'
1893
1962
                  NTRY=0
1894
1963
                ENDIF
1895
1964
              ENDIF
1898
1967
                WRITE(*,*) '##W02B WARNING Inconsistent helicity '
1899
1968
     $           ,HELPICKED
1900
1969
                IF(HELINITSTARTOVER) THEN
1901
 
                  WRITE(*,*) '##I01 INFO Initialization starting ove'
1902
 
     $             //'r because of inconsistency in the helicit'
1903
 
     $             //'y filter setup.'
 
1970
                  WRITE(*,*) '##I01 INFO Initialization starting over'
 
1971
     $             //' because of inconsistency in the helicity filter'
 
1972
     $             //' setup.'
1904
1973
                  NTRY=0
1905
1974
                ELSE
1906
1975
                  GOODHEL(HELPICKED)=.TRUE.
1920
1989
                  WRITE(*,*) '##W02 WARNING Inconsistent loop amp ',I
1921
1990
     $             ,' for helicity ',HELPICKED,'.'
1922
1991
                  IF(LOOPINITSTARTOVER) THEN
1923
 
                    WRITE(*,*) '##I01 INFO Initialization startin'
1924
 
     $               //'g over because of inconsistency in the loo'
1925
 
     $               //'p filter setup.'
 
1992
                    WRITE(*,*) '##I01 INFO Initialization starting'
 
1993
     $               //' over because of inconsistency in the loop'
 
1994
     $               //' filter setup.'
1926
1995
                    NTRY=0
1927
1996
                  ELSE
1928
1997
                    GOODAMP(I,HELPICKED)=.TRUE.
1933
2002
          ENDIF
1934
2003
        ELSEIF (.NOT.HELDOUBLECHECKED)THEN
1935
2004
          IF ((.NOT.GOODHEL(HELPICKED)).AND.(.NOT.ML5_0_ISZERO(ABS(ANS(
1936
 
     $     1))+ABS(ANS(2))+ABS(ANS(3)),REF/DBLE(NCOMB),-1))) THEN
1937
 
            WRITE(*,*) '##W15 Helicity filter could not be successfull'
1938
 
     $       //'y double checked.'
1939
 
            WRITE(*,*) '##One reason for this is that you have change'
1940
 
     $       //'d sensible parameters which affected what are the zer'
1941
 
     $       //'o helicity configurations.'
1942
 
            WRITE(*,*) '##MadLoop will try to reset the Helicit'
1943
 
     $       //'y filter with the next PS points it receives.'
 
2005
     $1))+ABS(ANS(2))+ABS(ANS(3)),REF/DBLE(NCOMB),-1))) THEN
 
2006
            WRITE(*,*) '##W15 Helicity filter could not be'
 
2007
     $       //' successfully double checked.'
 
2008
            WRITE(*,*) '##One reason for this is that you have changed'
 
2009
     $       //' sensible parameters which affected what are the zero'
 
2010
     $       //' helicity configurations.'
 
2011
            WRITE(*,*) '##MadLoop will try to reset the Helicity'
 
2012
     $       //' filter with the next PS points it receives.'
1944
2013
            NTRY=0
1945
2014
            OPEN(30,FILE=HELFILTERFN,ERR=349)
1946
2015
 349        CONTINUE
1947
2016
            CLOSE(30,STATUS='delete')
1948
2017
          ENDIF
1949
2018
C         SET HELDOUBLECHECKED TO .TRUE. WHEN DONE
1950
 
C         even if it failed we do not want to redo the check afterwards
1951
 
C          if HELINITSTARTOVER=.FALSE.
 
2019
C         even if it failed we do not want to redo the check
 
2020
C          afterwards if HELINITSTARTOVER=.FALSE.
1952
2021
          IF (HELPICKED.EQ.NCOMB.AND.(NTRY.NE.0.OR..NOT.HELINITSTARTOVE
1953
 
     $     R)) THEN
 
2022
     $R)) THEN
1954
2023
            DONEHELDOUBLECHECK=.TRUE.
1955
2024
          ENDIF
1956
2025
        ENDIF
1965
2034
          ANS(2)=BUFFR(2)
1966
2035
          ANS(3)=BUFFR(3)
1967
2036
C         We add one here to the number of PS points used for building
1968
 
C          the reference scale for comparison (used only for loop-induc
1969
 
C         ed processes).
 
2037
C          the reference scale for comparison (used only for
 
2038
C          loop-induced processes).
1970
2039
          NPSPOINTS = NPSPOINTS+1
1971
2040
          IF(NTRY.EQ.0) THEN
1972
2041
            NATTEMPTS=NATTEMPTS+1
1973
2042
            IF(NATTEMPTS.EQ.MAXATTEMPTS) THEN
1974
 
              WRITE(*,*) '##E01 ERROR Could not initialize the filter'
1975
 
     $         //'s in ',MAXATTEMPTS,' trials'
 
2043
              WRITE(*,*) '##E01 ERROR Could not initialize the filters'
 
2044
     $         //' in ',MAXATTEMPTS,' trials'
1976
2045
              STOP
1977
2046
            ENDIF
1978
2047
          ENDIF
1987
2056
        ENDIF
1988
2057
      ENDDO
1989
2058
 
1990
 
      IF(.NOT.CHECKPHASE.AND.HELDOUBLECHECKED.AND.(CTMODERUN.LE.
1991
 
     $ -1)) THEN
 
2059
      IF(.NOT.CHECKPHASE.AND.HELDOUBLECHECKED.AND.(CTMODERUN.LE.-1))
 
2060
     $  THEN
1992
2061
        STAB_INDEX=STAB_INDEX+1
1993
2062
        IF(DOING_QP_EVALS) THEN
1994
2063
          QP_RES(1,STAB_INDEX)=ANS(1)
2018
2087
        CTMODE=BASIC_CT_MODE
2019
2088
 
2020
2089
        IF(.NOT.EVAL_DONE(3).AND. ((DOING_QP_EVALS.AND.NROTATIONS_QP.GE
2021
 
     $   .1).OR.((.NOT.DOING_QP_EVALS).AND.NROTATIONS_DP.GE.1)) ) THEN
 
2090
     $.1).OR.((.NOT.DOING_QP_EVALS).AND.NROTATIONS_DP.GE.1)) ) THEN
2022
2091
          EVAL_DONE(3)=.TRUE.
2023
2092
          CALL ML5_0_ROTATE_PS(PS,P,1)
2024
2093
          IF (DOING_QP_EVALS) CALL ML5_0_MP_ROTATE_PS(MP_PS,MP_P,1)
2026
2095
        ENDIF
2027
2096
 
2028
2097
        IF(.NOT.EVAL_DONE(4).AND. ((DOING_QP_EVALS.AND.NROTATIONS_QP.GE
2029
 
     $   .2).OR.((.NOT.DOING_QP_EVALS).AND.NROTATIONS_DP.GE.2)) ) THEN
 
2098
     $.2).OR.((.NOT.DOING_QP_EVALS).AND.NROTATIONS_DP.GE.2)) ) THEN
2030
2099
          EVAL_DONE(4)=.TRUE.
2031
2100
          CALL ML5_0_ROTATE_PS(PS,P,2)
2032
2101
          IF (DOING_QP_EVALS) CALL ML5_0_MP_ROTATE_PS(MP_PS,MP_P,2)
2046
2115
            RET_CODE_H=4
2047
2116
            NEPS=NEPS+1
2048
2117
            CALL ML5_0_COMPUTE_ACCURACY(DP_RES,N_DP_EVAL,TEMP1,TEMP)
2049
 
            WRITE(*,*) '##W03 WARNING An unstable PS point was'
2050
 
     $       ,       ' detected.'
 
2118
            WRITE(*,*) '##W03 WARNING An unstable PS point was',      
 
2119
     $        ' detected.'
2051
2120
            WRITE(*,*) '##(DP,QP) accuracies : (',TEMP1,',',ACC,')'
2052
2121
            WRITE(*,*) '##Best estimate (fin,1eps,2eps) :',(ANS(I),I=1
2053
2122
     $       ,3)
2064
2133
              ENDDO
2065
2134
            ENDIF
2066
2135
            IF(NEPS.EQ.10) THEN
2067
 
              WRITE(*,*) '##Further output of the details of thes'
2068
 
     $         //'e unstable PS points will now be suppressed.'
 
2136
              WRITE(*,*) '##Further output of the details of these'
 
2137
     $         //' unstable PS points will now be suppressed.'
2069
2138
            ENDIF
2070
2139
          ENDIF
2071
2140
        ELSE
2115
2184
        CTMODEINIT=CTMODEINIT_BU
2116
2185
      ENDIF
2117
2186
 
 
2187
C     Reinitialize the Lorentz test if it had been disabled because
 
2188
C      spin-2 particles are in the external states.
 
2189
      NROTATIONS_DP = NROTATIONS_DP_BU
 
2190
      NROTATIONS_QP = NROTATIONS_QP_BU
 
2191
 
2118
2192
C     Conform to the returned synthax of split orders even though the
2119
2193
C      default output does not support it (this then done only for
2120
2194
C      compatibility purpose).
2140
2214
 
2141
2215
      END
2142
2216
 
2143
 
      SUBROUTINE ML5_0_COMPUTE_ACCURACY(FULLLIST, LENGTH, ACC
2144
 
     $ , ESTIMATE)
 
2217
      SUBROUTINE ML5_0_COMPUTE_ACCURACY(FULLLIST, LENGTH, ACC,
 
2218
     $  ESTIMATE)
2145
2219
      IMPLICIT NONE
2146
2220
C     
2147
2221
C     PARAMETERS 
2267
2341
C     ----------
2268
2342
C     BEGIN CODE
2269
2343
C     ----------
2270
 
      WRITE(*,*) '##WARNING:: Ignored, the possibility of selectin'
2271
 
     $ //'g specific squared order contributions is not available i'
2272
 
     $ //'n the default mode.'
 
2344
      WRITE(*,*) '##WARNING:: Ignored, the possibility of selecting'
 
2345
     $ //' specific squared order contributions is not available in'
 
2346
     $ //' the default mode.'
2273
2347
 
2274
2348
      END
2275
2349
 
2355
2429
C     ----------
2356
2430
      USER_STAB_PREC = PREC_ASKED
2357
2431
      CALL ML5_0_SLOOPMATRIXHEL(P,HEL,ANS)
2358
 
      IF(ALWAYS_TEST_STABILITY.AND.(H.EQ.1.OR.ACCURACY(0).LT.0.0D0)
2359
 
     $ ) THEN
 
2432
      IF(ALWAYS_TEST_STABILITY.AND.(H.EQ.1.OR.ACCURACY(0).LT.0.0D0))
 
2433
     $  THEN
2360
2434
        BYPASS_CHECK = .TRUE.
2361
2435
        CALL ML5_0_SLOOPMATRIXHEL(P,HEL,ANS)
2362
2436
        BYPASS_CHECK = .FALSE.
2378
2452
     $ ,RET_CODE)
2379
2453
C     
2380
2454
C     Inputs are:
2381
 
C     P(0:3, Nexternal)  double  :: Kinematic configuration (E,px,py,pz
2382
 
C     )
 
2455
C     P(0:3, Nexternal)  double  :: Kinematic configuration
 
2456
C      (E,px,py,pz)
2383
2457
C     PEC_ASKED          double  :: Target relative accuracy, -1 for
2384
2458
C      default
2385
2459
C     
2451
2525
C     ----------
2452
2526
      USER_STAB_PREC = PREC_ASKED
2453
2527
      CALL ML5_0_SLOOPMATRIX(P,ANS)
2454
 
      IF(ALWAYS_TEST_STABILITY.AND.(H.EQ.1.OR.ACCURACY(0).LT.0.0D0)
2455
 
     $ ) THEN
 
2528
      IF(ALWAYS_TEST_STABILITY.AND.(H.EQ.1.OR.ACCURACY(0).LT.0.0D0))
 
2529
     $  THEN
2456
2530
        BYPASS_CHECK = .TRUE.
2457
2531
        CALL ML5_0_SLOOPMATRIX(P,ANS)
2458
2532
        BYPASS_CHECK = .FALSE.