1% Tests of Cvitanovic Package. 2 3 4% COPYRIGHT (C) 1990, INSTITUTE OF NUCLEAR PHYSICS, MOSCOW STATE UNIV. 5% CVITBUBLE TEST OF CVITANOVIC PACKAGE 6% AUTHOR A. KRYUKOV, ARODIONOV, A.TARANOV 7% VERSION 1.1 8% RELEASE 18-SEP-90 9 10index j1,j2,j3,j4,j5,j6,j7,j8,j9,j0; 11 12 13vecdim n$ 14 15 16 17% Tests of the weels with buble 18% (Use notation from SIGSAM Bull, 1989, v.23, no.4, pp.15-24) 19 20g(l,j1,j2,j2,j1); 21 22 23 2 24n 25 26 27g(l,j1,j2)*g(l1,j3,j1,j2,j3); 28 29 30 2 31n 32 33g(l,j1,j2)*g(l1,j3,j1,j3,j2); 34 35 36n*( - n + 2) 37 38g(l,j1,j2)*g(l1,j3,j3,j2,j1); 39 40 41 2 42n 43 44 45g(l,j1,j2,j3,j4)*g(l1,j1,j2,j3,j4); 46 47 48n*(3*n - 2) 49 50g(l,j1,j2)*g(l1,j3,j4,j1,j2,j4,j3); 51 52 53 3 54n 55 56g(l,j1,j2,j3,j4)*g(l1,j1,j4,j2,j3); 57 58 59n*( - n + 2) 60 61g(l,j1,j2)*g(l1,j3,j4,j1,j4,j3,j2); 62 63 64 2 65n*(n - 4*n + 4) 66 67 68g(l,j1,j2)*g(l1,j3,j4,j5,j1,j2,j3,j4,j5); 69 70 71 2 2 72n *( - n + 6*n - 4) 73 74g(l,j1,j2,j3,j4)*g(l1,j5,j1,j2,j3,j5,j4); 75 76 77 2 78n*( - 3*n + 8*n - 4) 79 80g(l,j1,j2,j3,j4,j5,j1)*g(l1,j2,j5,j3,j4); 81 82 83 2 84n *( - n + 2) 85 86g(l,j1,j2,j3,j4,j5,j1,j2,j5)*g(l1,j4,j3); 87 88 89 2 2 90n *( - n + 6*n - 4) 91 92 93g(l,j1,j2)*g(l1,j3,j4,j5,j6,j1,j2,j3,j4,j5,j6); 94 95 96 2 3 2 97n *(n - 12*n + 28*n - 16) 98 99g(l,j1,j2,j3,j4)*g(l1,j5,j6,j1,j2,j3,j4,j6,j5); 100 101 102 3 103n *(3*n - 2) 104 105g(l,j1,j2,j3,j4,j5,j6)*g(l1,j1,j2,j4,j3,j6,j5); 106 107 108 2 109n*(7*n - 22*n + 16) 110 111g(l,j1,j2,j3,j4,j5,j6,j1,j2)*g(l1,j6,j3,j4,j5); 112 113 114 2 2 115n *( - 3*n + 8*n - 4) 116 117g(l,j1,j2,j3,j4,j5,j6,j7,j1,j2,j3,j4,j5)*g(l1,j6,j7); 118 119 120 2 4 3 2 121n *(n - 20*n + 100*n - 160*n + 80) 122 123g(l,j1,j2,j3,j4,j5,j6,j7,j1,j2,j3)*g(l1,j4,j5,j7,j6); 124 125 126 2 3 2 127n *(n - 8*n + 16*n - 8) 128 129g(l,j1,j2,j3,j4,j5,j6,j7,j2)*g(l1,j1,j3,j4,j5,j6,j7); 130 131 132 3 2 133n*( - 15*n + 60*n - 76*n + 32) 134 135 136 137% COPYRIGHT (C) 1988,1990, INSTITUTE OF NUCLEAR PHYSICS, MOSCOW STATE U. 138% CVITEST Test of CVITANOVIC PACKAGE 139% AUTHOR A. KRYUKOV, A.RODIONOV, A.TARANOV 140% VERSION 1.2 141% RELEASE 11-MAR-90 142 143% 144% Test for trace of Dirac matrices. 145% 146% All tests are the lattices with difference lines 147% (Use notation from SIGSAM Bull, 1989, v.4,no.23, pp.15-24) 148 149index m1,m2,m3,m4,m5,m6,m7,m8,m9,m0; 150 151 152index n1,n2,n3,n4,n5,n6,n7,n8,n9,n0; 153 154 155 156vecdim n; 157 158 159 160 g(l,n1,n1); 161 162 163n 164 165 g(l,n1,m1,n1,m1); 166 167 168n*( - n + 2) 169 170 171 g(l,n1,n2,n2,n1); 172 173 174 2 175n 176 177 g(l,n1,n2,m1,n2,n1,m1); 178 179 180 2 181n*(n - 4*n + 4) 182 183 g(l,n1,n2,m1,m2,n2,n1,m2,m1); 184 185 186 3 2 187n*(n - 8*n + 24*n - 16) 188 189 190 g(l,n1,n2,n3,n3,n2,n1); 191 192 193 3 194n 195 196 g(l,n1,n2,n3,m1,n3,n2,n1,m1); 197 198 199 3 2 200n*( - n + 6*n - 12*n + 8) 201 202 g(l,n1,n2,n3,m1,m2,n3,n2,n1,m2,m1); 203 204 205 4 3 2 206n*(n - 12*n + 60*n - 112*n + 64) 207 208 g(l,n1,n2,n3,m1,m2,m3,n3,n2,n1,m3,m2,m1); 209 210 211 5 4 3 2 212n*( - n + 18*n - 144*n + 528*n - 816*n + 416) 213 214 g(l,n1,n2,n3,m1,n3,n1,n2,m1); 215 216 217 3 2 218n*(n - 8*n + 16*n - 8) 219 220 g(l,n1,n2,n3,m1,m2,n3,n1,n2,m1,m2); 221 222 223 4 3 2 224n*(n - 16*n + 72*n - 120*n + 64) 225 226 g(l,n1,n2,n3,m1,m2,m3,n2,n3,n1,m3,m1,m2); 227 228 229 5 4 3 2 230n*( - n + 22*n - 172*n + 584*n - 848*n + 416) 231 232 233 234% COPYRIGHT (C) 1988,1990, INSTITUTE OF NUCLEAR PHYSICS, MOSCOW STATE U. 235% CVITWEEL TEST OF CVITANOVIC PACKAGE 236% AUTHOR A. KRYUKOV, ARODIONOV, A.TARANOV 237% VERSION 1.2 238% RELEASE 11-MAR-90 239 240index j1,j2,j3,j4,j5,j6,j7,j8,j9,j0; 241 242 243vecdim n$ 244 245 246 247% Test of CVITANOVIC PACKAGE 248% 249% All tests are the weels with defferent spoke 250% (Use notation from SIGSAM Bull, 1989, v.23, no.4, pp.15-24) 251 252g(l,j1,j2,j2,j1); 253 254 255 2 256n 257 258 259g(l,j1,j2,j3,j1,j2,j3); 260 261 262 2 263n*( - n + 6*n - 4) 264 265g(l,j1,j2,j3,j1,j3,j2); 266 267 268 2 269n*(n - 4*n + 4) 270 271g(l,j1,j2,j3,j3,j2,j1); 272 273 274 3 275n 276 277 278g(l,j1,j2,j3,j4,j1,j2,j3,j4); 279 280 281 3 2 282n*(n - 12*n + 28*n - 16) 283 284g(l,j1,j2,j3,j4,j1,j2,j4,j3); 285 286 287 3 2 288n*( - n + 10*n - 24*n + 16) 289 290g(l,j1,j2,j3,j4,j1,j4,j2,j3); 291 292 293 3 2 294n*(n - 8*n + 16*n - 8) 295 296g(l,j1,j2,j3,j4,j1,j4,j3,j2); 297 298 299 3 2 300n*( - n + 6*n - 12*n + 8) 301 302 303g(l,j1,j2,j3,j4,j5,j1,j2,j3,j4,j5); 304 305 306 4 3 2 307n*(n - 20*n + 100*n - 160*n + 80) 308 309g(l,j1,j2,j3,j4,j5,j1,j2,j3,j5,j4); 310 311 312 4 3 2 313n*( - n + 18*n - 88*n + 152*n - 80) 314 315g(l,j1,j2,j3,j4,j5,j1,j2,j5,j3,j4); 316 317 318 4 3 2 319n*(n - 16*n + 72*n - 120*n + 64) 320 321g(l,j1,j2,j3,j4,j5,j1,j2,j5,j4,j3); 322 323 324 4 3 2 325n*( - n + 14*n - 60*n + 112*n - 64) 326 327 328g(l,j1,j2,j3,j4,j5,j6,j1,j2,j3,j4,j5,j6); 329 330 331 5 4 3 2 332n*( - n + 30*n - 260*n + 840*n - 1120*n + 512) 333 334g(l,j1,j2,j3,j4,j5,j6,j1,j2,j3,j4,j6,j5); 335 336 337 5 4 3 2 338n*(n - 28*n + 236*n - 784*n + 1088*n - 512) 339 340g(l,j1,j2,j3,j4,j5,j6,j1,j2,j4,j3,j6,j5); 341 342 343 5 4 3 2 344n*( - n + 26*n - 216*n + 736*n - 1056*n + 512) 345 346g(l,j1,j2,j3,j4,j5,j6,j1,j2,j6,j3,j4,j5); 347 348 349 5 4 3 2 350n*(n - 24*n + 176*n - 536*n + 704*n - 320) 351 352g(l,j1,j2,j3,j4,j5,j6,j7,j1,j2,j3,j4,j5,j6,j7); 353 354 355 6 5 4 3 2 356n*( - n + 42*n - 560*n + 3080*n - 7840*n + 9184*n - 3904) 357 358g(l,j1,j2,j3,j4,j5,j6,j7,j1,j2,j3,j4,j5,j7,j6); 359 360 361 6 5 4 3 2 362n*(n - 40*n + 520*n - 2880*n + 7520*n - 9024*n + 3904) 363 364g(l,j1,j2,j3,j4,j5,j6,j7,j2,j1,j3,j4,j5,j6,j7); 365 366 367 6 5 4 3 2 368n*(n - 40*n + 520*n - 2880*n + 7520*n - 9024*n + 3904) 369 370 371% Test of example that calculated incorrectly in earlier package. 372 373index ix,iy,iz; 374 375 376mass p1=mm, p2=mm, p3=mm, p4=mm, k1=0; 377 378 379mshell p1,p2,p3,p4,k1; 380 381 382vector q1,q2; 383 384 385operator ga,gb; 386 387 388for all p let ga(p)=g(la,p) + mm, 389 gb(p)=g(lb,p) + mm; 390 391 392xx := g(la,ix)*g(la,iy)*(g(lb,ix)*gb(p1)*g(lb,iy)*gb(q2) + 393 gb(p3)*g(lb,ix)*g(lb,iy)); 394 395 396 2 397xx := - p1.q2*n + 2*p1.q2 + mm *n + mm*n 398 399let q1=p1-k1, q2=p3+k1; 400 401 402 403xx; 404 405 406 2 407 - k1.p1*n + 2*k1.p1 - p1.p3*n + 2*p1.p3 + mm *n + mm*n 408 409 410end; 411 412Tested on x86_64-pc-windows CSL 413Time (counter 1): 16 ms 414 415End of Lisp run after 0.01+0.06 seconds 416real 0.21 417user 0.00 418sys 0.06 419