auxiliary/all4ptJregs.frm

   1         Format 250;
   2         Off statistics;
   3
   4         #include- pss2-polylogs.h
   5
   6 #procedure IntegrateOrder(i)
   7
   8         id order('i') = Int('i')*KN('Npts');
   9         #call IntegrationOrdering()
  10 * According to Table 1 of 1609.XXXX, expanding the KN factor
  11 * up to weight 6 computes the regular 4pt integrals up to zeta7.
  12 * You can select up to weight 6 here. If you want more, you'll
  13 * need to generate z-removal identities at higher length.
  14         #call KobaNielsen('Npts',6)
  15         #call IntegrateG(z'i')
  16
  17 #endprocedure
  18
  19 #procedure Sliced(i)
  20
  21         transform Z replace(1,last)=(1,z1,2,z2,3,z3,4,z4,5,z5,6,z6,7,z7,8,z8);
  22
  23         argument Z;
  24         id z1 = 0;
  25         id z{'Npts'-1} = 1;
  26         endargument;
  27         .sort
  28
  29         #call IntegrateOrder('i')
  30
  31         .sort
  32         Polyfun;
  33         id gather(x1?) = x1;
  34
  35         bracket zeta2,zeta3,zeta5,zeta7;
  36         print +s -f;
  37         .sort
  38         drop;
  39         .sort
  40
  41 #endprocedure
  42
  43         #define Npts "4"
  44
  45         L [Z(1,2)] = order(2)*Z(1,2);
  46         L [Z(2,3)] = order(2)*Z(2,3);
  47         #call Sliced(2)
  48
  49         .end