1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
NextTest("Tensors");
TestMathPiper(TSimplify( TSum({j}) Delta(i,j)*v(j) ),v(i));
TestMathPiper(TSimplify( TSum({j,i}) Delta(i,j)*Delta(i,j) ), Ndim);
TestMathPiper(TSimplify( TSum({j,i}) Delta(i,j)*Delta(j,i) ), Ndim);
TestMathPiper(TSimplify( TSum({j}) Delta(i,j)*Delta(j,k) ), Delta(i,k));
TestMathPiper(TSimplify( TSum({i}) v(i)*v(i) ), TSum({i})(v(i)^2));
Retract("v",1);
Rulebase("v",{ii});
f(i,j):=v(i)*v(j);
TestMathPiper(f(i,i),v(i)^2);
TestMathPiper(TSimplify( TSum({i}) f(i,i) ),TSum({i})(v(i)^2));
TestMathPiper(TSimplify( TSum({j}) Delta(i,j)*f(j,k) ),v(i)*v(k));
TestMathPiper(TSimplify(TSum({i,j}) Delta(i,j)*f(i,j) ), TSum({j})v(j)^2);
TestMathPiper(TSimplify(TSum({i})X(j)*TD(i)X(i)), Ndim*X(j));
TestMathPiper(TSimplify(TSum({i}) TD(i)(X(i)*X(j)) ), Ndim*X(j)+X(j));
TestMathPiper(TSimplify(TSum({i}) X(i)*TD(i)X(j) ), X(j));
TestMathPiper(TSimplify(TSum({i})TD(i)v(i)), TSum({i})TD(i)v(i));
TestMathPiper(TSimplify(TSum({i,j})TD(i)TD(j)(X(i)*X(j))), Ndim+Ndim^2);
TestMathPiper(TSimplify(TSum({i})TD(i)(X(i)*X(j)*X(j))), Ndim*X(j)^2+2*X(j)^2);
TestMathPiper(TSimplify(TSum({i,j,k})TD(i)TD(j)TD(k)(X(i)*X(j)*X(k))), 3*Ndim^2+2*Ndim+Ndim^3);
|
