NptJregs/4ptJregs.h

   1
   2    id Jpt4reg1(?m1)*Jpt4reg1(?m2)*Jpt4reg1(?m3) =
   3
   4        + zeta2*zwgt^2 * (
   5           - KK(?m3,[c],?m2)
   6           )
   7
   8           #if 'w' >= 3
   9        + zeta3*zwgt^3 * (
  10           + KK(?m2,[c],?m1)*KK(?m3,[c],?m2)
  11           + KK(?m3,[c],?m2)^2
  12           )
  13
  14           #if 'w' >= 4
  15        + zeta2^2*zwgt^4 * (
  16           - 2/5*KK(?m2,[c],?m1)^2*KK(?m3,[c],?m2)
  17           - 1/10*KK(?m2,[c],?m1)*KK(?m3,[c],?m2)^2
  18           - 2/5*KK(?m3,[c],?m2)^3
  19           )
  20
  21           #if 'w' >= 5
  22        + zeta2*zeta3*zwgt^5 * (
  23           - KK(?m2,[c],?m1)^2*KK(?m3,[c],?m2)^2
  24           - KK(?m2,[c],?m1)*KK(?m3,[c],?m2)^3
  25           )
  26
  27        + zeta5*zwgt^5 * (
  28           + KK(?m2,[c],?m1)^3*KK(?m3,[c],?m2)
  29           + 2*KK(?m2,[c],?m1)^2*KK(?m3,[c],?m2)^2
  30           + 2*KK(?m2,[c],?m1)*KK(?m3,[c],?m2)^3
  31           + KK(?m3,[c],?m2)^4
  32           )
  33
  34           #if 'w' >= 6
  35        + zeta2^3*zwgt^6 * (
  36           - 8/35*KK(?m2,[c],?m1)^4*KK(?m3,[c],?m2)
  37           - 6/35*KK(?m2,[c],?m1)^3*KK(?m3,[c],?m2)^2
  38           - 23/70*KK(?m2,[c],?m1)^2*KK(?m3,[c],?m2)^3
  39           - 6/35*KK(?m2,[c],?m1)*KK(?m3,[c],?m2)^4
  40           - 8/35*KK(?m3,[c],?m2)^5
  41           )
  42
  43        + zeta3^2*zwgt^6 * (
  44           + 1/2*KK(?m2,[c],?m1)^3*KK(?m3,[c],?m2)^2
  45           + KK(?m2,[c],?m1)^2*KK(?m3,[c],?m2)^3
  46           + 1/2*KK(?m2,[c],?m1)*KK(?m3,[c],?m2)^4
  47           )
  48
  49           #if 'w' >= 7
  50        + zeta2^2*zeta3*zwgt^7 * (
  51           - 2/5*KK(?m2,[c],?m1)^4*KK(?m3,[c],?m2)^2
  52           - 1/2*KK(?m2,[c],?m1)^3*KK(?m3,[c],?m2)^3
  53           - 1/2*KK(?m2,[c],?m1)^2*KK(?m3,[c],?m2)^4
  54           - 2/5*KK(?m2,[c],?m1)*KK(?m3,[c],?m2)^5
  55           )
  56
  57        + zeta2*zeta5*zwgt^7 * (
  58           - KK(?m2,[c],?m1)^4*KK(?m3,[c],?m2)^2
  59           - 2*KK(?m2,[c],?m1)^3*KK(?m3,[c],?m2)^3
  60           - 2*KK(?m2,[c],?m1)^2*KK(?m3,[c],?m2)^4
  61           - KK(?m2,[c],?m1)*KK(?m3,[c],?m2)^5
  62           )
  63
  64        + zeta7*zwgt^7 * (
  65           + KK(?m2,[c],?m1)^5*KK(?m3,[c],?m2)
  66           + 3*KK(?m2,[c],?m1)^4*KK(?m3,[c],?m2)^2
  67           + 5*KK(?m2,[c],?m1)^3*KK(?m3,[c],?m2)^3
  68           + 5*KK(?m2,[c],?m1)^2*KK(?m3,[c],?m2)^4
  69           + 3*KK(?m2,[c],?m1)*KK(?m3,[c],?m2)^5
  70           + KK(?m3,[c],?m2)^6
  71           )
  72           #endif
  73           #endif
  74           #endif
  75           #endif
  76           #endif
  77           ;
  78
  79    id Jpt4reg2(?m1)*Jpt4reg2(?m2)*Jpt4reg2(?m3) =
  80
  81        + zeta2*zwgt^2 * (
  82           - KK(?m2,[c],?m1)
  83           )
  84
  85           #if 'w' >= 3
  86        + zeta3*zwgt^3 * (
  87           + KK(?m2,[c],?m1)^2
  88           + KK(?m2,[c],?m1)*KK(?m3,[c],?m2)
  89           )
  90
  91           #if 'w' >= 4
  92        + zeta2^2*zwgt^4 * (
  93           - 2/5*KK(?m2,[c],?m1)^3
  94           - 1/10*KK(?m2,[c],?m1)^2*KK(?m3,[c],?m2)
  95           - 2/5*KK(?m2,[c],?m1)*KK(?m3,[c],?m2)^2
  96           )
  97
  98           #if 'w' >= 5
  99        + zeta2*zeta3*zwgt^5 * (
 100           - KK(?m2,[c],?m1)^3*KK(?m3,[c],?m2)
 101           - KK(?m2,[c],?m1)^2*KK(?m3,[c],?m2)^2
 102           )
 103
 104        + zeta5*zwgt^5 * (
 105           + KK(?m2,[c],?m1)^4
 106           + 2*KK(?m2,[c],?m1)^3*KK(?m3,[c],?m2)
 107           + 2*KK(?m2,[c],?m1)^2*KK(?m3,[c],?m2)^2
 108           + KK(?m2,[c],?m1)*KK(?m3,[c],?m2)^3
 109           )
 110
 111           #if 'w' >= 6
 112        + zeta2^3*zwgt^6 * (
 113           - 8/35*KK(?m2,[c],?m1)^5
 114           - 6/35*KK(?m2,[c],?m1)^4*KK(?m3,[c],?m2)
 115           - 23/70*KK(?m2,[c],?m1)^3*KK(?m3,[c],?m2)^2
 116           - 6/35*KK(?m2,[c],?m1)^2*KK(?m3,[c],?m2)^3
 117           - 8/35*KK(?m2,[c],?m1)*KK(?m3,[c],?m2)^4
 118           )
 119
 120        + zeta3^2*zwgt^6 * (
 121           + 1/2*KK(?m2,[c],?m1)^4*KK(?m3,[c],?m2)
 122           + KK(?m2,[c],?m1)^3*KK(?m3,[c],?m2)^2
 123           + 1/2*KK(?m2,[c],?m1)^2*KK(?m3,[c],?m2)^3
 124           )
 125
 126           #if 'w' >= 7
 127        + zeta2^2*zeta3*zwgt^7 * (
 128           - 2/5*KK(?m2,[c],?m1)^5*KK(?m3,[c],?m2)
 129           - 1/2*KK(?m2,[c],?m1)^4*KK(?m3,[c],?m2)^2
 130           - 1/2*KK(?m2,[c],?m1)^3*KK(?m3,[c],?m2)^3
 131           - 2/5*KK(?m2,[c],?m1)^2*KK(?m3,[c],?m2)^4
 132           )
 133
 134        + zeta2*zeta5*zwgt^7 * (
 135           - KK(?m2,[c],?m1)^5*KK(?m3,[c],?m2)
 136           - 2*KK(?m2,[c],?m1)^4*KK(?m3,[c],?m2)^2
 137           - 2*KK(?m2,[c],?m1)^3*KK(?m3,[c],?m2)^3
 138           - KK(?m2,[c],?m1)^2*KK(?m3,[c],?m2)^4
 139           )
 140
 141        + zeta7*zwgt^7 * (
 142           + KK(?m2,[c],?m1)^6
 143           + 3*KK(?m2,[c],?m1)^5*KK(?m3,[c],?m2)
 144           + 5*KK(?m2,[c],?m1)^4*KK(?m3,[c],?m2)^2
 145           + 5*KK(?m2,[c],?m1)^3*KK(?m3,[c],?m2)^3
 146           + 3*KK(?m2,[c],?m1)^2*KK(?m3,[c],?m2)^4
 147           + KK(?m2,[c],?m1)*KK(?m3,[c],?m2)^5
 148           )
 149           #endif
 150           #endif
 151           #endif
 152           #endif
 153           #endif
 154           ;
 155
 156