← Back to branch summary

~maddevelopers/mg5amcnlo/2.7.1.3

« back to all changes in this revision

Viewing changes to madgraph/iolibs/template_files/loop_optimized/polynomial.inc

  • 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:
madgraph/iolibs/template_files/loop_optimized/polynomial.inc
3
3
C MULTIPLY BY THE BORN
4
4
 
5
5
      SUBROUTINE %(mp_prefix)s%(proc_prefix)sCREATE_LOOP_COEFS(LOOP_WF,RANK,LCUT_SIZE,LOOP_GROUP_NUMBER,SYMFACT,MULTIPLIER,COLOR_ID,HELCONFIG)
 
6
          USE %(proc_prefix)sPOLYNOMIAL_CONSTANTS         
6
7
          implicit none
7
8
C  
8
9
C CONSTANTS 
17
18
      PARAMETER (IMAG1=(ZERO,ONE))
18
19
          %(complex_format)s CMPLX_ZERO
19
20
          PARAMETER (CMPLX_ZERO=(ZERO,ZERO))
20
 
          INTEGER MAXLWFSIZE
21
 
          PARAMETER (MAXLWFSIZE=%(max_lwf_size)d)
22
 
          INTEGER LOOPMAXCOEFS
23
 
          PARAMETER (LOOPMAXCOEFS=%(loop_max_coefs)d)
24
21
      INTEGER    NCOLORROWS
25
22
          PARAMETER (NCOLORROWS=%(nloopamps)d)
26
23
          INTEGER    NLOOPGROUPS
99
96
C amplitude level so that no multiplication is performed.
100
97
 
101
98
      SUBROUTINE %(mp_prefix)s%(proc_prefix)sCREATE_LOOP_COEFS(LOOP_WF,RANK,LCUT_SIZE,LOOP_GROUP_NUMBER,SYMFACT,MULTIPLIER)
 
99
          USE %(proc_prefix)sPOLYNOMIAL_CONSTANTS         
102
100
          implicit none
103
101
C  
104
102
C CONSTANTS 
109
107
      PARAMETER (IMAG1=(ZERO,ONE))
110
108
          %(complex_format)s CMPLX_ZERO
111
109
          PARAMETER (CMPLX_ZERO=(ZERO,ZERO))
112
 
          INTEGER MAXLWFSIZE
113
 
          PARAMETER (MAXLWFSIZE=%(max_lwf_size)d)
114
 
          INTEGER LOOPMAXCOEFS
115
 
          PARAMETER (LOOPMAXCOEFS=%(loop_max_coefs)d)
116
110
          INTEGER    NLOOPGROUPS
117
111
      PARAMETER (NLOOPGROUPS=%(nloop_groups)d)
118
112
          INTEGER    NCOMB
141
135
 
142
136
          END
143
137
## }
 
138
 
 
139
## if(ninja_available){
 
140
      SUBROUTINE %(mp_prefix)s%(proc_prefix)sINVERT_MOMENTA_IN_POLYNOMIAL(NCOEFS,POLYNOMIAL)
 
141
C       Just a handy subroutine to modify the coefficients for the
 
142
C       tranformation q_loop -> -q_loop
 
143
C       It is only used for the NINJA interface
 
144
        USE %(proc_prefix)sPOLYNOMIAL_CONSTANTS
 
145
        IMPLICIT NONE
 
146
 
 
147
        INTEGER I, NCOEFS
 
148
 
 
149
                %(complex_format)s POLYNOMIAL(0:NCOEFS-1)
 
150
 
 
151
        DO I=0,NCOEFS-1
 
152
          IF (MOD(COEFTORANK_MAP(I),2).eq.1) then
 
153
            POLYNOMIAL(I)=-POLYNOMIAL(I)
 
154
          ENDIF
 
155
        ENDDO
 
156
 
 
157
      END
 
158
## }
 
159
 
 
160
C Now the routines to update the wavefunctions
 
161