← 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_optimized/dux_mumvmxg/mp_helas_calls_ampb_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/long_ML_SMQCD_optimized/dux_mumvmxg/mp_helas_calls_ampb_1.f
 
1
      SUBROUTINE ML5_0_MP_HELAS_CALLS_AMPB_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=5)
 
10
      INTEGER    NCOMB
 
11
      PARAMETER (NCOMB=32)
 
12
 
 
13
      INTEGER NBORNAMPS
 
14
      PARAMETER (NBORNAMPS=2)
 
15
      INTEGER    NLOOPS, NLOOPGROUPS, NCTAMPS
 
16
      PARAMETER (NLOOPS=11, NLOOPGROUPS=9, NCTAMPS=28)
 
17
      INTEGER    NLOOPAMPS
 
18
      PARAMETER (NLOOPAMPS=39)
 
19
      INTEGER    NWAVEFUNCS,NLOOPWAVEFUNCS
 
20
      PARAMETER (NWAVEFUNCS=10,NLOOPWAVEFUNCS=27)
 
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_CT_REQ_SO_DONE) THEN
 
79
        GOTO 1001
 
80
      ENDIF
 
81
 
 
82
      CALL MP_IXXXXX(P(0,1),ZERO,NHEL(1),+1*IC(1),W(1,1))
 
83
      CALL MP_OXXXXX(P(0,2),ZERO,NHEL(2),-1*IC(2),W(1,2))
 
84
      CALL MP_OXXXXX(P(0,3),ZERO,NHEL(3),+1*IC(3),W(1,3))
 
85
      CALL MP_IXXXXX(P(0,4),ZERO,NHEL(4),-1*IC(4),W(1,4))
 
86
      CALL MP_VXXXXX(P(0,5),ZERO,NHEL(5),+1*IC(5),W(1,5))
 
87
      CALL MP_FFV1_2(W(1,1),W(1,5),GC_5,ZERO,ZERO,W(1,6))
 
88
      CALL MP_FFV2_3(W(1,4),W(1,3),GC_47,MDL_MW,MDL_WW,W(1,7))
 
89
C     Amplitude(s) for born diagram with ID 1
 
90
      CALL MP_FFV2_0(W(1,6),W(1,2),W(1,7),GC_47,AMP(1))
 
91
      CALL MP_FFV1_1(W(1,2),W(1,5),GC_5,ZERO,ZERO,W(1,8))
 
92
C     Amplitude(s) for born diagram with ID 2
 
93
      CALL MP_FFV2_0(W(1,1),W(1,8),W(1,7),GC_47,AMP(2))
 
94
      CALL MP_FFV2_1(W(1,2),W(1,7),GC_47,ZERO,ZERO,W(1,9))
 
95
C     Counter-term amplitude(s) for loop diagram number 3
 
96
      CALL MP_R2_QQ_1_0(W(1,6),W(1,9),R2_QQQ,AMPL(1,1))
 
97
C     Counter-term amplitude(s) for loop diagram number 4
 
98
      CALL MP_FFV2_0(W(1,6),W(1,2),W(1,7),R2_SXCW,AMPL(1,2))
 
99
C     Counter-term amplitude(s) for loop diagram number 5
 
100
      CALL MP_FFV2_0(W(1,1),W(1,8),W(1,7),R2_SXCW,AMPL(1,3))
 
101
C     Counter-term amplitude(s) for loop diagram number 7
 
102
      CALL MP_FFV1_0(W(1,1),W(1,9),W(1,5),UV_GQQQ_1EPS,AMPL(2,4))
 
103
      CALL MP_FFV1_0(W(1,1),W(1,9),W(1,5),UV_GQQQ_1EPS,AMPL(2,5))
 
104
      CALL MP_FFV1_0(W(1,1),W(1,9),W(1,5),UV_GQQQ_1EPS,AMPL(2,6))
 
105
      CALL MP_FFV1_0(W(1,1),W(1,9),W(1,5),UV_GQQQ_1EPS,AMPL(2,7))
 
106
      CALL MP_FFV1_0(W(1,1),W(1,9),W(1,5),UV_GQQB,AMPL(1,8))
 
107
      CALL MP_FFV1_0(W(1,1),W(1,9),W(1,5),UV_GQQQ_1EPS,AMPL(2,9))
 
108
      CALL MP_FFV1_0(W(1,1),W(1,9),W(1,5),UV_GQQT,AMPL(1,10))
 
109
      CALL MP_FFV1_0(W(1,1),W(1,9),W(1,5),UV_GQQQ_1EPS,AMPL(2,11))
 
110
      CALL MP_FFV1_0(W(1,1),W(1,9),W(1,5),UV_GQQG_1EPS,AMPL(2,12))
 
111
      CALL MP_FFV1_0(W(1,1),W(1,9),W(1,5),R2_GQQ,AMPL(1,13))
 
112
      CALL MP_FFV2_2(W(1,1),W(1,7),GC_47,ZERO,ZERO,W(1,10))
 
113
C     Counter-term amplitude(s) for loop diagram number 11
 
114
      CALL MP_R2_QQ_1_0(W(1,10),W(1,8),R2_QQQ,AMPL(1,14))
 
115
C     Counter-term amplitude(s) for loop diagram number 12
 
116
      CALL MP_FFV1_0(W(1,10),W(1,2),W(1,5),UV_GQQQ_1EPS,AMPL(2,15))
 
117
      CALL MP_FFV1_0(W(1,10),W(1,2),W(1,5),UV_GQQQ_1EPS,AMPL(2,16))
 
118
      CALL MP_FFV1_0(W(1,10),W(1,2),W(1,5),UV_GQQQ_1EPS,AMPL(2,17))
 
119
      CALL MP_FFV1_0(W(1,10),W(1,2),W(1,5),UV_GQQQ_1EPS,AMPL(2,18))
 
120
      CALL MP_FFV1_0(W(1,10),W(1,2),W(1,5),UV_GQQB,AMPL(1,19))
 
121
      CALL MP_FFV1_0(W(1,10),W(1,2),W(1,5),UV_GQQQ_1EPS,AMPL(2,20))
 
122
      CALL MP_FFV1_0(W(1,10),W(1,2),W(1,5),UV_GQQT,AMPL(1,21))
 
123
      CALL MP_FFV1_0(W(1,10),W(1,2),W(1,5),UV_GQQQ_1EPS,AMPL(2,22))
 
124
      CALL MP_FFV1_0(W(1,10),W(1,2),W(1,5),UV_GQQG_1EPS,AMPL(2,23))
 
125
      CALL MP_FFV1_0(W(1,10),W(1,2),W(1,5),R2_GQQ,AMPL(1,24))
 
126
 
 
127
      GOTO 1001
 
128
 2000 CONTINUE
 
129
      MP_CT_REQ_SO_DONE=.TRUE.
 
130
 1001 CONTINUE
 
131
      END
 
132