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

  • Committer: olivier Mattelaer
  • Date: 2016-05-12 11:00:18 UTC
  • mfrom: (262.1.150 2.3.4)
  • Revision ID: olivier.mattelaer@uclouvain.be-20160512110018-sevb79f0wm4g8mpp
pass to 2.4.0

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/short_ML_SMQCD_optimized/gg_ttx/mp_helas_calls_uvct_1.f
 
1
      SUBROUTINE ML5_0_MP_HELAS_CALLS_UVCT_1(P,NHEL,H,IC)
 
2
C     
 
3
      USE ML5_0_POLYNOMIAL_CONSTANTS
 
4
      IMPLICIT NONE
 
5
C     
 
6
C     CONSTANTS
 
7
C     
 
8
      INTEGER    NEXTERNAL
 
9
      PARAMETER (NEXTERNAL=4)
 
10
      INTEGER    NCOMB
 
11
      PARAMETER (NCOMB=16)
 
12
 
 
13
      INTEGER NBORNAMPS
 
14
      PARAMETER (NBORNAMPS=3)
 
15
      INTEGER    NLOOPS, NLOOPGROUPS, NCTAMPS
 
16
      PARAMETER (NLOOPS=44, NLOOPGROUPS=26, NCTAMPS=85)
 
17
      INTEGER    NLOOPAMPS
 
18
      PARAMETER (NLOOPAMPS=129)
 
19
      INTEGER    NWAVEFUNCS,NLOOPWAVEFUNCS
 
20
      PARAMETER (NWAVEFUNCS=10,NLOOPWAVEFUNCS=93)
 
21
      REAL*16     ZERO
 
22
      PARAMETER (ZERO=0.0E0_16)
 
23
      COMPLEX*32     IZERO
 
24
      PARAMETER (IZERO=CMPLX(0.0E0_16,0.0E0_16,KIND=16))
 
25
C     These are constants related to the split orders
 
26
      INTEGER    NSO, NSQUAREDSO, NAMPSO
 
27
      PARAMETER (NSO=0, NSQUAREDSO=0, NAMPSO=0)
 
28
C     
 
29
C     ARGUMENTS
 
30
C     
 
31
      REAL*16 P(0:3,NEXTERNAL)
 
32
      INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
 
33
      INTEGER H
 
34
C     
 
35
C     LOCAL VARIABLES
 
36
C     
 
37
      INTEGER I,J,K
 
38
      COMPLEX*32 COEFS(MAXLWFSIZE,0:VERTEXMAXCOEFS-1,MAXLWFSIZE)
 
39
C     
 
40
C     GLOBAL VARIABLES
 
41
C     
 
42
      INCLUDE 'mp_coupl_same_name.inc'
 
43
 
 
44
      INTEGER GOODHEL(NCOMB)
 
45
      LOGICAL GOODAMP(NSQUAREDSO,NLOOPGROUPS)
 
46
      COMMON/ML5_0_FILTERS/GOODAMP,GOODHEL
 
47
 
 
48
      INTEGER SQSO_TARGET
 
49
      COMMON/ML5_0_SOCHOICE/SQSO_TARGET
 
50
 
 
51
      LOGICAL UVCT_REQ_SO_DONE,MP_UVCT_REQ_SO_DONE,CT_REQ_SO_DONE
 
52
     $ ,MP_CT_REQ_SO_DONE,LOOP_REQ_SO_DONE,MP_LOOP_REQ_SO_DONE
 
53
     $ ,CTCALL_REQ_SO_DONE,FILTER_SO
 
54
      COMMON/ML5_0_SO_REQS/UVCT_REQ_SO_DONE,MP_UVCT_REQ_SO_DONE
 
55
     $ ,CT_REQ_SO_DONE,MP_CT_REQ_SO_DONE,LOOP_REQ_SO_DONE
 
56
     $ ,MP_LOOP_REQ_SO_DONE,CTCALL_REQ_SO_DONE,FILTER_SO
 
57
 
 
58
      COMPLEX*32 AMP(NBORNAMPS)
 
59
      COMMON/ML5_0_MP_AMPS/AMP
 
60
      COMPLEX*32 W(20,NWAVEFUNCS)
 
61
      COMMON/ML5_0_MP_W/W
 
62
 
 
63
      COMPLEX*32 WL(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE
 
64
     $ ,0:NLOOPWAVEFUNCS)
 
65
      COMPLEX*32 PL(0:3,0:NLOOPWAVEFUNCS)
 
66
      COMMON/ML5_0_MP_WL/WL,PL
 
67
 
 
68
      COMPLEX*32 AMPL(3,NCTAMPS)
 
69
      COMMON/ML5_0_MP_AMPL/AMPL
 
70
 
 
71
C     
 
72
C     ----------
 
73
C     BEGIN CODE
 
74
C     ----------
 
75
 
 
76
C     The target squared split order contribution is already reached
 
77
C      if true.
 
78
      IF (FILTER_SO.AND.MP_UVCT_REQ_SO_DONE) THEN
 
79
        GOTO 1001
 
80
      ENDIF
 
81
 
 
82
C     Amplitude(s) for UVCT diagram with ID 40
 
83
      CALL MP_FFV1_0(W(1,4),W(1,3),W(1,5),GC_5,AMPL(1,80))
 
84
      AMPL(1,80)=AMPL(1,80)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1+2.0D0
 
85
     $ *UVWFCT_T_0)
 
86
C     Amplitude(s) for UVCT diagram with ID 41
 
87
      CALL MP_FFV1_0(W(1,4),W(1,3),W(1,5),GC_5,AMPL(2,81))
 
88
      AMPL(2,81)=AMPL(2,81)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0
 
89
     $ *UVWFCT_G_2_1EPS)
 
90
C     Amplitude(s) for UVCT diagram with ID 42
 
91
      CALL MP_FFV1_0(W(1,4),W(1,6),W(1,2),GC_5,AMPL(1,82))
 
92
      AMPL(1,82)=AMPL(1,82)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1+2.0D0
 
93
     $ *UVWFCT_T_0)
 
94
C     Amplitude(s) for UVCT diagram with ID 43
 
95
      CALL MP_FFV1_0(W(1,4),W(1,6),W(1,2),GC_5,AMPL(2,83))
 
96
      AMPL(2,83)=AMPL(2,83)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0
 
97
     $ *UVWFCT_G_2_1EPS)
 
98
C     Amplitude(s) for UVCT diagram with ID 44
 
99
      CALL MP_FFV1_0(W(1,7),W(1,3),W(1,2),GC_5,AMPL(1,84))
 
100
      AMPL(1,84)=AMPL(1,84)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1+2.0D0
 
101
     $ *UVWFCT_T_0)
 
102
C     Amplitude(s) for UVCT diagram with ID 45
 
103
      CALL MP_FFV1_0(W(1,7),W(1,3),W(1,2),GC_5,AMPL(2,85))
 
104
      AMPL(2,85)=AMPL(2,85)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0
 
105
     $ *UVWFCT_G_2_1EPS)
 
106
 
 
107
      GOTO 1001
 
108
 3000 CONTINUE
 
109
      MP_UVCT_REQ_SO_DONE=.TRUE.
 
110
 1001 CONTINUE
 
111
      END
 
112