← 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_optimized/gg_ttx/polynomial.f

mwrge with 2.7.3

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/short_ML_SMQCD_optimized/gg_ttx/polynomial.f
99
99
C     
100
100
C     FUNCTIONS
101
101
C     
102
 
      INTEGER ML5_0_ML5SOINDEX_FOR_BORN_AMP, ML5_0_ML5SOINDEX_FOR_LOOP_
103
 
     $AMP, ML5_0_ML5SQSOINDEX
 
102
      INTEGER ML5_0_ML5SOINDEX_FOR_BORN_AMP,
 
103
     $  ML5_0_ML5SOINDEX_FOR_LOOP_AMP, ML5_0_ML5SQSOINDEX
104
104
C     
105
105
C     GLOBAL VARIABLES
106
106
C     
134
134
        CFTOT=CMPLX(CF_N(COLOR_ID,I)/(ONE*ABS(CF_D(COLOR_ID,I))),ZERO
135
135
     $   ,KIND=8)
136
136
        IF(CF_D(COLOR_ID,I).LT.0) CFTOT=CFTOT*IMAG1
137
 
        CONST(ML5_0_ML5SOINDEX_FOR_BORN_AMP(I))=CONST(ML5_0_ML5SOINDEX_
138
 
     $FOR_BORN_AMP(I))+CFTOT*CONJG(AMP(I))
 
137
        CONST(ML5_0_ML5SOINDEX_FOR_BORN_AMP(I))
 
138
     $   =CONST(ML5_0_ML5SOINDEX_FOR_BORN_AMP(I))+CFTOT*CONJG(AMP(I))
139
139
      ENDDO
140
140
 
141
141
      DO I=1,NAMPSO
146
146
            CONST(I)=CONST(I)*GOODHEL(HELCONFIG)
147
147
          ENDIF
148
148
          CALL ML5_0_MERGE_WL(LOOP_WF,RANK,LCUT_SIZE,CONST(I)
149
 
     $     ,LOOPCOEFS(0,ML5_0_ML5SQSOINDEX(I,ML5_0_ML5SOINDEX_FOR_LOOP_
150
 
     $AMP(COLOR_ID)),LOOP_GROUP_NUMBER))
 
149
     $     ,LOOPCOEFS(0,ML5_0_ML5SQSOINDEX(I
 
150
     $     ,ML5_0_ML5SOINDEX_FOR_LOOP_AMP(COLOR_ID)),LOOP_GROUP_NUMBER)
 
151
     $     )
151
152
        ENDIF
152
153
      ENDDO
153
154
 
200
201
C     
201
202
C     FUNCTIONS
202
203
C     
203
 
      INTEGER ML5_0_ML5SOINDEX_FOR_BORN_AMP, ML5_0_ML5SOINDEX_FOR_LOOP_
204
 
     $AMP, ML5_0_ML5SQSOINDEX
 
204
      INTEGER ML5_0_ML5SOINDEX_FOR_BORN_AMP,
 
205
     $  ML5_0_ML5SOINDEX_FOR_LOOP_AMP, ML5_0_ML5SQSOINDEX
205
206
C     
206
207
C     GLOBAL VARIABLES
207
208
C     
235
236
        CFTOT=CMPLX(CF_N(COLOR_ID,I)/(ONE*ABS(CF_D(COLOR_ID,I))),ZERO
236
237
     $   ,KIND=16)
237
238
        IF(CF_D(COLOR_ID,I).LT.0) CFTOT=CFTOT*IMAG1
238
 
        CONST(ML5_0_ML5SOINDEX_FOR_BORN_AMP(I))=CONST(ML5_0_ML5SOINDEX_
239
 
     $FOR_BORN_AMP(I))+CFTOT*CONJG(AMP(I))
 
239
        CONST(ML5_0_ML5SOINDEX_FOR_BORN_AMP(I))
 
240
     $   =CONST(ML5_0_ML5SOINDEX_FOR_BORN_AMP(I))+CFTOT*CONJG(AMP(I))
240
241
      ENDDO
241
242
 
242
243
      DO I=1,NAMPSO
247
248
            CONST(I)=CONST(I)*GOODHEL(HELCONFIG)
248
249
          ENDIF
249
250
          CALL MP_ML5_0_MERGE_WL(LOOP_WF,RANK,LCUT_SIZE,CONST(I)
250
 
     $     ,LOOPCOEFS(0,ML5_0_ML5SQSOINDEX(I,ML5_0_ML5SOINDEX_FOR_LOOP_
251
 
     $AMP(COLOR_ID)),LOOP_GROUP_NUMBER))
 
251
     $     ,LOOPCOEFS(0,ML5_0_ML5SQSOINDEX(I
 
252
     $     ,ML5_0_ML5SOINDEX_FOR_LOOP_AMP(COLOR_ID)),LOOP_GROUP_NUMBER)
 
253
     $     )
252
254
        ENDIF
253
255
      ENDDO
254
256
 
775
777
            OUT(J,4,I)=OUT(J,4,I)+A(K,0,I)*B(J,4,K)+A(K,4,I)*B(J,0,K)
776
778
            OUT(J,5,I)=OUT(J,5,I)+A(K,1,I)*B(J,1,K)
777
779
            OUT(J,6,I)=OUT(J,6,I)+A(K,1,I)*B(J,2,K)+A(K,2,I)*B(J,1,K)
 
780
            OUT(J,8,I)=OUT(J,8,I)+A(K,1,I)*B(J,3,K)+A(K,3,I)*B(J,1,K)
 
781
            OUT(J,11,I)=OUT(J,11,I)+A(K,1,I)*B(J,4,K)+A(K,4,I)*B(J,1,K)
778
782
            OUT(J,7,I)=OUT(J,7,I)+A(K,2,I)*B(J,2,K)
779
 
            OUT(J,8,I)=OUT(J,8,I)+A(K,1,I)*B(J,3,K)+A(K,3,I)*B(J,1,K)
780
783
            OUT(J,9,I)=OUT(J,9,I)+A(K,2,I)*B(J,3,K)+A(K,3,I)*B(J,2,K)
 
784
            OUT(J,12,I)=OUT(J,12,I)+A(K,2,I)*B(J,4,K)+A(K,4,I)*B(J,2,K)
781
785
            OUT(J,10,I)=OUT(J,10,I)+A(K,3,I)*B(J,3,K)
782
 
            OUT(J,11,I)=OUT(J,11,I)+A(K,1,I)*B(J,4,K)+A(K,4,I)*B(J,1,K)
783
 
            OUT(J,12,I)=OUT(J,12,I)+A(K,2,I)*B(J,4,K)+A(K,4,I)*B(J,2,K)
784
786
            OUT(J,13,I)=OUT(J,13,I)+A(K,3,I)*B(J,4,K)+A(K,4,I)*B(J,3,K)
785
787
            OUT(J,14,I)=OUT(J,14,I)+A(K,4,I)*B(J,4,K)
786
788
          ENDDO
810
812
            OUT(J,4,I)=OUT(J,4,I)+A(K,0,I)*B(J,4,K)+A(K,4,I)*B(J,0,K)
811
813
            OUT(J,5,I)=OUT(J,5,I)+A(K,1,I)*B(J,1,K)
812
814
            OUT(J,6,I)=OUT(J,6,I)+A(K,1,I)*B(J,2,K)+A(K,2,I)*B(J,1,K)
 
815
            OUT(J,8,I)=OUT(J,8,I)+A(K,1,I)*B(J,3,K)+A(K,3,I)*B(J,1,K)
 
816
            OUT(J,11,I)=OUT(J,11,I)+A(K,1,I)*B(J,4,K)+A(K,4,I)*B(J,1,K)
813
817
            OUT(J,7,I)=OUT(J,7,I)+A(K,2,I)*B(J,2,K)
814
 
            OUT(J,8,I)=OUT(J,8,I)+A(K,1,I)*B(J,3,K)+A(K,3,I)*B(J,1,K)
815
818
            OUT(J,9,I)=OUT(J,9,I)+A(K,2,I)*B(J,3,K)+A(K,3,I)*B(J,2,K)
 
819
            OUT(J,12,I)=OUT(J,12,I)+A(K,2,I)*B(J,4,K)+A(K,4,I)*B(J,2,K)
816
820
            OUT(J,10,I)=OUT(J,10,I)+A(K,3,I)*B(J,3,K)
817
 
            OUT(J,11,I)=OUT(J,11,I)+A(K,1,I)*B(J,4,K)+A(K,4,I)*B(J,1,K)
818
 
            OUT(J,12,I)=OUT(J,12,I)+A(K,2,I)*B(J,4,K)+A(K,4,I)*B(J,2,K)
819
821
            OUT(J,13,I)=OUT(J,13,I)+A(K,3,I)*B(J,4,K)+A(K,4,I)*B(J,3,K)
820
822
            OUT(J,14,I)=OUT(J,14,I)+A(K,4,I)*B(J,4,K)
821
823
          ENDDO