← Back to branch summary

~maddevelopers/mg5amcnlo/MS_with_decayevents

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/short_ML_SMQCD_optimized/gg_ttx/mp_compute_loop_coefs.f

  • Committer: olivier Mattelaer
  • Date: 2014-03-31 16:58:22 UTC
  • mfrom: (252.1.59 2.1.1)
  • Revision ID: olivier.mattelaer@uclouvain.be-20140331165822-4eb1040586ic1y1v
pass to 2.1.1

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/short_ML_SMQCD_optimized/gg_ttx/mp_compute_loop_coefs.f
 
1
      SUBROUTINE MP_COMPUTE_LOOP_COEFS(PS,ANSDP)
 
2
C     
 
3
C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
 
4
C     By the MadGraph5_aMC@NLO Development Team
 
5
C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
 
6
C     
 
7
C     Returns amplitude squared summed/avg over colors
 
8
C     and helicities for the point in phase space P(0:3,NEXTERNAL)
 
9
C     and external lines W(0:6,NEXTERNAL)
 
10
C     
 
11
C     Process: g g > t t~ QED=0 QCD=2 [ virt = QCD ] WEIGHTED=6
 
12
C     
 
13
      IMPLICIT NONE
 
14
C     
 
15
C     CONSTANTS
 
16
C     
 
17
      CHARACTER*64 PARAMFILENAME
 
18
      PARAMETER ( PARAMFILENAME='MadLoopParams.dat')
 
19
      INTEGER NBORNAMPS
 
20
      PARAMETER (NBORNAMPS=3)
 
21
      INTEGER    NLOOPS, NLOOPGROUPS, NCTAMPS
 
22
      PARAMETER (NLOOPS=53, NLOOPGROUPS=26, NCTAMPS=85)
 
23
      INTEGER    NCOLORROWS
 
24
      PARAMETER (NCOLORROWS=138)
 
25
      INTEGER    NEXTERNAL
 
26
      PARAMETER (NEXTERNAL=4)
 
27
      INTEGER    NWAVEFUNCS,NLOOPWAVEFUNCS
 
28
      PARAMETER (NWAVEFUNCS=10,NLOOPWAVEFUNCS=117)
 
29
      INTEGER MAXLWFSIZE
 
30
      PARAMETER (MAXLWFSIZE=4)
 
31
      INTEGER LOOPMAXCOEFS, VERTEXMAXCOEFS
 
32
      PARAMETER (LOOPMAXCOEFS=35, VERTEXMAXCOEFS=5)
 
33
      INTEGER    NCOMB
 
34
      PARAMETER (NCOMB=16)
 
35
      REAL*16    ZERO
 
36
      PARAMETER (ZERO=0E0_16)
 
37
      COMPLEX*32 IMAG1
 
38
      PARAMETER (IMAG1=(0E0_16,1E0_16))
 
39
C     
 
40
C     ARGUMENTS 
 
41
C     
 
42
      REAL*16 PS(0:3,NEXTERNAL)
 
43
      REAL*8 ANSDP(3)
 
44
C     
 
45
C     LOCAL VARIABLES 
 
46
C     
 
47
      INTEGER I,J,K,H,DUMMY
 
48
      INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
 
49
      REAL*16 P(0:3,NEXTERNAL)
 
50
      DATA IC/NEXTERNAL*1/
 
51
      REAL*16 ANS(3)
 
52
      COMPLEX*32 COEFS(MAXLWFSIZE,0:VERTEXMAXCOEFS-1,MAXLWFSIZE)
 
53
      COMPLEX*32 CFTOT
 
54
C     
 
55
C     GLOBAL VARIABLES
 
56
C     
 
57
      INCLUDE 'mp_coupl_same_name.inc'
 
58
 
 
59
      LOGICAL CHECKPHASE, HELDOUBLECHECKED
 
60
      COMMON/INIT/CHECKPHASE, HELDOUBLECHECKED
 
61
 
 
62
      INTEGER HELOFFSET
 
63
      INTEGER GOODHEL(NCOMB)
 
64
      LOGICAL GOODAMP(NLOOPGROUPS)
 
65
      COMMON/FILTERS/GOODAMP,GOODHEL,HELOFFSET
 
66
 
 
67
      INTEGER HELPICKED
 
68
      COMMON/HELCHOICE/HELPICKED
 
69
 
 
70
      COMPLEX*32 AMP(NBORNAMPS)
 
71
      COMMON/MP_AMPS/AMP
 
72
      COMPLEX*32 W(20,NWAVEFUNCS)
 
73
      COMMON/MP_W/W
 
74
 
 
75
      COMPLEX*16 DPW(20,NWAVEFUNCS)
 
76
      COMMON/W/DPW
 
77
 
 
78
      COMPLEX*32 WL(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE,0:NLOOPWAVEF
 
79
     $ UNCS)
 
80
      COMPLEX*32 PL(0:3,0:NLOOPWAVEFUNCS)
 
81
      COMMON/MP_WL/WL,PL
 
82
 
 
83
      COMPLEX*32 LOOPCOEFS(0:LOOPMAXCOEFS-1,NLOOPS)
 
84
      COMMON/MP_LCOEFS/LOOPCOEFS
 
85
 
 
86
      COMPLEX*32 AMPL(3,NCTAMPS)
 
87
      COMMON/MP_AMPL/AMPL
 
88
 
 
89
      INTEGER CF_D(NCOLORROWS,NBORNAMPS)
 
90
      INTEGER CF_N(NCOLORROWS,NBORNAMPS)
 
91
      COMMON/CF/CF_D,CF_N
 
92
 
 
93
      INTEGER HELC(NEXTERNAL,NCOMB)
 
94
      COMMON/HELCONFIGS/HELC
 
95
 
 
96
      LOGICAL MP_DONE_ONCE
 
97
      COMMON/MP_DONE_ONCE/MP_DONE_ONCE
 
98
 
 
99
C     ----------
 
100
C     BEGIN CODE
 
101
C     ----------
 
102
 
 
103
C     To be on the safe side, we always update the MP params here.
 
104
C     It can be redundant as this routine can be called a couple of
 
105
C      times for the same PS point during the stability checks.
 
106
C     But it is really not time consuming and I would rather be safe.
 
107
      CALL MP_UPDATE_AS_PARAM()
 
108
 
 
109
      MP_DONE_ONCE = .TRUE.
 
110
 
 
111
C     AS A SAFETY MEASURE WE FIRST COPY HERE THE PS POINT
 
112
      DO I=1,NEXTERNAL
 
113
        DO J=0,3
 
114
          P(J,I)=PS(J,I)
 
115
        ENDDO
 
116
      ENDDO
 
117
 
 
118
      DO I=0,3
 
119
        PL(I,0)=(ZERO,ZERO)
 
120
      ENDDO
 
121
      DO I=1,MAXLWFSIZE
 
122
        DO J=0,LOOPMAXCOEFS-1
 
123
          DO K=1,MAXLWFSIZE
 
124
            IF(I.EQ.K.AND.J.EQ.0) THEN
 
125
              WL(I,J,K,0)=(1.0E0_16,ZERO)
 
126
            ELSE
 
127
              WL(I,J,K,0)=(ZERO,ZERO)
 
128
            ENDIF
 
129
          ENDDO
 
130
        ENDDO
 
131
      ENDDO
 
132
 
 
133
      DO K=1, 3
 
134
        DO I=1,NCTAMPS
 
135
          AMPL(K,I)=(ZERO,ZERO)
 
136
        ENDDO
 
137
      ENDDO
 
138
 
 
139
      DO I=1,NLOOPS
 
140
        DO J=0,LOOPMAXCOEFS-1
 
141
          LOOPCOEFS(J,I)=(ZERO,ZERO)
 
142
        ENDDO
 
143
      ENDDO
 
144
 
 
145
      DO K=1,3
 
146
        ANSDP(K)=0.0D0
 
147
        ANS(K)=ZERO
 
148
      ENDDO
 
149
 
 
150
      DO H=1,NCOMB
 
151
        IF ((HELPICKED.EQ.H).OR.((HELPICKED.EQ.-1).AND.(CHECKPHASE.OR.(
 
152
     $   .NOT.HELDOUBLECHECKED).OR.(GOODHEL(H).GT.-HELOFFSET.AND.GOODHE
 
153
     $   L(H).NE.0)))) THEN
 
154
          DO I=1,NEXTERNAL
 
155
            NHEL(I)=HELC(I,H)
 
156
          ENDDO
 
157
          CALL MP_VXXXXX(P(0,1),ZERO,NHEL(1),-1*IC(1),W(1,1))
 
158
          CALL MP_VXXXXX(P(0,2),ZERO,NHEL(2),-1*IC(2),W(1,2))
 
159
          CALL MP_OXXXXX(P(0,3),MDL_MT,NHEL(3),+1*IC(3),W(1,3))
 
160
          CALL MP_IXXXXX(P(0,4),MDL_MT,NHEL(4),-1*IC(4),W(1,4))
 
161
          CALL MP_VVV1P0_1(W(1,1),W(1,2),GC_4,ZERO,ZERO,W(1,5))
 
162
C         Amplitude(s) for born diagram with ID 1
 
163
          CALL MP_FFV1_0(W(1,4),W(1,3),W(1,5),GC_5,AMP(1))
 
164
          CALL MP_FFV1_1(W(1,3),W(1,1),GC_5,MDL_MT,MDL_WT,W(1,6))
 
165
C         Amplitude(s) for born diagram with ID 2
 
166
          CALL MP_FFV1_0(W(1,4),W(1,6),W(1,2),GC_5,AMP(2))
 
167
          CALL MP_FFV1_2(W(1,4),W(1,1),GC_5,MDL_MT,MDL_WT,W(1,7))
 
168
C         Amplitude(s) for born diagram with ID 3
 
169
          CALL MP_FFV1_0(W(1,7),W(1,3),W(1,2),GC_5,AMP(3))
 
170
          CALL MP_FFV1P0_3(W(1,4),W(1,3),GC_5,ZERO,ZERO,W(1,8))
 
171
C         Counter-term amplitude(s) for loop diagram number 4
 
172
          CALL MP_R2_GG_1_0(W(1,5),W(1,8),R2_GGQ,AMPL(1,1))
 
173
C         Counter-term amplitude(s) for loop diagram number 5
 
174
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),UV_3GB_1EPS,AMPL(2,2))
 
175
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),UV_3GB_1EPS,AMPL(2,3))
 
176
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),UV_3GB_1EPS,AMPL(2,4))
 
177
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),UV_3GB_1EPS,AMPL(2,5))
 
178
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),UV_3GB,AMPL(1,6))
 
179
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),UV_3GB_1EPS,AMPL(2,7))
 
180
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),UV_3GT,AMPL(1,8))
 
181
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),UV_3GB_1EPS,AMPL(2,9))
 
182
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),UV_3GG_1EPS,AMPL(2,10))
 
183
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),R2_3GQ,AMPL(1,11))
 
184
C         Counter-term amplitude(s) for loop diagram number 7
 
185
          CALL MP_R2_GG_1_0(W(1,5),W(1,8),R2_GGQ,AMPL(1,12))
 
186
C         Counter-term amplitude(s) for loop diagram number 8
 
187
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),R2_3GQ,AMPL(1,13))
 
188
C         Counter-term amplitude(s) for loop diagram number 10
 
189
          CALL MP_R2_GG_1_0(W(1,5),W(1,8),R2_GGQ,AMPL(1,14))
 
190
C         Counter-term amplitude(s) for loop diagram number 11
 
191
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),R2_3GQ,AMPL(1,15))
 
192
C         Counter-term amplitude(s) for loop diagram number 13
 
193
          CALL MP_R2_GG_1_0(W(1,5),W(1,8),R2_GGQ,AMPL(1,16))
 
194
C         Counter-term amplitude(s) for loop diagram number 14
 
195
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),R2_3GQ,AMPL(1,17))
 
196
C         Counter-term amplitude(s) for loop diagram number 16
 
197
          CALL MP_R2_GG_1_R2_GG_3_0(W(1,5),W(1,8),R2_GGQ,R2_GGB,AMPL(1
 
198
     $     ,18))
 
199
C         Counter-term amplitude(s) for loop diagram number 17
 
200
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),R2_3GQ,AMPL(1,19))
 
201
C         Counter-term amplitude(s) for loop diagram number 19
 
202
          CALL MP_R2_GG_1_R2_GG_3_0(W(1,5),W(1,8),R2_GGQ,R2_GGT,AMPL(1
 
203
     $     ,20))
 
204
C         Counter-term amplitude(s) for loop diagram number 20
 
205
          CALL MP_FFV1_0(W(1,4),W(1,3),W(1,5),UV_GQQQ_1EPS,AMPL(2,21))
 
206
          CALL MP_FFV1_0(W(1,4),W(1,3),W(1,5),UV_GQQQ_1EPS,AMPL(2,22))
 
207
          CALL MP_FFV1_0(W(1,4),W(1,3),W(1,5),UV_GQQQ_1EPS,AMPL(2,23))
 
208
          CALL MP_FFV1_0(W(1,4),W(1,3),W(1,5),UV_GQQQ_1EPS,AMPL(2,24))
 
209
          CALL MP_FFV1_0(W(1,4),W(1,3),W(1,5),UV_GQQB,AMPL(1,25))
 
210
          CALL MP_FFV1_0(W(1,4),W(1,3),W(1,5),UV_GQQQ_1EPS,AMPL(2,26))
 
211
          CALL MP_FFV1_0(W(1,4),W(1,3),W(1,5),UV_GQQT,AMPL(1,27))
 
212
          CALL MP_FFV1_0(W(1,4),W(1,3),W(1,5),UV_GQQQ_1EPS,AMPL(2,28))
 
213
          CALL MP_FFV1_0(W(1,4),W(1,3),W(1,5),UV_GQQG_1EPS,AMPL(2,29))
 
214
          CALL MP_FFV1_0(W(1,4),W(1,3),W(1,5),R2_GQQ,AMPL(1,30))
 
215
          CALL MP_FFV1_2(W(1,4),W(1,2),GC_5,MDL_MT,MDL_WT,W(1,9))
 
216
C         Counter-term amplitude(s) for loop diagram number 22
 
217
          CALL MP_R2_QQ_1_R2_QQ_2_0(W(1,9),W(1,6),R2_QQQ,R2_QQT,AMPL(1
 
218
     $     ,31))
 
219
          CALL MP_R2_QQ_2_0(W(1,9),W(1,6),UV_TMASS,AMPL(1,32))
 
220
          CALL MP_R2_QQ_2_0(W(1,9),W(1,6),UV_TMASS_1EPS,AMPL(2,33))
 
221
C         Counter-term amplitude(s) for loop diagram number 23
 
222
          CALL MP_FFV1_0(W(1,4),W(1,6),W(1,2),UV_GQQQ_1EPS,AMPL(2,34))
 
223
          CALL MP_FFV1_0(W(1,4),W(1,6),W(1,2),UV_GQQQ_1EPS,AMPL(2,35))
 
224
          CALL MP_FFV1_0(W(1,4),W(1,6),W(1,2),UV_GQQQ_1EPS,AMPL(2,36))
 
225
          CALL MP_FFV1_0(W(1,4),W(1,6),W(1,2),UV_GQQQ_1EPS,AMPL(2,37))
 
226
          CALL MP_FFV1_0(W(1,4),W(1,6),W(1,2),UV_GQQB,AMPL(1,38))
 
227
          CALL MP_FFV1_0(W(1,4),W(1,6),W(1,2),UV_GQQQ_1EPS,AMPL(2,39))
 
228
          CALL MP_FFV1_0(W(1,4),W(1,6),W(1,2),UV_GQQT,AMPL(1,40))
 
229
          CALL MP_FFV1_0(W(1,4),W(1,6),W(1,2),UV_GQQQ_1EPS,AMPL(2,41))
 
230
          CALL MP_FFV1_0(W(1,4),W(1,6),W(1,2),UV_GQQG_1EPS,AMPL(2,42))
 
231
          CALL MP_FFV1_0(W(1,4),W(1,6),W(1,2),R2_GQQ,AMPL(1,43))
 
232
          CALL MP_FFV1_1(W(1,3),W(1,2),GC_5,MDL_MT,MDL_WT,W(1,10))
 
233
C         Counter-term amplitude(s) for loop diagram number 25
 
234
          CALL MP_R2_QQ_1_R2_QQ_2_0(W(1,7),W(1,10),R2_QQQ,R2_QQT
 
235
     $     ,AMPL(1,44))
 
236
          CALL MP_R2_QQ_2_0(W(1,7),W(1,10),UV_TMASS,AMPL(1,45))
 
237
          CALL MP_R2_QQ_2_0(W(1,7),W(1,10),UV_TMASS_1EPS,AMPL(2,46))
 
238
C         Counter-term amplitude(s) for loop diagram number 26
 
239
          CALL MP_FFV1_0(W(1,7),W(1,3),W(1,2),UV_GQQQ_1EPS,AMPL(2,47))
 
240
          CALL MP_FFV1_0(W(1,7),W(1,3),W(1,2),UV_GQQQ_1EPS,AMPL(2,48))
 
241
          CALL MP_FFV1_0(W(1,7),W(1,3),W(1,2),UV_GQQQ_1EPS,AMPL(2,49))
 
242
          CALL MP_FFV1_0(W(1,7),W(1,3),W(1,2),UV_GQQQ_1EPS,AMPL(2,50))
 
243
          CALL MP_FFV1_0(W(1,7),W(1,3),W(1,2),UV_GQQB,AMPL(1,51))
 
244
          CALL MP_FFV1_0(W(1,7),W(1,3),W(1,2),UV_GQQQ_1EPS,AMPL(2,52))
 
245
          CALL MP_FFV1_0(W(1,7),W(1,3),W(1,2),UV_GQQT,AMPL(1,53))
 
246
          CALL MP_FFV1_0(W(1,7),W(1,3),W(1,2),UV_GQQQ_1EPS,AMPL(2,54))
 
247
          CALL MP_FFV1_0(W(1,7),W(1,3),W(1,2),UV_GQQG_1EPS,AMPL(2,55))
 
248
          CALL MP_FFV1_0(W(1,7),W(1,3),W(1,2),R2_GQQ,AMPL(1,56))
 
249
C         Counter-term amplitude(s) for loop diagram number 28
 
250
          CALL MP_FFV1_0(W(1,4),W(1,10),W(1,1),UV_GQQQ_1EPS,AMPL(2,57))
 
251
          CALL MP_FFV1_0(W(1,4),W(1,10),W(1,1),UV_GQQQ_1EPS,AMPL(2,58))
 
252
          CALL MP_FFV1_0(W(1,4),W(1,10),W(1,1),UV_GQQQ_1EPS,AMPL(2,59))
 
253
          CALL MP_FFV1_0(W(1,4),W(1,10),W(1,1),UV_GQQQ_1EPS,AMPL(2,60))
 
254
          CALL MP_FFV1_0(W(1,4),W(1,10),W(1,1),UV_GQQB,AMPL(1,61))
 
255
          CALL MP_FFV1_0(W(1,4),W(1,10),W(1,1),UV_GQQQ_1EPS,AMPL(2,62))
 
256
          CALL MP_FFV1_0(W(1,4),W(1,10),W(1,1),UV_GQQT,AMPL(1,63))
 
257
          CALL MP_FFV1_0(W(1,4),W(1,10),W(1,1),UV_GQQQ_1EPS,AMPL(2,64))
 
258
          CALL MP_FFV1_0(W(1,4),W(1,10),W(1,1),UV_GQQG_1EPS,AMPL(2,65))
 
259
          CALL MP_FFV1_0(W(1,4),W(1,10),W(1,1),R2_GQQ,AMPL(1,66))
 
260
C         Counter-term amplitude(s) for loop diagram number 29
 
261
          CALL MP_FFV1_0(W(1,9),W(1,3),W(1,1),UV_GQQQ_1EPS,AMPL(2,67))
 
262
          CALL MP_FFV1_0(W(1,9),W(1,3),W(1,1),UV_GQQQ_1EPS,AMPL(2,68))
 
263
          CALL MP_FFV1_0(W(1,9),W(1,3),W(1,1),UV_GQQQ_1EPS,AMPL(2,69))
 
264
          CALL MP_FFV1_0(W(1,9),W(1,3),W(1,1),UV_GQQQ_1EPS,AMPL(2,70))
 
265
          CALL MP_FFV1_0(W(1,9),W(1,3),W(1,1),UV_GQQB,AMPL(1,71))
 
266
          CALL MP_FFV1_0(W(1,9),W(1,3),W(1,1),UV_GQQQ_1EPS,AMPL(2,72))
 
267
          CALL MP_FFV1_0(W(1,9),W(1,3),W(1,1),UV_GQQT,AMPL(1,73))
 
268
          CALL MP_FFV1_0(W(1,9),W(1,3),W(1,1),UV_GQQQ_1EPS,AMPL(2,74))
 
269
          CALL MP_FFV1_0(W(1,9),W(1,3),W(1,1),UV_GQQG_1EPS,AMPL(2,75))
 
270
          CALL MP_FFV1_0(W(1,9),W(1,3),W(1,1),R2_GQQ,AMPL(1,76))
 
271
C         Counter-term amplitude(s) for loop diagram number 31
 
272
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),R2_3GQ,AMPL(1,77))
 
273
C         Counter-term amplitude(s) for loop diagram number 41
 
274
          CALL MP_R2_GG_1_R2_GG_2_0(W(1,5),W(1,8),R2_GGG_1,R2_GGG_2
 
275
     $     ,AMPL(1,78))
 
276
C         Counter-term amplitude(s) for loop diagram number 42
 
277
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),R2_3GG,AMPL(1,79))
 
278
C         Amplitude(s) for UVCT diagram with ID 49
 
279
          CALL MP_FFV1_0(W(1,4),W(1,3),W(1,5),GC_5,AMPL(1,80))
 
280
          AMPL(1,80)=AMPL(1,80)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
 
281
     $     +2.0D0*UVWFCT_T_0)
 
282
C         Amplitude(s) for UVCT diagram with ID 50
 
283
          CALL MP_FFV1_0(W(1,4),W(1,3),W(1,5),GC_5,AMPL(2,81))
 
284
          AMPL(2,81)=AMPL(2,81)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0*UVWFCT_G_2
 
285
     $     _1EPS)
 
286
C         Amplitude(s) for UVCT diagram with ID 51
 
287
          CALL MP_FFV1_0(W(1,4),W(1,6),W(1,2),GC_5,AMPL(1,82))
 
288
          AMPL(1,82)=AMPL(1,82)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
 
289
     $     +2.0D0*UVWFCT_T_0)
 
290
C         Amplitude(s) for UVCT diagram with ID 52
 
291
          CALL MP_FFV1_0(W(1,4),W(1,6),W(1,2),GC_5,AMPL(2,83))
 
292
          AMPL(2,83)=AMPL(2,83)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0*UVWFCT_G_2
 
293
     $     _1EPS)
 
294
C         Amplitude(s) for UVCT diagram with ID 53
 
295
          CALL MP_FFV1_0(W(1,7),W(1,3),W(1,2),GC_5,AMPL(1,84))
 
296
          AMPL(1,84)=AMPL(1,84)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
 
297
     $     +2.0D0*UVWFCT_T_0)
 
298
C         Amplitude(s) for UVCT diagram with ID 54
 
299
          CALL MP_FFV1_0(W(1,7),W(1,3),W(1,2),GC_5,AMPL(2,85))
 
300
          AMPL(2,85)=AMPL(2,85)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0*UVWFCT_G_2
 
301
     $     _1EPS)
 
302
          IF (.NOT.CHECKPHASE.AND.HELDOUBLECHECKED.AND.HELPICKED.EQ.
 
303
     $     -1) THEN
 
304
            DUMMY=GOODHEL(H)
 
305
          ELSE
 
306
            DUMMY=1
 
307
          ENDIF
 
308
          DO I=1,NCTAMPS
 
309
            DO J=1,NBORNAMPS
 
310
              CFTOT=CMPLX(CF_N(I,J)/REAL(ABS(CF_D(I,J)),KIND=16)
 
311
     $         ,0.0E0_16,KIND=16)
 
312
              IF(CF_D(I,J).LT.0) CFTOT=CFTOT*IMAG1
 
313
              DO K=1,3
 
314
                ANS(K)=ANS(K)+DUMMY*2.0E0_16*REAL(CFTOT*AMPL(K,I)
 
315
     $           *CONJG(AMP(J)),KIND=16)
 
316
              ENDDO
 
317
            ENDDO
 
318
          ENDDO
 
319
C         Coefficient construction for loop diagram with ID 4
 
320
          CALL MP_FFV1L2_1(PL(0,0),W(1,5),GC_5,ZERO,ZERO,PL(0,1),COEFS)
 
321
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,1))
 
322
          CALL MP_FFV1L2_1(PL(0,1),W(1,8),GC_5,ZERO,ZERO,PL(0,2),COEFS)
 
323
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,1),4,COEFS,4,4,WL(1,0,1,2))
 
324
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,2),2,4,LOOPCOEFS(0,1),1
 
325
     $     ,86,H)
 
326
C         Coefficient construction for loop diagram with ID 5
 
327
          CALL MP_FFV1L1_2(PL(0,0),W(1,1),GC_5,ZERO,ZERO,PL(0,3),COEFS)
 
328
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,3))
 
329
          CALL MP_FFV1L1_2(PL(0,3),W(1,2),GC_5,ZERO,ZERO,PL(0,4),COEFS)
 
330
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,3),4,COEFS,4,4,WL(1,0,1,4))
 
331
          CALL MP_FFV1L1_2(PL(0,4),W(1,8),GC_5,ZERO,ZERO,PL(0,5),COEFS)
 
332
          CALL MP_UPDATE_WL_2_1(WL(1,0,1,4),4,COEFS,4,4,WL(1,0,1,5))
 
333
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,5),3,4,LOOPCOEFS(0,2),1
 
334
     $     ,87,H)
 
335
C         Coefficient construction for loop diagram with ID 6
 
336
          CALL MP_FFV1L2_1(PL(0,0),W(1,1),GC_5,ZERO,ZERO,PL(0,6),COEFS)
 
337
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,6))
 
338
          CALL MP_FFV1L2_1(PL(0,6),W(1,2),GC_5,ZERO,ZERO,PL(0,7),COEFS)
 
339
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,6),4,COEFS,4,4,WL(1,0,1,7))
 
340
          CALL MP_FFV1L2_1(PL(0,7),W(1,8),GC_5,ZERO,ZERO,PL(0,8),COEFS)
 
341
          CALL MP_UPDATE_WL_2_1(WL(1,0,1,7),4,COEFS,4,4,WL(1,0,1,8))
 
342
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,8),3,4,LOOPCOEFS(0,3),1
 
343
     $     ,88,H)
 
344
C         Coefficient construction for loop diagram with ID 7
 
345
          CALL MP_FFV1L2_1(PL(0,0),W(1,5),GC_5,ZERO,ZERO,PL(0,9),COEFS)
 
346
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,9))
 
347
          CALL MP_FFV1L2_1(PL(0,9),W(1,8),GC_5,ZERO,ZERO,PL(0,10)
 
348
     $     ,COEFS)
 
349
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,9),4,COEFS,4,4,WL(1,0,1,10))
 
350
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,10),2,4,LOOPCOEFS(0,4),1
 
351
     $     ,89,H)
 
352
C         Coefficient construction for loop diagram with ID 8
 
353
          CALL MP_FFV1L1_2(PL(0,0),W(1,1),GC_5,ZERO,ZERO,PL(0,11)
 
354
     $     ,COEFS)
 
355
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,11))
 
356
          CALL MP_FFV1L1_2(PL(0,11),W(1,2),GC_5,ZERO,ZERO,PL(0,12)
 
357
     $     ,COEFS)
 
358
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,11),4,COEFS,4,4,WL(1,0,1,12))
 
359
          CALL MP_FFV1L1_2(PL(0,12),W(1,8),GC_5,ZERO,ZERO,PL(0,13)
 
360
     $     ,COEFS)
 
361
          CALL MP_UPDATE_WL_2_1(WL(1,0,1,12),4,COEFS,4,4,WL(1,0,1,13))
 
362
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,13),3,4,LOOPCOEFS(0,5),1
 
363
     $     ,90,H)
 
364
C         Coefficient construction for loop diagram with ID 9
 
365
          CALL MP_FFV1L2_1(PL(0,0),W(1,1),GC_5,ZERO,ZERO,PL(0,14)
 
366
     $     ,COEFS)
 
367
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,14))
 
368
          CALL MP_FFV1L2_1(PL(0,14),W(1,2),GC_5,ZERO,ZERO,PL(0,15)
 
369
     $     ,COEFS)
 
370
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,14),4,COEFS,4,4,WL(1,0,1,15))
 
371
          CALL MP_FFV1L2_1(PL(0,15),W(1,8),GC_5,ZERO,ZERO,PL(0,16)
 
372
     $     ,COEFS)
 
373
          CALL MP_UPDATE_WL_2_1(WL(1,0,1,15),4,COEFS,4,4,WL(1,0,1,16))
 
374
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,16),3,4,LOOPCOEFS(0,6),1
 
375
     $     ,91,H)
 
376
C         Coefficient construction for loop diagram with ID 10
 
377
          CALL MP_FFV1L2_1(PL(0,0),W(1,5),GC_5,ZERO,ZERO,PL(0,17)
 
378
     $     ,COEFS)
 
379
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,17))
 
380
          CALL MP_FFV1L2_1(PL(0,17),W(1,8),GC_5,ZERO,ZERO,PL(0,18)
 
381
     $     ,COEFS)
 
382
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,17),4,COEFS,4,4,WL(1,0,1,18))
 
383
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,18),2,4,LOOPCOEFS(0,7),1
 
384
     $     ,92,H)
 
385
C         Coefficient construction for loop diagram with ID 11
 
386
          CALL MP_FFV1L1_2(PL(0,0),W(1,1),GC_5,ZERO,ZERO,PL(0,19)
 
387
     $     ,COEFS)
 
388
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,19))
 
389
          CALL MP_FFV1L1_2(PL(0,19),W(1,2),GC_5,ZERO,ZERO,PL(0,20)
 
390
     $     ,COEFS)
 
391
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,19),4,COEFS,4,4,WL(1,0,1,20))
 
392
          CALL MP_FFV1L1_2(PL(0,20),W(1,8),GC_5,ZERO,ZERO,PL(0,21)
 
393
     $     ,COEFS)
 
394
          CALL MP_UPDATE_WL_2_1(WL(1,0,1,20),4,COEFS,4,4,WL(1,0,1,21))
 
395
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,21),3,4,LOOPCOEFS(0,8),1
 
396
     $     ,93,H)
 
397
C         Coefficient construction for loop diagram with ID 12
 
398
          CALL MP_FFV1L2_1(PL(0,0),W(1,1),GC_5,ZERO,ZERO,PL(0,22)
 
399
     $     ,COEFS)
 
400
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,22))
 
401
          CALL MP_FFV1L2_1(PL(0,22),W(1,2),GC_5,ZERO,ZERO,PL(0,23)
 
402
     $     ,COEFS)
 
403
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,22),4,COEFS,4,4,WL(1,0,1,23))
 
404
          CALL MP_FFV1L2_1(PL(0,23),W(1,8),GC_5,ZERO,ZERO,PL(0,24)
 
405
     $     ,COEFS)
 
406
          CALL MP_UPDATE_WL_2_1(WL(1,0,1,23),4,COEFS,4,4,WL(1,0,1,24))
 
407
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,24),3,4,LOOPCOEFS(0,9),1
 
408
     $     ,94,H)
 
409
C         Coefficient construction for loop diagram with ID 13
 
410
          CALL MP_FFV1L2_1(PL(0,0),W(1,5),GC_5,ZERO,ZERO,PL(0,25)
 
411
     $     ,COEFS)
 
412
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,25))
 
413
          CALL MP_FFV1L2_1(PL(0,25),W(1,8),GC_5,ZERO,ZERO,PL(0,26)
 
414
     $     ,COEFS)
 
415
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,25),4,COEFS,4,4,WL(1,0,1,26))
 
416
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,26),2,4,LOOPCOEFS(0,10),1
 
417
     $     ,95,H)
 
418
C         Coefficient construction for loop diagram with ID 14
 
419
          CALL MP_FFV1L1_2(PL(0,0),W(1,1),GC_5,ZERO,ZERO,PL(0,27)
 
420
     $     ,COEFS)
 
421
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,27))
 
422
          CALL MP_FFV1L1_2(PL(0,27),W(1,2),GC_5,ZERO,ZERO,PL(0,28)
 
423
     $     ,COEFS)
 
424
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,27),4,COEFS,4,4,WL(1,0,1,28))
 
425
          CALL MP_FFV1L1_2(PL(0,28),W(1,8),GC_5,ZERO,ZERO,PL(0,29)
 
426
     $     ,COEFS)
 
427
          CALL MP_UPDATE_WL_2_1(WL(1,0,1,28),4,COEFS,4,4,WL(1,0,1,29))
 
428
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,29),3,4,LOOPCOEFS(0,11),1
 
429
     $     ,96,H)
 
430
C         Coefficient construction for loop diagram with ID 15
 
431
          CALL MP_FFV1L2_1(PL(0,0),W(1,1),GC_5,ZERO,ZERO,PL(0,30)
 
432
     $     ,COEFS)
 
433
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,30))
 
434
          CALL MP_FFV1L2_1(PL(0,30),W(1,2),GC_5,ZERO,ZERO,PL(0,31)
 
435
     $     ,COEFS)
 
436
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,30),4,COEFS,4,4,WL(1,0,1,31))
 
437
          CALL MP_FFV1L2_1(PL(0,31),W(1,8),GC_5,ZERO,ZERO,PL(0,32)
 
438
     $     ,COEFS)
 
439
          CALL MP_UPDATE_WL_2_1(WL(1,0,1,31),4,COEFS,4,4,WL(1,0,1,32))
 
440
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,32),3,4,LOOPCOEFS(0,12),1
 
441
     $     ,97,H)
 
442
C         Coefficient construction for loop diagram with ID 16
 
443
          CALL MP_FFV1L2_1(PL(0,0),W(1,5),GC_5,MDL_MB,ZERO,PL(0,33)
 
444
     $     ,COEFS)
 
445
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,33))
 
446
          CALL MP_FFV1L2_1(PL(0,33),W(1,8),GC_5,MDL_MB,ZERO,PL(0,34)
 
447
     $     ,COEFS)
 
448
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,33),4,COEFS,4,4,WL(1,0,1,34))
 
449
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,34),2,4,LOOPCOEFS(0,13),1
 
450
     $     ,98,H)
 
451
C         Coefficient construction for loop diagram with ID 17
 
452
          CALL MP_FFV1L1_2(PL(0,0),W(1,1),GC_5,MDL_MB,ZERO,PL(0,35)
 
453
     $     ,COEFS)
 
454
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,35))
 
455
          CALL MP_FFV1L1_2(PL(0,35),W(1,2),GC_5,MDL_MB,ZERO,PL(0,36)
 
456
     $     ,COEFS)
 
457
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,35),4,COEFS,4,4,WL(1,0,1,36))
 
458
          CALL MP_FFV1L1_2(PL(0,36),W(1,8),GC_5,MDL_MB,ZERO,PL(0,37)
 
459
     $     ,COEFS)
 
460
          CALL MP_UPDATE_WL_2_1(WL(1,0,1,36),4,COEFS,4,4,WL(1,0,1,37))
 
461
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,37),3,4,LOOPCOEFS(0,14),1
 
462
     $     ,99,H)
 
463
C         Coefficient construction for loop diagram with ID 18
 
464
          CALL MP_FFV1L2_1(PL(0,0),W(1,1),GC_5,MDL_MB,ZERO,PL(0,38)
 
465
     $     ,COEFS)
 
466
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,38))
 
467
          CALL MP_FFV1L2_1(PL(0,38),W(1,2),GC_5,MDL_MB,ZERO,PL(0,39)
 
468
     $     ,COEFS)
 
469
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,38),4,COEFS,4,4,WL(1,0,1,39))
 
470
          CALL MP_FFV1L2_1(PL(0,39),W(1,8),GC_5,MDL_MB,ZERO,PL(0,40)
 
471
     $     ,COEFS)
 
472
          CALL MP_UPDATE_WL_2_1(WL(1,0,1,39),4,COEFS,4,4,WL(1,0,1,40))
 
473
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,40),3,4,LOOPCOEFS(0,15),1
 
474
     $     ,100,H)
 
475
C         Coefficient construction for loop diagram with ID 19
 
476
          CALL MP_FFV1L2_1(PL(0,0),W(1,5),GC_5,MDL_MT,MDL_WT,PL(0,41)
 
477
     $     ,COEFS)
 
478
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,41))
 
479
          CALL MP_FFV1L2_1(PL(0,41),W(1,8),GC_5,MDL_MT,MDL_WT,PL(0,42)
 
480
     $     ,COEFS)
 
481
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,41),4,COEFS,4,4,WL(1,0,1,42))
 
482
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,42),2,4,LOOPCOEFS(0,16),1
 
483
     $     ,101,H)
 
484
C         Coefficient construction for loop diagram with ID 20
 
485
          CALL MP_FFV1L1P0_3(PL(0,0),W(1,3),GC_5,ZERO,ZERO,PL(0,43)
 
486
     $     ,COEFS)
 
487
          CALL MP_UPDATE_WL_0_0(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,43))
 
488
          CALL MP_FFV1L3_2(PL(0,43),W(1,4),GC_5,MDL_MT,MDL_WT,PL(0,44)
 
489
     $     ,COEFS)
 
490
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,43),4,COEFS,4,4,WL(1,0,1,44))
 
491
          CALL MP_FFV1L1_2(PL(0,44),W(1,5),GC_5,MDL_MT,MDL_WT,PL(0,45)
 
492
     $     ,COEFS)
 
493
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,44),4,COEFS,4,4,WL(1,0,1,45))
 
494
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,45),2,4,LOOPCOEFS(0,17),1
 
495
     $     ,102,H)
 
496
C         Coefficient construction for loop diagram with ID 21
 
497
          CALL MP_FFV1L3_1(PL(0,0),W(1,3),GC_5,MDL_MT,MDL_WT,PL(0,46)
 
498
     $     ,COEFS)
 
499
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,46))
 
500
          CALL MP_FFV1L2P0_3(PL(0,46),W(1,4),GC_5,ZERO,ZERO,PL(0,47)
 
501
     $     ,COEFS)
 
502
          CALL MP_UPDATE_WL_1_0(WL(1,0,1,46),4,COEFS,4,4,WL(1,0,1,47))
 
503
          CALL MP_VVV1L2P0_1(PL(0,47),W(1,5),GC_4,ZERO,ZERO,PL(0,48)
 
504
     $     ,COEFS)
 
505
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,47),4,COEFS,4,4,WL(1,0,1,48))
 
506
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,48),2,4,LOOPCOEFS(0,18),1
 
507
     $     ,103,H)
 
508
C         Coefficient construction for loop diagram with ID 22
 
509
          CALL MP_FFV1L1P0_3(PL(0,0),W(1,6),GC_5,ZERO,ZERO,PL(0,49)
 
510
     $     ,COEFS)
 
511
          CALL MP_UPDATE_WL_0_0(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,49))
 
512
          CALL MP_FFV1L3_2(PL(0,49),W(1,9),GC_5,MDL_MT,MDL_WT,PL(0,50)
 
513
     $     ,COEFS)
 
514
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,49),4,COEFS,4,4,WL(1,0,1,50))
 
515
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,50),1,4,LOOPCOEFS(0,19),1
 
516
     $     ,104,H)
 
517
C         Coefficient construction for loop diagram with ID 23
 
518
          CALL MP_FFV1L2_1(PL(0,0),W(1,2),GC_5,MDL_MT,MDL_WT,PL(0,51)
 
519
     $     ,COEFS)
 
520
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,51))
 
521
          CALL MP_FFV1L2P0_3(PL(0,51),W(1,4),GC_5,ZERO,ZERO,PL(0,52)
 
522
     $     ,COEFS)
 
523
          CALL MP_UPDATE_WL_1_0(WL(1,0,1,51),4,COEFS,4,4,WL(1,0,1,52))
 
524
          CALL MP_FFV1L3_1(PL(0,52),W(1,6),GC_5,MDL_MT,MDL_WT,PL(0,53)
 
525
     $     ,COEFS)
 
526
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,52),4,COEFS,4,4,WL(1,0,1,53))
 
527
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,53),2,4,LOOPCOEFS(0,20),1
 
528
     $     ,105,H)
 
529
C         Coefficient construction for loop diagram with ID 24
 
530
          CALL MP_VVV1L2P0_1(PL(0,0),W(1,2),GC_4,ZERO,ZERO,PL(0,54)
 
531
     $     ,COEFS)
 
532
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,54))
 
533
          CALL MP_FFV1L3_2(PL(0,54),W(1,4),GC_5,MDL_MT,MDL_WT,PL(0,55)
 
534
     $     ,COEFS)
 
535
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,54),4,COEFS,4,4,WL(1,0,1,55))
 
536
          CALL MP_FFV1L1P0_3(PL(0,55),W(1,6),GC_5,ZERO,ZERO,PL(0,56)
 
537
     $     ,COEFS)
 
538
          CALL MP_UPDATE_WL_2_0(WL(1,0,1,55),4,COEFS,4,4,WL(1,0,1,56))
 
539
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,56),2,4,LOOPCOEFS(0,21),1
 
540
     $     ,106,H)
 
541
C         Coefficient construction for loop diagram with ID 25
 
542
          CALL MP_FFV1L1P0_3(PL(0,0),W(1,10),GC_5,ZERO,ZERO,PL(0,57)
 
543
     $     ,COEFS)
 
544
          CALL MP_UPDATE_WL_0_0(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,57))
 
545
          CALL MP_FFV1L3_2(PL(0,57),W(1,7),GC_5,MDL_MT,MDL_WT,PL(0,58)
 
546
     $     ,COEFS)
 
547
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,57),4,COEFS,4,4,WL(1,0,1,58))
 
548
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,58),1,4,LOOPCOEFS(0,22),1
 
549
     $     ,107,H)
 
550
C         Coefficient construction for loop diagram with ID 26
 
551
          CALL MP_FFV1L1_2(PL(0,0),W(1,2),GC_5,MDL_MT,MDL_WT,PL(0,59)
 
552
     $     ,COEFS)
 
553
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,59))
 
554
          CALL MP_FFV1L1P0_3(PL(0,59),W(1,3),GC_5,ZERO,ZERO,PL(0,60)
 
555
     $     ,COEFS)
 
556
          CALL MP_UPDATE_WL_1_0(WL(1,0,1,59),4,COEFS,4,4,WL(1,0,1,60))
 
557
          CALL MP_FFV1L3_2(PL(0,60),W(1,7),GC_5,MDL_MT,MDL_WT,PL(0,61)
 
558
     $     ,COEFS)
 
559
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,60),4,COEFS,4,4,WL(1,0,1,61))
 
560
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,61),2,4,LOOPCOEFS(0,23),1
 
561
     $     ,108,H)
 
562
C         Coefficient construction for loop diagram with ID 27
 
563
          CALL MP_FFV1L3_1(PL(0,54),W(1,3),GC_5,MDL_MT,MDL_WT,PL(0,62)
 
564
     $     ,COEFS)
 
565
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,54),4,COEFS,4,4,WL(1,0,1,62))
 
566
          CALL MP_FFV1L2P0_3(PL(0,62),W(1,7),GC_5,ZERO,ZERO,PL(0,63)
 
567
     $     ,COEFS)
 
568
          CALL MP_UPDATE_WL_2_0(WL(1,0,1,62),4,COEFS,4,4,WL(1,0,1,63))
 
569
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,63),2,4,LOOPCOEFS(0,24),1
 
570
     $     ,109,H)
 
571
C         Coefficient construction for loop diagram with ID 28
 
572
          CALL MP_FFV1L2_1(PL(0,0),W(1,1),GC_5,MDL_MT,MDL_WT,PL(0,64)
 
573
     $     ,COEFS)
 
574
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,64))
 
575
          CALL MP_FFV1L2P0_3(PL(0,64),W(1,4),GC_5,ZERO,ZERO,PL(0,65)
 
576
     $     ,COEFS)
 
577
          CALL MP_UPDATE_WL_1_0(WL(1,0,1,64),4,COEFS,4,4,WL(1,0,1,65))
 
578
          CALL MP_FFV1L3_1(PL(0,65),W(1,10),GC_5,MDL_MT,MDL_WT,PL(0
 
579
     $     ,66),COEFS)
 
580
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,65),4,COEFS,4,4,WL(1,0,1,66))
 
581
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,66),2,4,LOOPCOEFS(0,25),1
 
582
     $     ,110,H)
 
583
C         Coefficient construction for loop diagram with ID 29
 
584
          CALL MP_FFV1L1_2(PL(0,0),W(1,1),GC_5,MDL_MT,MDL_WT,PL(0,67)
 
585
     $     ,COEFS)
 
586
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,67))
 
587
          CALL MP_FFV1L1P0_3(PL(0,67),W(1,3),GC_5,ZERO,ZERO,PL(0,68)
 
588
     $     ,COEFS)
 
589
          CALL MP_UPDATE_WL_1_0(WL(1,0,1,67),4,COEFS,4,4,WL(1,0,1,68))
 
590
          CALL MP_FFV1L3_2(PL(0,68),W(1,9),GC_5,MDL_MT,MDL_WT,PL(0,69)
 
591
     $     ,COEFS)
 
592
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,68),4,COEFS,4,4,WL(1,0,1,69))
 
593
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,69),2,4,LOOPCOEFS(0,26),1
 
594
     $     ,111,H)
 
595
C         Coefficient construction for loop diagram with ID 30
 
596
          CALL MP_FFV1L1_2(PL(0,67),W(1,2),GC_5,MDL_MT,MDL_WT,PL(0,70)
 
597
     $     ,COEFS)
 
598
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,67),4,COEFS,4,4,WL(1,0,1,70))
 
599
          CALL MP_FFV1L1P0_3(PL(0,70),W(1,3),GC_5,ZERO,ZERO,PL(0,71)
 
600
     $     ,COEFS)
 
601
          CALL MP_UPDATE_WL_2_0(WL(1,0,1,70),4,COEFS,4,4,WL(1,0,1,71))
 
602
          CALL MP_FFV1L3_2(PL(0,71),W(1,4),GC_5,MDL_MT,MDL_WT,PL(0,72)
 
603
     $     ,COEFS)
 
604
          CALL MP_UPDATE_WL_2_1(WL(1,0,1,71),4,COEFS,4,4,WL(1,0,1,72))
 
605
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,72),3,4,LOOPCOEFS(0,27),1
 
606
     $     ,112,H)
 
607
C         Coefficient construction for loop diagram with ID 31
 
608
          CALL MP_FFV1L1_2(PL(0,70),W(1,8),GC_5,MDL_MT,MDL_WT,PL(0,73)
 
609
     $     ,COEFS)
 
610
          CALL MP_UPDATE_WL_2_1(WL(1,0,1,70),4,COEFS,4,4,WL(1,0,1,73))
 
611
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,73),3,4,LOOPCOEFS(0,28),1
 
612
     $     ,113,H)
 
613
C         Coefficient construction for loop diagram with ID 32
 
614
          CALL MP_FFV1L2_1(PL(0,64),W(1,2),GC_5,MDL_MT,MDL_WT,PL(0,74)
 
615
     $     ,COEFS)
 
616
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,64),4,COEFS,4,4,WL(1,0,1,74))
 
617
          CALL MP_FFV1L2_1(PL(0,74),W(1,8),GC_5,MDL_MT,MDL_WT,PL(0,75)
 
618
     $     ,COEFS)
 
619
          CALL MP_UPDATE_WL_2_1(WL(1,0,1,74),4,COEFS,4,4,WL(1,0,1,75))
 
620
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,75),3,4,LOOPCOEFS(0,29),1
 
621
     $     ,114,H)
 
622
C         Coefficient construction for loop diagram with ID 33
 
623
          CALL MP_FFV1L2P0_3(PL(0,74),W(1,4),GC_5,ZERO,ZERO,PL(0,76)
 
624
     $     ,COEFS)
 
625
          CALL MP_UPDATE_WL_2_0(WL(1,0,1,74),4,COEFS,4,4,WL(1,0,1,76))
 
626
          CALL MP_FFV1L3_1(PL(0,76),W(1,3),GC_5,MDL_MT,MDL_WT,PL(0,77)
 
627
     $     ,COEFS)
 
628
          CALL MP_UPDATE_WL_2_1(WL(1,0,1,76),4,COEFS,4,4,WL(1,0,1,77))
 
629
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,77),3,4,LOOPCOEFS(0,30),1
 
630
     $     ,115,H)
 
631
C         Coefficient construction for loop diagram with ID 34
 
632
          CALL MP_VVV1L2P0_1(PL(0,68),W(1,2),GC_4,ZERO,ZERO,PL(0,78)
 
633
     $     ,COEFS)
 
634
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,68),4,COEFS,4,4,WL(1,0,1,78))
 
635
          CALL MP_FFV1L3_2(PL(0,78),W(1,4),GC_5,MDL_MT,MDL_WT,PL(0,79)
 
636
     $     ,COEFS)
 
637
          CALL MP_UPDATE_WL_2_1(WL(1,0,1,78),4,COEFS,4,4,WL(1,0,1,79))
 
638
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,79),3,4,LOOPCOEFS(0,31),1
 
639
     $     ,116,H)
 
640
C         Coefficient construction for loop diagram with ID 35
 
641
          CALL MP_VVV1L2P0_1(PL(0,0),W(1,1),GC_4,ZERO,ZERO,PL(0,80)
 
642
     $     ,COEFS)
 
643
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,80))
 
644
          CALL MP_FFV1L3_2(PL(0,80),W(1,4),GC_5,MDL_MT,MDL_WT,PL(0,81)
 
645
     $     ,COEFS)
 
646
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,80),4,COEFS,4,4,WL(1,0,1,81))
 
647
          CALL MP_FFV1L1P0_3(PL(0,81),W(1,10),GC_5,ZERO,ZERO,PL(0,82)
 
648
     $     ,COEFS)
 
649
          CALL MP_UPDATE_WL_2_0(WL(1,0,1,81),4,COEFS,4,4,WL(1,0,1,82))
 
650
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,82),2,4,LOOPCOEFS(0,32),1
 
651
     $     ,117,H)
 
652
C         Coefficient construction for loop diagram with ID 36
 
653
          CALL MP_FFV1L3_1(PL(0,80),W(1,3),GC_5,MDL_MT,MDL_WT,PL(0,83)
 
654
     $     ,COEFS)
 
655
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,80),4,COEFS,4,4,WL(1,0,1,83))
 
656
          CALL MP_FFV1L2P0_3(PL(0,83),W(1,9),GC_5,ZERO,ZERO,PL(0,84)
 
657
     $     ,COEFS)
 
658
          CALL MP_UPDATE_WL_2_0(WL(1,0,1,83),4,COEFS,4,4,WL(1,0,1,84))
 
659
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,84),2,4,LOOPCOEFS(0,33),1
 
660
     $     ,118,H)
 
661
C         Coefficient construction for loop diagram with ID 37
 
662
          CALL MP_FFV1L2_1(PL(0,83),W(1,2),GC_5,MDL_MT,MDL_WT,PL(0,85)
 
663
     $     ,COEFS)
 
664
          CALL MP_UPDATE_WL_2_1(WL(1,0,1,83),4,COEFS,4,4,WL(1,0,1,85))
 
665
          CALL MP_FFV1L2P0_3(PL(0,85),W(1,4),GC_5,ZERO,ZERO,PL(0,86)
 
666
     $     ,COEFS)
 
667
          CALL MP_UPDATE_WL_3_0(WL(1,0,1,85),4,COEFS,4,4,WL(1,0,1,86))
 
668
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,86),3,4,LOOPCOEFS(0,34),1
 
669
     $     ,119,H)
 
670
C         Coefficient construction for loop diagram with ID 38
 
671
          CALL MP_VVVV1L2P0_1(PL(0,47),W(1,1),W(1,2),GC_6,ZERO,ZERO
 
672
     $     ,PL(0,87),COEFS)
 
673
          CALL MP_UPDATE_WL_1_0(WL(1,0,1,47),4,COEFS,4,4,WL(1,0,1,87))
 
674
          CALL MP_VVVV3L2P0_1(PL(0,47),W(1,1),W(1,2),GC_6,ZERO,ZERO
 
675
     $     ,PL(0,88),COEFS)
 
676
          CALL MP_UPDATE_WL_1_0(WL(1,0,1,47),4,COEFS,4,4,WL(1,0,1,88))
 
677
          CALL MP_VVVV4L2P0_1(PL(0,47),W(1,1),W(1,2),GC_6,ZERO,ZERO
 
678
     $     ,PL(0,89),COEFS)
 
679
          CALL MP_UPDATE_WL_1_0(WL(1,0,1,47),4,COEFS,4,4,WL(1,0,1,89))
 
680
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,87),1,4,LOOPCOEFS(0,35),1
 
681
     $     ,120,H)
 
682
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,88),1,4,LOOPCOEFS(0,36),1
 
683
     $     ,121,H)
 
684
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,89),1,4,LOOPCOEFS(0,37),1
 
685
     $     ,122,H)
 
686
C         Coefficient construction for loop diagram with ID 39
 
687
          CALL MP_VVV1L2P0_1(PL(0,80),W(1,2),GC_4,ZERO,ZERO,PL(0,90)
 
688
     $     ,COEFS)
 
689
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,80),4,COEFS,4,4,WL(1,0,1,90))
 
690
          CALL MP_FFV1L3_2(PL(0,90),W(1,4),GC_5,MDL_MT,MDL_WT,PL(0,91)
 
691
     $     ,COEFS)
 
692
          CALL MP_UPDATE_WL_2_1(WL(1,0,1,90),4,COEFS,4,4,WL(1,0,1,91))
 
693
          CALL MP_FFV1L1P0_3(PL(0,91),W(1,3),GC_5,ZERO,ZERO,PL(0,92)
 
694
     $     ,COEFS)
 
695
          CALL MP_UPDATE_WL_3_0(WL(1,0,1,91),4,COEFS,4,4,WL(1,0,1,92))
 
696
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,92),3,4,LOOPCOEFS(0,38),1
 
697
     $     ,123,H)
 
698
C         Coefficient construction for loop diagram with ID 40
 
699
          CALL MP_FFV1L3_1(PL(0,90),W(1,3),GC_5,MDL_MT,MDL_WT,PL(0,93)
 
700
     $     ,COEFS)
 
701
          CALL MP_UPDATE_WL_2_1(WL(1,0,1,90),4,COEFS,4,4,WL(1,0,1,93))
 
702
          CALL MP_FFV1L2P0_3(PL(0,93),W(1,4),GC_5,ZERO,ZERO,PL(0,94)
 
703
     $     ,COEFS)
 
704
          CALL MP_UPDATE_WL_3_0(WL(1,0,1,93),4,COEFS,4,4,WL(1,0,1,94))
 
705
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,94),3,4,LOOPCOEFS(0,39),1
 
706
     $     ,124,H)
 
707
C         Coefficient construction for loop diagram with ID 41
 
708
          CALL MP_VVV1L2P0_1(PL(0,0),W(1,5),GC_4,ZERO,ZERO,PL(0,95)
 
709
     $     ,COEFS)
 
710
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,95))
 
711
          CALL MP_VVV1L2P0_1(PL(0,95),W(1,8),GC_4,ZERO,ZERO,PL(0,96)
 
712
     $     ,COEFS)
 
713
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,95),4,COEFS,4,4,WL(1,0,1,96))
 
714
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,96),2,4,LOOPCOEFS(0,40),2
 
715
     $     ,125,H)
 
716
C         Coefficient construction for loop diagram with ID 42
 
717
          CALL MP_VVV1L2P0_1(PL(0,90),W(1,8),GC_4,ZERO,ZERO,PL(0,97)
 
718
     $     ,COEFS)
 
719
          CALL MP_UPDATE_WL_2_1(WL(1,0,1,90),4,COEFS,4,4,WL(1,0,1,97))
 
720
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,97),3,4,LOOPCOEFS(0,41),1
 
721
     $     ,126,H)
 
722
C         Coefficient construction for loop diagram with ID 43
 
723
          CALL MP_VVVV1L2P0_1(PL(0,80),W(1,8),W(1,2),GC_6,ZERO,ZERO
 
724
     $     ,PL(0,98),COEFS)
 
725
          CALL MP_UPDATE_WL_1_0(WL(1,0,1,80),4,COEFS,4,4,WL(1,0,1,98))
 
726
          CALL MP_VVVV3L2P0_1(PL(0,80),W(1,8),W(1,2),GC_6,ZERO,ZERO
 
727
     $     ,PL(0,99),COEFS)
 
728
          CALL MP_UPDATE_WL_1_0(WL(1,0,1,80),4,COEFS,4,4,WL(1,0,1,99))
 
729
          CALL MP_VVVV4L2P0_1(PL(0,80),W(1,8),W(1,2),GC_6,ZERO,ZERO
 
730
     $     ,PL(0,100),COEFS)
 
731
          CALL MP_UPDATE_WL_1_0(WL(1,0,1,80),4,COEFS,4,4,WL(1,0,1,100))
 
732
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,98),1,4,LOOPCOEFS(0,42),2
 
733
     $     ,127,H)
 
734
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,99),1,4,LOOPCOEFS(0,43),2
 
735
     $     ,128,H)
 
736
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,100),1,4,LOOPCOEFS(0,44)
 
737
     $     ,2,129,H)
 
738
C         Coefficient construction for loop diagram with ID 44
 
739
          CALL MP_VVVV1L2P0_1(PL(0,54),W(1,8),W(1,1),GC_6,ZERO,ZERO
 
740
     $     ,PL(0,101),COEFS)
 
741
          CALL MP_UPDATE_WL_1_0(WL(1,0,1,54),4,COEFS,4,4,WL(1,0,1,101))
 
742
          CALL MP_VVVV3L2P0_1(PL(0,54),W(1,8),W(1,1),GC_6,ZERO,ZERO
 
743
     $     ,PL(0,102),COEFS)
 
744
          CALL MP_UPDATE_WL_1_0(WL(1,0,1,54),4,COEFS,4,4,WL(1,0,1,102))
 
745
          CALL MP_VVVV4L2P0_1(PL(0,54),W(1,8),W(1,1),GC_6,ZERO,ZERO
 
746
     $     ,PL(0,103),COEFS)
 
747
          CALL MP_UPDATE_WL_1_0(WL(1,0,1,54),4,COEFS,4,4,WL(1,0,1,103))
 
748
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,101),1,4,LOOPCOEFS(0,45)
 
749
     $     ,2,130,H)
 
750
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,102),1,4,LOOPCOEFS(0,46)
 
751
     $     ,2,131,H)
 
752
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,103),1,4,LOOPCOEFS(0,47)
 
753
     $     ,2,132,H)
 
754
C         Coefficient construction for loop diagram with ID 45
 
755
          CALL MP_VVVV1L2P0_1(PL(0,0),W(1,1),W(1,2),GC_6,ZERO,ZERO
 
756
     $     ,PL(0,104),COEFS)
 
757
          CALL MP_UPDATE_WL_0_0(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,104))
 
758
          CALL MP_VVVV3L2P0_1(PL(0,0),W(1,1),W(1,2),GC_6,ZERO,ZERO
 
759
     $     ,PL(0,105),COEFS)
 
760
          CALL MP_UPDATE_WL_0_0(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,105))
 
761
          CALL MP_VVVV4L2P0_1(PL(0,0),W(1,1),W(1,2),GC_6,ZERO,ZERO
 
762
     $     ,PL(0,106),COEFS)
 
763
          CALL MP_UPDATE_WL_0_0(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1,106))
 
764
          CALL MP_VVV1L2P0_1(PL(0,104),W(1,8),GC_4,ZERO,ZERO,PL(0,107)
 
765
     $     ,COEFS)
 
766
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,104),4,COEFS,4,4,WL(1,0,1
 
767
     $     ,107))
 
768
          CALL MP_VVV1L2P0_1(PL(0,105),W(1,8),GC_4,ZERO,ZERO,PL(0,108)
 
769
     $     ,COEFS)
 
770
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,105),4,COEFS,4,4,WL(1,0,1
 
771
     $     ,108))
 
772
          CALL MP_VVV1L2P0_1(PL(0,106),W(1,8),GC_4,ZERO,ZERO,PL(0,109)
 
773
     $     ,COEFS)
 
774
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,106),4,COEFS,4,4,WL(1,0,1
 
775
     $     ,109))
 
776
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,107),1,4,LOOPCOEFS(0,48)
 
777
     $     ,2,133,H)
 
778
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,108),1,4,LOOPCOEFS(0,49)
 
779
     $     ,2,134,H)
 
780
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,109),1,4,LOOPCOEFS(0,50)
 
781
     $     ,2,135,H)
 
782
C         Coefficient construction for loop diagram with ID 46
 
783
          CALL MP_GHGHGL2_1(PL(0,0),W(1,5),GC_4,ZERO,ZERO,PL(0,110)
 
784
     $     ,COEFS)
 
785
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),1,COEFS,1,1,WL(1,0,1,110))
 
786
          CALL MP_GHGHGL2_1(PL(0,110),W(1,8),GC_4,ZERO,ZERO,PL(0,111)
 
787
     $     ,COEFS)
 
788
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,110),1,COEFS,1,1,WL(1,0,1
 
789
     $     ,111))
 
790
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,111),2,1,LOOPCOEFS(0,51)
 
791
     $     ,1,136,H)
 
792
C         Coefficient construction for loop diagram with ID 47
 
793
          CALL MP_GHGHGL1_2(PL(0,0),W(1,1),GC_4,ZERO,ZERO,PL(0,112)
 
794
     $     ,COEFS)
 
795
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),1,COEFS,1,1,WL(1,0,1,112))
 
796
          CALL MP_GHGHGL1_2(PL(0,112),W(1,2),GC_4,ZERO,ZERO,PL(0,113)
 
797
     $     ,COEFS)
 
798
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,112),1,COEFS,1,1,WL(1,0,1
 
799
     $     ,113))
 
800
          CALL MP_GHGHGL1_2(PL(0,113),W(1,8),GC_4,ZERO,ZERO,PL(0,114)
 
801
     $     ,COEFS)
 
802
          CALL MP_UPDATE_WL_2_1(WL(1,0,1,113),1,COEFS,1,1,WL(1,0,1
 
803
     $     ,114))
 
804
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,114),3,1,LOOPCOEFS(0,52)
 
805
     $     ,1,137,H)
 
806
C         Coefficient construction for loop diagram with ID 48
 
807
          CALL MP_GHGHGL2_1(PL(0,0),W(1,1),GC_4,ZERO,ZERO,PL(0,115)
 
808
     $     ,COEFS)
 
809
          CALL MP_UPDATE_WL_0_1(WL(1,0,1,0),1,COEFS,1,1,WL(1,0,1,115))
 
810
          CALL MP_GHGHGL2_1(PL(0,115),W(1,2),GC_4,ZERO,ZERO,PL(0,116)
 
811
     $     ,COEFS)
 
812
          CALL MP_UPDATE_WL_1_1(WL(1,0,1,115),1,COEFS,1,1,WL(1,0,1
 
813
     $     ,116))
 
814
          CALL MP_GHGHGL2_1(PL(0,116),W(1,8),GC_4,ZERO,ZERO,PL(0,117)
 
815
     $     ,COEFS)
 
816
          CALL MP_UPDATE_WL_2_1(WL(1,0,1,116),1,COEFS,1,1,WL(1,0,1
 
817
     $     ,117))
 
818
          CALL MP_CREATE_LOOP_COEFS(WL(1,0,1,117),3,1,LOOPCOEFS(0,53)
 
819
     $     ,1,138,H)
 
820
        ENDIF
 
821
      ENDDO
 
822
 
 
823
C     Copy the qp wfs to the dp ones as they are used to setup the CT
 
824
C      calls.
 
825
      DO I=1,NWAVEFUNCS
 
826
        DO J=1,MAXLWFSIZE+4
 
827
          DPW(J,I)=W(J,I)
 
828
        ENDDO
 
829
      ENDDO
 
830
 
 
831
      DO I=1,3
 
832
        ANSDP(I)=ANS(I)
 
833
      ENDDO
 
834
 
 
835
      CALL MP_ADD_COEFS(LOOPCOEFS(0,1),2,LOOPCOEFS(0,4),2)
 
836
      CALL MP_ADD_COEFS(LOOPCOEFS(0,1),2,LOOPCOEFS(0,7),2)
 
837
      CALL MP_ADD_COEFS(LOOPCOEFS(0,1),2,LOOPCOEFS(0,10),2)
 
838
      CALL MP_ADD_COEFS(LOOPCOEFS(0,1),2,LOOPCOEFS(0,40),2)
 
839
      CALL MP_ADD_COEFS(LOOPCOEFS(0,1),2,LOOPCOEFS(0,51),2)
 
840
      CALL MP_ADD_COEFS(LOOPCOEFS(0,1),2,LOOPCOEFS(0,48),1)
 
841
      CALL MP_ADD_COEFS(LOOPCOEFS(0,1),2,LOOPCOEFS(0,49),1)
 
842
      CALL MP_ADD_COEFS(LOOPCOEFS(0,1),2,LOOPCOEFS(0,50),1)
 
843
      CALL MP_ADD_COEFS(LOOPCOEFS(0,2),3,LOOPCOEFS(0,3),3)
 
844
      CALL MP_ADD_COEFS(LOOPCOEFS(0,2),3,LOOPCOEFS(0,5),3)
 
845
      CALL MP_ADD_COEFS(LOOPCOEFS(0,2),3,LOOPCOEFS(0,6),3)
 
846
      CALL MP_ADD_COEFS(LOOPCOEFS(0,2),3,LOOPCOEFS(0,8),3)
 
847
      CALL MP_ADD_COEFS(LOOPCOEFS(0,2),3,LOOPCOEFS(0,9),3)
 
848
      CALL MP_ADD_COEFS(LOOPCOEFS(0,2),3,LOOPCOEFS(0,11),3)
 
849
      CALL MP_ADD_COEFS(LOOPCOEFS(0,2),3,LOOPCOEFS(0,12),3)
 
850
      CALL MP_ADD_COEFS(LOOPCOEFS(0,2),3,LOOPCOEFS(0,41),3)
 
851
      CALL MP_ADD_COEFS(LOOPCOEFS(0,2),3,LOOPCOEFS(0,52),3)
 
852
      CALL MP_ADD_COEFS(LOOPCOEFS(0,2),3,LOOPCOEFS(0,53),3)
 
853
      CALL MP_ADD_COEFS(LOOPCOEFS(0,14),3,LOOPCOEFS(0,15),3)
 
854
      CALL MP_ADD_COEFS(LOOPCOEFS(0,18),2,LOOPCOEFS(0,35),1)
 
855
      CALL MP_ADD_COEFS(LOOPCOEFS(0,18),2,LOOPCOEFS(0,36),1)
 
856
      CALL MP_ADD_COEFS(LOOPCOEFS(0,18),2,LOOPCOEFS(0,37),1)
 
857
      CALL MP_ADD_COEFS(LOOPCOEFS(0,28),3,LOOPCOEFS(0,29),3)
 
858
      CALL MP_ADD_COEFS(LOOPCOEFS(0,42),1,LOOPCOEFS(0,43),1)
 
859
      CALL MP_ADD_COEFS(LOOPCOEFS(0,42),1,LOOPCOEFS(0,44),1)
 
860
      CALL MP_ADD_COEFS(LOOPCOEFS(0,45),1,LOOPCOEFS(0,46),1)
 
861
      CALL MP_ADD_COEFS(LOOPCOEFS(0,45),1,LOOPCOEFS(0,47),1)
 
862
 
 
863
      END
 
864