← Back to branch summary

~maddevelopers/mg5amcnlo/2.5.4_run_py8_at_evtgen

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/short_ML_SMQCD_optimized/ddx_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/ddx_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=1)
 
15
      INTEGER    NLOOPS, NLOOPGROUPS, NCTAMPS
 
16
      PARAMETER (NLOOPS=11, NLOOPGROUPS=8, NCTAMPS=29)
 
17
      INTEGER    NLOOPAMPS
 
18
      PARAMETER (NLOOPAMPS=40)
 
19
      INTEGER    NWAVEFUNCS,NLOOPWAVEFUNCS
 
20
      PARAMETER (NWAVEFUNCS=6,NLOOPWAVEFUNCS=26)
 
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 13
 
83
      CALL MP_FFV1_0(W(1,4),W(1,3),W(1,5),GC_5,AMPL(1,28))
 
84
      AMPL(1,28)=AMPL(1,28)*(2.0D0*UVWFCT_T_0)
 
85
C     Amplitude(s) for UVCT diagram with ID 14
 
86
      CALL MP_FFV1_0(W(1,4),W(1,3),W(1,5),GC_5,AMPL(2,29))
 
87
      AMPL(2,29)=AMPL(2,29)*(2.0D0*UVWFCT_B_0_1EPS)
 
88
 
 
89
      GOTO 1001
 
90
 3000 CONTINUE
 
91
      MP_UVCT_REQ_SO_DONE=.TRUE.
 
92
 1001 CONTINUE
 
93
      END
 
94