← Back to branch summary

~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

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

  • Committer: olivier Mattelaer
  • Date: 2015-03-05 00:14:16 UTC
  • mfrom: (258.1.9 2.3)
  • mto: (258.8.1 2.3)
  • mto: This revision was merged to the branch mainline in revision 259.
  • Revision ID: olivier.mattelaer@uclouvain.be-20150305001416-y9mzeykfzwnl9t0j
partial merge

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/short_ML_SMQCD_default/ddx_ttx/loop_matrix.f
28
28
C     CONSTANTS
29
29
C     
30
30
      INTEGER    NLOOPAMPS
31
 
      PARAMETER (NLOOPAMPS=43)
 
31
      PARAMETER (NLOOPAMPS=40)
32
32
C     
33
33
C     ARGUMENTS 
34
34
C     
93
93
      INTEGER NBORNAMPS
94
94
      PARAMETER (NBORNAMPS=1)
95
95
      INTEGER    NLOOPAMPS, NCTAMPS
96
 
      PARAMETER (NLOOPAMPS=43, NCTAMPS=29)
 
96
      PARAMETER (NLOOPAMPS=40, NCTAMPS=29)
97
97
      INTEGER    NEXTERNAL
98
98
      PARAMETER (NEXTERNAL=4)
99
99
      INTEGER    NWAVEFUNCS
313
313
        ML_INIT = .FALSE.
314
314
      ENDIF
315
315
 
 
316
C     Setup the file paths
 
317
      CALL JOINPATH(MLPATH,PARAMFNAME,PARAMFN)
 
318
      CALL JOINPATH(MLPATH,PROC_PREFIX,TMP)
 
319
      CALL JOINPATH(TMP,HELCONFIGFNAME,HELCONFIGFN)
 
320
      CALL JOINPATH(TMP,LOOPFILTERFNAME,LOOPFILTERFN)
 
321
      CALL JOINPATH(TMP,COLORNUMFNAME,COLORNUMFN)
 
322
      CALL JOINPATH(TMP,COLORDENOMFNAME,COLORDENOMFN)
 
323
      CALL JOINPATH(TMP,HELFILTERFNAME,HELFILTERFN)
 
324
 
 
325
      OPEN(1, FILE=COLORNUMFN, ERR=104, STATUS='OLD',          
 
326
     $  ACTION='READ')
 
327
      DO I=1,NLOOPAMPS
 
328
        READ(1,*,END=105) (CF_N(I,J),J=1,NBORNAMPS)
 
329
      ENDDO
 
330
      GOTO 105
 
331
 104  CONTINUE
 
332
      STOP 'Color factors could not be initialized from file ML5_0_Col'
 
333
     $ //'orNumFactors.dat. File not found'
 
334
 105  CONTINUE
 
335
      CLOSE(1)
 
336
      OPEN(1, FILE=COLORDENOMFN, ERR=106, STATUS='OLD',          
 
337
     $  ACTION='READ')
 
338
      DO I=1,NLOOPAMPS
 
339
        READ(1,*,END=107) (CF_D(I,J),J=1,NBORNAMPS)
 
340
      ENDDO
 
341
      GOTO 107
 
342
 106  CONTINUE
 
343
      STOP 'Color factors could not be initialized from file ML5_0_Col'
 
344
     $ //'orDenomFactors.dat. File not found'
 
345
 107  CONTINUE
 
346
      CLOSE(1)
 
347
      OPEN(1, FILE=HELCONFIGFN, ERR=108, STATUS='OLD',                
 
348
     $   ACTION='READ')
 
349
      DO H=1,NCOMB
 
350
        READ(1,*,END=109) (HELC(I,H),I=1,NEXTERNAL)
 
351
      ENDDO
 
352
      GOTO 109
 
353
 108  CONTINUE
 
354
      STOP 'Color helictiy configurations could not be initialize'
 
355
     $ //'d from file ML5_0_HelConfigs.dat. File not found'
 
356
 109  CONTINUE
 
357
      CLOSE(1)
 
358
      IF(BOOTANDSTOP) THEN
 
359
        WRITE(*,*) 'Stopped by user request.'
 
360
        STOP
 
361
      ENDIF
 
362
 
316
363
      IF(NTRY.EQ.0) THEN
317
 
 
318
 
C       Setup the file paths
319
 
        CALL JOINPATH(MLPATH,PARAMFNAME,PARAMFN)
320
 
        CALL JOINPATH(MLPATH,PROC_PREFIX,TMP)
321
 
        CALL JOINPATH(TMP,HELCONFIGFNAME,HELCONFIGFN)
322
 
        CALL JOINPATH(TMP,LOOPFILTERFNAME,LOOPFILTERFN)
323
 
        CALL JOINPATH(TMP,COLORNUMFNAME,COLORNUMFN)
324
 
        CALL JOINPATH(TMP,COLORDENOMFNAME,COLORDENOMFN)
325
 
        CALL JOINPATH(TMP,HELFILTERFNAME,HELFILTERFN)
326
 
 
327
364
        CALL ML5_0_SET_N_EVALS(N_DP_EVAL,N_QP_EVAL)
328
 
        HELDOUBLECHECKED=.NOT.DOUBLECHECKHELICITYFILTER
 
365
        HELDOUBLECHECKED=(.NOT.DOUBLECHECKHELICITYFILTER).OR.(HELICITYF
 
366
     $   ILTERLEVEL.EQ.0)
329
367
        DO J=1,NCOMB
330
368
          DO I=1,NCTAMPS
331
369
            GOODAMP(I,J)=.TRUE.
346
384
        ENDDO
347
385
 101    CONTINUE
348
386
        CLOSE(1)
 
387
        IF (HELICITYFILTERLEVEL.EQ.0) THEN
 
388
          FOUNDHELFILTER=.TRUE.
 
389
          DO J=1,NCOMB
 
390
            GOODHEL(J)=.TRUE.
 
391
          ENDDO
 
392
          GOTO 122
 
393
        ENDIF
349
394
        OPEN(1, FILE=HELFILTERFN, ERR=102, STATUS='OLD',          
350
395
     $    ACTION='READ')
351
396
        READ(1,*,END=103) (GOODHEL(I),I=1,NCOMB)
357
402
        ENDDO
358
403
 103    CONTINUE
359
404
        CLOSE(1)
360
 
        OPEN(1, FILE=COLORNUMFN, ERR=104, STATUS='OLD',          
361
 
     $    ACTION='READ')
362
 
        DO I=1,NLOOPAMPS
363
 
          READ(1,*,END=105) (CF_N(I,J),J=1,NBORNAMPS)
364
 
        ENDDO
365
 
        GOTO 105
366
 
 104    CONTINUE
367
 
        STOP 'Color factors could not be initialized from fil'
368
 
     $   //'e ML5_0_ColorNumFactors.dat. File not found'
369
 
 105    CONTINUE
370
 
        CLOSE(1)
371
 
        OPEN(1, FILE=COLORDENOMFN, ERR=106, STATUS='OLD',          
372
 
     $    ACTION='READ')
373
 
        DO I=1,NLOOPAMPS
374
 
          READ(1,*,END=107) (CF_D(I,J),J=1,NBORNAMPS)
375
 
        ENDDO
376
 
        GOTO 107
377
 
 106    CONTINUE
378
 
        STOP 'Color factors could not be initialized from fil'
379
 
     $   //'e ML5_0_ColorDenomFactors.dat. File not found'
380
 
 107    CONTINUE
381
 
        CLOSE(1)
382
 
        OPEN(1, FILE=HELCONFIGFN, ERR=108, STATUS='OLD',              
383
 
     $       ACTION='READ')
384
 
        DO H=1,NCOMB
385
 
          READ(1,*,END=109) (HELC(I,H),I=1,NEXTERNAL)
386
 
        ENDDO
387
 
        GOTO 109
388
 
 108    CONTINUE
389
 
        STOP 'Color helictiy configurations could not be initialize'
390
 
     $   //'d from file ML5_0_HelConfigs.dat. File not found'
391
 
 109    CONTINUE
392
 
        CLOSE(1)
393
 
        IF(BOOTANDSTOP) THEN
394
 
          WRITE(*,*) 'Stopped by user request.'
395
 
          STOP
396
 
        ENDIF
 
405
 122    CONTINUE
397
406
      ENDIF
398
407
 
399
408
      MP_DONE=.FALSE.
577
586
          CALL FFV1P0_3(W(1,4,H),W(1,3,H),GC_5,ZERO,ZERO,W(1,6,H))
578
587
C         Counter-term amplitude(s) for loop diagram number 2
579
588
          CALL R2_GG_1_0(W(1,5,H),W(1,6,H),R2_GGQ,AMPL(1,1))
 
589
          CALL R2_GG_1_0(W(1,5,H),W(1,6,H),R2_GGQ,AMPL(1,2))
 
590
          CALL R2_GG_1_0(W(1,5,H),W(1,6,H),R2_GGQ,AMPL(1,3))
 
591
          CALL R2_GG_1_0(W(1,5,H),W(1,6,H),R2_GGQ,AMPL(1,4))
580
592
C         Counter-term amplitude(s) for loop diagram number 5
581
593
          CALL FFV1_0(W(1,1,H),W(1,2,H),W(1,6,H),UV_GQQQ_1EPS,AMPL(2
582
 
     $     ,2))
583
 
          CALL FFV1_0(W(1,1,H),W(1,2,H),W(1,6,H),UV_GQQQ_1EPS,AMPL(2
584
 
     $     ,3))
585
 
          CALL FFV1_0(W(1,1,H),W(1,2,H),W(1,6,H),UV_GQQQ_1EPS,AMPL(2
586
 
     $     ,4))
587
 
          CALL FFV1_0(W(1,1,H),W(1,2,H),W(1,6,H),UV_GQQQ_1EPS,AMPL(2
588
594
     $     ,5))
589
 
          CALL FFV1_0(W(1,1,H),W(1,2,H),W(1,6,H),UV_GQQB,AMPL(1,6))
 
595
          CALL FFV1_0(W(1,1,H),W(1,2,H),W(1,6,H),UV_GQQQ_1EPS,AMPL(2
 
596
     $     ,6))
590
597
          CALL FFV1_0(W(1,1,H),W(1,2,H),W(1,6,H),UV_GQQQ_1EPS,AMPL(2
591
598
     $     ,7))
592
 
          CALL FFV1_0(W(1,1,H),W(1,2,H),W(1,6,H),UV_GQQT,AMPL(1,8))
593
 
          CALL FFV1_0(W(1,1,H),W(1,2,H),W(1,6,H),UV_GQQQ_1EPS,AMPL(2
594
 
     $     ,9))
 
599
          CALL FFV1_0(W(1,1,H),W(1,2,H),W(1,6,H),UV_GQQQ_1EPS,AMPL(2
 
600
     $     ,8))
 
601
          CALL FFV1_0(W(1,1,H),W(1,2,H),W(1,6,H),UV_GQQB,AMPL(1,9))
 
602
          CALL FFV1_0(W(1,1,H),W(1,2,H),W(1,6,H),UV_GQQQ_1EPS,AMPL(2
 
603
     $     ,10))
 
604
          CALL FFV1_0(W(1,1,H),W(1,2,H),W(1,6,H),UV_GQQT,AMPL(1,11))
 
605
          CALL FFV1_0(W(1,1,H),W(1,2,H),W(1,6,H),UV_GQQQ_1EPS,AMPL(2
 
606
     $     ,12))
595
607
          CALL FFV1_0(W(1,1,H),W(1,2,H),W(1,6,H),UV_GQQG_1EPS,AMPL(2
596
 
     $     ,10))
597
 
          CALL FFV1_0(W(1,1,H),W(1,2,H),W(1,6,H),R2_GQQ,AMPL(1,11))
 
608
     $     ,13))
 
609
          CALL FFV1_0(W(1,1,H),W(1,2,H),W(1,6,H),R2_GQQ,AMPL(1,14))
598
610
C         Counter-term amplitude(s) for loop diagram number 7
599
 
          CALL R2_GG_1_0(W(1,5,H),W(1,6,H),R2_GGQ,AMPL(1,12))
600
 
C         Counter-term amplitude(s) for loop diagram number 8
601
 
          CALL R2_GG_1_0(W(1,5,H),W(1,6,H),R2_GGQ,AMPL(1,13))
602
 
C         Counter-term amplitude(s) for loop diagram number 9
603
 
          CALL R2_GG_1_0(W(1,5,H),W(1,6,H),R2_GGQ,AMPL(1,14))
604
 
C         Counter-term amplitude(s) for loop diagram number 10
605
611
          CALL R2_GG_1_R2_GG_3_0(W(1,5,H),W(1,6,H),R2_GGQ,R2_GGB
606
612
     $     ,AMPL(1,15))
607
 
C         Counter-term amplitude(s) for loop diagram number 11
 
613
C         Counter-term amplitude(s) for loop diagram number 8
608
614
          CALL R2_GG_1_R2_GG_3_0(W(1,5,H),W(1,6,H),R2_GGQ,R2_GGT
609
615
     $     ,AMPL(1,16))
610
 
C         Counter-term amplitude(s) for loop diagram number 12
 
616
C         Counter-term amplitude(s) for loop diagram number 9
611
617
          CALL FFV1_0(W(1,4,H),W(1,3,H),W(1,5,H),UV_GQQQ_1EPS,AMPL(2
612
618
     $     ,17))
613
619
          CALL FFV1_0(W(1,4,H),W(1,3,H),W(1,5,H),UV_GQQQ_1EPS,AMPL(2
625
631
          CALL FFV1_0(W(1,4,H),W(1,3,H),W(1,5,H),UV_GQQG_1EPS,AMPL(2
626
632
     $     ,25))
627
633
          CALL FFV1_0(W(1,4,H),W(1,3,H),W(1,5,H),R2_GQQ,AMPL(1,26))
628
 
C         Counter-term amplitude(s) for loop diagram number 14
 
634
C         Counter-term amplitude(s) for loop diagram number 11
629
635
          CALL R2_GG_1_R2_GG_2_0(W(1,5,H),W(1,6,H),R2_GGG_1,R2_GGG_2
630
636
     $     ,AMPL(1,27))
631
 
C         Amplitude(s) for UVCT diagram with ID 16
 
637
C         Amplitude(s) for UVCT diagram with ID 13
632
638
          CALL FFV1_0(W(1,4,H),W(1,3,H),W(1,5,H),GC_5,AMPL(1,28))
633
639
          AMPL(1,28)=AMPL(1,28)*(2.0D0*UVWFCT_T_0)
634
 
C         Amplitude(s) for UVCT diagram with ID 17
 
640
C         Amplitude(s) for UVCT diagram with ID 14
635
641
          CALL FFV1_0(W(1,4,H),W(1,3,H),W(1,5,H),GC_5,AMPL(2,29))
636
642
          AMPL(2,29)=AMPL(2,29)*(2.0D0*UVWFCT_B_0_1EPS)
637
643
 300      CONTINUE
664
670
C     Loop amplitude for loop diagram with ID 2
665
671
      CALL ML5_0_LOOP_2_2(1,5,6,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
666
672
     $ ,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),GC_5,MP__GC_5,GC_5
667
 
     $ ,MP__GC_5,2,1,30,AMPL(1,30),S(30))
 
673
     $ ,MP__GC_5,2,1,4,30,AMPL(1,30),S(30))
668
674
C     Loop amplitude for loop diagram with ID 3
669
675
      CALL ML5_0_LOOP_4_4(2,1,2,4,3,DCMPLX(ZERO),CMPLX(MP__ZERO
670
676
     $ ,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MT)
671
677
     $ ,CMPLX(MP__MDL_MT,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
672
 
     $ ,GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5,MP__GC_5,2,1,31
673
 
     $ ,AMPL(1,31),S(31))
 
678
     $ ,GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5,MP__GC_5,2,1,1
 
679
     $ ,31,AMPL(1,31),S(31))
674
680
C     Loop amplitude for loop diagram with ID 4
675
681
      CALL ML5_0_LOOP_4_4(3,1,2,3,4,DCMPLX(ZERO),CMPLX(MP__ZERO
676
682
     $ ,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(MDL_MT)
677
683
     $ ,CMPLX(MP__MDL_MT,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
678
 
     $ ,GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5,MP__GC_5,2,1,32
679
 
     $ ,AMPL(1,32),S(32))
 
684
     $ ,GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5,MP__GC_5,2,1,1
 
685
     $ ,32,AMPL(1,32),S(32))
680
686
C     Loop amplitude for loop diagram with ID 5
681
687
      CALL ML5_0_LOOP_3_3(4,1,2,6,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
682
688
     $ ,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZER
683
 
     $ O,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5,GC_4,MP__GC_4,2,1,33
 
689
     $ O,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5,GC_4,MP__GC_4,2,1,1,33
684
690
     $ ,AMPL(1,33),S(33))
685
691
C     Loop amplitude for loop diagram with ID 6
686
692
      CALL ML5_0_LOOP_3_3(5,1,2,6,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
687
693
     $ ,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZER
688
 
     $ O,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5,MP__GC_5,2,1,34
 
694
     $ O,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5,MP__GC_5,2,1,1,34
689
695
     $ ,AMPL(1,34),S(34))
690
696
C     Loop amplitude for loop diagram with ID 7
691
 
      CALL ML5_0_LOOP_2_2(1,5,6,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
692
 
     $ ,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),GC_5,MP__GC_5,GC_5
693
 
     $ ,MP__GC_5,2,1,35,AMPL(1,35),S(35))
694
 
C     Loop amplitude for loop diagram with ID 8
695
 
      CALL ML5_0_LOOP_2_2(1,5,6,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
696
 
     $ ,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),GC_5,MP__GC_5,GC_5
697
 
     $ ,MP__GC_5,2,1,36,AMPL(1,36),S(36))
698
 
C     Loop amplitude for loop diagram with ID 9
699
 
      CALL ML5_0_LOOP_2_2(1,5,6,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
700
 
     $ ,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),GC_5,MP__GC_5,GC_5
701
 
     $ ,MP__GC_5,2,1,37,AMPL(1,37),S(37))
702
 
C     Loop amplitude for loop diagram with ID 10
703
697
      CALL ML5_0_LOOP_2_2(1,5,6,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB
704
698
     $ ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),GC_5
705
 
     $ ,MP__GC_5,GC_5,MP__GC_5,2,1,38,AMPL(1,38),S(38))
706
 
C     Loop amplitude for loop diagram with ID 11
 
699
     $ ,MP__GC_5,GC_5,MP__GC_5,2,1,1,35,AMPL(1,35),S(35))
 
700
C     Loop amplitude for loop diagram with ID 8
707
701
      CALL ML5_0_LOOP_2_2(1,5,6,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT
708
702
     $ ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),GC_5
709
 
     $ ,MP__GC_5,GC_5,MP__GC_5,2,1,39,AMPL(1,39),S(39))
710
 
C     Loop amplitude for loop diagram with ID 12
 
703
     $ ,MP__GC_5,GC_5,MP__GC_5,2,1,1,36,AMPL(1,36),S(36))
 
704
C     Loop amplitude for loop diagram with ID 9
711
705
      CALL ML5_0_LOOP_3_3(6,3,4,5,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
712
706
     $ ,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(MDL_MT)
713
707
     $ ,CMPLX(MP__MDL_MT,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5,GC_5
714
 
     $ ,MP__GC_5,2,1,40,AMPL(1,40),S(40))
715
 
C     Loop amplitude for loop diagram with ID 13
 
708
     $ ,MP__GC_5,2,1,1,37,AMPL(1,37),S(37))
 
709
C     Loop amplitude for loop diagram with ID 10
716
710
      CALL ML5_0_LOOP_3_3(7,3,4,5,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT
717
711
     $ ,KIND=16),DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),DCMPLX(ZERO)
718
712
     $ ,CMPLX(MP__ZERO,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5,GC_4
719
 
     $ ,MP__GC_4,2,1,41,AMPL(1,41),S(41))
720
 
C     Loop amplitude for loop diagram with ID 14
 
713
     $ ,MP__GC_4,2,1,1,38,AMPL(1,38),S(38))
 
714
C     Loop amplitude for loop diagram with ID 11
721
715
      CALL ML5_0_LOOP_2_2(8,5,6,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
722
716
     $ ,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),GC_4,MP__GC_4,GC_4
723
 
     $ ,MP__GC_4,2,2,42,AMPL(1,42),S(42))
724
 
C     Loop amplitude for loop diagram with ID 15
 
717
     $ ,MP__GC_4,2,2,1,39,AMPL(1,39),S(39))
 
718
C     Loop amplitude for loop diagram with ID 12
725
719
      CALL ML5_0_LOOP_2_2(9,5,6,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16)
726
720
     $ ,DCMPLX(ZERO),CMPLX(MP__ZERO,KIND=16),GC_4,MP__GC_4,GC_4
727
 
     $ ,MP__GC_4,2,1,43,AMPL(1,43),S(43))
 
721
     $ ,MP__GC_4,2,1,1,40,AMPL(1,40),S(40))
728
722
 
729
723
      DO I=NCTAMPS+1,NLOOPAMPS
730
724
        ANS(1)=ANS(1)+AMPL(1,I)