4 SUBROUTINE uux(P,I3,I4,H1,H2,KEY,ANS)
7 INTEGER I3,I4, H1,H2, KEY
8 real*8 p(0:3,6), ans, ans1,ans2,ans3,ans4, ans5,ans6
27 IF(i3.EQ.1 .AND. i4.EQ.-1)
CALL uux_ddx_h(p,h1,h2,ans)
28 IF(i3.EQ.2 .AND. i4.EQ.-2)
CALL uux_uux_h(p,h1,h2,ans)
29 IF(i3.EQ.3 .AND. i4.EQ.-3)
CAll uux_ssx_h(p,h1,h2,ans)
30 IF(i3.EQ.4 .AND. i4.EQ.-4)
CALL uux_ccx_h(p,h1,h2,ans)
31 IF(i3.EQ.1 .AND. i4.EQ.-3)
CALL uux_dsx_h(p,h1,h2,ans)
32 IF(i3.EQ.3 .AND. i4.EQ.-1)
CALL uux_sdx_h(p,h1,h2,ans)
34 IF(i3.EQ.-1 .AND. i4.EQ.1)
CALL uux_ddx_h(pp,h1,h2,ans)
35 IF(i3.EQ.-2 .AND. i4.EQ.2)
CALL uux_uux_h(pp,h1,h2,ans)
36 IF(i3.EQ.-3 .AND. i4.EQ.3)
CAll uux_ssx_h(pp,h1,h2,ans)
37 IF(i3.EQ.-4 .AND. i4.EQ.4)
CALL uux_ccx_h(pp,h1,h2,ans)
38 IF(i3.EQ.-3 .AND. i4.EQ.1)
CALL uux_dsx_h(pp,h1,h2,ans)
39 IF(i3.EQ.-1 .AND. i4.EQ.3)
CALL uux_sdx_h(pp,h1,h2,ans)
40 IF(i3.EQ.0 .AND. i4.EQ.0)
THEN
41 CALL uux_ddx_h(p,h1,h2,ans1)
42 CALL uux_uux_h(p,h1,h2,ans2)
43 CALL uux_ssx_h(p,h1,h2,ans3)
44 CALL uux_ccx_h(p,h1,h2,ans4)
45 CALL uux_dsx_h(p,h1,h2,ans5)
46 CALL uux_sdx_h(p,h1,h2,ans6)
47 ans=ans1+ans2+ans3+ans4+ans5+ans6
48 CALL uux_ddx_h(pp,h1,h2,ans1)
49 CALL uux_uux_h(pp,h1,h2,ans2)
50 CALL uux_ssx_h(pp,h1,h2,ans3)
51 CALL uux_ccx_h(pp,h1,h2,ans4)
52 CALL uux_dsx_h(pp,h1,h2,ans5)
53 CALL uux_sdx_h(pp,h1,h2,ans6)
54 ans=ans1+ans2+ans3+ans4+ans5+ans6+ans
56 ELSE IF(key.EQ.0)
THEN
57 IF(i3.EQ.1 .AND. i4.EQ.-1)
CALL uux_ddx_noh(p,h1,h2,ans)
58 IF(i3.EQ.2 .AND. i4.EQ.-2)
CALL uux_uux_noh(p,h1,h2,ans)
59 IF(i3.EQ.3 .AND. i4.EQ.-3)
CAll uux_ssx_noh(p,h1,h2,ans)
60 IF(i3.EQ.4 .AND. i4.EQ.-4)
CALL uux_ccx_noh(p,h1,h2,ans)
61 IF(i3.EQ.1 .AND. i4.EQ.-3)
CALL uux_dsx_noh(p,h1,h2,ans)
62 IF(i3.EQ.3 .AND. i4.EQ.-1)
CALL uux_sdx_noh(p,h1,h2,ans)
64 IF(i3.EQ.21 .and .i4.EQ.21)
CALL uux_gg_noh(p,h1,h2,ans)
66 IF(i3.EQ.-1 .AND. i4.EQ.1)
CALL uux_ddx_noh(pp,h1,h2,ans)
67 IF(i3.EQ.-2 .AND. i4.EQ.2)
CALL uux_uux_noh(pp,h1,h2,ans)
68 IF(i3.EQ.-3 .AND. i4.EQ.3)
CAll uux_ssx_noh(pp,h1,h2,ans)
69 IF(i3.EQ.-4 .AND. i4.EQ.4)
CALL uux_ccx_noh(pp,h1,h2,ans)
70 IF(i3.EQ.-3 .AND. i4.EQ.1)
CALL uux_dsx_noh(pp,h1,h2,ans)
71 IF(i3.EQ.-1 .AND. i4.EQ.3)
CALL uux_sdx_noh(pp,h1,h2,ans)
73 IF(i3.EQ.0 .AND. i4.EQ.0)
THEN
74 CALL uux_ddx_noh(p,h1,h2,ans1)
75 CALL uux_uux_noh(p,h1,h2,ans2)
76 CALL uux_ssx_noh(p,h1,h2,ans3)
77 CALL uux_ccx_noh(p,h1,h2,ans4)
78 CALL uux_dsx_noh(p,h1,h2,ans5)
79 CALL uux_sdx_noh(p,h1,h2,ans6)
80 ans=ans1+ans2+ans3+ans4+ans5+ans6
81 CALL uux_ddx_noh(pp,h1,h2,ans1)
82 CALL uux_uux_noh(pp,h1,h2,ans2)
83 CALL uux_ssx_noh(pp,h1,h2,ans3)
84 CALL uux_ccx_noh(pp,h1,h2,ans4)
85 CALL uux_dsx_noh(pp,h1,h2,ans5)
86 CALL uux_sdx_noh(pp,h1,h2,ans6)
87 ans=ans1+ans2+ans3+ans4+ans5+ans6+ans
88 CALL uux_gg_noh(p,h1,h2,ans1)
92 WRITE(*,*)
'NOT FINISHED'
101 SUBROUTINE uux_ddx_h(P,H1,H2,ANS)
121 parameter(nexternal=6)
127 real*8 p(0:3,nexternal),ans
132 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
134 real*8 matrix_uux_ddx_h
136 INTEGER JC(NEXTERNAL)
137 LOGICAL GOODHEL(NCOMB)
139 DATA goodhel/ncomb*.false./
140 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
141 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
142 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
143 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
144 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
145 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
146 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
147 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
148 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
149 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
150 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
151 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
152 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
153 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
154 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
155 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
156 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
157 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
158 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
159 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
160 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
161 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
162 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
163 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
164 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
165 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
166 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
167 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
168 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
169 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
170 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
171 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
172 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
173 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
174 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
175 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
176 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
177 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
178 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
179 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
180 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
181 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
182 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
183 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
184 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
185 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
186 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
187 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
188 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
189 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
190 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
191 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
192 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
193 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
194 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
195 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
196 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
197 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
198 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
199 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
200 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
201 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
202 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
203 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
213 t=matrix_uux_ddx_h(p ,h1,h2,nhel(1,ihel),jc(1))
220 REAL*8 FUNCTION matrix_uux_ddx_h(P,H1,H2,NHEL,IC)
239 parameter(nexternal=6)
240 INTEGER NWAVEFUNCS, NCOLOR
241 parameter(nwavefuncs=7, ncolor=2)
245 parameter(imag1=(0d0,1d0))
249 real*8 p(0:3,nexternal)
250 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
256 real*8 denom(ncolor), cf(ncolor,ncolor)
257 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
258 COMPLEX*16 W(18,NWAVEFUNCS)
260 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
269 DATA (cf(i, 1),i= 1, 2) / 9, 3/
272 DATA (cf(i, 2),i= 1, 2) / 3, 9/
277 matrix_uux_ddx_h=0.d0
278 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
279 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
284 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
285 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
286 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
287 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
288 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
289 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
290 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
291 CALL ffv2_5_3(w(1,1),w(1,2),gc_51,gc_58,mz,wz,w(1,6))
292 CALL ffv2_3_3(w(1,4),w(1,3),gc_50,gc_58,mz,wz,w(1,5))
294 CALL vvs1_0(w(1,6),w(1,5),w(1,7),gc_81,amp(1))
295 CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,5))
296 CALL ffv2_3(w(1,4),w(1,2),gc_100,mw,ww,w(1,3))
298 CALL vvs1_0(w(1,5),w(1,3),w(1,7),gc_72,amp(2))
306 ztemp = ztemp + cf(j,i)*jamp(j)
308 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
310 matrix_uux_ddx_h=matrix
317 SUBROUTINE uux_uux_h(P,H1,H2,ANS)
337 parameter(nexternal=6)
343 real*8 p(0:3,nexternal),ans
348 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
350 real*8 matrix_uux_uux_h
352 INTEGER JC(NEXTERNAL)
353 LOGICAL GOODHEL(NCOMB)
355 DATA goodhel/ncomb*.false./
356 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
357 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
358 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
359 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
360 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
361 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
362 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
363 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
364 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
365 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
366 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
367 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
368 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
369 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
370 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
371 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
372 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
373 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
374 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
375 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
376 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
377 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
378 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
379 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
380 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
381 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
382 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
383 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
384 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
385 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
386 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
387 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
388 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
389 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
390 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
391 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
392 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
393 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
394 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
395 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
396 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
397 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
398 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
399 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
400 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
401 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
402 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
403 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
404 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
405 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
406 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
407 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
408 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
409 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
410 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
411 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
412 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
413 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
414 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
415 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
416 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
417 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
418 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
419 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
429 t=matrix_uux_uux_h(p ,h1,h2,nhel(1,ihel),jc(1))
436 REAL*8 FUNCTION matrix_uux_uux_h(P,H1,H2,NHEL,IC)
455 parameter(nexternal=6)
456 INTEGER NWAVEFUNCS, NCOLOR
457 parameter(nwavefuncs=7, ncolor=2)
461 parameter(imag1=(0d0,1d0))
465 real*8 p(0:3,nexternal)
466 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
472 real*8 denom(ncolor), cf(ncolor,ncolor)
473 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
474 COMPLEX*16 W(18,NWAVEFUNCS)
476 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
485 DATA (cf(i, 1),i= 1, 2) / 9, 3/
488 DATA (cf(i, 2),i= 1, 2) / 3, 9/
493 matrix_uux_uux_h=0.d0
494 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
495 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
500 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
501 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
502 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
503 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
504 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
505 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
506 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
507 CALL ffv2_5_3(w(1,1),w(1,2),gc_51,gc_58,mz,wz,w(1,6))
508 CALL ffv2_5_3(w(1,4),w(1,3),gc_51,gc_58,mz,wz,w(1,5))
510 CALL vvs1_0(w(1,6),w(1,5),w(1,7),gc_81,amp(1))
511 CALL ffv2_5_3(w(1,1),w(1,3),gc_51,gc_58,mz,wz,w(1,5))
512 CALL ffv2_5_3(w(1,4),w(1,2),gc_51,gc_58,mz,wz,w(1,3))
514 CALL vvs1_0(w(1,5),w(1,3),w(1,7),gc_81,amp(2))
522 ztemp = ztemp + cf(j,i)*jamp(j)
524 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
526 matrix_uux_uux_h=matrix
533 SUBROUTINE uux_ssx_h(P,H1,H2,ANS)
553 parameter(nexternal=6)
559 real*8 p(0:3,nexternal),ans
564 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
566 real*8 matrix_uux_ssx_h
568 INTEGER JC(NEXTERNAL)
569 LOGICAL GOODHEL(NCOMB)
571 DATA goodhel/ncomb*.false./
572 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
573 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
574 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
575 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
576 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
577 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
578 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
579 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
580 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
581 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
582 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
583 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
584 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
585 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
586 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
587 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
588 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
589 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
590 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
591 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
592 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
593 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
594 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
595 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
596 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
597 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
598 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
599 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
600 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
601 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
602 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
603 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
604 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
605 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
606 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
607 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
608 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
609 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
610 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
611 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
612 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
613 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
614 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
615 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
616 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
617 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
618 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
619 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
620 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
621 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
622 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
623 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
624 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
625 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
626 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
627 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
628 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
629 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
630 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
631 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
632 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
633 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
634 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
635 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
645 t=matrix_uux_ssx_h(p ,h1,h2,nhel(1,ihel),jc(1))
652 REAL*8 FUNCTION matrix_uux_ssx_h(P,H1,H2,NHEL,IC)
671 parameter(nexternal=6)
672 INTEGER NWAVEFUNCS, NCOLOR
673 parameter(nwavefuncs=7, ncolor=2)
677 parameter(imag1=(0d0,1d0))
681 real*8 p(0:3,nexternal)
682 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
688 real*8 denom(ncolor), cf(ncolor,ncolor)
689 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
690 COMPLEX*16 W(18,NWAVEFUNCS)
692 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
701 DATA (cf(i, 1),i= 1, 2) / 9, 3/
704 DATA (cf(i, 2),i= 1, 2) / 3, 9/
709 matrix_uux_ssx_h=0.d0
710 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
711 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
716 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
717 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
718 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
719 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
720 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
721 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
722 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
723 CALL ffv2_5_3(w(1,1),w(1,2),gc_51,gc_58,mz,wz,w(1,6))
724 CALL ffv2_3_3(w(1,4),w(1,3),gc_50,gc_58,mz,wz,w(1,5))
726 CALL vvs1_0(w(1,6),w(1,5),w(1,7),gc_81,amp(1))
727 CALL ffv2_3(w(1,1),w(1,3),gc_101,mw,ww,w(1,5))
728 CALL ffv2_3(w(1,4),w(1,2),gc_101,mw,ww,w(1,3))
730 CALL vvs1_0(w(1,5),w(1,3),w(1,7),gc_72,amp(2))
738 ztemp = ztemp + cf(j,i)*jamp(j)
740 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
742 matrix_uux_ssx_h=matrix
749 SUBROUTINE uux_ccx_h(P,H1,H2,ANS)
769 parameter(nexternal=6)
775 real*8 p(0:3,nexternal),ans
780 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
782 real*8 matrix_uux_ccx_h
784 INTEGER JC(NEXTERNAL)
785 LOGICAL GOODHEL(NCOMB)
787 DATA goodhel/ncomb*.false./
788 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
789 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
790 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
791 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
792 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
793 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
794 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
795 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
796 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
797 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
798 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
799 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
800 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
801 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
802 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
803 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
804 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
805 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
806 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
807 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
808 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
809 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
810 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
811 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
812 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
813 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
814 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
815 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
816 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
817 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
818 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
819 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
820 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
821 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
822 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
823 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
824 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
825 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
826 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
827 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
828 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
829 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
830 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
831 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
832 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
833 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
834 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
835 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
836 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
837 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
838 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
839 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
840 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
841 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
842 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
843 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
844 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
845 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
846 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
847 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
848 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
849 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
850 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
851 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
861 t=matrix_uux_ccx_h(p ,h1,h2,nhel(1,ihel),jc(1))
868 REAL*8 FUNCTION matrix_uux_ccx_h(P,H1,H2,NHEL,IC)
887 parameter(nexternal=6)
888 INTEGER NWAVEFUNCS, NCOLOR
889 parameter(nwavefuncs=7, ncolor=1)
893 parameter(imag1=(0d0,1d0))
897 real*8 p(0:3,nexternal)
898 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
904 real*8 denom(ncolor), cf(ncolor,ncolor)
905 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
906 COMPLEX*16 W(18,NWAVEFUNCS)
908 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
917 DATA (cf(i, 1),i= 1, 1) / 9/
922 matrix_uux_ccx_h=0.d0
923 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
924 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
929 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
930 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
931 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
932 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
933 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
934 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
935 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
936 CALL ffv2_5_3(w(1,1),w(1,2),gc_51,gc_58,mz,wz,w(1,6))
937 CALL ffv2_5_3(w(1,4),w(1,3),gc_51,gc_58,mz,wz,w(1,2))
939 CALL vvs1_0(w(1,6),w(1,2),w(1,7),gc_81,amp(1))
946 ztemp = ztemp + cf(j,i)*jamp(j)
948 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
950 matrix_uux_ccx_h=matrix
957 SUBROUTINE uux_dsx_h(P,H1,H2,ANS)
977 parameter(nexternal=6)
983 real*8 p(0:3,nexternal),ans
988 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
990 real*8 matrix_uux_dsx_h
992 INTEGER JC(NEXTERNAL)
993 LOGICAL GOODHEL(NCOMB)
995 DATA goodhel/ncomb*.false./
996 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
997 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
998 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
999 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1000 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1001 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1002 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1003 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1004 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1005 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1006 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1007 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1008 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1009 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1010 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1011 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1012 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1013 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1014 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1015 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1016 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1017 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1018 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1019 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1020 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1021 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1022 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1023 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1024 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1025 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1026 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1027 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1028 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1029 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1030 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1031 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1032 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1033 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1034 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1035 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1036 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1037 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1038 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1039 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1040 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1041 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1042 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1043 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1044 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1045 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1046 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1047 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1048 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1049 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1050 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1051 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1052 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1053 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1054 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1055 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1056 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1057 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1058 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1059 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1069 t=matrix_uux_dsx_h(p ,h1,h2,nhel(1,ihel),jc(1))
1076 REAL*8 FUNCTION matrix_uux_dsx_h(P,H1,H2,NHEL,IC)
1093 parameter(ngraphs=1)
1095 parameter(nexternal=6)
1096 INTEGER NWAVEFUNCS, NCOLOR
1097 parameter(nwavefuncs=7, ncolor=1)
1101 parameter(imag1=(0d0,1d0))
1105 real*8 p(0:3,nexternal)
1106 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1112 real*8 denom(ncolor), cf(ncolor,ncolor)
1113 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1114 COMPLEX*16 W(18,NWAVEFUNCS)
1115 COMPLEX*16 DUM0,DUM1
1116 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1125 DATA (cf(i, 1),i= 1, 1) / 9/
1130 matrix_uux_dsx_h=0.d0
1131 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
1132 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
1137 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1138 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1139 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1140 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1141 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1142 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1143 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
1144 CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,6))
1145 CALL ffv2_3(w(1,4),w(1,2),gc_101,mw,ww,w(1,3))
1147 CALL vvs1_0(w(1,6),w(1,3),w(1,7),gc_72,amp(1))
1154 ztemp = ztemp + cf(j,i)*jamp(j)
1156 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1158 matrix_uux_dsx_h=matrix
1165 SUBROUTINE uux_sdx_h(P,H1,H2,ANS)
1185 parameter(nexternal=6)
1187 parameter( ncomb=64)
1191 real*8 p(0:3,nexternal),ans
1196 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1198 real*8 matrix_uux_sdx_h
1199 INTEGER IHEL,IDEN, I
1200 INTEGER JC(NEXTERNAL)
1201 LOGICAL GOODHEL(NCOMB)
1203 DATA goodhel/ncomb*.false./
1204 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1205 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1206 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1207 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1208 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1209 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1210 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1211 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1212 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1213 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1214 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1215 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1216 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1217 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1218 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1219 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1220 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1221 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1222 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1223 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1224 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1225 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1226 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1227 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1228 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1229 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1230 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1231 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1232 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1233 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1234 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1235 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1236 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1237 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1238 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1239 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1240 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1241 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1242 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1243 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1244 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1245 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1246 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1247 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1248 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1249 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1250 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1251 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1252 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1253 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1254 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1255 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1256 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1257 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1258 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1259 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1260 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1261 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1262 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1263 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1264 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1265 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1266 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1267 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1277 t=matrix_uux_sdx_h(p ,h1,h2,nhel(1,ihel),jc(1))
1284 REAL*8 FUNCTION matrix_uux_sdx_h(P,H1,H2,NHEL,IC)
1301 parameter(ngraphs=1)
1303 parameter(nexternal=6)
1304 INTEGER NWAVEFUNCS, NCOLOR
1305 parameter(nwavefuncs=7, ncolor=1)
1309 parameter(imag1=(0d0,1d0))
1313 real*8 p(0:3,nexternal)
1314 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1320 real*8 denom(ncolor), cf(ncolor,ncolor)
1321 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1322 COMPLEX*16 W(18,NWAVEFUNCS)
1323 COMPLEX*16 DUM0,DUM1
1324 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1333 DATA (cf(i, 1),i= 1, 1) / 9/
1338 matrix_uux_sdx_h=0.d0
1339 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
1340 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
1345 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1346 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1347 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1348 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1349 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1350 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1351 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
1352 CALL ffv2_3(w(1,1),w(1,3),gc_101,mw,ww,w(1,6))
1353 CALL ffv2_3(w(1,4),w(1,2),gc_100,mw,ww,w(1,3))
1355 CALL vvs1_0(w(1,6),w(1,3),w(1,7),gc_72,amp(1))
1362 ztemp = ztemp + cf(j,i)*jamp(j)
1364 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1366 matrix_uux_sdx_h=matrix
1376 SUBROUTINE uux_ddx_noh(P,H1,H2,ANS)
1395 parameter(nexternal=6)
1397 parameter( ncomb=64)
1401 real*8 p(0:3,nexternal),ans
1406 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1408 real*8 matrix_uux_ddx_noh
1409 INTEGER IHEL,IDEN, I
1410 INTEGER JC(NEXTERNAL)
1411 LOGICAL GOODHEL(NCOMB)
1413 DATA goodhel/ncomb*.false./
1414 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1415 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1416 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1417 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1418 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1419 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1420 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1421 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1422 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1423 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1424 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1425 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1426 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1427 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1428 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1429 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1430 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1431 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1432 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1433 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1434 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1435 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1436 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1437 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1438 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1439 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1440 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1441 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1442 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1443 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1444 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1445 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1446 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1447 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1448 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1449 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1450 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1451 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1452 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1453 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1454 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1455 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1456 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1457 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1458 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1459 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1460 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1461 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1462 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1463 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1464 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1465 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1466 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1467 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1468 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1469 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1470 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1471 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1472 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1473 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1474 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1475 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1476 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1477 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1487 t=matrix_uux_ddx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1494 REAL*8 FUNCTION matrix_uux_ddx_noh(P,H1,H2,NHEL,IC)
1510 parameter(ngraphs=43)
1512 parameter(nexternal=6)
1513 INTEGER NWAVEFUNCS, NCOLOR
1514 parameter(nwavefuncs=14, ncolor=2)
1518 parameter(imag1=(0d0,1d0))
1522 real*8 p(0:3,nexternal)
1523 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1529 real*8 denom(ncolor), cf(ncolor,ncolor)
1530 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1531 COMPLEX*16 W(18,NWAVEFUNCS)
1532 COMPLEX*16 DUM0,DUM1
1533 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1542 DATA (cf(i, 1),i= 1, 2) / 9, 3/
1545 DATA (cf(i, 2),i= 1, 2) / 3, 9/
1550 matrix_uux_ddx_noh=0.d0
1551 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
1552 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
1558 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1559 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1560 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1561 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1562 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1563 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1564 CALL ffv1p0_3(w(1,1),w(1,2),gc_2,zero,zero,w(1,7))
1565 CALL ffv1p0_3(w(1,4),w(1,3),gc_1,zero,zero,w(1,8))
1566 CALL ffv1_2(w(1,5),w(1,7),gc_3,mta,zero,w(1,9))
1568 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(1))
1569 CALL ffv1_1(w(1,6),w(1,7),gc_3,mta,zero,w(1,10))
1571 CALL ffv1_0(w(1,5),w(1,10),w(1,8),gc_3,amp(2))
1572 CALL ffv2_3_3(w(1,4),w(1,3),gc_50,gc_58,mz,wz,w(1,11))
1574 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(3))
1576 CALL ffv2_4_0(w(1,5),w(1,10),w(1,11),gc_50,gc_59,amp(4))
1577 CALL ffv2_5_3(w(1,1),w(1,2),gc_51,gc_58,mz,wz,w(1,10))
1578 CALL ffv2_4_2(w(1,5),w(1,10),gc_50,gc_59,mta,zero,w(1,9))
1580 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(5))
1581 CALL ffv2_4_1(w(1,6),w(1,10),gc_50,gc_59,mta,zero,w(1,12))
1583 CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
1585 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(7))
1587 CALL ffv2_4_0(w(1,5),w(1,12),w(1,11),gc_50,gc_59,amp(8))
1588 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,12))
1589 CALL ffv1_1(w(1,3),w(1,7),gc_1,zero,zero,w(1,9))
1591 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_1,amp(9))
1592 CALL ffv1_2(w(1,4),w(1,7),gc_1,zero,zero,w(1,13))
1594 CALL ffv1_0(w(1,13),w(1,3),w(1,12),gc_1,amp(10))
1595 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
1597 CALL ffv2_3_0(w(1,4),w(1,9),w(1,7),gc_50,gc_58,amp(11))
1599 CALL ffv2_3_0(w(1,13),w(1,3),w(1,7),gc_50,gc_58,amp(12))
1600 CALL ffv1p0_3(w(1,1),w(1,2),gc_11,zero,zero,w(1,13))
1601 CALL ffv1_1(w(1,3),w(1,13),gc_11,zero,zero,w(1,9))
1603 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_1,amp(13))
1604 CALL ffv1_2(w(1,4),w(1,13),gc_11,zero,zero,w(1,14))
1606 CALL ffv1_0(w(1,14),w(1,3),w(1,12),gc_1,amp(14))
1608 CALL ffv2_3_0(w(1,4),w(1,9),w(1,7),gc_50,gc_58,amp(15))
1610 CALL ffv2_3_0(w(1,14),w(1,3),w(1,7),gc_50,gc_58,amp(16))
1611 CALL ffv2_3_1(w(1,3),w(1,10),gc_50,gc_58,zero,zero,w(1,14))
1613 CALL ffv1_0(w(1,4),w(1,14),w(1,12),gc_1,amp(17))
1614 CALL ffv2_3_2(w(1,4),w(1,10),gc_50,gc_58,zero,zero,w(1,9))
1616 CALL ffv1_0(w(1,9),w(1,3),w(1,12),gc_1,amp(18))
1618 CALL ffv2_3_0(w(1,4),w(1,14),w(1,7),gc_50,gc_58,amp(19))
1620 CALL ffv2_3_0(w(1,9),w(1,3),w(1,7),gc_50,gc_58,amp(20))
1621 CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,9))
1622 CALL ffv2_3(w(1,4),w(1,2),gc_100,mw,ww,w(1,14))
1623 CALL ffv2_2(w(1,5),w(1,9),gc_108,zero,zero,w(1,10))
1625 CALL ffv2_0(w(1,10),w(1,6),w(1,14),gc_108,amp(21))
1627 CALL vvv1_0(w(1,12),w(1,9),w(1,14),gc_4,amp(22))
1629 CALL vvv1_0(w(1,9),w(1,14),w(1,7),gc_53,amp(23))
1630 CALL ffv2_1(w(1,2),w(1,9),gc_100,zero,zero,w(1,10))
1632 CALL ffv1_0(w(1,4),w(1,10),w(1,12),gc_1,amp(24))
1633 CALL ffv2_2(w(1,4),w(1,9),gc_100,zero,zero,w(1,6))
1635 CALL ffv1_0(w(1,6),w(1,2),w(1,12),gc_2,amp(25))
1637 CALL ffv2_3_0(w(1,4),w(1,10),w(1,7),gc_50,gc_58,amp(26))
1639 CALL ffv2_5_0(w(1,6),w(1,2),w(1,7),gc_51,gc_58,amp(27))
1640 CALL ffv2_2(w(1,1),w(1,14),gc_100,zero,zero,w(1,6))
1642 CALL ffv1_0(w(1,6),w(1,3),w(1,12),gc_1,amp(28))
1643 CALL ffv1_2(w(1,1),w(1,12),gc_2,zero,zero,w(1,10))
1645 CALL ffv2_0(w(1,10),w(1,3),w(1,14),gc_100,amp(29))
1647 CALL ffv2_3_0(w(1,6),w(1,3),w(1,7),gc_50,gc_58,amp(30))
1648 CALL ffv2_5_2(w(1,1),w(1,7),gc_51,gc_58,zero,zero,w(1,6))
1650 CALL ffv2_0(w(1,6),w(1,3),w(1,14),gc_100,amp(31))
1651 CALL ffv1_2(w(1,1),w(1,8),gc_2,zero,zero,w(1,14))
1653 CALL ffv1_0(w(1,14),w(1,2),w(1,12),gc_2,amp(32))
1655 CALL ffv1_0(w(1,10),w(1,2),w(1,8),gc_2,amp(33))
1657 CALL ffv2_5_0(w(1,14),w(1,2),w(1,7),gc_51,gc_58,amp(34))
1659 CALL ffv1_0(w(1,6),w(1,2),w(1,8),gc_2,amp(35))
1660 CALL ffv1p0_3(w(1,4),w(1,3),gc_11,zero,zero,w(1,8))
1661 CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,4))
1663 CALL ffv1_0(w(1,4),w(1,2),w(1,12),gc_2,amp(36))
1665 CALL ffv1_0(w(1,10),w(1,2),w(1,8),gc_11,amp(37))
1667 CALL ffv2_5_0(w(1,4),w(1,2),w(1,7),gc_51,gc_58,amp(38))
1669 CALL ffv1_0(w(1,6),w(1,2),w(1,8),gc_11,amp(39))
1670 CALL ffv2_5_2(w(1,1),w(1,11),gc_51,gc_58,zero,zero,w(1,8))
1672 CALL ffv1_0(w(1,8),w(1,2),w(1,12),gc_2,amp(40))
1674 CALL ffv2_5_0(w(1,10),w(1,2),w(1,11),gc_51,gc_58,amp(41))
1676 CALL ffv2_5_0(w(1,8),w(1,2),w(1,7),gc_51,gc_58,amp(42))
1678 CALL ffv2_5_0(w(1,6),w(1,2),w(1,11),gc_51,gc_58,amp(43))
1679 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
1680 $ -amp(9)-amp(10)-amp(11)-amp(12)+1./6.*amp(13)+1./6.*amp(14)
1681 $ +1./6.*amp(15)+1./6.*amp(16)-amp(17)-amp(18)-amp(19)-amp(20)
1682 $ -amp(32)-amp(33)-amp(34)-amp(35)+1./6.*amp(36)+1./6.*amp(37)
1683 $ +1./6.*amp(38)+1./6.*amp(39)-amp(40)-amp(41)-amp(42)-amp(43)
1684 jamp(2)=-1./2.*amp(13)-1./2.*amp(14)-1./2.*amp(15)-1./2.*amp(16)
1685 $ +amp(21)+amp(22)+amp(23)+amp(24)+amp(25)+amp(26)+amp(27)
1686 $ +amp(28)+amp(29)+amp(30)+amp(31)-1./2.*amp(36)-1./2.*amp(37)
1687 $ -1./2.*amp(38)-1./2.*amp(39)
1693 ztemp = ztemp + cf(j,i)*jamp(j)
1695 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1697 matrix_uux_ddx_noh=matrix
1704 SUBROUTINE uux_uux_noh(P,H1,H2,ANS)
1723 parameter(nexternal=6)
1725 parameter( ncomb=64)
1729 real*8 p(0:3,nexternal),ans
1734 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1736 real*8 matrix_uux_uux_noh
1737 INTEGER IHEL,IDEN, I
1738 INTEGER JC(NEXTERNAL)
1739 LOGICAL GOODHEL(NCOMB)
1741 DATA goodhel/ncomb*.false./
1742 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1743 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1744 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1745 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1746 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1747 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1748 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1749 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1750 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1751 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1752 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1753 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1754 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1755 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1756 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1757 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1758 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1759 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1760 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1761 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1762 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1763 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1764 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1765 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1766 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1767 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1768 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1769 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1770 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1771 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1772 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1773 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1774 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1775 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1776 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1777 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1778 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1779 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1780 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1781 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1782 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1783 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1784 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1785 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1786 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1787 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1788 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1789 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1790 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1791 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1792 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1793 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1794 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1795 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1796 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1797 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1798 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1799 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1800 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1801 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1802 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1803 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1804 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1805 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1815 t=matrix_uux_uux_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1822 REAL*8 FUNCTION matrix_uux_uux_noh(P,H1,H2,NHEL,IC)
1838 parameter(ngraphs=64)
1840 parameter(nexternal=6)
1841 INTEGER NWAVEFUNCS, NCOLOR
1842 parameter(nwavefuncs=16, ncolor=2)
1846 parameter(imag1=(0d0,1d0))
1850 real*8 p(0:3,nexternal)
1851 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1857 real*8 denom(ncolor), cf(ncolor,ncolor)
1858 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1859 COMPLEX*16 W(18,NWAVEFUNCS)
1860 COMPLEX*16 DUM0,DUM1
1861 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1870 DATA (cf(i, 1),i= 1, 2) / 9, 3/
1873 DATA (cf(i, 2),i= 1, 2) / 3, 9/
1878 matrix_uux_uux_noh=0.d0
1879 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
1880 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
1885 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1886 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1887 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1888 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1889 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1890 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1891 CALL ffv1p0_3(w(1,1),w(1,2),gc_2,zero,zero,w(1,7))
1892 CALL ffv1p0_3(w(1,4),w(1,3),gc_2,zero,zero,w(1,8))
1893 CALL ffv1_2(w(1,5),w(1,7),gc_3,mta,zero,w(1,9))
1895 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(1))
1896 CALL ffv1_1(w(1,6),w(1,7),gc_3,mta,zero,w(1,10))
1898 CALL ffv1_0(w(1,5),w(1,10),w(1,8),gc_3,amp(2))
1899 CALL ffv2_5_3(w(1,4),w(1,3),gc_51,gc_58,mz,wz,w(1,11))
1901 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(3))
1903 CALL ffv2_4_0(w(1,5),w(1,10),w(1,11),gc_50,gc_59,amp(4))
1904 CALL ffv2_5_3(w(1,1),w(1,2),gc_51,gc_58,mz,wz,w(1,10))
1905 CALL ffv2_4_2(w(1,5),w(1,10),gc_50,gc_59,mta,zero,w(1,9))
1907 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(5))
1908 CALL ffv2_4_1(w(1,6),w(1,10),gc_50,gc_59,mta,zero,w(1,12))
1910 CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
1912 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(7))
1914 CALL ffv2_4_0(w(1,5),w(1,12),w(1,11),gc_50,gc_59,amp(8))
1915 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,12))
1916 CALL ffv1_1(w(1,3),w(1,7),gc_2,zero,zero,w(1,9))
1918 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_2,amp(9))
1919 CALL ffv1_2(w(1,4),w(1,7),gc_2,zero,zero,w(1,13))
1921 CALL ffv1_0(w(1,13),w(1,3),w(1,12),gc_2,amp(10))
1922 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
1924 CALL ffv2_5_0(w(1,4),w(1,9),w(1,7),gc_51,gc_58,amp(11))
1926 CALL ffv2_5_0(w(1,13),w(1,3),w(1,7),gc_51,gc_58,amp(12))
1927 CALL ffv1p0_3(w(1,1),w(1,2),gc_11,zero,zero,w(1,13))
1928 CALL ffv1_1(w(1,3),w(1,13),gc_11,zero,zero,w(1,9))
1930 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_2,amp(13))
1931 CALL ffv1_2(w(1,4),w(1,13),gc_11,zero,zero,w(1,14))
1933 CALL ffv1_0(w(1,14),w(1,3),w(1,12),gc_2,amp(14))
1935 CALL ffv2_5_0(w(1,4),w(1,9),w(1,7),gc_51,gc_58,amp(15))
1937 CALL ffv2_5_0(w(1,14),w(1,3),w(1,7),gc_51,gc_58,amp(16))
1938 CALL ffv2_5_1(w(1,3),w(1,10),gc_51,gc_58,zero,zero,w(1,14))
1940 CALL ffv1_0(w(1,4),w(1,14),w(1,12),gc_2,amp(17))
1941 CALL ffv2_5_2(w(1,4),w(1,10),gc_51,gc_58,zero,zero,w(1,9))
1943 CALL ffv1_0(w(1,9),w(1,3),w(1,12),gc_2,amp(18))
1945 CALL ffv2_5_0(w(1,4),w(1,14),w(1,7),gc_51,gc_58,amp(19))
1947 CALL ffv2_5_0(w(1,9),w(1,3),w(1,7),gc_51,gc_58,amp(20))
1948 CALL ffv1p0_3(w(1,1),w(1,3),gc_2,zero,zero,w(1,9))
1949 CALL ffv1p0_3(w(1,4),w(1,2),gc_2,zero,zero,w(1,14))
1950 CALL ffv1_2(w(1,5),w(1,9),gc_3,mta,zero,w(1,10))
1952 CALL ffv1_0(w(1,10),w(1,6),w(1,14),gc_3,amp(21))
1953 CALL ffv1_1(w(1,6),w(1,9),gc_3,mta,zero,w(1,13))
1955 CALL ffv1_0(w(1,5),w(1,13),w(1,14),gc_3,amp(22))
1956 CALL ffv2_5_3(w(1,4),w(1,2),gc_51,gc_58,mz,wz,w(1,15))
1958 CALL ffv2_4_0(w(1,10),w(1,6),w(1,15),gc_50,gc_59,amp(23))
1960 CALL ffv2_4_0(w(1,5),w(1,13),w(1,15),gc_50,gc_59,amp(24))
1961 CALL ffv2_5_3(w(1,1),w(1,3),gc_51,gc_58,mz,wz,w(1,13))
1962 CALL ffv2_4_2(w(1,5),w(1,13),gc_50,gc_59,mta,zero,w(1,10))
1964 CALL ffv1_0(w(1,10),w(1,6),w(1,14),gc_3,amp(25))
1965 CALL ffv2_4_1(w(1,6),w(1,13),gc_50,gc_59,mta,zero,w(1,16))
1967 CALL ffv1_0(w(1,5),w(1,16),w(1,14),gc_3,amp(26))
1969 CALL ffv2_4_0(w(1,10),w(1,6),w(1,15),gc_50,gc_59,amp(27))
1971 CALL ffv2_4_0(w(1,5),w(1,16),w(1,15),gc_50,gc_59,amp(28))
1972 CALL ffv1_1(w(1,2),w(1,9),gc_2,zero,zero,w(1,16))
1974 CALL ffv1_0(w(1,4),w(1,16),w(1,12),gc_2,amp(29))
1975 CALL ffv1_2(w(1,4),w(1,9),gc_2,zero,zero,w(1,5))
1977 CALL ffv1_0(w(1,5),w(1,2),w(1,12),gc_2,amp(30))
1979 CALL ffv2_5_0(w(1,4),w(1,16),w(1,7),gc_51,gc_58,amp(31))
1981 CALL ffv2_5_0(w(1,5),w(1,2),w(1,7),gc_51,gc_58,amp(32))
1982 CALL ffv1p0_3(w(1,1),w(1,3),gc_11,zero,zero,w(1,5))
1983 CALL ffv1_1(w(1,2),w(1,5),gc_11,zero,zero,w(1,16))
1985 CALL ffv1_0(w(1,4),w(1,16),w(1,12),gc_2,amp(33))
1986 CALL ffv1_2(w(1,4),w(1,5),gc_11,zero,zero,w(1,9))
1988 CALL ffv1_0(w(1,9),w(1,2),w(1,12),gc_2,amp(34))
1990 CALL ffv2_5_0(w(1,4),w(1,16),w(1,7),gc_51,gc_58,amp(35))
1992 CALL ffv2_5_0(w(1,9),w(1,2),w(1,7),gc_51,gc_58,amp(36))
1993 CALL ffv2_5_1(w(1,2),w(1,13),gc_51,gc_58,zero,zero,w(1,9))
1995 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_2,amp(37))
1996 CALL ffv2_5_2(w(1,4),w(1,13),gc_51,gc_58,zero,zero,w(1,16))
1998 CALL ffv1_0(w(1,16),w(1,2),w(1,12),gc_2,amp(38))
2000 CALL ffv2_5_0(w(1,4),w(1,9),w(1,7),gc_51,gc_58,amp(39))
2002 CALL ffv2_5_0(w(1,16),w(1,2),w(1,7),gc_51,gc_58,amp(40))
2003 CALL ffv1_2(w(1,1),w(1,14),gc_2,zero,zero,w(1,16))
2005 CALL ffv1_0(w(1,16),w(1,3),w(1,12),gc_2,amp(41))
2006 CALL ffv1_2(w(1,1),w(1,12),gc_2,zero,zero,w(1,9))
2008 CALL ffv1_0(w(1,9),w(1,3),w(1,14),gc_2,amp(42))
2010 CALL ffv2_5_0(w(1,16),w(1,3),w(1,7),gc_51,gc_58,amp(43))
2011 CALL ffv2_5_2(w(1,1),w(1,7),gc_51,gc_58,zero,zero,w(1,16))
2013 CALL ffv1_0(w(1,16),w(1,3),w(1,14),gc_2,amp(44))
2014 CALL ffv1p0_3(w(1,4),w(1,2),gc_11,zero,zero,w(1,14))
2015 CALL ffv1_2(w(1,1),w(1,14),gc_11,zero,zero,w(1,13))
2017 CALL ffv1_0(w(1,13),w(1,3),w(1,12),gc_2,amp(45))
2019 CALL ffv1_0(w(1,9),w(1,3),w(1,14),gc_11,amp(46))
2021 CALL ffv2_5_0(w(1,13),w(1,3),w(1,7),gc_51,gc_58,amp(47))
2023 CALL ffv1_0(w(1,16),w(1,3),w(1,14),gc_11,amp(48))
2024 CALL ffv2_5_2(w(1,1),w(1,15),gc_51,gc_58,zero,zero,w(1,14))
2026 CALL ffv1_0(w(1,14),w(1,3),w(1,12),gc_2,amp(49))
2028 CALL ffv2_5_0(w(1,9),w(1,3),w(1,15),gc_51,gc_58,amp(50))
2030 CALL ffv2_5_0(w(1,14),w(1,3),w(1,7),gc_51,gc_58,amp(51))
2032 CALL ffv2_5_0(w(1,16),w(1,3),w(1,15),gc_51,gc_58,amp(52))
2033 CALL ffv1_2(w(1,1),w(1,8),gc_2,zero,zero,w(1,15))
2035 CALL ffv1_0(w(1,15),w(1,2),w(1,12),gc_2,amp(53))
2037 CALL ffv1_0(w(1,9),w(1,2),w(1,8),gc_2,amp(54))
2039 CALL ffv2_5_0(w(1,15),w(1,2),w(1,7),gc_51,gc_58,amp(55))
2041 CALL ffv1_0(w(1,16),w(1,2),w(1,8),gc_2,amp(56))
2042 CALL ffv1p0_3(w(1,4),w(1,3),gc_11,zero,zero,w(1,8))
2043 CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,4))
2045 CALL ffv1_0(w(1,4),w(1,2),w(1,12),gc_2,amp(57))
2047 CALL ffv1_0(w(1,9),w(1,2),w(1,8),gc_11,amp(58))
2049 CALL ffv2_5_0(w(1,4),w(1,2),w(1,7),gc_51,gc_58,amp(59))
2051 CALL ffv1_0(w(1,16),w(1,2),w(1,8),gc_11,amp(60))
2052 CALL ffv2_5_2(w(1,1),w(1,11),gc_51,gc_58,zero,zero,w(1,8))
2054 CALL ffv1_0(w(1,8),w(1,2),w(1,12),gc_2,amp(61))
2056 CALL ffv2_5_0(w(1,9),w(1,2),w(1,11),gc_51,gc_58,amp(62))
2058 CALL ffv2_5_0(w(1,8),w(1,2),w(1,7),gc_51,gc_58,amp(63))
2060 CALL ffv2_5_0(w(1,16),w(1,2),w(1,11),gc_51,gc_58,amp(64))
2061 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
2062 $ -amp(9)-amp(10)-amp(11)-amp(12)+1./6.*amp(13)+1./6.*amp(14)
2063 $ +1./6.*amp(15)+1./6.*amp(16)-amp(17)-amp(18)-amp(19)-amp(20)
2064 $ +1./2.*amp(33)+1./2.*amp(34)+1./2.*amp(35)+1./2.*amp(36)
2065 $ +1./2.*amp(45)+1./2.*amp(46)+1./2.*amp(47)+1./2.*amp(48)
2066 $ -amp(53)-amp(54)-amp(55)-amp(56)+1./6.*amp(57)+1./6.*amp(58)
2067 $ +1./6.*amp(59)+1./6.*amp(60)-amp(61)-amp(62)-amp(63)-amp(64)
2068 jamp(2)=-1./2.*amp(13)-1./2.*amp(14)-1./2.*amp(15)-1./2.*amp(16)
2069 $ +amp(21)+amp(22)+amp(23)+amp(24)+amp(25)+amp(26)+amp(27)
2070 $ +amp(28)+amp(29)+amp(30)+amp(31)+amp(32)-1./6.*amp(33)
2071 $ -1./6.*amp(34)-1./6.*amp(35)-1./6.*amp(36)+amp(37)+amp(38)
2072 $ +amp(39)+amp(40)+amp(41)+amp(42)+amp(43)+amp(44)-1./6.*amp(45)
2073 $ -1./6.*amp(46)-1./6.*amp(47)-1./6.*amp(48)+amp(49)+amp(50)
2074 $ +amp(51)+amp(52)-1./2.*amp(57)-1./2.*amp(58)-1./2.*amp(59)
2081 ztemp = ztemp + cf(j,i)*jamp(j)
2083 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
2085 matrix_uux_uux_noh=matrix
2092 SUBROUTINE uux_ssx_noh(P,H1,H2,ANS)
2111 parameter(nexternal=6)
2113 parameter( ncomb=64)
2117 real*8 p(0:3,nexternal),ans
2122 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
2124 real*8 matrix_uux_ssx_noh
2125 INTEGER IHEL,IDEN, I
2126 INTEGER JC(NEXTERNAL)
2127 LOGICAL GOODHEL(NCOMB)
2129 DATA goodhel/ncomb*.false./
2130 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
2131 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
2132 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
2133 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
2134 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
2135 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
2136 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
2137 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
2138 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
2139 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
2140 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
2141 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
2142 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
2143 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
2144 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
2145 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
2146 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
2147 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
2148 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
2149 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
2150 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
2151 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
2152 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
2153 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
2154 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
2155 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
2156 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
2157 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
2158 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
2159 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
2160 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
2161 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
2162 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
2163 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
2164 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
2165 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
2166 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
2167 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
2168 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
2169 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
2170 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
2171 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
2172 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
2173 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
2174 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
2175 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
2176 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
2177 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
2178 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
2179 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
2180 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
2181 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
2182 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
2183 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
2184 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
2185 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
2186 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
2187 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
2188 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
2189 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
2190 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
2191 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
2192 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
2193 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
2203 t=matrix_uux_ssx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
2210 REAL*8 FUNCTION matrix_uux_ssx_noh(P,H1,H2,NHEL,IC)
2226 parameter(ngraphs=43)
2228 parameter(nexternal=6)
2229 INTEGER NWAVEFUNCS, NCOLOR
2230 parameter(nwavefuncs=14, ncolor=2)
2234 parameter(imag1=(0d0,1d0))
2238 real*8 p(0:3,nexternal)
2239 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
2245 real*8 denom(ncolor), cf(ncolor,ncolor)
2246 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
2247 COMPLEX*16 W(18,NWAVEFUNCS)
2248 COMPLEX*16 DUM0,DUM1
2249 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
2258 DATA (cf(i, 1),i= 1, 2) / 9, 3/
2261 DATA (cf(i, 2),i= 1, 2) / 3, 9/
2266 matrix_uux_ssx_noh=0.d0
2267 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
2268 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
2273 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
2274 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
2275 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
2276 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
2277 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
2278 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
2279 CALL ffv1p0_3(w(1,1),w(1,2),gc_2,zero,zero,w(1,7))
2280 CALL ffv1p0_3(w(1,4),w(1,3),gc_1,zero,zero,w(1,8))
2281 CALL ffv1_2(w(1,5),w(1,7),gc_3,mta,zero,w(1,9))
2283 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(1))
2284 CALL ffv1_1(w(1,6),w(1,7),gc_3,mta,zero,w(1,10))
2286 CALL ffv1_0(w(1,5),w(1,10),w(1,8),gc_3,amp(2))
2287 CALL ffv2_3_3(w(1,4),w(1,3),gc_50,gc_58,mz,wz,w(1,11))
2289 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(3))
2291 CALL ffv2_4_0(w(1,5),w(1,10),w(1,11),gc_50,gc_59,amp(4))
2292 CALL ffv2_5_3(w(1,1),w(1,2),gc_51,gc_58,mz,wz,w(1,10))
2293 CALL ffv2_4_2(w(1,5),w(1,10),gc_50,gc_59,mta,zero,w(1,9))
2295 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(5))
2296 CALL ffv2_4_1(w(1,6),w(1,10),gc_50,gc_59,mta,zero,w(1,12))
2298 CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
2300 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(7))
2302 CALL ffv2_4_0(w(1,5),w(1,12),w(1,11),gc_50,gc_59,amp(8))
2303 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,12))
2304 CALL ffv1_1(w(1,3),w(1,7),gc_1,zero,zero,w(1,9))
2306 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_1,amp(9))
2307 CALL ffv1_2(w(1,4),w(1,7),gc_1,zero,zero,w(1,13))
2309 CALL ffv1_0(w(1,13),w(1,3),w(1,12),gc_1,amp(10))
2310 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
2312 CALL ffv2_3_0(w(1,4),w(1,9),w(1,7),gc_50,gc_58,amp(11))
2314 CALL ffv2_3_0(w(1,13),w(1,3),w(1,7),gc_50,gc_58,amp(12))
2315 CALL ffv1p0_3(w(1,1),w(1,2),gc_11,zero,zero,w(1,13))
2316 CALL ffv1_1(w(1,3),w(1,13),gc_11,zero,zero,w(1,9))
2318 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_1,amp(13))
2319 CALL ffv1_2(w(1,4),w(1,13),gc_11,zero,zero,w(1,14))
2321 CALL ffv1_0(w(1,14),w(1,3),w(1,12),gc_1,amp(14))
2323 CALL ffv2_3_0(w(1,4),w(1,9),w(1,7),gc_50,gc_58,amp(15))
2325 CALL ffv2_3_0(w(1,14),w(1,3),w(1,7),gc_50,gc_58,amp(16))
2326 CALL ffv2_3_1(w(1,3),w(1,10),gc_50,gc_58,zero,zero,w(1,14))
2328 CALL ffv1_0(w(1,4),w(1,14),w(1,12),gc_1,amp(17))
2329 CALL ffv2_3_2(w(1,4),w(1,10),gc_50,gc_58,zero,zero,w(1,9))
2331 CALL ffv1_0(w(1,9),w(1,3),w(1,12),gc_1,amp(18))
2333 CALL ffv2_3_0(w(1,4),w(1,14),w(1,7),gc_50,gc_58,amp(19))
2335 CALL ffv2_3_0(w(1,9),w(1,3),w(1,7),gc_50,gc_58,amp(20))
2336 CALL ffv2_3(w(1,1),w(1,3),gc_101,mw,ww,w(1,9))
2337 CALL ffv2_3(w(1,4),w(1,2),gc_101,mw,ww,w(1,14))
2338 CALL ffv2_2(w(1,5),w(1,9),gc_108,zero,zero,w(1,10))
2340 CALL ffv2_0(w(1,10),w(1,6),w(1,14),gc_108,amp(21))
2342 CALL vvv1_0(w(1,12),w(1,9),w(1,14),gc_4,amp(22))
2344 CALL vvv1_0(w(1,9),w(1,14),w(1,7),gc_53,amp(23))
2345 CALL ffv2_1(w(1,2),w(1,9),gc_101,zero,zero,w(1,10))
2347 CALL ffv1_0(w(1,4),w(1,10),w(1,12),gc_1,amp(24))
2348 CALL ffv2_2(w(1,4),w(1,9),gc_101,zero,zero,w(1,6))
2350 CALL ffv1_0(w(1,6),w(1,2),w(1,12),gc_2,amp(25))
2352 CALL ffv2_3_0(w(1,4),w(1,10),w(1,7),gc_50,gc_58,amp(26))
2354 CALL ffv2_5_0(w(1,6),w(1,2),w(1,7),gc_51,gc_58,amp(27))
2355 CALL ffv2_2(w(1,1),w(1,14),gc_101,zero,zero,w(1,6))
2357 CALL ffv1_0(w(1,6),w(1,3),w(1,12),gc_1,amp(28))
2358 CALL ffv1_2(w(1,1),w(1,12),gc_2,zero,zero,w(1,10))
2360 CALL ffv2_0(w(1,10),w(1,3),w(1,14),gc_101,amp(29))
2362 CALL ffv2_3_0(w(1,6),w(1,3),w(1,7),gc_50,gc_58,amp(30))
2363 CALL ffv2_5_2(w(1,1),w(1,7),gc_51,gc_58,zero,zero,w(1,6))
2365 CALL ffv2_0(w(1,6),w(1,3),w(1,14),gc_101,amp(31))
2366 CALL ffv1_2(w(1,1),w(1,8),gc_2,zero,zero,w(1,14))
2368 CALL ffv1_0(w(1,14),w(1,2),w(1,12),gc_2,amp(32))
2370 CALL ffv1_0(w(1,10),w(1,2),w(1,8),gc_2,amp(33))
2372 CALL ffv2_5_0(w(1,14),w(1,2),w(1,7),gc_51,gc_58,amp(34))
2374 CALL ffv1_0(w(1,6),w(1,2),w(1,8),gc_2,amp(35))
2375 CALL ffv1p0_3(w(1,4),w(1,3),gc_11,zero,zero,w(1,8))
2376 CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,4))
2378 CALL ffv1_0(w(1,4),w(1,2),w(1,12),gc_2,amp(36))
2380 CALL ffv1_0(w(1,10),w(1,2),w(1,8),gc_11,amp(37))
2382 CALL ffv2_5_0(w(1,4),w(1,2),w(1,7),gc_51,gc_58,amp(38))
2384 CALL ffv1_0(w(1,6),w(1,2),w(1,8),gc_11,amp(39))
2385 CALL ffv2_5_2(w(1,1),w(1,11),gc_51,gc_58,zero,zero,w(1,8))
2387 CALL ffv1_0(w(1,8),w(1,2),w(1,12),gc_2,amp(40))
2389 CALL ffv2_5_0(w(1,10),w(1,2),w(1,11),gc_51,gc_58,amp(41))
2391 CALL ffv2_5_0(w(1,8),w(1,2),w(1,7),gc_51,gc_58,amp(42))
2393 CALL ffv2_5_0(w(1,6),w(1,2),w(1,11),gc_51,gc_58,amp(43))
2394 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
2395 $ -amp(9)-amp(10)-amp(11)-amp(12)+1./6.*amp(13)+1./6.*amp(14)
2396 $ +1./6.*amp(15)+1./6.*amp(16)-amp(17)-amp(18)-amp(19)-amp(20)
2397 $ -amp(32)-amp(33)-amp(34)-amp(35)+1./6.*amp(36)+1./6.*amp(37)
2398 $ +1./6.*amp(38)+1./6.*amp(39)-amp(40)-amp(41)-amp(42)-amp(43)
2399 jamp(2)=-1./2.*amp(13)-1./2.*amp(14)-1./2.*amp(15)-1./2.*amp(16)
2400 $ +amp(21)+amp(22)+amp(23)+amp(24)+amp(25)+amp(26)+amp(27)
2401 $ +amp(28)+amp(29)+amp(30)+amp(31)-1./2.*amp(36)-1./2.*amp(37)
2402 $ -1./2.*amp(38)-1./2.*amp(39)
2408 ztemp = ztemp + cf(j,i)*jamp(j)
2410 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
2412 matrix_uux_ssx_noh=matrix
2419 SUBROUTINE uux_ccx_noh(P,H1,H2,ANS)
2438 parameter(nexternal=6)
2440 parameter( ncomb=64)
2444 real*8 p(0:3,nexternal),ans
2449 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
2451 real*8 matrix_uux_ccx_noh
2452 INTEGER IHEL,IDEN, I
2453 INTEGER JC(NEXTERNAL)
2454 LOGICAL GOODHEL(NCOMB)
2456 DATA goodhel/ncomb*.false./
2457 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
2458 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
2459 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
2460 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
2461 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
2462 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
2463 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
2464 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
2465 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
2466 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
2467 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
2468 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
2469 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
2470 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
2471 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
2472 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
2473 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
2474 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
2475 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
2476 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
2477 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
2478 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
2479 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
2480 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
2481 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
2482 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
2483 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
2484 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
2485 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
2486 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
2487 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
2488 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
2489 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
2490 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
2491 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
2492 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
2493 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
2494 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
2495 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
2496 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
2497 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
2498 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
2499 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
2500 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
2501 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
2502 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
2503 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
2504 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
2505 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
2506 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
2507 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
2508 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
2509 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
2510 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
2511 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
2512 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
2513 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
2514 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
2515 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
2516 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
2517 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
2518 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
2519 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
2520 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
2530 t=matrix_uux_ccx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
2537 REAL*8 FUNCTION matrix_uux_ccx_noh(P,H1,H2,NHEL,IC)
2553 parameter(ngraphs=32)
2555 parameter(nexternal=6)
2556 INTEGER NWAVEFUNCS, NCOLOR
2557 parameter(nwavefuncs=13, ncolor=2)
2561 parameter(imag1=(0d0,1d0))
2565 real*8 p(0:3,nexternal)
2566 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
2572 real*8 denom(ncolor), cf(ncolor,ncolor)
2573 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
2574 COMPLEX*16 W(18,NWAVEFUNCS)
2575 COMPLEX*16 DUM0,DUM1
2576 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
2585 DATA (cf(i, 1),i= 1, 2) / 9, 3/
2588 DATA (cf(i, 2),i= 1, 2) / 3, 9/
2593 matrix_uux_ccx_noh=0.d0
2594 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
2595 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
2600 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
2601 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
2602 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
2603 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
2604 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
2605 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
2606 CALL ffv1p0_3(w(1,1),w(1,2),gc_2,zero,zero,w(1,7))
2607 CALL ffv1p0_3(w(1,4),w(1,3),gc_2,zero,zero,w(1,8))
2608 CALL ffv1_2(w(1,5),w(1,7),gc_3,mta,zero,w(1,9))
2610 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(1))
2611 CALL ffv1_1(w(1,6),w(1,7),gc_3,mta,zero,w(1,10))
2613 CALL ffv1_0(w(1,5),w(1,10),w(1,8),gc_3,amp(2))
2614 CALL ffv2_5_3(w(1,4),w(1,3),gc_51,gc_58,mz,wz,w(1,11))
2616 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(3))
2618 CALL ffv2_4_0(w(1,5),w(1,10),w(1,11),gc_50,gc_59,amp(4))
2619 CALL ffv2_5_3(w(1,1),w(1,2),gc_51,gc_58,mz,wz,w(1,10))
2620 CALL ffv2_4_2(w(1,5),w(1,10),gc_50,gc_59,mta,zero,w(1,9))
2622 CALL ffv1_0(w(1,9),w(1,6),w(1,8),gc_3,amp(5))
2623 CALL ffv2_4_1(w(1,6),w(1,10),gc_50,gc_59,mta,zero,w(1,12))
2625 CALL ffv1_0(w(1,5),w(1,12),w(1,8),gc_3,amp(6))
2627 CALL ffv2_4_0(w(1,9),w(1,6),w(1,11),gc_50,gc_59,amp(7))
2629 CALL ffv2_4_0(w(1,5),w(1,12),w(1,11),gc_50,gc_59,amp(8))
2630 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,12))
2631 CALL ffv1_1(w(1,3),w(1,7),gc_2,zero,zero,w(1,9))
2633 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_2,amp(9))
2634 CALL ffv1_2(w(1,4),w(1,7),gc_2,zero,zero,w(1,13))
2636 CALL ffv1_0(w(1,13),w(1,3),w(1,12),gc_2,amp(10))
2637 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
2639 CALL ffv2_5_0(w(1,4),w(1,9),w(1,7),gc_51,gc_58,amp(11))
2641 CALL ffv2_5_0(w(1,13),w(1,3),w(1,7),gc_51,gc_58,amp(12))
2642 CALL ffv1p0_3(w(1,1),w(1,2),gc_11,zero,zero,w(1,13))
2643 CALL ffv1_1(w(1,3),w(1,13),gc_11,zero,zero,w(1,9))
2645 CALL ffv1_0(w(1,4),w(1,9),w(1,12),gc_2,amp(13))
2646 CALL ffv1_2(w(1,4),w(1,13),gc_11,zero,zero,w(1,6))
2648 CALL ffv1_0(w(1,6),w(1,3),w(1,12),gc_2,amp(14))
2650 CALL ffv2_5_0(w(1,4),w(1,9),w(1,7),gc_51,gc_58,amp(15))
2652 CALL ffv2_5_0(w(1,6),w(1,3),w(1,7),gc_51,gc_58,amp(16))
2653 CALL ffv2_5_1(w(1,3),w(1,10),gc_51,gc_58,zero,zero,w(1,6))
2655 CALL ffv1_0(w(1,4),w(1,6),w(1,12),gc_2,amp(17))
2656 CALL ffv2_5_2(w(1,4),w(1,10),gc_51,gc_58,zero,zero,w(1,9))
2658 CALL ffv1_0(w(1,9),w(1,3),w(1,12),gc_2,amp(18))
2660 CALL ffv2_5_0(w(1,4),w(1,6),w(1,7),gc_51,gc_58,amp(19))
2662 CALL ffv2_5_0(w(1,9),w(1,3),w(1,7),gc_51,gc_58,amp(20))
2663 CALL ffv1_2(w(1,1),w(1,8),gc_2,zero,zero,w(1,9))
2665 CALL ffv1_0(w(1,9),w(1,2),w(1,12),gc_2,amp(21))
2666 CALL ffv1_2(w(1,1),w(1,12),gc_2,zero,zero,w(1,6))
2668 CALL ffv1_0(w(1,6),w(1,2),w(1,8),gc_2,amp(22))
2670 CALL ffv2_5_0(w(1,9),w(1,2),w(1,7),gc_51,gc_58,amp(23))
2671 CALL ffv2_5_2(w(1,1),w(1,7),gc_51,gc_58,zero,zero,w(1,9))
2673 CALL ffv1_0(w(1,9),w(1,2),w(1,8),gc_2,amp(24))
2674 CALL ffv1p0_3(w(1,4),w(1,3),gc_11,zero,zero,w(1,8))
2675 CALL ffv1_2(w(1,1),w(1,8),gc_11,zero,zero,w(1,4))
2677 CALL ffv1_0(w(1,4),w(1,2),w(1,12),gc_2,amp(25))
2679 CALL ffv1_0(w(1,6),w(1,2),w(1,8),gc_11,amp(26))
2681 CALL ffv2_5_0(w(1,4),w(1,2),w(1,7),gc_51,gc_58,amp(27))
2683 CALL ffv1_0(w(1,9),w(1,2),w(1,8),gc_11,amp(28))
2684 CALL ffv2_5_2(w(1,1),w(1,11),gc_51,gc_58,zero,zero,w(1,8))
2686 CALL ffv1_0(w(1,8),w(1,2),w(1,12),gc_2,amp(29))
2688 CALL ffv2_5_0(w(1,6),w(1,2),w(1,11),gc_51,gc_58,amp(30))
2690 CALL ffv2_5_0(w(1,8),w(1,2),w(1,7),gc_51,gc_58,amp(31))
2692 CALL ffv2_5_0(w(1,9),w(1,2),w(1,11),gc_51,gc_58,amp(32))
2693 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
2694 $ -amp(9)-amp(10)-amp(11)-amp(12)+1./6.*amp(13)+1./6.*amp(14)
2695 $ +1./6.*amp(15)+1./6.*amp(16)-amp(17)-amp(18)-amp(19)-amp(20)
2696 $ -amp(21)-amp(22)-amp(23)-amp(24)+1./6.*amp(25)+1./6.*amp(26)
2697 $ +1./6.*amp(27)+1./6.*amp(28)-amp(29)-amp(30)-amp(31)-amp(32)
2698 jamp(2)=+1./2.*(-amp(13)-amp(14)-amp(15)-amp(16)-amp(25)-amp(26)
2705 ztemp = ztemp + cf(j,i)*jamp(j)
2707 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
2709 matrix_uux_ccx_noh=matrix
2716 SUBROUTINE uux_dsx_noh(P,H1,H2,ANS)
2735 parameter(nexternal=6)
2737 parameter( ncomb=64)
2741 real*8 p(0:3,nexternal),ans
2746 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
2748 real*8 matrix_uux_dsx_noh
2749 INTEGER IHEL,IDEN, I
2750 INTEGER JC(NEXTERNAL)
2751 LOGICAL GOODHEL(NCOMB)
2753 DATA goodhel/ncomb*.false./
2754 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
2755 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
2756 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
2757 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
2758 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
2759 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
2760 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
2761 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
2762 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
2763 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
2764 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
2765 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
2766 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
2767 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
2768 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
2769 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
2770 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
2771 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
2772 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
2773 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
2774 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
2775 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
2776 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
2777 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
2778 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
2779 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
2780 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
2781 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
2782 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
2783 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
2784 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
2785 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
2786 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
2787 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
2788 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
2789 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
2790 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
2791 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
2792 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
2793 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
2794 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
2795 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
2796 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
2797 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
2798 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
2799 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
2800 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
2801 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
2802 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
2803 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
2804 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
2805 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
2806 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
2807 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
2808 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
2809 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
2810 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
2811 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
2812 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
2813 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
2814 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
2815 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
2816 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
2817 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
2827 t=matrix_uux_dsx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
2834 REAL*8 FUNCTION matrix_uux_dsx_noh(P,H1,H2,NHEL,IC)
2850 parameter(ngraphs=11)
2852 parameter(nexternal=6)
2853 INTEGER NWAVEFUNCS, NCOLOR
2854 parameter(nwavefuncs=10, ncolor=1)
2858 parameter(imag1=(0d0,1d0))
2862 real*8 p(0:3,nexternal)
2863 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
2869 real*8 denom(ncolor), cf(ncolor,ncolor)
2870 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
2871 COMPLEX*16 W(18,NWAVEFUNCS)
2872 COMPLEX*16 DUM0,DUM1
2873 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
2882 DATA (cf(i, 1),i= 1, 1) / 9/
2887 matrix_uux_dsx_noh=0.d0
2888 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
2889 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
2894 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
2895 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
2896 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
2897 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
2898 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
2899 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
2900 CALL ffv2_3(w(1,1),w(1,3),gc_100,mw,ww,w(1,7))
2901 CALL ffv2_3(w(1,4),w(1,2),gc_101,mw,ww,w(1,8))
2902 CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
2904 CALL ffv2_0(w(1,9),w(1,6),w(1,8),gc_108,amp(1))
2905 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
2907 CALL vvv1_0(w(1,9),w(1,7),w(1,8),gc_4,amp(2))
2908 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
2910 CALL vvv1_0(w(1,7),w(1,8),w(1,10),gc_53,amp(3))
2911 CALL ffv2_1(w(1,2),w(1,7),gc_101,zero,zero,w(1,6))
2913 CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_1,amp(4))
2914 CALL ffv2_2(w(1,4),w(1,7),gc_101,zero,zero,w(1,5))
2916 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_2,amp(5))
2918 CALL ffv2_3_0(w(1,4),w(1,6),w(1,10),gc_50,gc_58,amp(6))
2920 CALL ffv2_5_0(w(1,5),w(1,2),w(1,10),gc_51,gc_58,amp(7))
2921 CALL ffv2_2(w(1,1),w(1,8),gc_100,zero,zero,w(1,5))
2923 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_1,amp(8))
2924 CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,2))
2926 CALL ffv2_0(w(1,2),w(1,3),w(1,8),gc_100,amp(9))
2928 CALL ffv2_3_0(w(1,5),w(1,3),w(1,10),gc_50,gc_58,amp(10))
2929 CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,5))
2931 CALL ffv2_0(w(1,5),w(1,3),w(1,8),gc_100,amp(11))
2932 jamp(1)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
2933 $ +amp(9)+amp(10)+amp(11)
2939 ztemp = ztemp + cf(j,i)*jamp(j)
2941 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
2943 matrix_uux_dsx_noh=matrix
2950 SUBROUTINE uux_sdx_noh(P,H1,H2,ANS)
2969 parameter(nexternal=6)
2971 parameter( ncomb=64)
2975 real*8 p(0:3,nexternal),ans
2980 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
2982 real*8 matrix_uux_sdx_noh
2983 INTEGER IHEL,IDEN, I
2984 INTEGER JC(NEXTERNAL)
2985 LOGICAL GOODHEL(NCOMB)
2987 DATA goodhel/ncomb*.false./
2988 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
2989 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
2990 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
2991 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
2992 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
2993 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
2994 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
2995 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
2996 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
2997 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
2998 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
2999 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
3000 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
3001 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
3002 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
3003 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
3004 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
3005 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
3006 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
3007 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
3008 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
3009 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
3010 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
3011 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
3012 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
3013 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
3014 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
3015 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
3016 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
3017 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
3018 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
3019 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
3020 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
3021 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
3022 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
3023 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
3024 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
3025 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
3026 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
3027 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
3028 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
3029 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
3030 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
3031 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
3032 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
3033 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
3034 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
3035 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
3036 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
3037 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
3038 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
3039 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
3040 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
3041 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
3042 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
3043 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
3044 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
3045 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
3046 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
3047 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
3048 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
3049 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
3050 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
3051 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
3061 t=matrix_uux_sdx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
3068 REAL*8 FUNCTION matrix_uux_sdx_noh(P,H1,H2,NHEL,IC)
3084 parameter(ngraphs=11)
3086 parameter(nexternal=6)
3087 INTEGER NWAVEFUNCS, NCOLOR
3088 parameter(nwavefuncs=10, ncolor=1)
3092 parameter(imag1=(0d0,1d0))
3096 real*8 p(0:3,nexternal)
3097 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
3103 real*8 denom(ncolor), cf(ncolor,ncolor)
3104 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
3105 COMPLEX*16 W(18,NWAVEFUNCS)
3106 COMPLEX*16 DUM0,DUM1
3107 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
3116 DATA (cf(i, 1),i= 1, 1) / 9/
3121 matrix_uux_sdx_noh=0.d0
3122 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
3123 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
3128 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
3129 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
3130 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
3131 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
3132 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
3133 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
3134 CALL ffv2_3(w(1,1),w(1,3),gc_101,mw,ww,w(1,7))
3135 CALL ffv2_3(w(1,4),w(1,2),gc_100,mw,ww,w(1,8))
3136 CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
3138 CALL ffv2_0(w(1,9),w(1,6),w(1,8),gc_108,amp(1))
3139 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
3141 CALL vvv1_0(w(1,9),w(1,7),w(1,8),gc_4,amp(2))
3142 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
3144 CALL vvv1_0(w(1,7),w(1,8),w(1,10),gc_53,amp(3))
3145 CALL ffv2_1(w(1,2),w(1,7),gc_100,zero,zero,w(1,6))
3147 CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_1,amp(4))
3148 CALL ffv2_2(w(1,4),w(1,7),gc_100,zero,zero,w(1,5))
3150 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_2,amp(5))
3152 CALL ffv2_3_0(w(1,4),w(1,6),w(1,10),gc_50,gc_58,amp(6))
3154 CALL ffv2_5_0(w(1,5),w(1,2),w(1,10),gc_51,gc_58,amp(7))
3155 CALL ffv2_2(w(1,1),w(1,8),gc_101,zero,zero,w(1,5))
3157 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_1,amp(8))
3158 CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,2))
3160 CALL ffv2_0(w(1,2),w(1,3),w(1,8),gc_101,amp(9))
3162 CALL ffv2_3_0(w(1,5),w(1,3),w(1,10),gc_50,gc_58,amp(10))
3163 CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,5))
3165 CALL ffv2_0(w(1,5),w(1,3),w(1,8),gc_101,amp(11))
3166 jamp(1)=+amp(1)+amp(2)+amp(3)+amp(4)+amp(5)+amp(6)+amp(7)+amp(8)
3167 $ +amp(9)+amp(10)+amp(11)
3173 ztemp = ztemp + cf(j,i)*jamp(j)
3175 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
3177 matrix_uux_sdx_noh=matrix
3189 SUBROUTINE uux_gg_noh(P,H1,H2,ANS)
3209 parameter(nexternal=6)
3211 parameter( ncomb=64)
3215 real*8 p(0:3,nexternal),ans
3220 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
3222 real*8 matrix_uux_gg_noh
3223 INTEGER IHEL,IDEN, I
3224 INTEGER JC(NEXTERNAL)
3225 LOGICAL GOODHEL(NCOMB)
3227 DATA goodhel/ncomb*.false./
3228 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
3229 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
3230 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
3231 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
3232 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
3233 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
3234 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
3235 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
3236 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
3237 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
3238 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
3239 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
3240 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
3241 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
3242 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
3243 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
3244 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
3245 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
3246 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
3247 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
3248 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
3249 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
3250 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
3251 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
3252 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
3253 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
3254 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
3255 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
3256 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
3257 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
3258 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
3259 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
3260 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
3261 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
3262 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
3263 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
3264 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
3265 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
3266 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
3267 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
3268 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
3269 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
3270 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
3271 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
3272 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
3273 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
3274 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
3275 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
3276 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
3277 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
3278 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
3279 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
3280 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
3281 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
3282 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
3283 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
3284 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
3285 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
3286 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
3287 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
3288 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
3289 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
3290 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
3291 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
3302 t=matrix_uux_gg_noh(p ,h1,h2,nhel(1,ihel),jc(1))
3309 REAL*8 FUNCTION matrix_uux_gg_noh(P,H1,H2,NHEL,IC)
3326 parameter(ngraphs=16)
3328 parameter(nexternal=6)
3329 INTEGER NWAVEFUNCS, NCOLOR
3330 parameter(nwavefuncs=10, ncolor=2)
3334 parameter(imag1=(0d0,1d0))
3338 real*8 p(0:3,nexternal)
3339 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
3345 real*8 denom(ncolor), cf(ncolor,ncolor)
3346 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
3347 COMPLEX*16 W(18,NWAVEFUNCS)
3348 COMPLEX*16 DUM0,DUM1
3349 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
3358 DATA (cf(i, 1),i= 1, 2) / 16, -2/
3361 DATA (cf(i, 2),i= 1, 2) / -2, 16/
3366 matrix_uux_gg_noh=0.d0
3367 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
3368 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
3373 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
3374 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
3375 CALL vxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
3376 CALL vxxxxx(p(0,4),zero,nhel(4),+1*ic(4),w(1,4))
3377 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
3378 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
3379 CALL vvv1p0_1(w(1,3),w(1,4),gc_10,zero,zero,w(1,7))
3380 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,8))
3381 CALL ffv1_1(w(1,2),w(1,7),gc_11,zero,zero,w(1,9))
3383 CALL ffv1_0(w(1,1),w(1,9),w(1,8),gc_2,amp(1))
3384 CALL ffv1_2(w(1,1),w(1,7),gc_11,zero,zero,w(1,10))
3386 CALL ffv1_0(w(1,10),w(1,2),w(1,8),gc_2,amp(2))
3387 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,7))
3389 CALL ffv2_5_0(w(1,1),w(1,9),w(1,7),gc_51,gc_58,amp(3))
3391 CALL ffv2_5_0(w(1,10),w(1,2),w(1,7),gc_51,gc_58,amp(4))
3392 CALL ffv1_1(w(1,2),w(1,3),gc_11,zero,zero,w(1,10))
3393 CALL ffv1_2(w(1,1),w(1,4),gc_11,zero,zero,w(1,9))
3395 CALL ffv1_0(w(1,9),w(1,10),w(1,8),gc_2,amp(5))
3397 CALL ffv2_5_0(w(1,9),w(1,10),w(1,7),gc_51,gc_58,amp(6))
3398 CALL ffv1_1(w(1,10),w(1,4),gc_11,zero,zero,w(1,6))
3400 CALL ffv1_0(w(1,1),w(1,6),w(1,8),gc_2,amp(7))
3402 CALL ffv2_5_0(w(1,1),w(1,6),w(1,7),gc_51,gc_58,amp(8))
3403 CALL ffv1_2(w(1,1),w(1,3),gc_11,zero,zero,w(1,6))
3404 CALL ffv1_1(w(1,2),w(1,4),gc_11,zero,zero,w(1,10))
3406 CALL ffv1_0(w(1,6),w(1,10),w(1,8),gc_2,amp(9))
3408 CALL ffv2_5_0(w(1,6),w(1,10),w(1,7),gc_51,gc_58,amp(10))
3409 CALL ffv1_2(w(1,6),w(1,4),gc_11,zero,zero,w(1,5))
3411 CALL ffv1_0(w(1,5),w(1,2),w(1,8),gc_2,amp(11))
3413 CALL ffv2_5_0(w(1,5),w(1,2),w(1,7),gc_51,gc_58,amp(12))
3414 CALL ffv1_1(w(1,10),w(1,3),gc_11,zero,zero,w(1,5))
3416 CALL ffv1_0(w(1,1),w(1,5),w(1,8),gc_2,amp(13))
3418 CALL ffv2_5_0(w(1,1),w(1,5),w(1,7),gc_51,gc_58,amp(14))
3419 CALL ffv1_2(w(1,9),w(1,3),gc_11,zero,zero,w(1,5))
3421 CALL ffv1_0(w(1,5),w(1,2),w(1,8),gc_2,amp(15))
3423 CALL ffv2_5_0(w(1,5),w(1,2),w(1,7),gc_51,gc_58,amp(16))
3424 jamp(1)=-imag1*amp(1)-imag1*amp(2)-imag1*amp(3)-imag1*amp(4)
3425 $ +amp(5)+amp(6)+amp(7)+amp(8)+amp(15)+amp(16)
3426 jamp(2)=+imag1*amp(1)+imag1*amp(2)+imag1*amp(3)+imag1*amp(4)
3427 $ +amp(9)+amp(10)+amp(11)+amp(12)+amp(13)+amp(14)
3433 ztemp = ztemp + cf(j,i)*jamp(j)
3435 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
3437 matrix_uux_gg_noh=matrix