auxiliary/zremovals/zremovalCte-length4.h

   1
   2         bracket G;
   3         .sort:beg1cte4;
   4         keep brackets;
   5
   6 * shuffle ids to eliminate z1 from the first slot
   7         id G('z1','z1','z1',w2?!{,'z1'}) = 1/3*G('z1',w2)*G('z1','z1',w2);
   8
   9         id G('z1','z1',a3?!{,'z1'},w2?!{,'z1'}) =
  10         + 1/2*G('z1',w2)*G('z1',a3,w2)
  11         - 1/2*G('z1',a3,'z1',w2)
  12         ;
  13
  14         id G('z1',a2?!{,'z1'},'z1',w2?!{,'z1'}) =
  15         + G('z1',w2)*G(a2,'z1',w2)
  16         - 2*G(a2,'z1','z1',w2)
  17         ;
  18
  19         id G('z1',a2?!{,'z1'},a3?!{,'z1'},w2?!{,'z1'}) =
  20         + G('z1',w2)*G(a2,a3,w2)
  21         - G(a2,'z1',a3,w2)
  22         - G(a2,a3,'z1',w2)
  23         ;
  24
  25 * end shuffles
  26
  27    id G(a1?!{,'z1'},a2?!{,'z1'},'z1',w2?!{,'z1'}) =
  28        + G(a1,a2,w2)*G(a2,'z1')
  29        + G(a1,a2,0,w2)
  30        - G(a1,w2)*G(a2,a1,'z1')
  31        + G(a1,w2)*G(a2,0,'z1')
  32        - G(a1,0,w2)*G(a2,'z1')
  33        + G(a2,a1,w2,'z1')
  34        - G(a2,a1,0,'z1')
  35        + 2*G(a2,'z1')*zeta2*delta(a1,0)
  36        + G(0,w2)*G(a2,a1,'z1')
  37        - G(0,'z1')*G(a1,a2,w2)
  38        - 2*G(a1,w2)*zeta2*delta(a2,0)
  39 * erik
  40         + (
  41         - Pi*sign(w2,'z1')*G(a2,a1,'z1')*i_
  42         - Pi*sign(w2,'z1')*G(a1,a2,w2)*i_
  43         + Pi*sign(w2,'z1')*G(a2,'z1')*G(a1,w2)*i_
  44         )
  45        ;
  46
  47         id G(a1?!{,'z1'},'z1',a2?!{,'z1'},w2?!{,'z1'}) =
  48        + G(a1,a2,w2)*G(a1,'z1')
  49        - G(a1,a2,w2)*G(a2,'z1')
  50        - G(a1,a2,w2,'z1')
  51        + G(a1,a2,'z1')*G(a2,w2)
  52        + G(a1,a2,0,'z1')
  53        + G(a1,w2)*G(a2,a1,'z1')
  54        - G(a1,w2)*G(a2,0,'z1')
  55        - G(a1,w2,'z1')*G(a2,w2)
  56        - 2*G(a1,'z1')*zeta2*delta(a2,0)
  57        + G(a1,0,a2,w2)
  58        + G(a1,0,w2)*G(a2,'z1')
  59        - G(a2,a1,w2,'z1')
  60        + G(a2,a1,0,'z1')
  61        - 2*G(a2,'z1')*zeta2*delta(a1,0)
  62        - G(0,a2,w2)*G(a1,'z1')
  63        - G(0,w2)*G(a1,a2,'z1')
  64        - G(0,w2)*G(a2,a1,'z1')
  65         + 2*G(a1,w2)*zeta2*delta(a2,0)
  66
  67 * erik
  68         + (
  69         + Pi*sign(w2,'z1')*G(a1,a2,'z1')*i_
  70         - Pi*sign(w2,'z1')*G(a2,'z1')*G(a1,w2)*i_
  71         + Pi*sign(w2,'z1')*G(a2,a1,'z1')*i_
  72         )
  73       ;
  74
  75         id G(a1?!{,'z1'},'z1','z1',w2?!{,'z1'}) =
  76        - G(a1,a1,w2,'z1')
  77        + G(a1,a1,'z1')*G(a1,w2)
  78        + G(a1,a1,0,'z1')
  79        - G(a1,w2,w2,'z1')
  80        + G(a1,w2,0,'z1')
  81        + G(a1,0,w2)*G(a1,'z1')
  82        + G(a1,0,w2,'z1')
  83        - G(a1,0,'z1')*G(a1,w2)
  84        + G(a1,0,0,w2)
  85        - G(a1,0,0,'z1')
  86        + G(0,a1,w2,'z1')
  87        - G(0,a1,'z1')*G(a1,w2)
  88        + G(0,a1,'z1')*G(0,w2)
  89        - G(0,a1,0,'z1')
  90        - G(0,w2)*G(a1,a1,'z1')
  91        - G(0,w2)*G(a1,w2,'z1')
  92        + G(0,w2)*G(a1,0,'z1')
  93        - G(0,'z1')*G(a1,0,w2)
  94        - G(0,0,w2)*G(a1,'z1')
  95         + G(0,0,'z1')*G(a1,w2)
  96         + zeta3*delta(a1,0)
  97 * erik
  98         + (
  99         + Pi*sign(w2,'z1')*G(a1,'z1')*G(0,w2)*i_
 100         + Pi*sign(w2,'z1')*G(0,'z1')*G(a1,w2)*i_
 101         - Pi*sign(w2,'z1')*G(a1,'z1')*G(a1,w2)*i_
 102         - Pi*sign(w2,'z1')*G(0,a1,'z1')*i_
 103         + Pi*sign(w2,'z1')*G(a1,w2,'z1')*i_
 104         - Pi*sign(w2,'z1')*G(a1,0,'z1')*i_
 105         + Pi*sign(w2,'z1')*G(a1,a1,'z1')*i_
 106         - Pi*sign(w2,'z1')*G(a1,0,w2)*i_
 107         - 1/2*Pi^2*sign(w2,'z1')^2*G(a1,w2)
 108         + 1/2*Pi^2*sign(w2,'z1')^2*G(a1,'z1')
 109 * erik again
 110         + Pi*sign(w2,'z1')*zeta2*i_*delta(a1,0)
 111         )
 112       ;
 113
 114
 115         id Pi^2 = 6*zeta2;
 116         id sign(i?,j?)^2 = 1;
 117         id delta(0,0) = 1;
 118         id once delta(0,1) = 0;
 119         #call ContourSign()
 120         .sort:Ctelen4;