3 SUBROUTINE udx(P,I3,I4,H1,H2,KEY,ANS)
6 INTEGER I3,I4, H1,H2, KEY
7 real*8 p(0:3,6), ans, ans1,ans2,ans3,ans4
25 IF(i3.EQ.2 .AND. i4.EQ.-1)
CALL udx_udx_h(p,h1,h2,ans)
26 IF(i3.EQ.2 .AND. i4.EQ.-3)
CALL udx_usx_h(p,h1,h2,ans)
27 IF(i3.EQ.4 .AND. i4.EQ.-1)
CAll udx_cdx_h(p,h1,h2,ans)
28 IF(i3.EQ.4 .AND. i4.EQ.-3)
CALL udx_csx_h(p,h1,h2,ans)
29 IF(i3.EQ.-1 .AND. i4.EQ.2)
CALL udx_udx_h(pp,h1,h2,ans)
30 IF(i3.EQ.-3 .AND. i4.EQ.2)
CALL udx_usx_h(pp,h1,h2,ans)
31 IF(i3.EQ.-1 .AND. i4.EQ.4)
CAll udx_cdx_h(pp,h1,h2,ans)
32 IF(i3.EQ.-3 .AND. i4.EQ.4)
CALL udx_csx_h(pp,h1,h2,ans)
33 IF(i3.EQ.0 .AND. i4.EQ.0)
THEN
34 CALL udx_udx_h(p,h1,h2,ans1)
35 CALL udx_usx_h(p,h1,h2,ans2)
36 CALL udx_cdx_h(p,h1,h2,ans3)
37 CALL udx_csx_h(p,h1,h2,ans4)
38 ans=ans1+ans2+ans3+ans4
39 CALL udx_udx_h(pp,h1,h2,ans1)
40 CALL udx_usx_h(pp,h1,h2,ans2)
41 CALL udx_cdx_h(pp,h1,h2,ans3)
42 CALL udx_csx_h(pp,h1,h2,ans4)
43 ans=ans1+ans2+ans3+ans4+ans
45 ELSE IF(key.EQ.0)
THEN
46 IF(i3.EQ.2 .AND. i4.EQ.-1)
CALL udx_udx_noh(p,h1,h2,ans)
47 IF(i3.EQ.2 .AND. i4.EQ.-3)
CALL udx_usx_noh(p,h1,h2,ans)
48 IF(i3.EQ.4 .AND. i4.EQ.-1)
CAll udx_cdx_noh(p,h1,h2,ans)
49 IF(i3.EQ.4 .AND. i4.EQ.-3)
CALL udx_csx_noh(p,h1,h2,ans)
51 IF(i3.EQ.-1 .AND. i4.EQ.2)
CALL udx_udx_noh(pp,h1,h2,ans)
52 IF(i3.EQ.-3 .AND. i4.EQ.2)
CALL udx_usx_noh(pp,h1,h2,ans)
53 IF(i3.EQ.-1 .AND. i4.EQ.4)
CAll udx_cdx_noh(pp,h1,h2,ans)
54 IF(i3.EQ.-3 .AND. i4.EQ.4)
CALL udx_csx_noh(pp,h1,h2,ans)
55 IF(i3.EQ.0 .AND. i4.EQ.0)
THEN
56 CALL udx_udx_noh(p,h1,h2,ans1)
57 CALL udx_usx_noh(p,h1,h2,ans2)
58 CALL udx_cdx_noh(p,h1,h2,ans3)
59 CALL udx_csx_noh(p,h1,h2,ans4)
60 ans=ans1+ans2+ans3+ans4
61 CALL udx_udx_noh(pp,h1,h2,ans1)
62 CALL udx_usx_noh(pp,h1,h2,ans2)
63 CALL udx_cdx_noh(pp,h1,h2,ans3)
64 CALL udx_csx_noh(pp,h1,h2,ans4)
65 ans=ans1+ans2+ans3+ans4+ans
68 WRITE(*,*)
'NOT FINISHED'
78 SUBROUTINE udx_udx_h(P,H1,H2,ANS)
98 parameter(nexternal=6)
104 real*8 p(0:3,nexternal),ans
109 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
111 real*8 matrix_udx_udx_h
113 INTEGER JC(NEXTERNAL)
114 LOGICAL GOODHEL(NCOMB)
116 DATA goodhel/ncomb*.false./
117 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
118 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
119 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
120 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
121 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
122 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
123 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
124 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
125 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
126 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
127 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
128 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
129 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
130 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
131 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
132 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
133 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
134 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
135 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
136 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
137 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
138 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
139 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
140 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
141 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
142 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
143 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
144 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
145 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
146 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
147 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
148 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
149 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
150 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
151 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
152 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
153 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
154 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
155 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
156 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
157 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
158 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
159 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
160 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
161 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
162 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
163 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
164 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
165 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
166 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
167 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
168 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
169 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
170 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
171 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
172 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
173 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
174 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
175 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
176 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
177 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
178 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
179 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
180 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
190 t=matrix_udx_udx_h(p ,h1,h2,nhel(1,ihel),jc(1))
197 REAL*8 FUNCTION matrix_udx_udx_h(P,H1,H2,NHEL,IC)
216 parameter(nexternal=6)
217 INTEGER NWAVEFUNCS, NCOLOR
218 parameter(nwavefuncs=7, ncolor=2)
222 parameter(imag1=(0d0,1d0))
226 real*8 p(0:3,nexternal)
227 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
233 real*8 denom(ncolor), cf(ncolor,ncolor)
234 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
235 COMPLEX*16 W(18,NWAVEFUNCS)
237 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
246 DATA (cf(i, 1),i= 1, 2) / 9, 3/
249 DATA (cf(i, 2),i= 1, 2) / 3, 9/
254 matrix_udx_udx_h=0.d0
255 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
256 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
261 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
262 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
263 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
264 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
265 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
266 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
267 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
268 CALL ffv2_3(w(1,1),w(1,2),gc_100,mw,ww,w(1,6))
269 CALL ffv2_3(w(1,4),w(1,3),gc_100,mw,ww,w(1,5))
271 CALL vvs1_0(w(1,6),w(1,5),w(1,7),gc_72,amp(1))
272 CALL ffv2_5_3(w(1,1),w(1,3),gc_51,gc_58,mz,wz,w(1,5))
273 CALL ffv2_3_3(w(1,4),w(1,2),gc_50,gc_58,mz,wz,w(1,3))
275 CALL vvs1_0(w(1,5),w(1,3),w(1,7),gc_81,amp(2))
283 ztemp = ztemp + cf(j,i)*jamp(j)
285 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
287 matrix_udx_udx_h=matrix
294 SUBROUTINE udx_usx_h(P,H1,H2,ANS)
314 parameter(nexternal=6)
320 real*8 p(0:3,nexternal),ans
325 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
327 real*8 matrix_udx_usx_h
329 INTEGER JC(NEXTERNAL)
330 LOGICAL GOODHEL(NCOMB)
332 DATA goodhel/ncomb*.false./
333 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
334 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
335 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
336 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
337 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
338 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
339 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
340 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
341 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
342 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
343 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
344 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
345 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
346 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
347 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
348 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
349 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
350 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
351 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
352 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
353 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
354 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
355 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
356 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
357 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
358 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
359 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
360 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
361 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
362 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
363 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
364 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
365 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
366 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
367 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
368 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
369 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
370 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
371 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
372 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
373 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
374 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
375 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
376 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
377 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
378 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
379 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
380 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
381 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
382 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
383 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
384 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
385 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
386 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
387 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
388 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
389 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
390 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
391 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
392 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
393 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
394 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
395 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
396 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
406 t=matrix_udx_usx_h(p ,h1,h2,nhel(1,ihel),jc(1))
413 REAL*8 FUNCTION matrix_udx_usx_h(P,H1,H2,NHEL,IC)
432 parameter(nexternal=6)
433 INTEGER NWAVEFUNCS, NCOLOR
434 parameter(nwavefuncs=7, ncolor=1)
438 parameter(imag1=(0d0,1d0))
442 real*8 p(0:3,nexternal)
443 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
449 real*8 denom(ncolor), cf(ncolor,ncolor)
450 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
451 COMPLEX*16 W(18,NWAVEFUNCS)
453 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
462 DATA (cf(i, 1),i= 1, 1) / 9/
467 matrix_udx_usx_h=0.d0
468 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
469 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
475 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
476 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
477 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
478 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
479 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
480 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
481 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
482 CALL ffv2_3(w(1,1),w(1,2),gc_100,mw,ww,w(1,6))
483 CALL ffv2_3(w(1,4),w(1,3),gc_101,mw,ww,w(1,2))
485 CALL vvs1_0(w(1,6),w(1,2),w(1,7),gc_72,amp(1))
492 ztemp = ztemp + cf(j,i)*jamp(j)
494 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
496 matrix_udx_usx_h=matrix
503 SUBROUTINE udx_cdx_h(P,H1,H2,ANS)
523 parameter(nexternal=6)
529 real*8 p(0:3,nexternal),ans
534 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
536 real*8 matrix_udx_cdx_h
538 INTEGER JC(NEXTERNAL)
539 LOGICAL GOODHEL(NCOMB)
541 DATA goodhel/ncomb*.false./
542 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
543 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
544 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
545 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
546 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
547 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
548 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
549 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
550 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
551 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
552 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
553 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
554 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
555 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
556 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
557 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
558 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
559 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
560 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
561 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
562 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
563 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
564 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
565 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
566 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
567 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
568 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
569 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
570 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
571 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
572 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
573 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
574 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
575 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
576 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
577 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
578 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
579 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
580 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
581 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
582 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
583 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
584 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
585 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
586 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
587 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
588 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
589 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
590 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
591 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
592 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
593 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
594 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
595 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
596 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
597 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
598 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
599 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
600 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
601 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
602 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
603 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
604 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
605 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
615 t=matrix_udx_cdx_h(p ,h1,h2,nhel(1,ihel),jc(1))
622 REAL*8 FUNCTION matrix_udx_cdx_h(P,H1,H2,NHEL,IC)
641 parameter(nexternal=6)
642 INTEGER NWAVEFUNCS, NCOLOR
643 parameter(nwavefuncs=7, ncolor=1)
647 parameter(imag1=(0d0,1d0))
651 real*8 p(0:3,nexternal)
652 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
658 real*8 denom(ncolor), cf(ncolor,ncolor)
659 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
660 COMPLEX*16 W(18,NWAVEFUNCS)
662 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
671 DATA (cf(i, 1),i= 1, 1) / 9/
676 matrix_udx_cdx_h=0.d0
677 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
678 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
683 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
684 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
685 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
686 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
687 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
688 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
689 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
690 CALL ffv2_3(w(1,1),w(1,2),gc_100,mw,ww,w(1,6))
691 CALL ffv2_3(w(1,4),w(1,3),gc_44,mw,ww,w(1,2))
693 CALL vvs1_0(w(1,6),w(1,2),w(1,7),gc_72,amp(1))
700 ztemp = ztemp + cf(j,i)*jamp(j)
702 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
704 matrix_udx_cdx_h=matrix
711 SUBROUTINE udx_csx_h(P,H1,H2,ANS)
731 parameter(nexternal=6)
737 real*8 p(0:3,nexternal),ans
742 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
744 real*8 matrix_udx_csx_h
746 INTEGER JC(NEXTERNAL)
747 LOGICAL GOODHEL(NCOMB)
749 DATA goodhel/ncomb*.false./
750 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
751 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
752 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
753 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
754 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
755 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
756 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
757 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
758 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
759 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
760 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
761 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
762 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
763 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
764 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
765 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
766 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
767 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
768 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
769 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
770 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
771 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
772 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
773 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
774 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
775 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
776 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
777 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
778 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
779 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
780 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
781 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
782 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
783 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
784 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
785 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
786 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
787 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
788 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
789 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
790 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
791 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
792 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
793 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
794 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
795 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
796 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
797 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
798 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
799 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
800 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
801 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
802 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
803 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
804 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
805 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
806 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
807 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
808 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
809 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
810 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
811 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
812 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
813 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
823 t=matrix_udx_csx_h(p ,h1,h2,nhel(1,ihel),jc(1))
830 REAL*8 FUNCTION matrix_udx_csx_h(P,H1,H2,NHEL,IC)
849 parameter(nexternal=6)
850 INTEGER NWAVEFUNCS, NCOLOR
851 parameter(nwavefuncs=7, ncolor=1)
855 parameter(imag1=(0d0,1d0))
859 real*8 p(0:3,nexternal)
860 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
866 real*8 denom(ncolor), cf(ncolor,ncolor)
867 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
868 COMPLEX*16 W(18,NWAVEFUNCS)
870 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
879 DATA (cf(i, 1),i= 1, 1) / 9/
884 matrix_udx_csx_h=0.d0
885 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
886 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
891 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
892 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
893 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
894 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
895 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
896 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
897 CALL ffs4_3(w(1,5),w(1,6),gc_99,mh,wh,w(1,7))
898 CALL ffv2_3(w(1,1),w(1,2),gc_100,mw,ww,w(1,6))
899 CALL ffv2_3(w(1,4),w(1,3),gc_100,mw,ww,w(1,2))
901 CALL vvs1_0(w(1,6),w(1,2),w(1,7),gc_72,amp(1))
908 ztemp = ztemp + cf(j,i)*jamp(j)
910 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
912 matrix_udx_csx_h=matrix
921 SUBROUTINE udx_udx_noh(P,H1,H2,ANS)
940 parameter(nexternal=6)
946 real*8 p(0:3,nexternal),ans
951 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
953 real*8 matrix_udx_udx_noh
955 INTEGER JC(NEXTERNAL)
956 LOGICAL GOODHEL(NCOMB)
958 DATA goodhel/ncomb*.false./
959 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
960 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
961 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
962 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
963 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
964 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
965 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
966 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
967 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
968 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
969 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
970 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
971 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
972 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
973 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
974 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
975 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
976 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
977 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
978 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
979 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
980 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
981 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
982 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
983 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
984 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
985 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
986 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
987 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
988 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
989 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
990 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
991 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
992 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
993 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
994 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
995 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
996 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
997 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
998 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
999 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1000 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1001 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1002 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1003 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1004 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1005 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1006 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1007 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1008 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1009 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1010 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1011 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1012 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1013 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1014 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1015 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1016 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1017 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1018 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1019 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1020 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1021 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1022 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1032 t=matrix_udx_udx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1039 REAL*8 FUNCTION matrix_udx_udx_noh(P,H1,H2,NHEL,IC)
1055 parameter(ngraphs=43)
1057 parameter(nexternal=6)
1058 INTEGER NWAVEFUNCS, NCOLOR
1059 parameter(nwavefuncs=15, ncolor=2)
1063 parameter(imag1=(0d0,1d0))
1067 real*8 p(0:3,nexternal)
1068 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1074 real*8 denom(ncolor), cf(ncolor,ncolor)
1075 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1076 COMPLEX*16 W(18,NWAVEFUNCS)
1077 COMPLEX*16 DUM0,DUM1
1078 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1087 DATA (cf(i, 1),i= 1, 2) / 9, 3/
1090 DATA (cf(i, 2),i= 1, 2) / 3, 9/
1095 matrix_udx_udx_noh=0.d0
1096 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
1097 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
1102 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1103 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1104 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1105 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1106 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1107 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1108 CALL ffv2_3(w(1,1),w(1,2),gc_100,mw,ww,w(1,7))
1109 CALL ffv2_3(w(1,4),w(1,3),gc_100,mw,ww,w(1,8))
1110 CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
1112 CALL ffv2_0(w(1,9),w(1,6),w(1,8),gc_108,amp(1))
1113 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1115 CALL vvv1_0(w(1,9),w(1,7),w(1,8),gc_4,amp(2))
1116 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1118 CALL vvv1_0(w(1,7),w(1,8),w(1,10),gc_53,amp(3))
1119 CALL ffv2_1(w(1,3),w(1,7),gc_100,zero,zero,w(1,11))
1121 CALL ffv1_0(w(1,4),w(1,11),w(1,9),gc_1,amp(4))
1122 CALL ffv2_2(w(1,4),w(1,7),gc_100,zero,zero,w(1,12))
1124 CALL ffv1_0(w(1,12),w(1,3),w(1,9),gc_2,amp(5))
1126 CALL ffv2_3_0(w(1,4),w(1,11),w(1,10),gc_50,gc_58,amp(6))
1128 CALL ffv2_5_0(w(1,12),w(1,3),w(1,10),gc_51,gc_58,amp(7))
1129 CALL ffv1p0_3(w(1,1),w(1,3),gc_2,zero,zero,w(1,12))
1130 CALL ffv1p0_3(w(1,4),w(1,2),gc_1,zero,zero,w(1,11))
1131 CALL ffv1_2(w(1,5),w(1,12),gc_3,mta,zero,w(1,7))
1133 CALL ffv1_0(w(1,7),w(1,6),w(1,11),gc_3,amp(8))
1134 CALL ffv1_1(w(1,6),w(1,12),gc_3,mta,zero,w(1,13))
1136 CALL ffv1_0(w(1,5),w(1,13),w(1,11),gc_3,amp(9))
1137 CALL ffv2_3_3(w(1,4),w(1,2),gc_50,gc_58,mz,wz,w(1,14))
1139 CALL ffv2_4_0(w(1,7),w(1,6),w(1,14),gc_50,gc_59,amp(10))
1141 CALL ffv2_4_0(w(1,5),w(1,13),w(1,14),gc_50,gc_59,amp(11))
1142 CALL ffv2_5_3(w(1,1),w(1,3),gc_51,gc_58,mz,wz,w(1,13))
1143 CALL ffv2_4_2(w(1,5),w(1,13),gc_50,gc_59,mta,zero,w(1,7))
1145 CALL ffv1_0(w(1,7),w(1,6),w(1,11),gc_3,amp(12))
1146 CALL ffv2_4_1(w(1,6),w(1,13),gc_50,gc_59,mta,zero,w(1,15))
1148 CALL ffv1_0(w(1,5),w(1,15),w(1,11),gc_3,amp(13))
1150 CALL ffv2_4_0(w(1,7),w(1,6),w(1,14),gc_50,gc_59,amp(14))
1152 CALL ffv2_4_0(w(1,5),w(1,15),w(1,14),gc_50,gc_59,amp(15))
1153 CALL ffv1_1(w(1,2),w(1,12),gc_1,zero,zero,w(1,15))
1155 CALL ffv1_0(w(1,4),w(1,15),w(1,9),gc_1,amp(16))
1156 CALL ffv1_2(w(1,4),w(1,12),gc_1,zero,zero,w(1,5))
1158 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(17))
1160 CALL ffv2_3_0(w(1,4),w(1,15),w(1,10),gc_50,gc_58,amp(18))
1162 CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(19))
1163 CALL ffv1p0_3(w(1,1),w(1,3),gc_11,zero,zero,w(1,5))
1164 CALL ffv1_1(w(1,2),w(1,5),gc_11,zero,zero,w(1,15))
1166 CALL ffv1_0(w(1,4),w(1,15),w(1,9),gc_1,amp(20))
1167 CALL ffv1_2(w(1,4),w(1,5),gc_11,zero,zero,w(1,12))
1169 CALL ffv1_0(w(1,12),w(1,2),w(1,9),gc_1,amp(21))
1171 CALL ffv2_3_0(w(1,4),w(1,15),w(1,10),gc_50,gc_58,amp(22))
1173 CALL ffv2_3_0(w(1,12),w(1,2),w(1,10),gc_50,gc_58,amp(23))
1174 CALL ffv2_3_1(w(1,2),w(1,13),gc_50,gc_58,zero,zero,w(1,12))
1176 CALL ffv1_0(w(1,4),w(1,12),w(1,9),gc_1,amp(24))
1177 CALL ffv2_3_2(w(1,4),w(1,13),gc_50,gc_58,zero,zero,w(1,15))
1179 CALL ffv1_0(w(1,15),w(1,2),w(1,9),gc_1,amp(25))
1181 CALL ffv2_3_0(w(1,4),w(1,12),w(1,10),gc_50,gc_58,amp(26))
1183 CALL ffv2_3_0(w(1,15),w(1,2),w(1,10),gc_50,gc_58,amp(27))
1184 CALL ffv1_2(w(1,1),w(1,11),gc_2,zero,zero,w(1,15))
1186 CALL ffv1_0(w(1,15),w(1,3),w(1,9),gc_2,amp(28))
1187 CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,12))
1189 CALL ffv1_0(w(1,12),w(1,3),w(1,11),gc_2,amp(29))
1191 CALL ffv2_5_0(w(1,15),w(1,3),w(1,10),gc_51,gc_58,amp(30))
1192 CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,15))
1194 CALL ffv1_0(w(1,15),w(1,3),w(1,11),gc_2,amp(31))
1195 CALL ffv1p0_3(w(1,4),w(1,2),gc_11,zero,zero,w(1,11))
1196 CALL ffv1_2(w(1,1),w(1,11),gc_11,zero,zero,w(1,4))
1198 CALL ffv1_0(w(1,4),w(1,3),w(1,9),gc_2,amp(32))
1200 CALL ffv1_0(w(1,12),w(1,3),w(1,11),gc_11,amp(33))
1202 CALL ffv2_5_0(w(1,4),w(1,3),w(1,10),gc_51,gc_58,amp(34))
1204 CALL ffv1_0(w(1,15),w(1,3),w(1,11),gc_11,amp(35))
1205 CALL ffv2_5_2(w(1,1),w(1,14),gc_51,gc_58,zero,zero,w(1,11))
1207 CALL ffv1_0(w(1,11),w(1,3),w(1,9),gc_2,amp(36))
1209 CALL ffv2_5_0(w(1,12),w(1,3),w(1,14),gc_51,gc_58,amp(37))
1211 CALL ffv2_5_0(w(1,11),w(1,3),w(1,10),gc_51,gc_58,amp(38))
1213 CALL ffv2_5_0(w(1,15),w(1,3),w(1,14),gc_51,gc_58,amp(39))
1214 CALL ffv2_2(w(1,1),w(1,8),gc_100,zero,zero,w(1,14))
1216 CALL ffv1_0(w(1,14),w(1,2),w(1,9),gc_1,amp(40))
1218 CALL ffv2_0(w(1,12),w(1,2),w(1,8),gc_100,amp(41))
1220 CALL ffv2_3_0(w(1,14),w(1,2),w(1,10),gc_50,gc_58,amp(42))
1222 CALL ffv2_0(w(1,15),w(1,2),w(1,8),gc_100,amp(43))
1223 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)
1224 $ +1./2.*amp(20)+1./2.*amp(21)+1./2.*amp(22)+1./2.*amp(23)
1225 $ +1./2.*amp(32)+1./2.*amp(33)+1./2.*amp(34)+1./2.*amp(35)
1226 $ -amp(40)-amp(41)-amp(42)-amp(43)
1227 jamp(2)=+amp(8)+amp(9)+amp(10)+amp(11)+amp(12)+amp(13)+amp(14)
1228 $ +amp(15)+amp(16)+amp(17)+amp(18)+amp(19)-1./6.*amp(20)
1229 $ -1./6.*amp(21)-1./6.*amp(22)-1./6.*amp(23)+amp(24)+amp(25)
1230 $ +amp(26)+amp(27)+amp(28)+amp(29)+amp(30)+amp(31)-1./6.*amp(32)
1231 $ -1./6.*amp(33)-1./6.*amp(34)-1./6.*amp(35)+amp(36)+amp(37)
1238 ztemp = ztemp + cf(j,i)*jamp(j)
1240 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1242 matrix_udx_udx_noh=matrix
1249 SUBROUTINE udx_usx_noh(P,H1,H2,ANS)
1268 parameter(nexternal=6)
1270 parameter( ncomb=64)
1274 real*8 p(0:3,nexternal),ans
1279 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1281 real*8 matrix_udx_usx_noh
1282 INTEGER IHEL,IDEN, I
1283 INTEGER JC(NEXTERNAL)
1284 LOGICAL GOODHEL(NCOMB)
1286 DATA goodhel/ncomb*.false./
1287 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1288 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1289 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1290 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1291 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1292 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1293 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1294 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1295 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1296 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1297 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1298 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1299 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1300 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1301 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1302 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1303 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1304 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1305 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1306 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1307 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1308 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1309 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1310 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1311 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1312 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1313 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1314 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1315 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1316 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1317 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1318 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1319 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1320 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1321 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1322 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1323 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1324 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1325 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1326 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1327 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1328 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1329 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1330 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1331 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1332 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1333 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1334 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1335 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1336 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1337 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1338 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1339 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1340 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1341 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1342 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1343 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1344 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1345 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1346 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1347 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1348 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1349 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1350 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1360 t=matrix_udx_usx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1367 REAL*8 FUNCTION matrix_udx_usx_noh(P,H1,H2,NHEL,IC)
1383 parameter(ngraphs=11)
1385 parameter(nexternal=6)
1386 INTEGER NWAVEFUNCS, NCOLOR
1387 parameter(nwavefuncs=10, ncolor=1)
1391 parameter(imag1=(0d0,1d0))
1395 real*8 p(0:3,nexternal)
1396 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1402 real*8 denom(ncolor), cf(ncolor,ncolor)
1403 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1404 COMPLEX*16 W(18,NWAVEFUNCS)
1405 COMPLEX*16 DUM0,DUM1
1406 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1415 DATA (cf(i, 1),i= 1, 1) / 9/
1420 matrix_udx_usx_noh=0.d0
1421 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
1422 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
1428 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1429 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1430 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1431 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1432 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1433 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1434 CALL ffv2_3(w(1,1),w(1,2),gc_100,mw,ww,w(1,7))
1435 CALL ffv2_3(w(1,4),w(1,3),gc_101,mw,ww,w(1,8))
1436 CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
1438 CALL ffv2_0(w(1,9),w(1,6),w(1,8),gc_108,amp(1))
1439 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1441 CALL vvv1_0(w(1,9),w(1,7),w(1,8),gc_4,amp(2))
1442 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1444 CALL vvv1_0(w(1,7),w(1,8),w(1,10),gc_53,amp(3))
1445 CALL ffv2_1(w(1,3),w(1,7),gc_101,zero,zero,w(1,6))
1447 CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_1,amp(4))
1448 CALL ffv2_2(w(1,4),w(1,7),gc_101,zero,zero,w(1,5))
1450 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(5))
1452 CALL ffv2_3_0(w(1,4),w(1,6),w(1,10),gc_50,gc_58,amp(6))
1454 CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(7))
1455 CALL ffv2_2(w(1,1),w(1,8),gc_100,zero,zero,w(1,5))
1457 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(8))
1458 CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,3))
1460 CALL ffv2_0(w(1,3),w(1,2),w(1,8),gc_100,amp(9))
1462 CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(10))
1463 CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,5))
1465 CALL ffv2_0(w(1,5),w(1,2),w(1,8),gc_100,amp(11))
1466 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
1467 $ -amp(9)-amp(10)-amp(11)
1473 ztemp = ztemp + cf(j,i)*jamp(j)
1475 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1477 matrix_udx_usx_noh=matrix
1484 SUBROUTINE udx_cdx_noh(P,H1,H2,ANS)
1503 parameter(nexternal=6)
1505 parameter( ncomb=64)
1509 real*8 p(0:3,nexternal),ans
1514 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1516 real*8 matrix_udx_cdx_noh
1517 INTEGER IHEL,IDEN, I
1518 INTEGER JC(NEXTERNAL)
1519 LOGICAL GOODHEL(NCOMB)
1521 DATA goodhel/ncomb*.false./
1522 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1523 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1524 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1525 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1526 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1527 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1528 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1529 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1530 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1531 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1532 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1533 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1534 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1535 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1536 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1537 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1538 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1539 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1540 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1541 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1542 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1543 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1544 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1545 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1546 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1547 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1548 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1549 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1550 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1551 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1552 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1553 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1554 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1555 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1556 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1557 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1558 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1559 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1560 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1561 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1562 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1563 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1564 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1565 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1566 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1567 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1568 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1569 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1570 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1571 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1572 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1573 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1574 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1575 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1576 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1577 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1578 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1579 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1580 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1581 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1582 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1583 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1584 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1585 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1595 t=matrix_udx_cdx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1602 REAL*8 FUNCTION matrix_udx_cdx_noh(P,H1,H2,NHEL,IC)
1618 parameter(ngraphs=11)
1620 parameter(nexternal=6)
1621 INTEGER NWAVEFUNCS, NCOLOR
1622 parameter(nwavefuncs=10, ncolor=1)
1626 parameter(imag1=(0d0,1d0))
1630 real*8 p(0:3,nexternal)
1631 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1637 real*8 denom(ncolor), cf(ncolor,ncolor)
1638 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1639 COMPLEX*16 W(18,NWAVEFUNCS)
1640 COMPLEX*16 DUM0,DUM1
1641 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1650 DATA (cf(i, 1),i= 1, 1) / 9/
1655 matrix_udx_cdx_noh=0.d0
1656 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
1657 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
1664 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1665 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1666 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1667 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1668 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1669 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1670 CALL ffv2_3(w(1,1),w(1,2),gc_100,mw,ww,w(1,7))
1671 CALL ffv2_3(w(1,4),w(1,3),gc_44,mw,ww,w(1,8))
1672 CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
1674 CALL ffv2_0(w(1,9),w(1,6),w(1,8),gc_108,amp(1))
1675 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1677 CALL vvv1_0(w(1,9),w(1,7),w(1,8),gc_4,amp(2))
1678 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1680 CALL vvv1_0(w(1,7),w(1,8),w(1,10),gc_53,amp(3))
1681 CALL ffv2_1(w(1,3),w(1,7),gc_44,zero,zero,w(1,6))
1683 CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_1,amp(4))
1684 CALL ffv2_2(w(1,4),w(1,7),gc_44,zero,zero,w(1,5))
1686 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(5))
1688 CALL ffv2_3_0(w(1,4),w(1,6),w(1,10),gc_50,gc_58,amp(6))
1690 CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(7))
1691 CALL ffv2_2(w(1,1),w(1,8),gc_100,zero,zero,w(1,5))
1693 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(8))
1694 CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,3))
1696 CALL ffv2_0(w(1,3),w(1,2),w(1,8),gc_100,amp(9))
1698 CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(10))
1699 CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,5))
1701 CALL ffv2_0(w(1,5),w(1,2),w(1,8),gc_100,amp(11))
1702 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
1703 $ -amp(9)-amp(10)-amp(11)
1709 ztemp = ztemp + cf(j,i)*jamp(j)
1711 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1713 matrix_udx_cdx_noh=matrix
1720 SUBROUTINE udx_csx_noh(P,H1,H2,ANS)
1739 parameter(nexternal=6)
1741 parameter( ncomb=64)
1745 real*8 p(0:3,nexternal),ans
1750 INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
1752 real*8 matrix_udx_csx_noh
1753 INTEGER IHEL,IDEN, I
1754 INTEGER JC(NEXTERNAL)
1755 LOGICAL GOODHEL(NCOMB)
1757 DATA goodhel/ncomb*.false./
1758 DATA (nhel(i, 1),i=1,6) /-1,-1,-1,-1,-1,-1/
1759 DATA (nhel(i, 2),i=1,6) /-1,-1,-1,-1,-1, 1/
1760 DATA (nhel(i, 3),i=1,6) /-1,-1,-1,-1, 1,-1/
1761 DATA (nhel(i, 4),i=1,6) /-1,-1,-1,-1, 1, 1/
1762 DATA (nhel(i, 5),i=1,6) /-1,-1,-1, 1,-1,-1/
1763 DATA (nhel(i, 6),i=1,6) /-1,-1,-1, 1,-1, 1/
1764 DATA (nhel(i, 7),i=1,6) /-1,-1,-1, 1, 1,-1/
1765 DATA (nhel(i, 8),i=1,6) /-1,-1,-1, 1, 1, 1/
1766 DATA (nhel(i, 9),i=1,6) /-1,-1, 1,-1,-1,-1/
1767 DATA (nhel(i, 10),i=1,6) /-1,-1, 1,-1,-1, 1/
1768 DATA (nhel(i, 11),i=1,6) /-1,-1, 1,-1, 1,-1/
1769 DATA (nhel(i, 12),i=1,6) /-1,-1, 1,-1, 1, 1/
1770 DATA (nhel(i, 13),i=1,6) /-1,-1, 1, 1,-1,-1/
1771 DATA (nhel(i, 14),i=1,6) /-1,-1, 1, 1,-1, 1/
1772 DATA (nhel(i, 15),i=1,6) /-1,-1, 1, 1, 1,-1/
1773 DATA (nhel(i, 16),i=1,6) /-1,-1, 1, 1, 1, 1/
1774 DATA (nhel(i, 17),i=1,6) /-1, 1,-1,-1,-1,-1/
1775 DATA (nhel(i, 18),i=1,6) /-1, 1,-1,-1,-1, 1/
1776 DATA (nhel(i, 19),i=1,6) /-1, 1,-1,-1, 1,-1/
1777 DATA (nhel(i, 20),i=1,6) /-1, 1,-1,-1, 1, 1/
1778 DATA (nhel(i, 21),i=1,6) /-1, 1,-1, 1,-1,-1/
1779 DATA (nhel(i, 22),i=1,6) /-1, 1,-1, 1,-1, 1/
1780 DATA (nhel(i, 23),i=1,6) /-1, 1,-1, 1, 1,-1/
1781 DATA (nhel(i, 24),i=1,6) /-1, 1,-1, 1, 1, 1/
1782 DATA (nhel(i, 25),i=1,6) /-1, 1, 1,-1,-1,-1/
1783 DATA (nhel(i, 26),i=1,6) /-1, 1, 1,-1,-1, 1/
1784 DATA (nhel(i, 27),i=1,6) /-1, 1, 1,-1, 1,-1/
1785 DATA (nhel(i, 28),i=1,6) /-1, 1, 1,-1, 1, 1/
1786 DATA (nhel(i, 29),i=1,6) /-1, 1, 1, 1,-1,-1/
1787 DATA (nhel(i, 30),i=1,6) /-1, 1, 1, 1,-1, 1/
1788 DATA (nhel(i, 31),i=1,6) /-1, 1, 1, 1, 1,-1/
1789 DATA (nhel(i, 32),i=1,6) /-1, 1, 1, 1, 1, 1/
1790 DATA (nhel(i, 33),i=1,6) / 1,-1,-1,-1,-1,-1/
1791 DATA (nhel(i, 34),i=1,6) / 1,-1,-1,-1,-1, 1/
1792 DATA (nhel(i, 35),i=1,6) / 1,-1,-1,-1, 1,-1/
1793 DATA (nhel(i, 36),i=1,6) / 1,-1,-1,-1, 1, 1/
1794 DATA (nhel(i, 37),i=1,6) / 1,-1,-1, 1,-1,-1/
1795 DATA (nhel(i, 38),i=1,6) / 1,-1,-1, 1,-1, 1/
1796 DATA (nhel(i, 39),i=1,6) / 1,-1,-1, 1, 1,-1/
1797 DATA (nhel(i, 40),i=1,6) / 1,-1,-1, 1, 1, 1/
1798 DATA (nhel(i, 41),i=1,6) / 1,-1, 1,-1,-1,-1/
1799 DATA (nhel(i, 42),i=1,6) / 1,-1, 1,-1,-1, 1/
1800 DATA (nhel(i, 43),i=1,6) / 1,-1, 1,-1, 1,-1/
1801 DATA (nhel(i, 44),i=1,6) / 1,-1, 1,-1, 1, 1/
1802 DATA (nhel(i, 45),i=1,6) / 1,-1, 1, 1,-1,-1/
1803 DATA (nhel(i, 46),i=1,6) / 1,-1, 1, 1,-1, 1/
1804 DATA (nhel(i, 47),i=1,6) / 1,-1, 1, 1, 1,-1/
1805 DATA (nhel(i, 48),i=1,6) / 1,-1, 1, 1, 1, 1/
1806 DATA (nhel(i, 49),i=1,6) / 1, 1,-1,-1,-1,-1/
1807 DATA (nhel(i, 50),i=1,6) / 1, 1,-1,-1,-1, 1/
1808 DATA (nhel(i, 51),i=1,6) / 1, 1,-1,-1, 1,-1/
1809 DATA (nhel(i, 52),i=1,6) / 1, 1,-1,-1, 1, 1/
1810 DATA (nhel(i, 53),i=1,6) / 1, 1,-1, 1,-1,-1/
1811 DATA (nhel(i, 54),i=1,6) / 1, 1,-1, 1,-1, 1/
1812 DATA (nhel(i, 55),i=1,6) / 1, 1,-1, 1, 1,-1/
1813 DATA (nhel(i, 56),i=1,6) / 1, 1,-1, 1, 1, 1/
1814 DATA (nhel(i, 57),i=1,6) / 1, 1, 1,-1,-1,-1/
1815 DATA (nhel(i, 58),i=1,6) / 1, 1, 1,-1,-1, 1/
1816 DATA (nhel(i, 59),i=1,6) / 1, 1, 1,-1, 1,-1/
1817 DATA (nhel(i, 60),i=1,6) / 1, 1, 1,-1, 1, 1/
1818 DATA (nhel(i, 61),i=1,6) / 1, 1, 1, 1,-1,-1/
1819 DATA (nhel(i, 62),i=1,6) / 1, 1, 1, 1,-1, 1/
1820 DATA (nhel(i, 63),i=1,6) / 1, 1, 1, 1, 1,-1/
1821 DATA (nhel(i, 64),i=1,6) / 1, 1, 1, 1, 1, 1/
1831 t=matrix_udx_csx_noh(p ,h1,h2,nhel(1,ihel),jc(1))
1838 REAL*8 FUNCTION matrix_udx_csx_noh(P,H1,H2,NHEL,IC)
1854 parameter(ngraphs=11)
1856 parameter(nexternal=6)
1857 INTEGER NWAVEFUNCS, NCOLOR
1858 parameter(nwavefuncs=10, ncolor=1)
1862 parameter(imag1=(0d0,1d0))
1866 real*8 p(0:3,nexternal)
1867 INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
1873 real*8 denom(ncolor), cf(ncolor,ncolor)
1874 COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
1875 COMPLEX*16 W(18,NWAVEFUNCS)
1876 COMPLEX*16 DUM0,DUM1
1877 DATA dum0, dum1/(0d0, 0d0), (1d0, 0d0)/
1886 DATA (cf(i, 1),i= 1, 1) / 9/
1891 matrix_udx_csx_noh=0.d0
1892 IF(h1.EQ.0. or .h1.EQ.nhel(5))
THEN
1893 IF(h2.EQ.0. or .h2.EQ.nhel(6))
THEN
1898 CALL ixxxxx(p(0,1),zero,nhel(1),+1*ic(1),w(1,1))
1899 CALL oxxxxx(p(0,2),zero,nhel(2),-1*ic(2),w(1,2))
1900 CALL oxxxxx(p(0,3),zero,nhel(3),+1*ic(3),w(1,3))
1901 CALL ixxxxx(p(0,4),zero,nhel(4),-1*ic(4),w(1,4))
1902 CALL ixxxxx(p(0,5),mta,nhel(5),-1*ic(5),w(1,5))
1903 CALL oxxxxx(p(0,6),mta,nhel(6),+1*ic(6),w(1,6))
1904 CALL ffv2_3(w(1,1),w(1,2),gc_100,mw,ww,w(1,7))
1905 CALL ffv2_3(w(1,4),w(1,3),gc_100,mw,ww,w(1,8))
1906 CALL ffv2_2(w(1,5),w(1,7),gc_108,zero,zero,w(1,9))
1908 CALL ffv2_0(w(1,9),w(1,6),w(1,8),gc_108,amp(1))
1909 CALL ffv1p0_3(w(1,5),w(1,6),gc_3,zero,zero,w(1,9))
1911 CALL vvv1_0(w(1,9),w(1,7),w(1,8),gc_4,amp(2))
1912 CALL ffv2_4_3(w(1,5),w(1,6),gc_50,gc_59,mz,wz,w(1,10))
1914 CALL vvv1_0(w(1,7),w(1,8),w(1,10),gc_53,amp(3))
1915 CALL ffv2_1(w(1,3),w(1,7),gc_100,zero,zero,w(1,6))
1917 CALL ffv1_0(w(1,4),w(1,6),w(1,9),gc_1,amp(4))
1918 CALL ffv2_2(w(1,4),w(1,7),gc_100,zero,zero,w(1,5))
1920 CALL ffv1_0(w(1,5),w(1,3),w(1,9),gc_2,amp(5))
1922 CALL ffv2_3_0(w(1,4),w(1,6),w(1,10),gc_50,gc_58,amp(6))
1924 CALL ffv2_5_0(w(1,5),w(1,3),w(1,10),gc_51,gc_58,amp(7))
1925 CALL ffv2_2(w(1,1),w(1,8),gc_100,zero,zero,w(1,5))
1927 CALL ffv1_0(w(1,5),w(1,2),w(1,9),gc_1,amp(8))
1928 CALL ffv1_2(w(1,1),w(1,9),gc_2,zero,zero,w(1,3))
1930 CALL ffv2_0(w(1,3),w(1,2),w(1,8),gc_100,amp(9))
1932 CALL ffv2_3_0(w(1,5),w(1,2),w(1,10),gc_50,gc_58,amp(10))
1933 CALL ffv2_5_2(w(1,1),w(1,10),gc_51,gc_58,zero,zero,w(1,5))
1935 CALL ffv2_0(w(1,5),w(1,2),w(1,8),gc_100,amp(11))
1936 jamp(1)=-amp(1)-amp(2)-amp(3)-amp(4)-amp(5)-amp(6)-amp(7)-amp(8)
1937 $ -amp(9)-amp(10)-amp(11)
1943 ztemp = ztemp + cf(j,i)*jamp(j)
1945 matrix = matrix+ztemp*dconjg(jamp(i))/denom(i)
1947 matrix_udx_csx_noh=matrix