← Back to branch summary

~maddevelopers/mg5amcnlo/3.0.1_4FS_merged3.0.2_latest

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOExportFKSTest/test_wprod_fksew/%SubProcesses%P0_udx_veep%V0_udx_veep%polynomial.f

  • Committer: marco zaro
  • Date: 2020-06-26 09:54:58 UTC
  • mfrom: (940.3.7 3.0.3)
  • Revision ID: marco.zaro@gmail.com-20200626095458-spujrl0ifwno0jgg
merged with 3.0.3

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/IOExportFKSTest/test_wprod_fksew/%SubProcesses%P0_udx_veep%V0_udx_veep%polynomial.f
113
113
        CFTOT=CMPLX(CF_N(COLOR_ID,I)/(ONE*ABS(CF_D(COLOR_ID,I))),ZERO
114
114
     $   ,KIND=8)
115
115
        IF(CF_D(COLOR_ID,I).LT.0) CFTOT=CFTOT*IMAG1
116
 
        CONST(ML5SOINDEX_FOR_BORN_AMP(I))=CONST(ML5SOINDEX_FOR_BORN_AMP
117
 
     $(I))+CFTOT*CONJG(AMP(I))
 
116
        CONST(ML5SOINDEX_FOR_BORN_AMP(I))
 
117
     $   =CONST(ML5SOINDEX_FOR_BORN_AMP(I))+CFTOT*CONJG(AMP(I))
118
118
      ENDDO
119
119
 
120
120
      DO I=1,NAMPSO
132
132
 
133
133
      END
134
134
 
 
135
      SUBROUTINE INVERT_MOMENTA_IN_POLYNOMIAL(NCOEFS,POLYNOMIAL)
 
136
C     Just a handy subroutine to modify the coefficients for the
 
137
C     tranformation q_loop -> -q_loop
 
138
C     It is only used for the NINJA interface
 
139
      USE POLYNOMIAL_CONSTANTS
 
140
      IMPLICIT NONE
 
141
 
 
142
      INTEGER I, NCOEFS
 
143
 
 
144
      COMPLEX*16 POLYNOMIAL(0:NCOEFS-1)
 
145
 
 
146
      DO I=0,NCOEFS-1
 
147
        IF (MOD(COEFTORANK_MAP(I),2).EQ.1) THEN
 
148
          POLYNOMIAL(I)=-POLYNOMIAL(I)
 
149
        ENDIF
 
150
      ENDDO
 
151
 
 
152
      END
135
153
 
136
154
C     Now the routines to update the wavefunctions
137
155
 
214
232
        CFTOT=CMPLX(CF_N(COLOR_ID,I)/(ONE*ABS(CF_D(COLOR_ID,I))),ZERO
215
233
     $   ,KIND=16)
216
234
        IF(CF_D(COLOR_ID,I).LT.0) CFTOT=CFTOT*IMAG1
217
 
        CONST(ML5SOINDEX_FOR_BORN_AMP(I))=CONST(ML5SOINDEX_FOR_BORN_AMP
218
 
     $(I))+CFTOT*CONJG(AMP(I))
 
235
        CONST(ML5SOINDEX_FOR_BORN_AMP(I))
 
236
     $   =CONST(ML5SOINDEX_FOR_BORN_AMP(I))+CFTOT*CONJG(AMP(I))
219
237
      ENDDO
220
238
 
221
239
      DO I=1,NAMPSO
233
251
 
234
252
      END
235
253
 
 
254
      SUBROUTINE MP_INVERT_MOMENTA_IN_POLYNOMIAL(NCOEFS,POLYNOMIAL)
 
255
C     Just a handy subroutine to modify the coefficients for the
 
256
C     tranformation q_loop -> -q_loop
 
257
C     It is only used for the NINJA interface
 
258
      USE POLYNOMIAL_CONSTANTS
 
259
      IMPLICIT NONE
 
260
 
 
261
      INTEGER I, NCOEFS
 
262
 
 
263
      COMPLEX*32 POLYNOMIAL(0:NCOEFS-1)
 
264
 
 
265
      DO I=0,NCOEFS-1
 
266
        IF (MOD(COEFTORANK_MAP(I),2).EQ.1) THEN
 
267
          POLYNOMIAL(I)=-POLYNOMIAL(I)
 
268
        ENDIF
 
269
      ENDDO
 
270
 
 
271
      END
236
272
 
237
273
C     Now the routines to update the wavefunctions
238
274