REDUCE 3.4.1, 15-Jul-92 ...
1:
(CVIT)
% Tests of Cvitanovic Package.
% COPYRIGHT (C) 1990, INSTITUTE OF NUCLEAR PHYSICS, MOSCOW STATE UNIV.
% CVITBUBLE TEST OF CVITANOVIC PACKAGE
% AUTHOR A. KRYUKOV, ARODIONOV, A.TARANOV
% VERSION 1.1
% RELEASE 18-SEP-90
index j1,j2,j3,j4,j5,j6,j7,j8,j9,j0;
vecdim n$
% Tests of the weels with buble
% (Use notation from SIGSAM Bull, 1989, v.23, no.4, pp.15-24)
g(l,j1,j2,j2,j1);
2
N
g(l,j1,j2)*g(l1,j3,j1,j2,j3);
2
N
g(l,j1,j2)*g(l1,j3,j1,j3,j2);
N*( - N + 2)
g(l,j1,j2)*g(l1,j3,j3,j2,j1);
2
N
g(l,j1,j2,j3,j4)*g(l1,j1,j2,j3,j4);
N*(3*N - 2)
g(l,j1,j2)*g(l1,j3,j4,j1,j2,j4,j3);
3
N
g(l,j1,j2,j3,j4)*g(l1,j1,j4,j2,j3);
N*( - N + 2)
g(l,j1,j2)*g(l1,j3,j4,j1,j4,j3,j2);
2
N*(N - 4*N + 4)
g(l,j1,j2)*g(l1,j3,j4,j5,j1,j2,j3,j4,j5);
2 2
N *( - N + 6*N - 4)
g(l,j1,j2,j3,j4)*g(l1,j5,j1,j2,j3,j5,j4);
2
N*( - 3*N + 8*N - 4)
g(l,j1,j2,j3,j4,j5,j1)*g(l1,j2,j5,j3,j4);
2
N *( - N + 2)
g(l,j1,j2,j3,j4,j5,j1,j2,j5)*g(l1,j4,j3);
2 2
N *( - N + 6*N - 4)
g(l,j1,j2)*g(l1,j3,j4,j5,j6,j1,j2,j3,j4,j5,j6);
2 3 2
N *(N - 12*N + 28*N - 16)
g(l,j1,j2,j3,j4)*g(l1,j5,j6,j1,j2,j3,j4,j6,j5);
3
N *(3*N - 2)
g(l,j1,j2,j3,j4,j5,j6)*g(l1,j1,j2,j4,j3,j6,j5);
2
N*(7*N - 22*N + 16)
g(l,j1,j2,j3,j4,j5,j6,j1,j2)*g(l1,j6,j3,j4,j5);
2 2
N *( - 3*N + 8*N - 4)
g(l,j1,j2,j3,j4,j5,j6,j7,j1,j2,j3,j4,j5)*g(l1,j6,j7);
2 4 3 2
N *(N - 20*N + 100*N - 160*N + 80)
g(l,j1,j2,j3,j4,j5,j6,j7,j1,j2,j3)*g(l1,j4,j5,j7,j6);
2 3 2
N *(N - 8*N + 16*N - 8)
g(l,j1,j2,j3,j4,j5,j6,j7,j2)*g(l1,j1,j3,j4,j5,j6,j7);
3 2
N*( - 15*N + 60*N - 76*N + 32)
% COPYRIGHT (C) 1988,1990, INSTITUTE OF NUCLEAR PHYSICS, MOSCOW STATE U.
% CVITEST Test of CVITANOVIC PACKAGE
% AUTHOR A. KRYUKOV, A.RODIONOV, A.TARANOV
% VERSION 1.2
% RELEASE 11-MAR-90
%
% Test for trace of Dirac matrices.
%
% All tests are the lattices with difference lines
% (Use notation from SIGSAM Bull, 1989, v.4,no.23, pp.15-24)
index m1,m2,m3,m4,m5,m6,m7,m8,m9,m0;
index n1,n2,n3,n4,n5,n6,n7,n8,n9,n0;
vecdim n;
g(l,n1,n1);
N
g(l,n1,m1,n1,m1);
N*( - N + 2)
g(l,n1,n2,n2,n1);
2
N
g(l,n1,n2,m1,n2,n1,m1);
2
N*(N - 4*N + 4)
g(l,n1,n2,m1,m2,n2,n1,m2,m1);
3 2
N*(N - 8*N + 24*N - 16)
g(l,n1,n2,n3,n3,n2,n1);
3
N
g(l,n1,n2,n3,m1,n3,n2,n1,m1);
3 2
N*( - N + 6*N - 12*N + 8)
g(l,n1,n2,n3,m1,m2,n3,n2,n1,m2,m1);
4 3 2
N*(N - 12*N + 60*N - 112*N + 64)
g(l,n1,n2,n3,m1,m2,m3,n3,n2,n1,m3,m2,m1);
5 4 3 2
N*( - N + 18*N - 144*N + 528*N - 816*N + 416)
g(l,n1,n2,n3,m1,n3,n1,n2,m1);
3 2
N*(N - 8*N + 16*N - 8)
g(l,n1,n2,n3,m1,m2,n3,n1,n2,m1,m2);
4 3 2
N*(N - 16*N + 72*N - 120*N + 64)
g(l,n1,n2,n3,m1,m2,m3,n2,n3,n1,m3,m1,m2);
5 4 3 2
N*( - N + 22*N - 172*N + 584*N - 848*N + 416)
% COPYRIGHT (C) 1988,1990, INSTITUTE OF NUCLEAR PHYSICS, MOSCOW STATE U.
% CVITWEEL TEST OF CVITANOVIC PACKAGE
% AUTHOR A. KRYUKOV, ARODIONOV, A.TARANOV
% VERSION 1.2
% RELEASE 11-MAR-90
index j1,j2,j3,j4,j5,j6,j7,j8,j9,j0;
vecdim n$
% Test of CVITANOVIC PACKAGE
%
% All tests are the weels with defferent spoke
% (Use notation from SIGSAM Bull, 1989, v.23, no.4, pp.15-24)
g(l,j1,j2,j2,j1);
2
N
g(l,j1,j2,j3,j1,j2,j3);
2
N*( - N + 6*N - 4)
g(l,j1,j2,j3,j1,j3,j2);
2
N*(N - 4*N + 4)
g(l,j1,j2,j3,j3,j2,j1);
3
N
g(l,j1,j2,j3,j4,j1,j2,j3,j4);
3 2
N*(N - 12*N + 28*N - 16)
g(l,j1,j2,j3,j4,j1,j2,j4,j3);
3 2
N*( - N + 10*N - 24*N + 16)
g(l,j1,j2,j3,j4,j1,j4,j2,j3);
3 2
N*(N - 8*N + 16*N - 8)
g(l,j1,j2,j3,j4,j1,j4,j3,j2);
3 2
N*( - N + 6*N - 12*N + 8)
g(l,j1,j2,j3,j4,j5,j1,j2,j3,j4,j5);
4 3 2
N*(N - 20*N + 100*N - 160*N + 80)
g(l,j1,j2,j3,j4,j5,j1,j2,j3,j5,j4);
4 3 2
N*( - N + 18*N - 88*N + 152*N - 80)
g(l,j1,j2,j3,j4,j5,j1,j2,j5,j3,j4);
4 3 2
N*(N - 16*N + 72*N - 120*N + 64)
g(l,j1,j2,j3,j4,j5,j1,j2,j5,j4,j3);
4 3 2
N*( - N + 14*N - 60*N + 112*N - 64)
g(l,j1,j2,j3,j4,j5,j6,j1,j2,j3,j4,j5,j6);
5 4 3 2
N*( - N + 30*N - 260*N + 840*N - 1120*N + 512)
g(l,j1,j2,j3,j4,j5,j6,j1,j2,j3,j4,j6,j5);
5 4 3 2
N*(N - 28*N + 236*N - 784*N + 1088*N - 512)
g(l,j1,j2,j3,j4,j5,j6,j1,j2,j4,j3,j6,j5);
5 4 3 2
N*( - N + 26*N - 216*N + 736*N - 1056*N + 512)
g(l,j1,j2,j3,j4,j5,j6,j1,j2,j6,j3,j4,j5);
5 4 3 2
N*(N - 24*N + 176*N - 536*N + 704*N - 320)
g(l,j1,j2,j3,j4,j5,j6,j7,j1,j2,j3,j4,j5,j6,j7);
6 5 4 3 2
N*( - N + 42*N - 560*N + 3080*N - 7840*N + 9184*N - 3904)
g(l,j1,j2,j3,j4,j5,j6,j7,j1,j2,j3,j4,j5,j7,j6);
6 5 4 3 2
N*(N - 40*N + 520*N - 2880*N + 7520*N - 9024*N + 3904)
g(l,j1,j2,j3,j4,j5,j6,j7,j2,j1,j3,j4,j5,j6,j7);
6 5 4 3 2
N*(N - 40*N + 520*N - 2880*N + 7520*N - 9024*N + 3904)
% Test of example that calculated incorrectly in earlier package.
index ix,iy,iz;
mass p1=mm, p2=mm, p3=mm, p4=mm, k1=0;
mshell p1,p2,p3,p4,k1;
vector q1,q2;
operator ga,gb;
for all p let ga(p)=g(la,p) + mm,
gb(p)=g(lb,p) + mm;
xx := g(la,ix)*g(la,iy)*(g(lb,ix)*gb(p1)*g(lb,iy)*gb(q2) +
gb(p3)*g(lb,ix)*g(lb,iy));
2
XX := - P1.Q2*N + 2*P1.Q2 + N*MM + N*MM
let q1=p1-k1, q2=p3+k1;
xx;
2
- P1.P3*N + 2*P1.P3 - P1.K1*N + 2*P1.K1 + N*MM + N*MM
end;
Time: 3536 ms plus GC time: 136 ms
Quitting