← Back to branch summary

~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/SquaredOrder_IOTest/sqso_uux_uuxuuxx/matrix_QCDsq_le_6.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/SquaredOrder_IOTest/sqso_uux_uuxuuxx/matrix_QCDsq_le_6.f
98
98
      LOGICAL GOODHEL(NCOMB)
99
99
      DATA NTRY/0/
100
100
      DATA GOODHEL/NCOMB*.FALSE./
101
 
      DATA (NHEL(I,   1),I=1,6) /-1,-1,-1,-1,-1,-1/
102
 
      DATA (NHEL(I,   2),I=1,6) /-1,-1,-1,-1,-1, 1/
103
 
      DATA (NHEL(I,   3),I=1,6) /-1,-1,-1,-1, 1,-1/
104
 
      DATA (NHEL(I,   4),I=1,6) /-1,-1,-1,-1, 1, 1/
105
 
      DATA (NHEL(I,   5),I=1,6) /-1,-1,-1, 1,-1,-1/
106
 
      DATA (NHEL(I,   6),I=1,6) /-1,-1,-1, 1,-1, 1/
107
 
      DATA (NHEL(I,   7),I=1,6) /-1,-1,-1, 1, 1,-1/
108
 
      DATA (NHEL(I,   8),I=1,6) /-1,-1,-1, 1, 1, 1/
109
 
      DATA (NHEL(I,   9),I=1,6) /-1,-1, 1,-1,-1,-1/
110
 
      DATA (NHEL(I,  10),I=1,6) /-1,-1, 1,-1,-1, 1/
111
 
      DATA (NHEL(I,  11),I=1,6) /-1,-1, 1,-1, 1,-1/
112
 
      DATA (NHEL(I,  12),I=1,6) /-1,-1, 1,-1, 1, 1/
113
 
      DATA (NHEL(I,  13),I=1,6) /-1,-1, 1, 1,-1,-1/
114
 
      DATA (NHEL(I,  14),I=1,6) /-1,-1, 1, 1,-1, 1/
115
 
      DATA (NHEL(I,  15),I=1,6) /-1,-1, 1, 1, 1,-1/
116
 
      DATA (NHEL(I,  16),I=1,6) /-1,-1, 1, 1, 1, 1/
117
 
      DATA (NHEL(I,  17),I=1,6) /-1, 1,-1,-1,-1,-1/
118
 
      DATA (NHEL(I,  18),I=1,6) /-1, 1,-1,-1,-1, 1/
119
 
      DATA (NHEL(I,  19),I=1,6) /-1, 1,-1,-1, 1,-1/
120
 
      DATA (NHEL(I,  20),I=1,6) /-1, 1,-1,-1, 1, 1/
121
 
      DATA (NHEL(I,  21),I=1,6) /-1, 1,-1, 1,-1,-1/
122
 
      DATA (NHEL(I,  22),I=1,6) /-1, 1,-1, 1,-1, 1/
123
 
      DATA (NHEL(I,  23),I=1,6) /-1, 1,-1, 1, 1,-1/
124
 
      DATA (NHEL(I,  24),I=1,6) /-1, 1,-1, 1, 1, 1/
125
 
      DATA (NHEL(I,  25),I=1,6) /-1, 1, 1,-1,-1,-1/
126
 
      DATA (NHEL(I,  26),I=1,6) /-1, 1, 1,-1,-1, 1/
127
 
      DATA (NHEL(I,  27),I=1,6) /-1, 1, 1,-1, 1,-1/
128
 
      DATA (NHEL(I,  28),I=1,6) /-1, 1, 1,-1, 1, 1/
129
 
      DATA (NHEL(I,  29),I=1,6) /-1, 1, 1, 1,-1,-1/
130
 
      DATA (NHEL(I,  30),I=1,6) /-1, 1, 1, 1,-1, 1/
131
 
      DATA (NHEL(I,  31),I=1,6) /-1, 1, 1, 1, 1,-1/
132
 
      DATA (NHEL(I,  32),I=1,6) /-1, 1, 1, 1, 1, 1/
133
 
      DATA (NHEL(I,  33),I=1,6) / 1,-1,-1,-1,-1,-1/
134
 
      DATA (NHEL(I,  34),I=1,6) / 1,-1,-1,-1,-1, 1/
135
 
      DATA (NHEL(I,  35),I=1,6) / 1,-1,-1,-1, 1,-1/
136
 
      DATA (NHEL(I,  36),I=1,6) / 1,-1,-1,-1, 1, 1/
137
 
      DATA (NHEL(I,  37),I=1,6) / 1,-1,-1, 1,-1,-1/
138
 
      DATA (NHEL(I,  38),I=1,6) / 1,-1,-1, 1,-1, 1/
139
 
      DATA (NHEL(I,  39),I=1,6) / 1,-1,-1, 1, 1,-1/
140
 
      DATA (NHEL(I,  40),I=1,6) / 1,-1,-1, 1, 1, 1/
141
 
      DATA (NHEL(I,  41),I=1,6) / 1,-1, 1,-1,-1,-1/
142
 
      DATA (NHEL(I,  42),I=1,6) / 1,-1, 1,-1,-1, 1/
143
 
      DATA (NHEL(I,  43),I=1,6) / 1,-1, 1,-1, 1,-1/
144
 
      DATA (NHEL(I,  44),I=1,6) / 1,-1, 1,-1, 1, 1/
145
 
      DATA (NHEL(I,  45),I=1,6) / 1,-1, 1, 1,-1,-1/
146
 
      DATA (NHEL(I,  46),I=1,6) / 1,-1, 1, 1,-1, 1/
147
 
      DATA (NHEL(I,  47),I=1,6) / 1,-1, 1, 1, 1,-1/
148
 
      DATA (NHEL(I,  48),I=1,6) / 1,-1, 1, 1, 1, 1/
149
 
      DATA (NHEL(I,  49),I=1,6) / 1, 1,-1,-1,-1,-1/
150
 
      DATA (NHEL(I,  50),I=1,6) / 1, 1,-1,-1,-1, 1/
151
 
      DATA (NHEL(I,  51),I=1,6) / 1, 1,-1,-1, 1,-1/
152
 
      DATA (NHEL(I,  52),I=1,6) / 1, 1,-1,-1, 1, 1/
153
 
      DATA (NHEL(I,  53),I=1,6) / 1, 1,-1, 1,-1,-1/
154
 
      DATA (NHEL(I,  54),I=1,6) / 1, 1,-1, 1,-1, 1/
155
 
      DATA (NHEL(I,  55),I=1,6) / 1, 1,-1, 1, 1,-1/
156
 
      DATA (NHEL(I,  56),I=1,6) / 1, 1,-1, 1, 1, 1/
157
 
      DATA (NHEL(I,  57),I=1,6) / 1, 1, 1,-1,-1,-1/
158
 
      DATA (NHEL(I,  58),I=1,6) / 1, 1, 1,-1,-1, 1/
159
 
      DATA (NHEL(I,  59),I=1,6) / 1, 1, 1,-1, 1,-1/
160
 
      DATA (NHEL(I,  60),I=1,6) / 1, 1, 1,-1, 1, 1/
161
 
      DATA (NHEL(I,  61),I=1,6) / 1, 1, 1, 1,-1,-1/
162
 
      DATA (NHEL(I,  62),I=1,6) / 1, 1, 1, 1,-1, 1/
163
 
      DATA (NHEL(I,  63),I=1,6) / 1, 1, 1, 1, 1,-1/
164
 
      DATA (NHEL(I,  64),I=1,6) / 1, 1, 1, 1, 1, 1/
 
101
      DATA (NHEL(I,   1),I=1,6) / 1,-1,-1, 1,-1, 1/
 
102
      DATA (NHEL(I,   2),I=1,6) / 1,-1,-1, 1,-1,-1/
 
103
      DATA (NHEL(I,   3),I=1,6) / 1,-1,-1, 1, 1, 1/
 
104
      DATA (NHEL(I,   4),I=1,6) / 1,-1,-1, 1, 1,-1/
 
105
      DATA (NHEL(I,   5),I=1,6) / 1,-1,-1,-1,-1, 1/
 
106
      DATA (NHEL(I,   6),I=1,6) / 1,-1,-1,-1,-1,-1/
 
107
      DATA (NHEL(I,   7),I=1,6) / 1,-1,-1,-1, 1, 1/
 
108
      DATA (NHEL(I,   8),I=1,6) / 1,-1,-1,-1, 1,-1/
 
109
      DATA (NHEL(I,   9),I=1,6) / 1,-1, 1, 1,-1, 1/
 
110
      DATA (NHEL(I,  10),I=1,6) / 1,-1, 1, 1,-1,-1/
 
111
      DATA (NHEL(I,  11),I=1,6) / 1,-1, 1, 1, 1, 1/
 
112
      DATA (NHEL(I,  12),I=1,6) / 1,-1, 1, 1, 1,-1/
 
113
      DATA (NHEL(I,  13),I=1,6) / 1,-1, 1,-1,-1, 1/
 
114
      DATA (NHEL(I,  14),I=1,6) / 1,-1, 1,-1,-1,-1/
 
115
      DATA (NHEL(I,  15),I=1,6) / 1,-1, 1,-1, 1, 1/
 
116
      DATA (NHEL(I,  16),I=1,6) / 1,-1, 1,-1, 1,-1/
 
117
      DATA (NHEL(I,  17),I=1,6) / 1, 1,-1, 1,-1, 1/
 
118
      DATA (NHEL(I,  18),I=1,6) / 1, 1,-1, 1,-1,-1/
 
119
      DATA (NHEL(I,  19),I=1,6) / 1, 1,-1, 1, 1, 1/
 
120
      DATA (NHEL(I,  20),I=1,6) / 1, 1,-1, 1, 1,-1/
 
121
      DATA (NHEL(I,  21),I=1,6) / 1, 1,-1,-1,-1, 1/
 
122
      DATA (NHEL(I,  22),I=1,6) / 1, 1,-1,-1,-1,-1/
 
123
      DATA (NHEL(I,  23),I=1,6) / 1, 1,-1,-1, 1, 1/
 
124
      DATA (NHEL(I,  24),I=1,6) / 1, 1,-1,-1, 1,-1/
 
125
      DATA (NHEL(I,  25),I=1,6) / 1, 1, 1, 1,-1, 1/
 
126
      DATA (NHEL(I,  26),I=1,6) / 1, 1, 1, 1,-1,-1/
 
127
      DATA (NHEL(I,  27),I=1,6) / 1, 1, 1, 1, 1, 1/
 
128
      DATA (NHEL(I,  28),I=1,6) / 1, 1, 1, 1, 1,-1/
 
129
      DATA (NHEL(I,  29),I=1,6) / 1, 1, 1,-1,-1, 1/
 
130
      DATA (NHEL(I,  30),I=1,6) / 1, 1, 1,-1,-1,-1/
 
131
      DATA (NHEL(I,  31),I=1,6) / 1, 1, 1,-1, 1, 1/
 
132
      DATA (NHEL(I,  32),I=1,6) / 1, 1, 1,-1, 1,-1/
 
133
      DATA (NHEL(I,  33),I=1,6) /-1,-1,-1, 1,-1, 1/
 
134
      DATA (NHEL(I,  34),I=1,6) /-1,-1,-1, 1,-1,-1/
 
135
      DATA (NHEL(I,  35),I=1,6) /-1,-1,-1, 1, 1, 1/
 
136
      DATA (NHEL(I,  36),I=1,6) /-1,-1,-1, 1, 1,-1/
 
137
      DATA (NHEL(I,  37),I=1,6) /-1,-1,-1,-1,-1, 1/
 
138
      DATA (NHEL(I,  38),I=1,6) /-1,-1,-1,-1,-1,-1/
 
139
      DATA (NHEL(I,  39),I=1,6) /-1,-1,-1,-1, 1, 1/
 
140
      DATA (NHEL(I,  40),I=1,6) /-1,-1,-1,-1, 1,-1/
 
141
      DATA (NHEL(I,  41),I=1,6) /-1,-1, 1, 1,-1, 1/
 
142
      DATA (NHEL(I,  42),I=1,6) /-1,-1, 1, 1,-1,-1/
 
143
      DATA (NHEL(I,  43),I=1,6) /-1,-1, 1, 1, 1, 1/
 
144
      DATA (NHEL(I,  44),I=1,6) /-1,-1, 1, 1, 1,-1/
 
145
      DATA (NHEL(I,  45),I=1,6) /-1,-1, 1,-1,-1, 1/
 
146
      DATA (NHEL(I,  46),I=1,6) /-1,-1, 1,-1,-1,-1/
 
147
      DATA (NHEL(I,  47),I=1,6) /-1,-1, 1,-1, 1, 1/
 
148
      DATA (NHEL(I,  48),I=1,6) /-1,-1, 1,-1, 1,-1/
 
149
      DATA (NHEL(I,  49),I=1,6) /-1, 1,-1, 1,-1, 1/
 
150
      DATA (NHEL(I,  50),I=1,6) /-1, 1,-1, 1,-1,-1/
 
151
      DATA (NHEL(I,  51),I=1,6) /-1, 1,-1, 1, 1, 1/
 
152
      DATA (NHEL(I,  52),I=1,6) /-1, 1,-1, 1, 1,-1/
 
153
      DATA (NHEL(I,  53),I=1,6) /-1, 1,-1,-1,-1, 1/
 
154
      DATA (NHEL(I,  54),I=1,6) /-1, 1,-1,-1,-1,-1/
 
155
      DATA (NHEL(I,  55),I=1,6) /-1, 1,-1,-1, 1, 1/
 
156
      DATA (NHEL(I,  56),I=1,6) /-1, 1,-1,-1, 1,-1/
 
157
      DATA (NHEL(I,  57),I=1,6) /-1, 1, 1, 1,-1, 1/
 
158
      DATA (NHEL(I,  58),I=1,6) /-1, 1, 1, 1,-1,-1/
 
159
      DATA (NHEL(I,  59),I=1,6) /-1, 1, 1, 1, 1, 1/
 
160
      DATA (NHEL(I,  60),I=1,6) /-1, 1, 1, 1, 1,-1/
 
161
      DATA (NHEL(I,  61),I=1,6) /-1, 1, 1,-1,-1, 1/
 
162
      DATA (NHEL(I,  62),I=1,6) /-1, 1, 1,-1,-1,-1/
 
163
      DATA (NHEL(I,  63),I=1,6) /-1, 1, 1,-1, 1, 1/
 
164
      DATA (NHEL(I,  64),I=1,6) /-1, 1, 1,-1, 1,-1/
165
165
      DATA IDEN/144/
166
166
C     
167
167
C     GLOBAL VARIABLES
521
521
      INTEGER I, SQORDERS(NSO)
522
522
      INTEGER AMPSPLITORDERS(NAMPSO,NSO)
523
523
      DATA (AMPSPLITORDERS(  1,I),I=  1,  2) /    4,    0/
 
524
      COMMON/AMPSPLITORDERS/AMPSPLITORDERS
524
525
C     
525
526
C     FUNCTION
526
527
C     
556
557
      INTEGER I,J
557
558
      INTEGER SQSPLITORDERS(NSQSO,NSO)
558
559
      DATA (SQSPLITORDERS(  1,I),I=  1,  2) /    8,    0/
 
560
      COMMON/SQPLITORDERS/SQPLITORDERS
559
561
C     
560
562
C     BEGIN CODE
561
563
C     
568
570
 1009   CONTINUE
569
571
      ENDDO
570
572
 
571
 
      WRITE(*,*) 'ERROR:: Stopping function SOINDEX_FOR_SQUARED_ORDERS'
 
573
      WRITE(*,*) 'ERROR:: Stopping in function'
 
574
      WRITE(*,*) 'SOINDEX_FOR_SQUARED_ORDERS'
572
575
      WRITE(*,*) 'Could not find squared orders ',(ORDERS(I),I=1,NSO)
573
576
      STOP
574
577
 
589
592
 
590
593
      END
591
594
 
 
595
C     This is the inverse subroutine of SOINDEX_FOR_SQUARED_ORDERS.
 
596
C      Not directly useful, but provided nonetheless.
 
597
      SUBROUTINE GET_SQUARED_ORDERS_FOR_SOINDEX(SOINDEX,ORDERS)
 
598
C     
 
599
C     This functions returns the orders identified by the squared
 
600
C      split order index in argument. Order values correspond to
 
601
C      following list of couplings (and in this order):
 
602
C     ['QCD', 'QED']
 
603
C     
 
604
C     CONSTANTS
 
605
C     
 
606
      INTEGER    NSO, NSQSO
 
607
      PARAMETER (NSO=2, NSQSO=1)
 
608
C     
 
609
C     ARGUMENTS
 
610
C     
 
611
      INTEGER SOINDEX, ORDERS(NSO)
 
612
C     
 
613
C     LOCAL VARIABLES
 
614
C     
 
615
      INTEGER I
 
616
      INTEGER SQPLITORDERS(NSQSO,NSO)
 
617
      COMMON/SQPLITORDERS/SQPLITORDERS
 
618
C     
 
619
C     BEGIN CODE
 
620
C     
 
621
      IF (SOINDEX.GT.0.AND.SOINDEX.LE.NSQSO) THEN
 
622
        DO I=1,NSO
 
623
          ORDERS(I) =  SQPLITORDERS(SOINDEX,I)
 
624
        ENDDO
 
625
        RETURN
 
626
      ENDIF
 
627
 
 
628
      WRITE(*,*) 'ERROR:: Stopping function GET_SQUARED_ORDERS_FOR_SOI'
 
629
     $ //'NDEX'
 
630
      WRITE(*,*) 'Could not find squared orders index ',SOINDEX
 
631
      STOP
 
632
 
 
633
      END SUBROUTINE
 
634
 
 
635
C     This is the inverse subroutine of getting amplitude SO orders.
 
636
C      Not directly useful, but provided nonetheless.
 
637
      SUBROUTINE GET_ORDERS_FOR_AMPSOINDEX(SOINDEX,ORDERS)
 
638
C     
 
639
C     This functions returns the orders identified by the split order
 
640
C      index in argument. Order values correspond to following list of
 
641
C      couplings (and in this order):
 
642
C     ['QCD', 'QED']
 
643
C     
 
644
C     CONSTANTS
 
645
C     
 
646
      INTEGER    NSO, NAMPSO
 
647
      PARAMETER (NSO=2, NAMPSO=1)
 
648
C     
 
649
C     ARGUMENTS
 
650
C     
 
651
      INTEGER SOINDEX, ORDERS(NSO)
 
652
C     
 
653
C     LOCAL VARIABLES
 
654
C     
 
655
      INTEGER I
 
656
      INTEGER AMPSPLITORDERS(NAMPSO,NSO)
 
657
      COMMON/AMPSPLITORDERS/AMPSPLITORDERS
 
658
C     
 
659
C     BEGIN CODE
 
660
C     
 
661
      IF (SOINDEX.GT.0.AND.SOINDEX.LE.NAMPSO) THEN
 
662
        DO I=1,NSO
 
663
          ORDERS(I) =  AMPSPLITORDERS(SOINDEX,I)
 
664
        ENDDO
 
665
        RETURN
 
666
      ENDIF
 
667
 
 
668
      WRITE(*,*) 'ERROR:: Stopping function GET_ORDERS_FOR_AMPSOINDEX'
 
669
      WRITE(*,*) 'Could not find amplitude split orders index ',SOINDEX
 
670
      STOP
 
671
 
 
672
      END SUBROUTINE
 
673
 
 
674
C     This function is not directly useful, but included for completene
 
675
C     ss
 
676
      INTEGER FUNCTION SOINDEX_FOR_AMPORDERS(ORDERS)
 
677
C     
 
678
C     This functions returns the integer index identifying the
 
679
C      amplitude split orders passed in argument which correspond to
 
680
C      the values of the following list of couplings (and in this
 
681
C      order):
 
682
C     ['QCD', 'QED']
 
683
C     
 
684
C     CONSTANTS
 
685
C     
 
686
      INTEGER    NSO, NAMPSO
 
687
      PARAMETER (NSO=2, NAMPSO=1)
 
688
C     
 
689
C     ARGUMENTS
 
690
C     
 
691
      INTEGER ORDERS(NSO)
 
692
C     
 
693
C     LOCAL VARIABLES
 
694
C     
 
695
      INTEGER I,J
 
696
      INTEGER AMPSPLITORDERS(NAMPSO,NSO)
 
697
      COMMON/AMPSPLITORDERS/AMPSPLITORDERS
 
698
C     
 
699
C     BEGIN CODE
 
700
C     
 
701
      DO I=1,NAMPSO
 
702
        DO J=1,NSO
 
703
          IF (ORDERS(J).NE.AMPSPLITORDERS(I,J)) GOTO 1009
 
704
        ENDDO
 
705
        SOINDEX_FOR_AMPORDERS = I
 
706
        RETURN
 
707
 1009   CONTINUE
 
708
      ENDDO
 
709
 
 
710
      WRITE(*,*) 'ERROR:: Stopping function SOINDEX_FOR_AMPORDERS'
 
711
      WRITE(*,*) 'Could not find squared orders ',(ORDERS(I),I=1,NSO)
 
712
      STOP
 
713
 
 
714
      END
 
715