← 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/CT_interface.f

  • Committer: olivier-mattelaer
  • Date: 2020-08-21 09:16:56 UTC
  • mfrom: (284.2.24 python3)
  • Revision ID: olivier-mattelaer-20200821091656-iizux2mj94tkssuo
pass to 2.8.0 (and move to python3 branch)

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/long_ML_SMQCD_default/gg_wmtbx/CT_interface.f
547
547
 
548
548
      END
549
549
 
550
 
      SUBROUTINE ML5_0_LOOP_4_5_4( LID, P1, P2, P3, P4, W1, W2, W3, W4
551
 
     $ , W5, M1,MP_M1, M2,MP_M2, M3,MP_M3, M4,MP_M4, C1,MP_C1, C2
552
 
     $ ,MP_C2, C3,MP_C3, C4,MP_C4,  RANK, LSYMFACT, LMULTIPLIER, AMPLN
553
 
     $ , RES, STABLE)
554
 
 
555
 
      INTEGER    NEXTERNAL
556
 
      PARAMETER (NEXTERNAL=5)
557
 
      INTEGER    MAXLCOUPLINGS
558
 
      PARAMETER (MAXLCOUPLINGS=5)
559
 
      INTEGER    NLOOPLINE
560
 
      PARAMETER (NLOOPLINE=4)
561
 
      INTEGER    NWAVEFUNCS
562
 
      PARAMETER (NWAVEFUNCS=28)
563
 
      INTEGER    NCOMB
564
 
      PARAMETER (NCOMB=48)
565
 
C     
566
 
C     ARGUMENTS 
567
 
C     
568
 
      INTEGER W1, W2, W3, W4, W5
569
 
      COMPLEX*16 M1, M2, M3, M4
570
 
      COMPLEX*32 MP_M1, MP_M2, MP_M3, MP_M4
571
 
      COMPLEX*16 C1, C2, C3, C4
572
 
      COMPLEX*32 MP_C1, MP_C2, MP_C3, MP_C4
573
 
      INTEGER P1, P2, P3, P4
574
 
      COMPLEX*16 RES(3)
575
 
      INTEGER LID, RANK, LSYMFACT, LMULTIPLIER
576
 
      INTEGER AMPLN
577
 
      LOGICAL STABLE
578
 
C     
579
 
C     LOCAL VARIABLES 
580
 
C     
581
 
      REAL*8 PL(0:3,NLOOPLINE)
582
 
      COMPLEX*16 M2L(NLOOPLINE)
583
 
      INTEGER PAIRING(NLOOPLINE)
584
 
      INTEGER I, J, K, TEMP
585
 
C     
586
 
C     GLOBAL VARIABLES
587
 
C     
588
 
      INTEGER WE(NEXTERNAL)
589
 
      INTEGER ID, SYMFACT, MULTIPLIER, AMPLNUM
590
 
      COMMON/ML5_0_LOOP/WE,ID, SYMFACT, MULTIPLIER,AMPLNUM
591
 
 
592
 
      COMPLEX*16 LC(MAXLCOUPLINGS)
593
 
      COMPLEX*16 ML(NEXTERNAL+2)
594
 
      COMMON/ML5_0_DP_LOOP/LC,ML
595
 
 
596
 
      COMPLEX*32 MP_LC(MAXLCOUPLINGS)
597
 
      COMPLEX*32 MP_ML(NEXTERNAL+2)
598
 
      COMMON/ML5_0_MP_LOOP/MP_LC,MP_ML
599
 
 
600
 
      COMPLEX*16 W(20,NWAVEFUNCS,NCOMB)
601
 
      INTEGER VALIDH
602
 
      COMMON/ML5_0_WFCTS/W
603
 
      COMMON/ML5_0_VALIDH/VALIDH
604
 
 
605
 
C     ----------
606
 
C     BEGIN CODE
607
 
C     ----------
608
 
 
609
 
      WE(1)=W1
610
 
      WE(2)=W2
611
 
      WE(3)=W3
612
 
      WE(4)=W4
613
 
      WE(5)=W5
614
 
      M2L(1)=M4**2
615
 
      M2L(2)=M1**2
616
 
      M2L(3)=M2**2
617
 
      M2L(4)=M3**2
618
 
      ML(1)=M4
619
 
      ML(2)=M4
620
 
      MP_ML(1)=MP_M4
621
 
      MP_ML(2)=MP_M4
622
 
      ML(3)=M1
623
 
      MP_ML(3)=MP_M1
624
 
      ML(4)=M2
625
 
      MP_ML(4)=MP_M2
626
 
      ML(5)=M3
627
 
      MP_ML(5)=MP_M3
628
 
      ML(6)=M4
629
 
      MP_ML(6)=MP_M4
630
 
      PAIRING(1)=P1
631
 
      PAIRING(2)=P2
632
 
      PAIRING(3)=P3
633
 
      PAIRING(4)=P4
634
 
      LC(1)=C1
635
 
      MP_LC(1)=MP_C1
636
 
      LC(2)=C2
637
 
      MP_LC(2)=MP_C2
638
 
      LC(3)=C3
639
 
      MP_LC(3)=MP_C3
640
 
      LC(4)=C4
641
 
      MP_LC(4)=MP_C4
642
 
      AMPLNUM=AMPLN
643
 
      ID=LID
644
 
      SYMFACT=LSYMFACT
645
 
      MULTIPLIER=LMULTIPLIER
646
 
      DO I=0,3
647
 
        TEMP=1
648
 
        DO J=1,NLOOPLINE
649
 
          PL(I,J)=0.D0
650
 
          DO K=TEMP,(TEMP+PAIRING(J)-1)
651
 
            PL(I,J)=PL(I,J)-DBLE(W(1+I,WE(K),VALIDH))
652
 
          ENDDO
653
 
          TEMP=TEMP+PAIRING(J)
654
 
        ENDDO
655
 
      ENDDO
656
 
      CALL ML5_0_CTLOOP(NLOOPLINE,PL,M2L,RANK,RES,STABLE)
657
 
 
658
 
      END
659
 
 
660
 
      SUBROUTINE ML5_0_LOOP_4_4( LID, W1, W2, W3, W4, M1,MP_M1, M2
661
 
     $ ,MP_M2, M3,MP_M3, M4,MP_M4, C1,MP_C1, C2,MP_C2, C3,MP_C3, C4
662
 
     $ ,MP_C4,  RANK, LSYMFACT, LMULTIPLIER, AMPLN, RES, STABLE)
663
 
 
664
 
      INTEGER    NEXTERNAL
665
 
      PARAMETER (NEXTERNAL=5)
666
 
      INTEGER    MAXLCOUPLINGS
667
 
      PARAMETER (MAXLCOUPLINGS=5)
668
 
      INTEGER    NLOOPLINE
669
 
      PARAMETER (NLOOPLINE=4)
670
 
      INTEGER    NWAVEFUNCS
671
 
      PARAMETER (NWAVEFUNCS=28)
672
 
      INTEGER    NCOMB
673
 
      PARAMETER (NCOMB=48)
674
 
C     
675
 
C     ARGUMENTS 
676
 
C     
677
 
      INTEGER W1, W2, W3, W4
678
 
      COMPLEX*16 M1, M2, M3, M4
679
 
      COMPLEX*32 MP_M1, MP_M2, MP_M3, MP_M4
680
 
      COMPLEX*16 C1, C2, C3, C4
681
 
      COMPLEX*32 MP_C1, MP_C2, MP_C3, MP_C4
682
 
 
683
 
      COMPLEX*16 RES(3)
684
 
      INTEGER LID, RANK, LSYMFACT, LMULTIPLIER
685
 
      INTEGER AMPLN
686
 
      LOGICAL STABLE
687
 
C     
688
 
C     LOCAL VARIABLES 
689
 
C     
690
 
      REAL*8 PL(0:3,NLOOPLINE)
691
 
      COMPLEX*16 M2L(NLOOPLINE)
692
 
      INTEGER PAIRING(NLOOPLINE)
693
 
      INTEGER I, J, K, TEMP
694
 
C     
695
 
C     GLOBAL VARIABLES
696
 
C     
697
 
      INTEGER WE(NEXTERNAL)
698
 
      INTEGER ID, SYMFACT, MULTIPLIER, AMPLNUM
699
 
      COMMON/ML5_0_LOOP/WE,ID, SYMFACT, MULTIPLIER,AMPLNUM
700
 
 
701
 
      COMPLEX*16 LC(MAXLCOUPLINGS)
702
 
      COMPLEX*16 ML(NEXTERNAL+2)
703
 
      COMMON/ML5_0_DP_LOOP/LC,ML
704
 
 
705
 
      COMPLEX*32 MP_LC(MAXLCOUPLINGS)
706
 
      COMPLEX*32 MP_ML(NEXTERNAL+2)
707
 
      COMMON/ML5_0_MP_LOOP/MP_LC,MP_ML
708
 
 
709
 
      COMPLEX*16 W(20,NWAVEFUNCS,NCOMB)
710
 
      INTEGER VALIDH
711
 
      COMMON/ML5_0_WFCTS/W
712
 
      COMMON/ML5_0_VALIDH/VALIDH
713
 
 
714
 
C     ----------
715
 
C     BEGIN CODE
716
 
C     ----------
717
 
 
718
 
      WE(1)=W1
719
 
      WE(2)=W2
720
 
      WE(3)=W3
721
 
      WE(4)=W4
722
 
      M2L(1)=M4**2
723
 
      M2L(2)=M1**2
724
 
      M2L(3)=M2**2
725
 
      M2L(4)=M3**2
726
 
      ML(1)=M4
727
 
      ML(2)=M4
728
 
      MP_ML(1)=MP_M4
729
 
      MP_ML(2)=MP_M4
730
 
      ML(3)=M1
731
 
      MP_ML(3)=MP_M1
732
 
      ML(4)=M2
733
 
      MP_ML(4)=MP_M2
734
 
      ML(5)=M3
735
 
      MP_ML(5)=MP_M3
736
 
      ML(6)=M4
737
 
      MP_ML(6)=MP_M4
738
 
      DO I=1,NLOOPLINE
739
 
        PAIRING(I)=1
740
 
      ENDDO
741
 
 
742
 
      LC(1)=C1
743
 
      MP_LC(1)=MP_C1
744
 
      LC(2)=C2
745
 
      MP_LC(2)=MP_C2
746
 
      LC(3)=C3
747
 
      MP_LC(3)=MP_C3
748
 
      LC(4)=C4
749
 
      MP_LC(4)=MP_C4
750
 
      AMPLNUM=AMPLN
751
 
      ID=LID
752
 
      SYMFACT=LSYMFACT
753
 
      MULTIPLIER=LMULTIPLIER
754
 
      DO I=0,3
755
 
        TEMP=1
756
 
        DO J=1,NLOOPLINE
757
 
          PL(I,J)=0.D0
758
 
          DO K=TEMP,(TEMP+PAIRING(J)-1)
759
 
            PL(I,J)=PL(I,J)-DBLE(W(1+I,WE(K),VALIDH))
760
 
          ENDDO
761
 
          TEMP=TEMP+PAIRING(J)
762
 
        ENDDO
763
 
      ENDDO
764
 
      CALL ML5_0_CTLOOP(NLOOPLINE,PL,M2L,RANK,RES,STABLE)
765
 
 
766
 
      END
767
 
 
768
550
      SUBROUTINE ML5_0_LOOP_5_5( LID, W1, W2, W3, W4, W5, M1,MP_M1, M2
769
551
     $ ,MP_M2, M3,MP_M3, M4,MP_M4, M5,MP_M5, C1,MP_C1, C2,MP_C2, C3
770
552
     $ ,MP_C3, C4,MP_C4, C5,MP_C5,  RANK, LSYMFACT, LMULTIPLIER, AMPLN
880
662
 
881
663
      END
882
664
 
 
665
      SUBROUTINE ML5_0_LOOP_4_4( LID, W1, W2, W3, W4, M1,MP_M1, M2
 
666
     $ ,MP_M2, M3,MP_M3, M4,MP_M4, C1,MP_C1, C2,MP_C2, C3,MP_C3, C4
 
667
     $ ,MP_C4,  RANK, LSYMFACT, LMULTIPLIER, AMPLN, RES, STABLE)
 
668
 
 
669
      INTEGER    NEXTERNAL
 
670
      PARAMETER (NEXTERNAL=5)
 
671
      INTEGER    MAXLCOUPLINGS
 
672
      PARAMETER (MAXLCOUPLINGS=5)
 
673
      INTEGER    NLOOPLINE
 
674
      PARAMETER (NLOOPLINE=4)
 
675
      INTEGER    NWAVEFUNCS
 
676
      PARAMETER (NWAVEFUNCS=28)
 
677
      INTEGER    NCOMB
 
678
      PARAMETER (NCOMB=48)
 
679
C     
 
680
C     ARGUMENTS 
 
681
C     
 
682
      INTEGER W1, W2, W3, W4
 
683
      COMPLEX*16 M1, M2, M3, M4
 
684
      COMPLEX*32 MP_M1, MP_M2, MP_M3, MP_M4
 
685
      COMPLEX*16 C1, C2, C3, C4
 
686
      COMPLEX*32 MP_C1, MP_C2, MP_C3, MP_C4
 
687
 
 
688
      COMPLEX*16 RES(3)
 
689
      INTEGER LID, RANK, LSYMFACT, LMULTIPLIER
 
690
      INTEGER AMPLN
 
691
      LOGICAL STABLE
 
692
C     
 
693
C     LOCAL VARIABLES 
 
694
C     
 
695
      REAL*8 PL(0:3,NLOOPLINE)
 
696
      COMPLEX*16 M2L(NLOOPLINE)
 
697
      INTEGER PAIRING(NLOOPLINE)
 
698
      INTEGER I, J, K, TEMP
 
699
C     
 
700
C     GLOBAL VARIABLES
 
701
C     
 
702
      INTEGER WE(NEXTERNAL)
 
703
      INTEGER ID, SYMFACT, MULTIPLIER, AMPLNUM
 
704
      COMMON/ML5_0_LOOP/WE,ID, SYMFACT, MULTIPLIER,AMPLNUM
 
705
 
 
706
      COMPLEX*16 LC(MAXLCOUPLINGS)
 
707
      COMPLEX*16 ML(NEXTERNAL+2)
 
708
      COMMON/ML5_0_DP_LOOP/LC,ML
 
709
 
 
710
      COMPLEX*32 MP_LC(MAXLCOUPLINGS)
 
711
      COMPLEX*32 MP_ML(NEXTERNAL+2)
 
712
      COMMON/ML5_0_MP_LOOP/MP_LC,MP_ML
 
713
 
 
714
      COMPLEX*16 W(20,NWAVEFUNCS,NCOMB)
 
715
      INTEGER VALIDH
 
716
      COMMON/ML5_0_WFCTS/W
 
717
      COMMON/ML5_0_VALIDH/VALIDH
 
718
 
 
719
C     ----------
 
720
C     BEGIN CODE
 
721
C     ----------
 
722
 
 
723
      WE(1)=W1
 
724
      WE(2)=W2
 
725
      WE(3)=W3
 
726
      WE(4)=W4
 
727
      M2L(1)=M4**2
 
728
      M2L(2)=M1**2
 
729
      M2L(3)=M2**2
 
730
      M2L(4)=M3**2
 
731
      ML(1)=M4
 
732
      ML(2)=M4
 
733
      MP_ML(1)=MP_M4
 
734
      MP_ML(2)=MP_M4
 
735
      ML(3)=M1
 
736
      MP_ML(3)=MP_M1
 
737
      ML(4)=M2
 
738
      MP_ML(4)=MP_M2
 
739
      ML(5)=M3
 
740
      MP_ML(5)=MP_M3
 
741
      ML(6)=M4
 
742
      MP_ML(6)=MP_M4
 
743
      DO I=1,NLOOPLINE
 
744
        PAIRING(I)=1
 
745
      ENDDO
 
746
 
 
747
      LC(1)=C1
 
748
      MP_LC(1)=MP_C1
 
749
      LC(2)=C2
 
750
      MP_LC(2)=MP_C2
 
751
      LC(3)=C3
 
752
      MP_LC(3)=MP_C3
 
753
      LC(4)=C4
 
754
      MP_LC(4)=MP_C4
 
755
      AMPLNUM=AMPLN
 
756
      ID=LID
 
757
      SYMFACT=LSYMFACT
 
758
      MULTIPLIER=LMULTIPLIER
 
759
      DO I=0,3
 
760
        TEMP=1
 
761
        DO J=1,NLOOPLINE
 
762
          PL(I,J)=0.D0
 
763
          DO K=TEMP,(TEMP+PAIRING(J)-1)
 
764
            PL(I,J)=PL(I,J)-DBLE(W(1+I,WE(K),VALIDH))
 
765
          ENDDO
 
766
          TEMP=TEMP+PAIRING(J)
 
767
        ENDDO
 
768
      ENDDO
 
769
      CALL ML5_0_CTLOOP(NLOOPLINE,PL,M2L,RANK,RES,STABLE)
 
770
 
 
771
      END
 
772
 
 
773
      SUBROUTINE ML5_0_LOOP_4_5_4( LID, P1, P2, P3, P4, W1, W2, W3, W4
 
774
     $ , W5, M1,MP_M1, M2,MP_M2, M3,MP_M3, M4,MP_M4, C1,MP_C1, C2
 
775
     $ ,MP_C2, C3,MP_C3, C4,MP_C4,  RANK, LSYMFACT, LMULTIPLIER, AMPLN
 
776
     $ , RES, STABLE)
 
777
 
 
778
      INTEGER    NEXTERNAL
 
779
      PARAMETER (NEXTERNAL=5)
 
780
      INTEGER    MAXLCOUPLINGS
 
781
      PARAMETER (MAXLCOUPLINGS=5)
 
782
      INTEGER    NLOOPLINE
 
783
      PARAMETER (NLOOPLINE=4)
 
784
      INTEGER    NWAVEFUNCS
 
785
      PARAMETER (NWAVEFUNCS=28)
 
786
      INTEGER    NCOMB
 
787
      PARAMETER (NCOMB=48)
 
788
C     
 
789
C     ARGUMENTS 
 
790
C     
 
791
      INTEGER W1, W2, W3, W4, W5
 
792
      COMPLEX*16 M1, M2, M3, M4
 
793
      COMPLEX*32 MP_M1, MP_M2, MP_M3, MP_M4
 
794
      COMPLEX*16 C1, C2, C3, C4
 
795
      COMPLEX*32 MP_C1, MP_C2, MP_C3, MP_C4
 
796
      INTEGER P1, P2, P3, P4
 
797
      COMPLEX*16 RES(3)
 
798
      INTEGER LID, RANK, LSYMFACT, LMULTIPLIER
 
799
      INTEGER AMPLN
 
800
      LOGICAL STABLE
 
801
C     
 
802
C     LOCAL VARIABLES 
 
803
C     
 
804
      REAL*8 PL(0:3,NLOOPLINE)
 
805
      COMPLEX*16 M2L(NLOOPLINE)
 
806
      INTEGER PAIRING(NLOOPLINE)
 
807
      INTEGER I, J, K, TEMP
 
808
C     
 
809
C     GLOBAL VARIABLES
 
810
C     
 
811
      INTEGER WE(NEXTERNAL)
 
812
      INTEGER ID, SYMFACT, MULTIPLIER, AMPLNUM
 
813
      COMMON/ML5_0_LOOP/WE,ID, SYMFACT, MULTIPLIER,AMPLNUM
 
814
 
 
815
      COMPLEX*16 LC(MAXLCOUPLINGS)
 
816
      COMPLEX*16 ML(NEXTERNAL+2)
 
817
      COMMON/ML5_0_DP_LOOP/LC,ML
 
818
 
 
819
      COMPLEX*32 MP_LC(MAXLCOUPLINGS)
 
820
      COMPLEX*32 MP_ML(NEXTERNAL+2)
 
821
      COMMON/ML5_0_MP_LOOP/MP_LC,MP_ML
 
822
 
 
823
      COMPLEX*16 W(20,NWAVEFUNCS,NCOMB)
 
824
      INTEGER VALIDH
 
825
      COMMON/ML5_0_WFCTS/W
 
826
      COMMON/ML5_0_VALIDH/VALIDH
 
827
 
 
828
C     ----------
 
829
C     BEGIN CODE
 
830
C     ----------
 
831
 
 
832
      WE(1)=W1
 
833
      WE(2)=W2
 
834
      WE(3)=W3
 
835
      WE(4)=W4
 
836
      WE(5)=W5
 
837
      M2L(1)=M4**2
 
838
      M2L(2)=M1**2
 
839
      M2L(3)=M2**2
 
840
      M2L(4)=M3**2
 
841
      ML(1)=M4
 
842
      ML(2)=M4
 
843
      MP_ML(1)=MP_M4
 
844
      MP_ML(2)=MP_M4
 
845
      ML(3)=M1
 
846
      MP_ML(3)=MP_M1
 
847
      ML(4)=M2
 
848
      MP_ML(4)=MP_M2
 
849
      ML(5)=M3
 
850
      MP_ML(5)=MP_M3
 
851
      ML(6)=M4
 
852
      MP_ML(6)=MP_M4
 
853
      PAIRING(1)=P1
 
854
      PAIRING(2)=P2
 
855
      PAIRING(3)=P3
 
856
      PAIRING(4)=P4
 
857
      LC(1)=C1
 
858
      MP_LC(1)=MP_C1
 
859
      LC(2)=C2
 
860
      MP_LC(2)=MP_C2
 
861
      LC(3)=C3
 
862
      MP_LC(3)=MP_C3
 
863
      LC(4)=C4
 
864
      MP_LC(4)=MP_C4
 
865
      AMPLNUM=AMPLN
 
866
      ID=LID
 
867
      SYMFACT=LSYMFACT
 
868
      MULTIPLIER=LMULTIPLIER
 
869
      DO I=0,3
 
870
        TEMP=1
 
871
        DO J=1,NLOOPLINE
 
872
          PL(I,J)=0.D0
 
873
          DO K=TEMP,(TEMP+PAIRING(J)-1)
 
874
            PL(I,J)=PL(I,J)-DBLE(W(1+I,WE(K),VALIDH))
 
875
          ENDDO
 
876
          TEMP=TEMP+PAIRING(J)
 
877
        ENDDO
 
878
      ENDDO
 
879
      CALL ML5_0_CTLOOP(NLOOPLINE,PL,M2L,RANK,RES,STABLE)
 
880
 
 
881
      END
 
882