← 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_compute_loop_coefs.f

  • Committer: olivier Mattelaer
  • Date: 2015-03-05 00:14:16 UTC
  • mfrom: (258.1.9 2.3)
  • mto: (258.8.1 2.3)
  • mto: This revision was merged to the branch mainline in revision 259.
  • Revision ID: olivier.mattelaer@uclouvain.be-20150305001416-y9mzeykfzwnl9t0j
partial merge

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
19
19
      INTEGER NBORNAMPS
20
20
      PARAMETER (NBORNAMPS=3)
21
21
      INTEGER    NLOOPS, NLOOPGROUPS, NCTAMPS
22
 
      PARAMETER (NLOOPS=53, NLOOPGROUPS=26, NCTAMPS=85)
 
22
      PARAMETER (NLOOPS=44, NLOOPGROUPS=26, NCTAMPS=85)
 
23
      INTEGER    NLOOPAMPS
 
24
      PARAMETER (NLOOPAMPS=129)
23
25
      INTEGER    NCOLORROWS
24
 
      PARAMETER (NCOLORROWS=138)
 
26
      PARAMETER (NCOLORROWS=NLOOPAMPS)
25
27
      INTEGER    NEXTERNAL
26
28
      PARAMETER (NEXTERNAL=4)
27
29
      INTEGER    NWAVEFUNCS,NLOOPWAVEFUNCS
28
 
      PARAMETER (NWAVEFUNCS=10,NLOOPWAVEFUNCS=117)
 
30
      PARAMETER (NWAVEFUNCS=10,NLOOPWAVEFUNCS=93)
29
31
      INTEGER MAXLWFSIZE
30
32
      PARAMETER (MAXLWFSIZE=4)
31
33
      INTEGER LOOPMAXCOEFS, VERTEXMAXCOEFS
39
41
C     These are constants related to the split orders
40
42
      INTEGER    NSO, NSQUAREDSO, NAMPSO
41
43
      PARAMETER (NSO=0, NSQUAREDSO=0, NAMPSO=0)
 
44
 
 
45
 
42
46
C     
43
47
C     ARGUMENTS 
44
48
C     
47
51
C     
48
52
C     LOCAL VARIABLES 
49
53
C     
50
 
      INTEGER I,J,K,H,DUMMY,ITEMP
 
54
      LOGICAL DPW_COPIED
 
55
      INTEGER I,J,K,H,HEL_MULT,ITEMP
51
56
      REAL*16 TEMP2
52
57
      INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
53
58
      REAL*16 P(0:3,NEXTERNAL)
77
82
      INTEGER HELPICKED
78
83
      COMMON/ML5_0_HELCHOICE/HELPICKED
79
84
 
 
85
      INTEGER USERHEL
 
86
      COMMON/ML5_0_USERCHOICE/USERHEL
 
87
 
80
88
      INTEGER SQSO_TARGET
81
89
      COMMON/ML5_0_SOCHOICE/SQSO_TARGET
82
90
 
89
97
 
90
98
      COMPLEX*32 AMP(NBORNAMPS)
91
99
      COMMON/ML5_0_MP_AMPS/AMP
 
100
      COMPLEX*16 DP_AMP(NBORNAMPS)
 
101
      COMMON/ML5_0_AMPS/DP_AMP
92
102
      COMPLEX*32 W(20,NWAVEFUNCS)
93
103
      COMMON/ML5_0_MP_W/W
94
104
 
103
113
      COMPLEX*32 LOOPCOEFS(0:LOOPMAXCOEFS-1,NSQUAREDSO,NLOOPGROUPS)
104
114
      COMMON/ML5_0_MP_LCOEFS/LOOPCOEFS
105
115
 
 
116
 
106
117
      COMPLEX*32 AMPL(3,NCTAMPS)
107
118
      COMMON/ML5_0_MP_AMPL/AMPL
108
119
 
 
120
 
109
121
      INTEGER CF_D(NCOLORROWS,NBORNAMPS)
110
122
      INTEGER CF_N(NCOLORROWS,NBORNAMPS)
111
123
      COMMON/ML5_0_CF/CF_D,CF_N
156
168
        ENDDO
157
169
      ENDDO
158
170
 
 
171
 
 
172
      DO I=1, NBORNAMPS
 
173
        DP_AMP(I) = (0.0D0,0.0D0)
 
174
        AMP(I) = (ZERO, ZERO)
 
175
      ENDDO
 
176
 
159
177
      DO I=1,NLOOPGROUPS
160
178
        DO J=0,LOOPMAXCOEFS-1
161
179
          DO K=1,NSQUAREDSO
171
189
        ENDDO
172
190
      ENDDO
173
191
 
 
192
      DPW_COPIED = .FALSE.
174
193
      DO H=1,NCOMB
175
194
        IF ((HELPICKED.EQ.H).OR.((HELPICKED.EQ.-1).AND.(CHECKPHASE.OR.(
176
195
     $   .NOT.HELDOUBLECHECKED).OR.(GOODHEL(H).GT.-HELOFFSET.AND.GOODHE
183
202
          MP_CT_REQ_SO_DONE=.FALSE.
184
203
          MP_LOOP_REQ_SO_DONE=.FALSE.
185
204
 
 
205
          IF (.NOT.CHECKPHASE.AND.HELDOUBLECHECKED.AND.HELPICKED.EQ.
 
206
     $     -1) THEN
 
207
            HEL_MULT=GOODHEL(H)
 
208
          ELSE
 
209
            HEL_MULT=1
 
210
          ENDIF
 
211
 
 
212
 
186
213
          CALL MP_VXXXXX(P(0,1),ZERO,NHEL(1),-1*IC(1),W(1,1))
187
214
          CALL MP_VXXXXX(P(0,2),ZERO,NHEL(2),-1*IC(2),W(1,2))
188
215
          CALL MP_OXXXXX(P(0,3),MDL_MT,NHEL(3),+1*IC(3),W(1,3))
199
226
          CALL MP_FFV1P0_3(W(1,4),W(1,3),GC_5,ZERO,ZERO,W(1,8))
200
227
C         Counter-term amplitude(s) for loop diagram number 4
201
228
          CALL MP_R2_GG_1_0(W(1,5),W(1,8),R2_GGQ,AMPL(1,1))
 
229
          CALL MP_R2_GG_1_0(W(1,5),W(1,8),R2_GGQ,AMPL(1,2))
 
230
          CALL MP_R2_GG_1_0(W(1,5),W(1,8),R2_GGQ,AMPL(1,3))
 
231
          CALL MP_R2_GG_1_0(W(1,5),W(1,8),R2_GGQ,AMPL(1,4))
202
232
C         Counter-term amplitude(s) for loop diagram number 5
203
 
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),UV_3GB_1EPS,AMPL(2,2))
204
 
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),UV_3GB_1EPS,AMPL(2,3))
205
 
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),UV_3GB_1EPS,AMPL(2,4))
206
233
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),UV_3GB_1EPS,AMPL(2,5))
207
 
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),UV_3GB,AMPL(1,6))
 
234
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),UV_3GB_1EPS,AMPL(2,6))
208
235
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),UV_3GB_1EPS,AMPL(2,7))
209
 
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),UV_3GT,AMPL(1,8))
210
 
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),UV_3GB_1EPS,AMPL(2,9))
211
 
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),UV_3GG_1EPS,AMPL(2,10))
212
 
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),R2_3GQ,AMPL(1,11))
213
 
C         Counter-term amplitude(s) for loop diagram number 7
214
 
          CALL MP_R2_GG_1_0(W(1,5),W(1,8),R2_GGQ,AMPL(1,12))
215
 
C         Counter-term amplitude(s) for loop diagram number 8
216
 
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),R2_3GQ,AMPL(1,13))
217
 
C         Counter-term amplitude(s) for loop diagram number 10
218
 
          CALL MP_R2_GG_1_0(W(1,5),W(1,8),R2_GGQ,AMPL(1,14))
219
 
C         Counter-term amplitude(s) for loop diagram number 11
 
236
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),UV_3GB_1EPS,AMPL(2,8))
 
237
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),UV_3GB,AMPL(1,9))
 
238
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),UV_3GB_1EPS,AMPL(2,10))
 
239
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),UV_3GT,AMPL(1,11))
 
240
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),UV_3GB_1EPS,AMPL(2,12))
 
241
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),UV_3GG_1EPS,AMPL(2,13))
 
242
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),R2_3GQ,AMPL(1,14))
220
243
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),R2_3GQ,AMPL(1,15))
221
 
C         Counter-term amplitude(s) for loop diagram number 13
222
 
          CALL MP_R2_GG_1_0(W(1,5),W(1,8),R2_GGQ,AMPL(1,16))
223
 
C         Counter-term amplitude(s) for loop diagram number 14
 
244
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),R2_3GQ,AMPL(1,16))
224
245
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),R2_3GQ,AMPL(1,17))
225
 
C         Counter-term amplitude(s) for loop diagram number 16
 
246
C         Counter-term amplitude(s) for loop diagram number 7
226
247
          CALL MP_R2_GG_1_R2_GG_3_0(W(1,5),W(1,8),R2_GGQ,R2_GGB,AMPL(1
227
248
     $     ,18))
228
 
C         Counter-term amplitude(s) for loop diagram number 17
 
249
C         Counter-term amplitude(s) for loop diagram number 8
229
250
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),R2_3GQ,AMPL(1,19))
230
 
C         Counter-term amplitude(s) for loop diagram number 19
 
251
C         Counter-term amplitude(s) for loop diagram number 10
231
252
          CALL MP_R2_GG_1_R2_GG_3_0(W(1,5),W(1,8),R2_GGQ,R2_GGT,AMPL(1
232
253
     $     ,20))
233
 
C         Counter-term amplitude(s) for loop diagram number 20
 
254
C         Counter-term amplitude(s) for loop diagram number 11
234
255
          CALL MP_FFV1_0(W(1,4),W(1,3),W(1,5),UV_GQQQ_1EPS,AMPL(2,21))
235
256
          CALL MP_FFV1_0(W(1,4),W(1,3),W(1,5),UV_GQQQ_1EPS,AMPL(2,22))
236
257
          CALL MP_FFV1_0(W(1,4),W(1,3),W(1,5),UV_GQQQ_1EPS,AMPL(2,23))
242
263
          CALL MP_FFV1_0(W(1,4),W(1,3),W(1,5),UV_GQQG_1EPS,AMPL(2,29))
243
264
          CALL MP_FFV1_0(W(1,4),W(1,3),W(1,5),R2_GQQ,AMPL(1,30))
244
265
          CALL MP_FFV1_2(W(1,4),W(1,2),GC_5,MDL_MT,MDL_WT,W(1,9))
245
 
C         Counter-term amplitude(s) for loop diagram number 22
 
266
C         Counter-term amplitude(s) for loop diagram number 13
246
267
          CALL MP_R2_QQ_1_R2_QQ_2_0(W(1,9),W(1,6),R2_QQQ,R2_QQT,AMPL(1
247
268
     $     ,31))
248
269
          CALL MP_R2_QQ_2_0(W(1,9),W(1,6),UV_TMASS,AMPL(1,32))
249
270
          CALL MP_R2_QQ_2_0(W(1,9),W(1,6),UV_TMASS_1EPS,AMPL(2,33))
250
 
C         Counter-term amplitude(s) for loop diagram number 23
 
271
C         Counter-term amplitude(s) for loop diagram number 14
251
272
          CALL MP_FFV1_0(W(1,4),W(1,6),W(1,2),UV_GQQQ_1EPS,AMPL(2,34))
252
273
          CALL MP_FFV1_0(W(1,4),W(1,6),W(1,2),UV_GQQQ_1EPS,AMPL(2,35))
253
274
          CALL MP_FFV1_0(W(1,4),W(1,6),W(1,2),UV_GQQQ_1EPS,AMPL(2,36))
259
280
          CALL MP_FFV1_0(W(1,4),W(1,6),W(1,2),UV_GQQG_1EPS,AMPL(2,42))
260
281
          CALL MP_FFV1_0(W(1,4),W(1,6),W(1,2),R2_GQQ,AMPL(1,43))
261
282
          CALL MP_FFV1_1(W(1,3),W(1,2),GC_5,MDL_MT,MDL_WT,W(1,10))
262
 
C         Counter-term amplitude(s) for loop diagram number 25
 
283
C         Counter-term amplitude(s) for loop diagram number 16
263
284
          CALL MP_R2_QQ_1_R2_QQ_2_0(W(1,7),W(1,10),R2_QQQ,R2_QQT
264
285
     $     ,AMPL(1,44))
265
286
          CALL MP_R2_QQ_2_0(W(1,7),W(1,10),UV_TMASS,AMPL(1,45))
266
287
          CALL MP_R2_QQ_2_0(W(1,7),W(1,10),UV_TMASS_1EPS,AMPL(2,46))
267
 
C         Counter-term amplitude(s) for loop diagram number 26
 
288
C         Counter-term amplitude(s) for loop diagram number 17
268
289
          CALL MP_FFV1_0(W(1,7),W(1,3),W(1,2),UV_GQQQ_1EPS,AMPL(2,47))
269
290
          CALL MP_FFV1_0(W(1,7),W(1,3),W(1,2),UV_GQQQ_1EPS,AMPL(2,48))
270
291
          CALL MP_FFV1_0(W(1,7),W(1,3),W(1,2),UV_GQQQ_1EPS,AMPL(2,49))
275
296
          CALL MP_FFV1_0(W(1,7),W(1,3),W(1,2),UV_GQQQ_1EPS,AMPL(2,54))
276
297
          CALL MP_FFV1_0(W(1,7),W(1,3),W(1,2),UV_GQQG_1EPS,AMPL(2,55))
277
298
          CALL MP_FFV1_0(W(1,7),W(1,3),W(1,2),R2_GQQ,AMPL(1,56))
278
 
C         Counter-term amplitude(s) for loop diagram number 28
 
299
C         Counter-term amplitude(s) for loop diagram number 19
279
300
          CALL MP_FFV1_0(W(1,4),W(1,10),W(1,1),UV_GQQQ_1EPS,AMPL(2,57))
280
301
          CALL MP_FFV1_0(W(1,4),W(1,10),W(1,1),UV_GQQQ_1EPS,AMPL(2,58))
281
302
          CALL MP_FFV1_0(W(1,4),W(1,10),W(1,1),UV_GQQQ_1EPS,AMPL(2,59))
286
307
          CALL MP_FFV1_0(W(1,4),W(1,10),W(1,1),UV_GQQQ_1EPS,AMPL(2,64))
287
308
          CALL MP_FFV1_0(W(1,4),W(1,10),W(1,1),UV_GQQG_1EPS,AMPL(2,65))
288
309
          CALL MP_FFV1_0(W(1,4),W(1,10),W(1,1),R2_GQQ,AMPL(1,66))
289
 
C         Counter-term amplitude(s) for loop diagram number 29
 
310
C         Counter-term amplitude(s) for loop diagram number 20
290
311
          CALL MP_FFV1_0(W(1,9),W(1,3),W(1,1),UV_GQQQ_1EPS,AMPL(2,67))
291
312
          CALL MP_FFV1_0(W(1,9),W(1,3),W(1,1),UV_GQQQ_1EPS,AMPL(2,68))
292
313
          CALL MP_FFV1_0(W(1,9),W(1,3),W(1,1),UV_GQQQ_1EPS,AMPL(2,69))
297
318
          CALL MP_FFV1_0(W(1,9),W(1,3),W(1,1),UV_GQQQ_1EPS,AMPL(2,74))
298
319
          CALL MP_FFV1_0(W(1,9),W(1,3),W(1,1),UV_GQQG_1EPS,AMPL(2,75))
299
320
          CALL MP_FFV1_0(W(1,9),W(1,3),W(1,1),R2_GQQ,AMPL(1,76))
300
 
C         Counter-term amplitude(s) for loop diagram number 31
 
321
C         Counter-term amplitude(s) for loop diagram number 22
301
322
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),R2_3GQ,AMPL(1,77))
302
 
C         Counter-term amplitude(s) for loop diagram number 41
 
323
C         Counter-term amplitude(s) for loop diagram number 32
303
324
          CALL MP_R2_GG_1_R2_GG_2_0(W(1,5),W(1,8),R2_GGG_1,R2_GGG_2
304
325
     $     ,AMPL(1,78))
305
 
C         Counter-term amplitude(s) for loop diagram number 42
 
326
C         Counter-term amplitude(s) for loop diagram number 33
306
327
          CALL MP_VVV1_0(W(1,1),W(1,2),W(1,8),R2_3GG,AMPL(1,79))
307
328
 2000     CONTINUE
308
329
          MP_CT_REQ_SO_DONE=.TRUE.
309
330
 
310
 
C         Amplitude(s) for UVCT diagram with ID 49
 
331
C         Amplitude(s) for UVCT diagram with ID 40
311
332
          CALL MP_FFV1_0(W(1,4),W(1,3),W(1,5),GC_5,AMPL(1,80))
312
333
          AMPL(1,80)=AMPL(1,80)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
313
334
     $     +2.0D0*UVWFCT_T_0)
314
 
C         Amplitude(s) for UVCT diagram with ID 50
 
335
C         Amplitude(s) for UVCT diagram with ID 41
315
336
          CALL MP_FFV1_0(W(1,4),W(1,3),W(1,5),GC_5,AMPL(2,81))
316
337
          AMPL(2,81)=AMPL(2,81)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0*UVWFCT_G_2
317
338
     $     _1EPS)
318
 
C         Amplitude(s) for UVCT diagram with ID 51
 
339
C         Amplitude(s) for UVCT diagram with ID 42
319
340
          CALL MP_FFV1_0(W(1,4),W(1,6),W(1,2),GC_5,AMPL(1,82))
320
341
          AMPL(1,82)=AMPL(1,82)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
321
342
     $     +2.0D0*UVWFCT_T_0)
322
 
C         Amplitude(s) for UVCT diagram with ID 52
 
343
C         Amplitude(s) for UVCT diagram with ID 43
323
344
          CALL MP_FFV1_0(W(1,4),W(1,6),W(1,2),GC_5,AMPL(2,83))
324
345
          AMPL(2,83)=AMPL(2,83)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0*UVWFCT_G_2
325
346
     $     _1EPS)
326
 
C         Amplitude(s) for UVCT diagram with ID 53
 
347
C         Amplitude(s) for UVCT diagram with ID 44
327
348
          CALL MP_FFV1_0(W(1,7),W(1,3),W(1,2),GC_5,AMPL(1,84))
328
349
          AMPL(1,84)=AMPL(1,84)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
329
350
     $     +2.0D0*UVWFCT_T_0)
330
 
C         Amplitude(s) for UVCT diagram with ID 54
 
351
C         Amplitude(s) for UVCT diagram with ID 45
331
352
          CALL MP_FFV1_0(W(1,7),W(1,3),W(1,2),GC_5,AMPL(2,85))
332
353
          AMPL(2,85)=AMPL(2,85)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0*UVWFCT_G_2
333
354
     $     _1EPS)
334
355
 3000     CONTINUE
335
356
          MP_UVCT_REQ_SO_DONE=.TRUE.
336
357
 
337
 
          IF (.NOT.CHECKPHASE.AND.HELDOUBLECHECKED.AND.HELPICKED.EQ.
338
 
     $     -1) THEN
339
 
            DUMMY=GOODHEL(H)
340
 
          ELSE
341
 
            DUMMY=1
342
 
          ENDIF
343
358
          DO I=1,NCTAMPS
344
359
            DO J=1,NBORNAMPS
345
360
              CFTOT=CMPLX(CF_N(I,J)/REAL(ABS(CF_D(I,J)),KIND=16)
349
364
     $         I),ML5_0_ML5SOINDEX_FOR_BORN_AMP(J))
350
365
              IF (.NOT.FILTER_SO.OR.SQSO_TARGET.EQ.ITEMP) THEN
351
366
                DO K=1,3
352
 
                  TEMP2 = DUMMY*2.0E0_16*REAL(CFTOT*AMPL(K,I)
 
367
                  TEMP2 = HEL_MULT*2.0E0_16*REAL(CFTOT*AMPL(K,I)
353
368
     $             *CONJG(AMP(J)),KIND=16)
354
369
                  ANS(K,ITEMP)=ANS(K,ITEMP)+TEMP2
355
370
                  ANS(K,0)=ANS(K,0)+TEMP2
365
380
          CALL MP_FFV1L2_1(PL(0,1),W(1,8),GC_5,ZERO,ZERO,PL(0,2),COEFS)
366
381
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,1),4,COEFS,4,4,WL(1,0,1
367
382
     $     ,2))
368
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,2),2,4,1,1,86,H)
 
383
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,2),2,4,1,1,4,86,H)
369
384
C         Coefficient construction for loop diagram with ID 5
370
385
          CALL MP_FFV1L1_2(PL(0,0),W(1,1),GC_5,ZERO,ZERO,PL(0,3),COEFS)
371
386
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
376
391
          CALL MP_FFV1L1_2(PL(0,4),W(1,8),GC_5,ZERO,ZERO,PL(0,5),COEFS)
377
392
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,4),4,COEFS,4,4,WL(1,0,1
378
393
     $     ,5))
379
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,5),3,4,2,1,87,H)
 
394
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,5),3,4,2,1,4,87,H)
380
395
C         Coefficient construction for loop diagram with ID 6
381
396
          CALL MP_FFV1L2_1(PL(0,0),W(1,1),GC_5,ZERO,ZERO,PL(0,6),COEFS)
382
397
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
387
402
          CALL MP_FFV1L2_1(PL(0,7),W(1,8),GC_5,ZERO,ZERO,PL(0,8),COEFS)
388
403
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,7),4,COEFS,4,4,WL(1,0,1
389
404
     $     ,8))
390
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,8),3,4,2,1,88,H)
 
405
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,8),3,4,2,1,4,88,H)
391
406
C         Coefficient construction for loop diagram with ID 7
392
 
          CALL MP_FFV1L2_1(PL(0,0),W(1,5),GC_5,ZERO,ZERO,PL(0,9),COEFS)
 
407
          CALL MP_FFV1L2_1(PL(0,0),W(1,5),GC_5,MDL_MB,ZERO,PL(0,9)
 
408
     $     ,COEFS)
393
409
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
394
410
     $     ,9))
395
 
          CALL MP_FFV1L2_1(PL(0,9),W(1,8),GC_5,ZERO,ZERO,PL(0,10)
 
411
          CALL MP_FFV1L2_1(PL(0,9),W(1,8),GC_5,MDL_MB,ZERO,PL(0,10)
396
412
     $     ,COEFS)
397
413
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,9),4,COEFS,4,4,WL(1,0,1
398
414
     $     ,10))
399
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,10),2,4,1,1,89,H)
 
415
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,10),2,4,3,1,1,89,H)
400
416
C         Coefficient construction for loop diagram with ID 8
401
 
          CALL MP_FFV1L1_2(PL(0,0),W(1,1),GC_5,ZERO,ZERO,PL(0,11)
 
417
          CALL MP_FFV1L1_2(PL(0,0),W(1,1),GC_5,MDL_MB,ZERO,PL(0,11)
402
418
     $     ,COEFS)
403
419
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
404
420
     $     ,11))
405
 
          CALL MP_FFV1L1_2(PL(0,11),W(1,2),GC_5,ZERO,ZERO,PL(0,12)
 
421
          CALL MP_FFV1L1_2(PL(0,11),W(1,2),GC_5,MDL_MB,ZERO,PL(0,12)
406
422
     $     ,COEFS)
407
423
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,11),4,COEFS,4,4,WL(1,0
408
424
     $     ,1,12))
409
 
          CALL MP_FFV1L1_2(PL(0,12),W(1,8),GC_5,ZERO,ZERO,PL(0,13)
 
425
          CALL MP_FFV1L1_2(PL(0,12),W(1,8),GC_5,MDL_MB,ZERO,PL(0,13)
410
426
     $     ,COEFS)
411
427
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,12),4,COEFS,4,4,WL(1,0
412
428
     $     ,1,13))
413
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,13),3,4,2,1,90,H)
 
429
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,13),3,4,4,1,1,90,H)
414
430
C         Coefficient construction for loop diagram with ID 9
415
 
          CALL MP_FFV1L2_1(PL(0,0),W(1,1),GC_5,ZERO,ZERO,PL(0,14)
 
431
          CALL MP_FFV1L2_1(PL(0,0),W(1,1),GC_5,MDL_MB,ZERO,PL(0,14)
416
432
     $     ,COEFS)
417
433
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
418
434
     $     ,14))
419
 
          CALL MP_FFV1L2_1(PL(0,14),W(1,2),GC_5,ZERO,ZERO,PL(0,15)
 
435
          CALL MP_FFV1L2_1(PL(0,14),W(1,2),GC_5,MDL_MB,ZERO,PL(0,15)
420
436
     $     ,COEFS)
421
437
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,14),4,COEFS,4,4,WL(1,0
422
438
     $     ,1,15))
423
 
          CALL MP_FFV1L2_1(PL(0,15),W(1,8),GC_5,ZERO,ZERO,PL(0,16)
 
439
          CALL MP_FFV1L2_1(PL(0,15),W(1,8),GC_5,MDL_MB,ZERO,PL(0,16)
424
440
     $     ,COEFS)
425
441
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,15),4,COEFS,4,4,WL(1,0
426
442
     $     ,1,16))
427
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,16),3,4,2,1,91,H)
 
443
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,16),3,4,4,1,1,91,H)
428
444
C         Coefficient construction for loop diagram with ID 10
429
 
          CALL MP_FFV1L2_1(PL(0,0),W(1,5),GC_5,ZERO,ZERO,PL(0,17)
 
445
          CALL MP_FFV1L2_1(PL(0,0),W(1,5),GC_5,MDL_MT,MDL_WT,PL(0,17)
430
446
     $     ,COEFS)
431
447
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
432
448
     $     ,17))
433
 
          CALL MP_FFV1L2_1(PL(0,17),W(1,8),GC_5,ZERO,ZERO,PL(0,18)
 
449
          CALL MP_FFV1L2_1(PL(0,17),W(1,8),GC_5,MDL_MT,MDL_WT,PL(0,18)
434
450
     $     ,COEFS)
435
451
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,17),4,COEFS,4,4,WL(1,0
436
452
     $     ,1,18))
437
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,18),2,4,1,1,92,H)
 
453
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,18),2,4,5,1,1,92,H)
438
454
C         Coefficient construction for loop diagram with ID 11
439
 
          CALL MP_FFV1L1_2(PL(0,0),W(1,1),GC_5,ZERO,ZERO,PL(0,19)
 
455
          CALL MP_FFV1L1P0_3(PL(0,0),W(1,3),GC_5,ZERO,ZERO,PL(0,19)
440
456
     $     ,COEFS)
441
 
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
 
457
          CALL MP_ML5_0_UPDATE_WL_0_0(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
442
458
     $     ,19))
443
 
          CALL MP_FFV1L1_2(PL(0,19),W(1,2),GC_5,ZERO,ZERO,PL(0,20)
 
459
          CALL MP_FFV1L3_2(PL(0,19),W(1,4),GC_5,MDL_MT,MDL_WT,PL(0,20)
444
460
     $     ,COEFS)
445
 
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,19),4,COEFS,4,4,WL(1,0
 
461
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,19),4,COEFS,4,4,WL(1,0
446
462
     $     ,1,20))
447
 
          CALL MP_FFV1L1_2(PL(0,20),W(1,8),GC_5,ZERO,ZERO,PL(0,21)
 
463
          CALL MP_FFV1L1_2(PL(0,20),W(1,5),GC_5,MDL_MT,MDL_WT,PL(0,21)
448
464
     $     ,COEFS)
449
 
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,20),4,COEFS,4,4,WL(1,0
 
465
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,20),4,COEFS,4,4,WL(1,0
450
466
     $     ,1,21))
451
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,21),3,4,2,1,93,H)
 
467
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,21),2,4,6,1,1,93,H)
452
468
C         Coefficient construction for loop diagram with ID 12
453
 
          CALL MP_FFV1L2_1(PL(0,0),W(1,1),GC_5,ZERO,ZERO,PL(0,22)
 
469
          CALL MP_FFV1L3_1(PL(0,0),W(1,3),GC_5,MDL_MT,MDL_WT,PL(0,22)
454
470
     $     ,COEFS)
455
471
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
456
472
     $     ,22))
457
 
          CALL MP_FFV1L2_1(PL(0,22),W(1,2),GC_5,ZERO,ZERO,PL(0,23)
 
473
          CALL MP_FFV1L2P0_3(PL(0,22),W(1,4),GC_5,ZERO,ZERO,PL(0,23)
458
474
     $     ,COEFS)
459
 
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,22),4,COEFS,4,4,WL(1,0
 
475
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,22),4,COEFS,4,4,WL(1,0
460
476
     $     ,1,23))
461
 
          CALL MP_FFV1L2_1(PL(0,23),W(1,8),GC_5,ZERO,ZERO,PL(0,24)
 
477
          CALL MP_VVV1L2P0_1(PL(0,23),W(1,5),GC_4,ZERO,ZERO,PL(0,24)
462
478
     $     ,COEFS)
463
 
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,23),4,COEFS,4,4,WL(1,0
 
479
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,23),4,COEFS,4,4,WL(1,0
464
480
     $     ,1,24))
465
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,24),3,4,2,1,94,H)
 
481
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,24),2,4,7,1,1,94,H)
466
482
C         Coefficient construction for loop diagram with ID 13
467
 
          CALL MP_FFV1L2_1(PL(0,0),W(1,5),GC_5,ZERO,ZERO,PL(0,25)
 
483
          CALL MP_FFV1L1P0_3(PL(0,0),W(1,6),GC_5,ZERO,ZERO,PL(0,25)
468
484
     $     ,COEFS)
469
 
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
 
485
          CALL MP_ML5_0_UPDATE_WL_0_0(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
470
486
     $     ,25))
471
 
          CALL MP_FFV1L2_1(PL(0,25),W(1,8),GC_5,ZERO,ZERO,PL(0,26)
 
487
          CALL MP_FFV1L3_2(PL(0,25),W(1,9),GC_5,MDL_MT,MDL_WT,PL(0,26)
472
488
     $     ,COEFS)
473
 
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,25),4,COEFS,4,4,WL(1,0
 
489
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,25),4,COEFS,4,4,WL(1,0
474
490
     $     ,1,26))
475
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,26),2,4,1,1,95,H)
 
491
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,26),1,4,8,1,1,95,H)
476
492
C         Coefficient construction for loop diagram with ID 14
477
 
          CALL MP_FFV1L1_2(PL(0,0),W(1,1),GC_5,ZERO,ZERO,PL(0,27)
 
493
          CALL MP_FFV1L2_1(PL(0,0),W(1,2),GC_5,MDL_MT,MDL_WT,PL(0,27)
478
494
     $     ,COEFS)
479
495
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
480
496
     $     ,27))
481
 
          CALL MP_FFV1L1_2(PL(0,27),W(1,2),GC_5,ZERO,ZERO,PL(0,28)
 
497
          CALL MP_FFV1L2P0_3(PL(0,27),W(1,4),GC_5,ZERO,ZERO,PL(0,28)
482
498
     $     ,COEFS)
483
 
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,27),4,COEFS,4,4,WL(1,0
 
499
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,27),4,COEFS,4,4,WL(1,0
484
500
     $     ,1,28))
485
 
          CALL MP_FFV1L1_2(PL(0,28),W(1,8),GC_5,ZERO,ZERO,PL(0,29)
 
501
          CALL MP_FFV1L3_1(PL(0,28),W(1,6),GC_5,MDL_MT,MDL_WT,PL(0,29)
486
502
     $     ,COEFS)
487
 
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,28),4,COEFS,4,4,WL(1,0
 
503
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,28),4,COEFS,4,4,WL(1,0
488
504
     $     ,1,29))
489
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,29),3,4,2,1,96,H)
 
505
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,29),2,4,9,1,1,96,H)
490
506
C         Coefficient construction for loop diagram with ID 15
491
 
          CALL MP_FFV1L2_1(PL(0,0),W(1,1),GC_5,ZERO,ZERO,PL(0,30)
 
507
          CALL MP_VVV1L2P0_1(PL(0,0),W(1,2),GC_4,ZERO,ZERO,PL(0,30)
492
508
     $     ,COEFS)
493
509
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
494
510
     $     ,30))
495
 
          CALL MP_FFV1L2_1(PL(0,30),W(1,2),GC_5,ZERO,ZERO,PL(0,31)
 
511
          CALL MP_FFV1L3_2(PL(0,30),W(1,4),GC_5,MDL_MT,MDL_WT,PL(0,31)
496
512
     $     ,COEFS)
497
513
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,30),4,COEFS,4,4,WL(1,0
498
514
     $     ,1,31))
499
 
          CALL MP_FFV1L2_1(PL(0,31),W(1,8),GC_5,ZERO,ZERO,PL(0,32)
 
515
          CALL MP_FFV1L1P0_3(PL(0,31),W(1,6),GC_5,ZERO,ZERO,PL(0,32)
500
516
     $     ,COEFS)
501
 
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,31),4,COEFS,4,4,WL(1,0
 
517
          CALL MP_ML5_0_UPDATE_WL_2_0(WL(1,0,1,31),4,COEFS,4,4,WL(1,0
502
518
     $     ,1,32))
503
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,32),3,4,2,1,97,H)
 
519
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,32),2,4,10,1,1,97,H)
504
520
C         Coefficient construction for loop diagram with ID 16
505
 
          CALL MP_FFV1L2_1(PL(0,0),W(1,5),GC_5,MDL_MB,ZERO,PL(0,33)
 
521
          CALL MP_FFV1L1P0_3(PL(0,0),W(1,10),GC_5,ZERO,ZERO,PL(0,33)
506
522
     $     ,COEFS)
507
 
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
 
523
          CALL MP_ML5_0_UPDATE_WL_0_0(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
508
524
     $     ,33))
509
 
          CALL MP_FFV1L2_1(PL(0,33),W(1,8),GC_5,MDL_MB,ZERO,PL(0,34)
 
525
          CALL MP_FFV1L3_2(PL(0,33),W(1,7),GC_5,MDL_MT,MDL_WT,PL(0,34)
510
526
     $     ,COEFS)
511
 
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,33),4,COEFS,4,4,WL(1,0
 
527
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,33),4,COEFS,4,4,WL(1,0
512
528
     $     ,1,34))
513
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,34),2,4,3,1,98,H)
 
529
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,34),1,4,11,1,1,98,H)
514
530
C         Coefficient construction for loop diagram with ID 17
515
 
          CALL MP_FFV1L1_2(PL(0,0),W(1,1),GC_5,MDL_MB,ZERO,PL(0,35)
 
531
          CALL MP_FFV1L1_2(PL(0,0),W(1,2),GC_5,MDL_MT,MDL_WT,PL(0,35)
516
532
     $     ,COEFS)
517
533
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
518
534
     $     ,35))
519
 
          CALL MP_FFV1L1_2(PL(0,35),W(1,2),GC_5,MDL_MB,ZERO,PL(0,36)
 
535
          CALL MP_FFV1L1P0_3(PL(0,35),W(1,3),GC_5,ZERO,ZERO,PL(0,36)
520
536
     $     ,COEFS)
521
 
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,35),4,COEFS,4,4,WL(1,0
 
537
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,35),4,COEFS,4,4,WL(1,0
522
538
     $     ,1,36))
523
 
          CALL MP_FFV1L1_2(PL(0,36),W(1,8),GC_5,MDL_MB,ZERO,PL(0,37)
 
539
          CALL MP_FFV1L3_2(PL(0,36),W(1,7),GC_5,MDL_MT,MDL_WT,PL(0,37)
524
540
     $     ,COEFS)
525
 
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,36),4,COEFS,4,4,WL(1,0
 
541
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,36),4,COEFS,4,4,WL(1,0
526
542
     $     ,1,37))
527
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,37),3,4,4,1,99,H)
 
543
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,37),2,4,12,1,1,99,H)
528
544
C         Coefficient construction for loop diagram with ID 18
529
 
          CALL MP_FFV1L2_1(PL(0,0),W(1,1),GC_5,MDL_MB,ZERO,PL(0,38)
530
 
     $     ,COEFS)
531
 
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
532
 
     $     ,38))
533
 
          CALL MP_FFV1L2_1(PL(0,38),W(1,2),GC_5,MDL_MB,ZERO,PL(0,39)
534
 
     $     ,COEFS)
535
 
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,38),4,COEFS,4,4,WL(1,0
 
545
          CALL MP_FFV1L3_1(PL(0,30),W(1,3),GC_5,MDL_MT,MDL_WT,PL(0,38)
 
546
     $     ,COEFS)
 
547
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,30),4,COEFS,4,4,WL(1,0
 
548
     $     ,1,38))
 
549
          CALL MP_FFV1L2P0_3(PL(0,38),W(1,7),GC_5,ZERO,ZERO,PL(0,39)
 
550
     $     ,COEFS)
 
551
          CALL MP_ML5_0_UPDATE_WL_2_0(WL(1,0,1,38),4,COEFS,4,4,WL(1,0
536
552
     $     ,1,39))
537
 
          CALL MP_FFV1L2_1(PL(0,39),W(1,8),GC_5,MDL_MB,ZERO,PL(0,40)
538
 
     $     ,COEFS)
539
 
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,39),4,COEFS,4,4,WL(1,0
540
 
     $     ,1,40))
541
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,40),3,4,4,1,100,H)
 
553
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,39),2,4,13,1,1,100
 
554
     $     ,H)
542
555
C         Coefficient construction for loop diagram with ID 19
543
 
          CALL MP_FFV1L2_1(PL(0,0),W(1,5),GC_5,MDL_MT,MDL_WT,PL(0,41)
 
556
          CALL MP_FFV1L2_1(PL(0,0),W(1,1),GC_5,MDL_MT,MDL_WT,PL(0,40)
544
557
     $     ,COEFS)
545
558
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
546
 
     $     ,41))
547
 
          CALL MP_FFV1L2_1(PL(0,41),W(1,8),GC_5,MDL_MT,MDL_WT,PL(0,42)
 
559
     $     ,40))
 
560
          CALL MP_FFV1L2P0_3(PL(0,40),W(1,4),GC_5,ZERO,ZERO,PL(0,41)
548
561
     $     ,COEFS)
 
562
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,40),4,COEFS,4,4,WL(1,0
 
563
     $     ,1,41))
 
564
          CALL MP_FFV1L3_1(PL(0,41),W(1,10),GC_5,MDL_MT,MDL_WT,PL(0
 
565
     $     ,42),COEFS)
549
566
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,41),4,COEFS,4,4,WL(1,0
550
567
     $     ,1,42))
551
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,42),2,4,5,1,101,H)
 
568
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,42),2,4,14,1,1,101
 
569
     $     ,H)
552
570
C         Coefficient construction for loop diagram with ID 20
553
 
          CALL MP_FFV1L1P0_3(PL(0,0),W(1,3),GC_5,ZERO,ZERO,PL(0,43)
 
571
          CALL MP_FFV1L1_2(PL(0,0),W(1,1),GC_5,MDL_MT,MDL_WT,PL(0,43)
554
572
     $     ,COEFS)
555
 
          CALL MP_ML5_0_UPDATE_WL_0_0(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
 
573
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
556
574
     $     ,43))
557
 
          CALL MP_FFV1L3_2(PL(0,43),W(1,4),GC_5,MDL_MT,MDL_WT,PL(0,44)
 
575
          CALL MP_FFV1L1P0_3(PL(0,43),W(1,3),GC_5,ZERO,ZERO,PL(0,44)
558
576
     $     ,COEFS)
559
 
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,43),4,COEFS,4,4,WL(1,0
 
577
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,43),4,COEFS,4,4,WL(1,0
560
578
     $     ,1,44))
561
 
          CALL MP_FFV1L1_2(PL(0,44),W(1,5),GC_5,MDL_MT,MDL_WT,PL(0,45)
 
579
          CALL MP_FFV1L3_2(PL(0,44),W(1,9),GC_5,MDL_MT,MDL_WT,PL(0,45)
562
580
     $     ,COEFS)
563
581
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,44),4,COEFS,4,4,WL(1,0
564
582
     $     ,1,45))
565
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,45),2,4,6,1,102,H)
 
583
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,45),2,4,15,1,1,102
 
584
     $     ,H)
566
585
C         Coefficient construction for loop diagram with ID 21
567
 
          CALL MP_FFV1L3_1(PL(0,0),W(1,3),GC_5,MDL_MT,MDL_WT,PL(0,46)
568
 
     $     ,COEFS)
569
 
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
570
 
     $     ,46))
571
 
          CALL MP_FFV1L2P0_3(PL(0,46),W(1,4),GC_5,ZERO,ZERO,PL(0,47)
572
 
     $     ,COEFS)
573
 
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,46),4,COEFS,4,4,WL(1,0
 
586
          CALL MP_FFV1L1_2(PL(0,43),W(1,2),GC_5,MDL_MT,MDL_WT,PL(0,46)
 
587
     $     ,COEFS)
 
588
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,43),4,COEFS,4,4,WL(1,0
 
589
     $     ,1,46))
 
590
          CALL MP_FFV1L1P0_3(PL(0,46),W(1,3),GC_5,ZERO,ZERO,PL(0,47)
 
591
     $     ,COEFS)
 
592
          CALL MP_ML5_0_UPDATE_WL_2_0(WL(1,0,1,46),4,COEFS,4,4,WL(1,0
574
593
     $     ,1,47))
575
 
          CALL MP_VVV1L2P0_1(PL(0,47),W(1,5),GC_4,ZERO,ZERO,PL(0,48)
 
594
          CALL MP_FFV1L3_2(PL(0,47),W(1,4),GC_5,MDL_MT,MDL_WT,PL(0,48)
576
595
     $     ,COEFS)
577
 
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,47),4,COEFS,4,4,WL(1,0
 
596
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,47),4,COEFS,4,4,WL(1,0
578
597
     $     ,1,48))
579
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,48),2,4,7,1,103,H)
 
598
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,48),3,4,16,1,1,103
 
599
     $     ,H)
580
600
C         Coefficient construction for loop diagram with ID 22
581
 
          CALL MP_FFV1L1P0_3(PL(0,0),W(1,6),GC_5,ZERO,ZERO,PL(0,49)
582
 
     $     ,COEFS)
583
 
          CALL MP_ML5_0_UPDATE_WL_0_0(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
584
 
     $     ,49))
585
 
          CALL MP_FFV1L3_2(PL(0,49),W(1,9),GC_5,MDL_MT,MDL_WT,PL(0,50)
586
 
     $     ,COEFS)
587
 
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,49),4,COEFS,4,4,WL(1,0
 
601
          CALL MP_FFV1L1_2(PL(0,46),W(1,8),GC_5,MDL_MT,MDL_WT,PL(0,49)
 
602
     $     ,COEFS)
 
603
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,46),4,COEFS,4,4,WL(1,0
 
604
     $     ,1,49))
 
605
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,49),3,4,17,1,1,104
 
606
     $     ,H)
 
607
C         Coefficient construction for loop diagram with ID 23
 
608
          CALL MP_FFV1L2_1(PL(0,40),W(1,2),GC_5,MDL_MT,MDL_WT,PL(0,50)
 
609
     $     ,COEFS)
 
610
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,40),4,COEFS,4,4,WL(1,0
588
611
     $     ,1,50))
589
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,50),1,4,8,1,104,H)
590
 
C         Coefficient construction for loop diagram with ID 23
591
 
          CALL MP_FFV1L2_1(PL(0,0),W(1,2),GC_5,MDL_MT,MDL_WT,PL(0,51)
592
 
     $     ,COEFS)
593
 
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
594
 
     $     ,51))
595
 
          CALL MP_FFV1L2P0_3(PL(0,51),W(1,4),GC_5,ZERO,ZERO,PL(0,52)
596
 
     $     ,COEFS)
597
 
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,51),4,COEFS,4,4,WL(1,0
 
612
          CALL MP_FFV1L2_1(PL(0,50),W(1,8),GC_5,MDL_MT,MDL_WT,PL(0,51)
 
613
     $     ,COEFS)
 
614
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,50),4,COEFS,4,4,WL(1,0
 
615
     $     ,1,51))
 
616
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,51),3,4,17,1,1,105
 
617
     $     ,H)
 
618
C         Coefficient construction for loop diagram with ID 24
 
619
          CALL MP_FFV1L2P0_3(PL(0,50),W(1,4),GC_5,ZERO,ZERO,PL(0,52)
 
620
     $     ,COEFS)
 
621
          CALL MP_ML5_0_UPDATE_WL_2_0(WL(1,0,1,50),4,COEFS,4,4,WL(1,0
598
622
     $     ,1,52))
599
 
          CALL MP_FFV1L3_1(PL(0,52),W(1,6),GC_5,MDL_MT,MDL_WT,PL(0,53)
 
623
          CALL MP_FFV1L3_1(PL(0,52),W(1,3),GC_5,MDL_MT,MDL_WT,PL(0,53)
600
624
     $     ,COEFS)
601
 
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,52),4,COEFS,4,4,WL(1,0
 
625
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,52),4,COEFS,4,4,WL(1,0
602
626
     $     ,1,53))
603
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,53),2,4,9,1,105,H)
604
 
C         Coefficient construction for loop diagram with ID 24
605
 
          CALL MP_VVV1L2P0_1(PL(0,0),W(1,2),GC_4,ZERO,ZERO,PL(0,54)
 
627
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,53),3,4,18,1,1,106
 
628
     $     ,H)
 
629
C         Coefficient construction for loop diagram with ID 25
 
630
          CALL MP_VVV1L2P0_1(PL(0,44),W(1,2),GC_4,ZERO,ZERO,PL(0,54)
606
631
     $     ,COEFS)
607
 
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
608
 
     $     ,54))
 
632
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,44),4,COEFS,4,4,WL(1,0
 
633
     $     ,1,54))
609
634
          CALL MP_FFV1L3_2(PL(0,54),W(1,4),GC_5,MDL_MT,MDL_WT,PL(0,55)
610
635
     $     ,COEFS)
611
 
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,54),4,COEFS,4,4,WL(1,0
 
636
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,54),4,COEFS,4,4,WL(1,0
612
637
     $     ,1,55))
613
 
          CALL MP_FFV1L1P0_3(PL(0,55),W(1,6),GC_5,ZERO,ZERO,PL(0,56)
614
 
     $     ,COEFS)
615
 
          CALL MP_ML5_0_UPDATE_WL_2_0(WL(1,0,1,55),4,COEFS,4,4,WL(1,0
616
 
     $     ,1,56))
617
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,56),2,4,10,1,106,H)
618
 
C         Coefficient construction for loop diagram with ID 25
619
 
          CALL MP_FFV1L1P0_3(PL(0,0),W(1,10),GC_5,ZERO,ZERO,PL(0,57)
620
 
     $     ,COEFS)
621
 
          CALL MP_ML5_0_UPDATE_WL_0_0(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
622
 
     $     ,57))
623
 
          CALL MP_FFV1L3_2(PL(0,57),W(1,7),GC_5,MDL_MT,MDL_WT,PL(0,58)
624
 
     $     ,COEFS)
625
 
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,57),4,COEFS,4,4,WL(1,0
 
638
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,55),3,4,19,1,1,107
 
639
     $     ,H)
 
640
C         Coefficient construction for loop diagram with ID 26
 
641
          CALL MP_VVV1L2P0_1(PL(0,0),W(1,1),GC_4,ZERO,ZERO,PL(0,56)
 
642
     $     ,COEFS)
 
643
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
 
644
     $     ,56))
 
645
          CALL MP_FFV1L3_2(PL(0,56),W(1,4),GC_5,MDL_MT,MDL_WT,PL(0,57)
 
646
     $     ,COEFS)
 
647
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,56),4,COEFS,4,4,WL(1,0
 
648
     $     ,1,57))
 
649
          CALL MP_FFV1L1P0_3(PL(0,57),W(1,10),GC_5,ZERO,ZERO,PL(0,58)
 
650
     $     ,COEFS)
 
651
          CALL MP_ML5_0_UPDATE_WL_2_0(WL(1,0,1,57),4,COEFS,4,4,WL(1,0
626
652
     $     ,1,58))
627
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,58),1,4,11,1,107,H)
628
 
C         Coefficient construction for loop diagram with ID 26
629
 
          CALL MP_FFV1L1_2(PL(0,0),W(1,2),GC_5,MDL_MT,MDL_WT,PL(0,59)
630
 
     $     ,COEFS)
631
 
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
632
 
     $     ,59))
633
 
          CALL MP_FFV1L1P0_3(PL(0,59),W(1,3),GC_5,ZERO,ZERO,PL(0,60)
634
 
     $     ,COEFS)
635
 
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,59),4,COEFS,4,4,WL(1,0
 
653
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,58),2,4,20,1,1,108
 
654
     $     ,H)
 
655
C         Coefficient construction for loop diagram with ID 27
 
656
          CALL MP_FFV1L3_1(PL(0,56),W(1,3),GC_5,MDL_MT,MDL_WT,PL(0,59)
 
657
     $     ,COEFS)
 
658
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,56),4,COEFS,4,4,WL(1,0
 
659
     $     ,1,59))
 
660
          CALL MP_FFV1L2P0_3(PL(0,59),W(1,9),GC_5,ZERO,ZERO,PL(0,60)
 
661
     $     ,COEFS)
 
662
          CALL MP_ML5_0_UPDATE_WL_2_0(WL(1,0,1,59),4,COEFS,4,4,WL(1,0
636
663
     $     ,1,60))
637
 
          CALL MP_FFV1L3_2(PL(0,60),W(1,7),GC_5,MDL_MT,MDL_WT,PL(0,61)
 
664
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,60),2,4,21,1,1,109
 
665
     $     ,H)
 
666
C         Coefficient construction for loop diagram with ID 28
 
667
          CALL MP_FFV1L2_1(PL(0,59),W(1,2),GC_5,MDL_MT,MDL_WT,PL(0,61)
638
668
     $     ,COEFS)
639
 
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,60),4,COEFS,4,4,WL(1,0
 
669
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,59),4,COEFS,4,4,WL(1,0
640
670
     $     ,1,61))
641
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,61),2,4,12,1,108,H)
642
 
C         Coefficient construction for loop diagram with ID 27
643
 
          CALL MP_FFV1L3_1(PL(0,54),W(1,3),GC_5,MDL_MT,MDL_WT,PL(0,62)
 
671
          CALL MP_FFV1L2P0_3(PL(0,61),W(1,4),GC_5,ZERO,ZERO,PL(0,62)
644
672
     $     ,COEFS)
645
 
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,54),4,COEFS,4,4,WL(1,0
 
673
          CALL MP_ML5_0_UPDATE_WL_3_0(WL(1,0,1,61),4,COEFS,4,4,WL(1,0
646
674
     $     ,1,62))
647
 
          CALL MP_FFV1L2P0_3(PL(0,62),W(1,7),GC_5,ZERO,ZERO,PL(0,63)
648
 
     $     ,COEFS)
649
 
          CALL MP_ML5_0_UPDATE_WL_2_0(WL(1,0,1,62),4,COEFS,4,4,WL(1,0
 
675
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,62),3,4,22,1,1,110
 
676
     $     ,H)
 
677
C         Coefficient construction for loop diagram with ID 29
 
678
          CALL MP_VVVV1L2P0_1(PL(0,23),W(1,1),W(1,2),GC_6,ZERO,ZERO
 
679
     $     ,PL(0,63),COEFS)
 
680
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,23),4,COEFS,4,4,WL(1,0
650
681
     $     ,1,63))
651
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,63),2,4,13,1,109,H)
652
 
C         Coefficient construction for loop diagram with ID 28
653
 
          CALL MP_FFV1L2_1(PL(0,0),W(1,1),GC_5,MDL_MT,MDL_WT,PL(0,64)
654
 
     $     ,COEFS)
655
 
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
656
 
     $     ,64))
657
 
          CALL MP_FFV1L2P0_3(PL(0,64),W(1,4),GC_5,ZERO,ZERO,PL(0,65)
658
 
     $     ,COEFS)
659
 
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,64),4,COEFS,4,4,WL(1,0
 
682
          CALL MP_VVVV3L2P0_1(PL(0,23),W(1,1),W(1,2),GC_6,ZERO,ZERO
 
683
     $     ,PL(0,64),COEFS)
 
684
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,23),4,COEFS,4,4,WL(1,0
 
685
     $     ,1,64))
 
686
          CALL MP_VVVV4L2P0_1(PL(0,23),W(1,1),W(1,2),GC_6,ZERO,ZERO
 
687
     $     ,PL(0,65),COEFS)
 
688
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,23),4,COEFS,4,4,WL(1,0
660
689
     $     ,1,65))
661
 
          CALL MP_FFV1L3_1(PL(0,65),W(1,10),GC_5,MDL_MT,MDL_WT,PL(0
662
 
     $     ,66),COEFS)
663
 
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,65),4,COEFS,4,4,WL(1,0
 
690
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,63),1,4,7,1,1,111,H)
 
691
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,64),1,4,7,1,1,112,H)
 
692
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,65),1,4,7,1,1,113,H)
 
693
C         Coefficient construction for loop diagram with ID 30
 
694
          CALL MP_VVV1L2P0_1(PL(0,56),W(1,2),GC_4,ZERO,ZERO,PL(0,66)
 
695
     $     ,COEFS)
 
696
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,56),4,COEFS,4,4,WL(1,0
664
697
     $     ,1,66))
665
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,66),2,4,14,1,110,H)
666
 
C         Coefficient construction for loop diagram with ID 29
667
 
          CALL MP_FFV1L1_2(PL(0,0),W(1,1),GC_5,MDL_MT,MDL_WT,PL(0,67)
 
698
          CALL MP_FFV1L3_2(PL(0,66),W(1,4),GC_5,MDL_MT,MDL_WT,PL(0,67)
668
699
     $     ,COEFS)
669
 
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
670
 
     $     ,67))
 
700
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,66),4,COEFS,4,4,WL(1,0
 
701
     $     ,1,67))
671
702
          CALL MP_FFV1L1P0_3(PL(0,67),W(1,3),GC_5,ZERO,ZERO,PL(0,68)
672
703
     $     ,COEFS)
673
 
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,67),4,COEFS,4,4,WL(1,0
 
704
          CALL MP_ML5_0_UPDATE_WL_3_0(WL(1,0,1,67),4,COEFS,4,4,WL(1,0
674
705
     $     ,1,68))
675
 
          CALL MP_FFV1L3_2(PL(0,68),W(1,9),GC_5,MDL_MT,MDL_WT,PL(0,69)
 
706
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,68),3,4,23,1,1,114
 
707
     $     ,H)
 
708
C         Coefficient construction for loop diagram with ID 31
 
709
          CALL MP_FFV1L3_1(PL(0,66),W(1,3),GC_5,MDL_MT,MDL_WT,PL(0,69)
676
710
     $     ,COEFS)
677
 
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,68),4,COEFS,4,4,WL(1,0
 
711
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,66),4,COEFS,4,4,WL(1,0
678
712
     $     ,1,69))
679
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,69),2,4,15,1,111,H)
680
 
C         Coefficient construction for loop diagram with ID 30
681
 
          CALL MP_FFV1L1_2(PL(0,67),W(1,2),GC_5,MDL_MT,MDL_WT,PL(0,70)
 
713
          CALL MP_FFV1L2P0_3(PL(0,69),W(1,4),GC_5,ZERO,ZERO,PL(0,70)
682
714
     $     ,COEFS)
683
 
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,67),4,COEFS,4,4,WL(1,0
 
715
          CALL MP_ML5_0_UPDATE_WL_3_0(WL(1,0,1,69),4,COEFS,4,4,WL(1,0
684
716
     $     ,1,70))
685
 
          CALL MP_FFV1L1P0_3(PL(0,70),W(1,3),GC_5,ZERO,ZERO,PL(0,71)
686
 
     $     ,COEFS)
687
 
          CALL MP_ML5_0_UPDATE_WL_2_0(WL(1,0,1,70),4,COEFS,4,4,WL(1,0
688
 
     $     ,1,71))
689
 
          CALL MP_FFV1L3_2(PL(0,71),W(1,4),GC_5,MDL_MT,MDL_WT,PL(0,72)
690
 
     $     ,COEFS)
691
 
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,71),4,COEFS,4,4,WL(1,0
 
717
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,70),3,4,24,1,1,115
 
718
     $     ,H)
 
719
C         Coefficient construction for loop diagram with ID 32
 
720
          CALL MP_VVV1L2P0_1(PL(0,0),W(1,5),GC_4,ZERO,ZERO,PL(0,71)
 
721
     $     ,COEFS)
 
722
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
 
723
     $     ,71))
 
724
          CALL MP_VVV1L2P0_1(PL(0,71),W(1,8),GC_4,ZERO,ZERO,PL(0,72)
 
725
     $     ,COEFS)
 
726
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,71),4,COEFS,4,4,WL(1,0
692
727
     $     ,1,72))
693
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,72),3,4,16,1,112,H)
694
 
C         Coefficient construction for loop diagram with ID 31
695
 
          CALL MP_FFV1L1_2(PL(0,70),W(1,8),GC_5,MDL_MT,MDL_WT,PL(0,73)
 
728
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,72),2,4,1,2,1,116,H)
 
729
C         Coefficient construction for loop diagram with ID 33
 
730
          CALL MP_VVV1L2P0_1(PL(0,66),W(1,8),GC_4,ZERO,ZERO,PL(0,73)
696
731
     $     ,COEFS)
697
 
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,70),4,COEFS,4,4,WL(1,0
 
732
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,66),4,COEFS,4,4,WL(1,0
698
733
     $     ,1,73))
699
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,73),3,4,17,1,113,H)
700
 
C         Coefficient construction for loop diagram with ID 32
701
 
          CALL MP_FFV1L2_1(PL(0,64),W(1,2),GC_5,MDL_MT,MDL_WT,PL(0,74)
702
 
     $     ,COEFS)
703
 
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,64),4,COEFS,4,4,WL(1,0
 
734
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,73),3,4,2,1,1,117,H)
 
735
C         Coefficient construction for loop diagram with ID 34
 
736
          CALL MP_VVVV1L2P0_1(PL(0,56),W(1,8),W(1,2),GC_6,ZERO,ZERO
 
737
     $     ,PL(0,74),COEFS)
 
738
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,56),4,COEFS,4,4,WL(1,0
704
739
     $     ,1,74))
705
 
          CALL MP_FFV1L2_1(PL(0,74),W(1,8),GC_5,MDL_MT,MDL_WT,PL(0,75)
706
 
     $     ,COEFS)
707
 
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,74),4,COEFS,4,4,WL(1,0
 
740
          CALL MP_VVVV3L2P0_1(PL(0,56),W(1,8),W(1,2),GC_6,ZERO,ZERO
 
741
     $     ,PL(0,75),COEFS)
 
742
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,56),4,COEFS,4,4,WL(1,0
708
743
     $     ,1,75))
709
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,75),3,4,17,1,114,H)
710
 
C         Coefficient construction for loop diagram with ID 33
711
 
          CALL MP_FFV1L2P0_3(PL(0,74),W(1,4),GC_5,ZERO,ZERO,PL(0,76)
712
 
     $     ,COEFS)
713
 
          CALL MP_ML5_0_UPDATE_WL_2_0(WL(1,0,1,74),4,COEFS,4,4,WL(1,0
 
744
          CALL MP_VVVV4L2P0_1(PL(0,56),W(1,8),W(1,2),GC_6,ZERO,ZERO
 
745
     $     ,PL(0,76),COEFS)
 
746
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,56),4,COEFS,4,4,WL(1,0
714
747
     $     ,1,76))
715
 
          CALL MP_FFV1L3_1(PL(0,76),W(1,3),GC_5,MDL_MT,MDL_WT,PL(0,77)
716
 
     $     ,COEFS)
717
 
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,76),4,COEFS,4,4,WL(1,0
 
748
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,74),1,4,25,2,1,118
 
749
     $     ,H)
 
750
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,75),1,4,25,2,1,119
 
751
     $     ,H)
 
752
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,76),1,4,25,2,1,120
 
753
     $     ,H)
 
754
C         Coefficient construction for loop diagram with ID 35
 
755
          CALL MP_VVVV1L2P0_1(PL(0,30),W(1,8),W(1,1),GC_6,ZERO,ZERO
 
756
     $     ,PL(0,77),COEFS)
 
757
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,30),4,COEFS,4,4,WL(1,0
718
758
     $     ,1,77))
719
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,77),3,4,18,1,115,H)
720
 
C         Coefficient construction for loop diagram with ID 34
721
 
          CALL MP_VVV1L2P0_1(PL(0,68),W(1,2),GC_4,ZERO,ZERO,PL(0,78)
722
 
     $     ,COEFS)
723
 
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,68),4,COEFS,4,4,WL(1,0
 
759
          CALL MP_VVVV3L2P0_1(PL(0,30),W(1,8),W(1,1),GC_6,ZERO,ZERO
 
760
     $     ,PL(0,78),COEFS)
 
761
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,30),4,COEFS,4,4,WL(1,0
724
762
     $     ,1,78))
725
 
          CALL MP_FFV1L3_2(PL(0,78),W(1,4),GC_5,MDL_MT,MDL_WT,PL(0,79)
726
 
     $     ,COEFS)
727
 
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,78),4,COEFS,4,4,WL(1,0
 
763
          CALL MP_VVVV4L2P0_1(PL(0,30),W(1,8),W(1,1),GC_6,ZERO,ZERO
 
764
     $     ,PL(0,79),COEFS)
 
765
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,30),4,COEFS,4,4,WL(1,0
728
766
     $     ,1,79))
729
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,79),3,4,19,1,116,H)
730
 
C         Coefficient construction for loop diagram with ID 35
731
 
          CALL MP_VVV1L2P0_1(PL(0,0),W(1,1),GC_4,ZERO,ZERO,PL(0,80)
732
 
     $     ,COEFS)
733
 
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
 
767
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,77),1,4,26,2,1,121
 
768
     $     ,H)
 
769
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,78),1,4,26,2,1,122
 
770
     $     ,H)
 
771
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,79),1,4,26,2,1,123
 
772
     $     ,H)
 
773
C         Coefficient construction for loop diagram with ID 36
 
774
          CALL MP_VVVV1L2P0_1(PL(0,0),W(1,1),W(1,2),GC_6,ZERO,ZERO
 
775
     $     ,PL(0,80),COEFS)
 
776
          CALL MP_ML5_0_UPDATE_WL_0_0(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
734
777
     $     ,80))
735
 
          CALL MP_FFV1L3_2(PL(0,80),W(1,4),GC_5,MDL_MT,MDL_WT,PL(0,81)
736
 
     $     ,COEFS)
737
 
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,80),4,COEFS,4,4,WL(1,0
738
 
     $     ,1,81))
739
 
          CALL MP_FFV1L1P0_3(PL(0,81),W(1,10),GC_5,ZERO,ZERO,PL(0,82)
740
 
     $     ,COEFS)
741
 
          CALL MP_ML5_0_UPDATE_WL_2_0(WL(1,0,1,81),4,COEFS,4,4,WL(1,0
742
 
     $     ,1,82))
743
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,82),2,4,20,1,117,H)
744
 
C         Coefficient construction for loop diagram with ID 36
745
 
          CALL MP_FFV1L3_1(PL(0,80),W(1,3),GC_5,MDL_MT,MDL_WT,PL(0,83)
746
 
     $     ,COEFS)
747
 
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,80),4,COEFS,4,4,WL(1,0
 
778
          CALL MP_VVVV3L2P0_1(PL(0,0),W(1,1),W(1,2),GC_6,ZERO,ZERO
 
779
     $     ,PL(0,81),COEFS)
 
780
          CALL MP_ML5_0_UPDATE_WL_0_0(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
 
781
     $     ,81))
 
782
          CALL MP_VVVV4L2P0_1(PL(0,0),W(1,1),W(1,2),GC_6,ZERO,ZERO
 
783
     $     ,PL(0,82),COEFS)
 
784
          CALL MP_ML5_0_UPDATE_WL_0_0(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
 
785
     $     ,82))
 
786
          CALL MP_VVV1L2P0_1(PL(0,80),W(1,8),GC_4,ZERO,ZERO,PL(0,83)
 
787
     $     ,COEFS)
 
788
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,80),4,COEFS,4,4,WL(1,0
748
789
     $     ,1,83))
749
 
          CALL MP_FFV1L2P0_3(PL(0,83),W(1,9),GC_5,ZERO,ZERO,PL(0,84)
 
790
          CALL MP_VVV1L2P0_1(PL(0,81),W(1,8),GC_4,ZERO,ZERO,PL(0,84)
750
791
     $     ,COEFS)
751
 
          CALL MP_ML5_0_UPDATE_WL_2_0(WL(1,0,1,83),4,COEFS,4,4,WL(1,0
 
792
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,81),4,COEFS,4,4,WL(1,0
752
793
     $     ,1,84))
753
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,84),2,4,21,1,118,H)
 
794
          CALL MP_VVV1L2P0_1(PL(0,82),W(1,8),GC_4,ZERO,ZERO,PL(0,85)
 
795
     $     ,COEFS)
 
796
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,82),4,COEFS,4,4,WL(1,0
 
797
     $     ,1,85))
 
798
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,83),1,4,1,2,1,124,H)
 
799
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,84),1,4,1,2,1,125,H)
 
800
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,85),1,4,1,2,1,126,H)
754
801
C         Coefficient construction for loop diagram with ID 37
755
 
          CALL MP_FFV1L2_1(PL(0,83),W(1,2),GC_5,MDL_MT,MDL_WT,PL(0,85)
756
 
     $     ,COEFS)
757
 
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,83),4,COEFS,4,4,WL(1,0
758
 
     $     ,1,85))
759
 
          CALL MP_FFV1L2P0_3(PL(0,85),W(1,4),GC_5,ZERO,ZERO,PL(0,86)
760
 
     $     ,COEFS)
761
 
          CALL MP_ML5_0_UPDATE_WL_3_0(WL(1,0,1,85),4,COEFS,4,4,WL(1,0
762
 
     $     ,1,86))
763
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,86),3,4,22,1,119,H)
 
802
          CALL MP_GHGHGL2_1(PL(0,0),W(1,5),GC_4,ZERO,ZERO,PL(0,86)
 
803
     $     ,COEFS)
 
804
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),1,COEFS,1,1,WL(1,0,1
 
805
     $     ,86))
 
806
          CALL MP_GHGHGL2_1(PL(0,86),W(1,8),GC_4,ZERO,ZERO,PL(0,87)
 
807
     $     ,COEFS)
 
808
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,86),1,COEFS,1,1,WL(1,0
 
809
     $     ,1,87))
 
810
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,87),2,1,1,1,1,127,H)
764
811
C         Coefficient construction for loop diagram with ID 38
765
 
          CALL MP_VVVV1L2P0_1(PL(0,47),W(1,1),W(1,2),GC_6,ZERO,ZERO
766
 
     $     ,PL(0,87),COEFS)
767
 
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,47),4,COEFS,4,4,WL(1,0
768
 
     $     ,1,87))
769
 
          CALL MP_VVVV3L2P0_1(PL(0,47),W(1,1),W(1,2),GC_6,ZERO,ZERO
770
 
     $     ,PL(0,88),COEFS)
771
 
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,47),4,COEFS,4,4,WL(1,0
772
 
     $     ,1,88))
773
 
          CALL MP_VVVV4L2P0_1(PL(0,47),W(1,1),W(1,2),GC_6,ZERO,ZERO
774
 
     $     ,PL(0,89),COEFS)
775
 
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,47),4,COEFS,4,4,WL(1,0
 
812
          CALL MP_GHGHGL1_2(PL(0,0),W(1,1),GC_4,ZERO,ZERO,PL(0,88)
 
813
     $     ,COEFS)
 
814
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),1,COEFS,1,1,WL(1,0,1
 
815
     $     ,88))
 
816
          CALL MP_GHGHGL1_2(PL(0,88),W(1,2),GC_4,ZERO,ZERO,PL(0,89)
 
817
     $     ,COEFS)
 
818
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,88),1,COEFS,1,1,WL(1,0
776
819
     $     ,1,89))
777
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,87),1,4,7,1,120,H)
778
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,88),1,4,7,1,121,H)
779
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,89),1,4,7,1,122,H)
 
820
          CALL MP_GHGHGL1_2(PL(0,89),W(1,8),GC_4,ZERO,ZERO,PL(0,90)
 
821
     $     ,COEFS)
 
822
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,89),1,COEFS,1,1,WL(1,0
 
823
     $     ,1,90))
 
824
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,90),3,1,2,1,1,128,H)
780
825
C         Coefficient construction for loop diagram with ID 39
781
 
          CALL MP_VVV1L2P0_1(PL(0,80),W(1,2),GC_4,ZERO,ZERO,PL(0,90)
782
 
     $     ,COEFS)
783
 
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,80),4,COEFS,4,4,WL(1,0
784
 
     $     ,1,90))
785
 
          CALL MP_FFV1L3_2(PL(0,90),W(1,4),GC_5,MDL_MT,MDL_WT,PL(0,91)
786
 
     $     ,COEFS)
787
 
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,90),4,COEFS,4,4,WL(1,0
788
 
     $     ,1,91))
789
 
          CALL MP_FFV1L1P0_3(PL(0,91),W(1,3),GC_5,ZERO,ZERO,PL(0,92)
790
 
     $     ,COEFS)
791
 
          CALL MP_ML5_0_UPDATE_WL_3_0(WL(1,0,1,91),4,COEFS,4,4,WL(1,0
 
826
          CALL MP_GHGHGL2_1(PL(0,0),W(1,1),GC_4,ZERO,ZERO,PL(0,91)
 
827
     $     ,COEFS)
 
828
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),1,COEFS,1,1,WL(1,0,1
 
829
     $     ,91))
 
830
          CALL MP_GHGHGL2_1(PL(0,91),W(1,2),GC_4,ZERO,ZERO,PL(0,92)
 
831
     $     ,COEFS)
 
832
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,91),1,COEFS,1,1,WL(1,0
792
833
     $     ,1,92))
793
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,92),3,4,23,1,123,H)
794
 
C         Coefficient construction for loop diagram with ID 40
795
 
          CALL MP_FFV1L3_1(PL(0,90),W(1,3),GC_5,MDL_MT,MDL_WT,PL(0,93)
 
834
          CALL MP_GHGHGL2_1(PL(0,92),W(1,8),GC_4,ZERO,ZERO,PL(0,93)
796
835
     $     ,COEFS)
797
 
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,90),4,COEFS,4,4,WL(1,0
 
836
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,92),1,COEFS,1,1,WL(1,0
798
837
     $     ,1,93))
799
 
          CALL MP_FFV1L2P0_3(PL(0,93),W(1,4),GC_5,ZERO,ZERO,PL(0,94)
800
 
     $     ,COEFS)
801
 
          CALL MP_ML5_0_UPDATE_WL_3_0(WL(1,0,1,93),4,COEFS,4,4,WL(1,0
802
 
     $     ,1,94))
803
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,94),3,4,24,1,124,H)
804
 
C         Coefficient construction for loop diagram with ID 41
805
 
          CALL MP_VVV1L2P0_1(PL(0,0),W(1,5),GC_4,ZERO,ZERO,PL(0,95)
806
 
     $     ,COEFS)
807
 
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
808
 
     $     ,95))
809
 
          CALL MP_VVV1L2P0_1(PL(0,95),W(1,8),GC_4,ZERO,ZERO,PL(0,96)
810
 
     $     ,COEFS)
811
 
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,95),4,COEFS,4,4,WL(1,0
812
 
     $     ,1,96))
813
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,96),2,4,1,2,125,H)
814
 
C         Coefficient construction for loop diagram with ID 42
815
 
          CALL MP_VVV1L2P0_1(PL(0,90),W(1,8),GC_4,ZERO,ZERO,PL(0,97)
816
 
     $     ,COEFS)
817
 
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,90),4,COEFS,4,4,WL(1,0
818
 
     $     ,1,97))
819
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,97),3,4,2,1,126,H)
820
 
C         Coefficient construction for loop diagram with ID 43
821
 
          CALL MP_VVVV1L2P0_1(PL(0,80),W(1,8),W(1,2),GC_6,ZERO,ZERO
822
 
     $     ,PL(0,98),COEFS)
823
 
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,80),4,COEFS,4,4,WL(1,0
824
 
     $     ,1,98))
825
 
          CALL MP_VVVV3L2P0_1(PL(0,80),W(1,8),W(1,2),GC_6,ZERO,ZERO
826
 
     $     ,PL(0,99),COEFS)
827
 
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,80),4,COEFS,4,4,WL(1,0
828
 
     $     ,1,99))
829
 
          CALL MP_VVVV4L2P0_1(PL(0,80),W(1,8),W(1,2),GC_6,ZERO,ZERO
830
 
     $     ,PL(0,100),COEFS)
831
 
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,80),4,COEFS,4,4,WL(1,0
832
 
     $     ,1,100))
833
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,98),1,4,25,2,127,H)
834
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,99),1,4,25,2,128,H)
835
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,100),1,4,25,2,129,H)
836
 
C         Coefficient construction for loop diagram with ID 44
837
 
          CALL MP_VVVV1L2P0_1(PL(0,54),W(1,8),W(1,1),GC_6,ZERO,ZERO
838
 
     $     ,PL(0,101),COEFS)
839
 
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,54),4,COEFS,4,4,WL(1,0
840
 
     $     ,1,101))
841
 
          CALL MP_VVVV3L2P0_1(PL(0,54),W(1,8),W(1,1),GC_6,ZERO,ZERO
842
 
     $     ,PL(0,102),COEFS)
843
 
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,54),4,COEFS,4,4,WL(1,0
844
 
     $     ,1,102))
845
 
          CALL MP_VVVV4L2P0_1(PL(0,54),W(1,8),W(1,1),GC_6,ZERO,ZERO
846
 
     $     ,PL(0,103),COEFS)
847
 
          CALL MP_ML5_0_UPDATE_WL_1_0(WL(1,0,1,54),4,COEFS,4,4,WL(1,0
848
 
     $     ,1,103))
849
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,101),1,4,26,2,130,H)
850
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,102),1,4,26,2,131,H)
851
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,103),1,4,26,2,132,H)
852
 
C         Coefficient construction for loop diagram with ID 45
853
 
          CALL MP_VVVV1L2P0_1(PL(0,0),W(1,1),W(1,2),GC_6,ZERO,ZERO
854
 
     $     ,PL(0,104),COEFS)
855
 
          CALL MP_ML5_0_UPDATE_WL_0_0(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
856
 
     $     ,104))
857
 
          CALL MP_VVVV3L2P0_1(PL(0,0),W(1,1),W(1,2),GC_6,ZERO,ZERO
858
 
     $     ,PL(0,105),COEFS)
859
 
          CALL MP_ML5_0_UPDATE_WL_0_0(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
860
 
     $     ,105))
861
 
          CALL MP_VVVV4L2P0_1(PL(0,0),W(1,1),W(1,2),GC_6,ZERO,ZERO
862
 
     $     ,PL(0,106),COEFS)
863
 
          CALL MP_ML5_0_UPDATE_WL_0_0(WL(1,0,1,0),4,COEFS,4,4,WL(1,0,1
864
 
     $     ,106))
865
 
          CALL MP_VVV1L2P0_1(PL(0,104),W(1,8),GC_4,ZERO,ZERO,PL(0,107)
866
 
     $     ,COEFS)
867
 
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,104),4,COEFS,4,4,WL(1,0
868
 
     $     ,1,107))
869
 
          CALL MP_VVV1L2P0_1(PL(0,105),W(1,8),GC_4,ZERO,ZERO,PL(0,108)
870
 
     $     ,COEFS)
871
 
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,105),4,COEFS,4,4,WL(1,0
872
 
     $     ,1,108))
873
 
          CALL MP_VVV1L2P0_1(PL(0,106),W(1,8),GC_4,ZERO,ZERO,PL(0,109)
874
 
     $     ,COEFS)
875
 
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,106),4,COEFS,4,4,WL(1,0
876
 
     $     ,1,109))
877
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,107),1,4,1,2,133,H)
878
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,108),1,4,1,2,134,H)
879
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,109),1,4,1,2,135,H)
880
 
C         Coefficient construction for loop diagram with ID 46
881
 
          CALL MP_GHGHGL2_1(PL(0,0),W(1,5),GC_4,ZERO,ZERO,PL(0,110)
882
 
     $     ,COEFS)
883
 
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),1,COEFS,1,1,WL(1,0,1
884
 
     $     ,110))
885
 
          CALL MP_GHGHGL2_1(PL(0,110),W(1,8),GC_4,ZERO,ZERO,PL(0,111)
886
 
     $     ,COEFS)
887
 
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,110),1,COEFS,1,1,WL(1,0
888
 
     $     ,1,111))
889
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,111),2,1,1,1,136,H)
890
 
C         Coefficient construction for loop diagram with ID 47
891
 
          CALL MP_GHGHGL1_2(PL(0,0),W(1,1),GC_4,ZERO,ZERO,PL(0,112)
892
 
     $     ,COEFS)
893
 
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),1,COEFS,1,1,WL(1,0,1
894
 
     $     ,112))
895
 
          CALL MP_GHGHGL1_2(PL(0,112),W(1,2),GC_4,ZERO,ZERO,PL(0,113)
896
 
     $     ,COEFS)
897
 
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,112),1,COEFS,1,1,WL(1,0
898
 
     $     ,1,113))
899
 
          CALL MP_GHGHGL1_2(PL(0,113),W(1,8),GC_4,ZERO,ZERO,PL(0,114)
900
 
     $     ,COEFS)
901
 
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,113),1,COEFS,1,1,WL(1,0
902
 
     $     ,1,114))
903
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,114),3,1,2,1,137,H)
904
 
C         Coefficient construction for loop diagram with ID 48
905
 
          CALL MP_GHGHGL2_1(PL(0,0),W(1,1),GC_4,ZERO,ZERO,PL(0,115)
906
 
     $     ,COEFS)
907
 
          CALL MP_ML5_0_UPDATE_WL_0_1(WL(1,0,1,0),1,COEFS,1,1,WL(1,0,1
908
 
     $     ,115))
909
 
          CALL MP_GHGHGL2_1(PL(0,115),W(1,2),GC_4,ZERO,ZERO,PL(0,116)
910
 
     $     ,COEFS)
911
 
          CALL MP_ML5_0_UPDATE_WL_1_1(WL(1,0,1,115),1,COEFS,1,1,WL(1,0
912
 
     $     ,1,116))
913
 
          CALL MP_GHGHGL2_1(PL(0,116),W(1,8),GC_4,ZERO,ZERO,PL(0,117)
914
 
     $     ,COEFS)
915
 
          CALL MP_ML5_0_UPDATE_WL_2_1(WL(1,0,1,116),1,COEFS,1,1,WL(1,0
916
 
     $     ,1,117))
917
 
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,117),3,1,2,1,138,H)
 
838
          CALL MP_ML5_0_CREATE_LOOP_COEFS(WL(1,0,1,93),3,1,2,1,1,129,H)
918
839
 4000     CONTINUE
919
840
          MP_LOOP_REQ_SO_DONE=.TRUE.
920
841
 
 
842
 
 
843
C         Copy the qp wfs to the dp ones as they are used to setup the
 
844
C          CT calls.
 
845
C         This needs to be done once since only the momenta of these
 
846
C          WF matters.
 
847
          IF(.NOT.DPW_COPIED) THEN
 
848
            DO I=1,NWAVEFUNCS
 
849
              DO J=1,MAXLWFSIZE+4
 
850
                DPW(J,I)=DBLE(W(J,I))
 
851
              ENDDO
 
852
            ENDDO
 
853
            DPW_COPIED=.TRUE.
 
854
          ENDIF
 
855
 
 
856
 
 
857
 
 
858
 
 
859
 
921
860
        ENDIF
922
861
      ENDDO
923
862
 
924
 
C     Copy the qp wfs to the dp ones as they are used to setup the CT
925
 
C      calls.
926
 
      DO I=1,NWAVEFUNCS
927
 
        DO J=1,MAXLWFSIZE+4
928
 
          DPW(J,I)=W(J,I)
929
 
        ENDDO
930
 
      ENDDO
931
 
 
932
863
      DO I=1,3
933
864
        DO J=0,NSQUAREDSO
934
865
          ANSDP(I,J)=ANS(I,J)