← Back to branch summary

~ubuntu-branches/ubuntu/karmic/maxima/karmic

« back to all changes in this revision

Viewing changes to share/tensor/anew.mac

  • Committer: Bazaar Package Importer
  • Author(s): Camm Maguire
  • Date: 2004-11-13 18:39:14 UTC
  • mto: (2.1.2 hoary) (3.2.1 sid) (1.1.5 upstream)
  • mto: This revision was merged to the branch mainline in revision 3.
  • Revision ID: james.westby@ubuntu.com-20041113183914-ttig0evwuatnqosl
Tags: upstream-5.9.1
Import upstream version 5.9.1

Show diffs side-by-side

added added

removed removed

Lines of Context:
share/tensor/anew.mac
1
 
LOADFILE(ITENSR,FASL,MIKE)$DECLARE(E,CONSTANT)$DEFCON(E)$DEFCON(E,E,DELTA)$
2
 
METRIC:G$
3
 
G(L1,L2):=BLOCK([A,B,AL,BL,CL,DL],DUMMY(AL,BL,CL,DL),
4
 
IF L2=[] THEN [A:L1[1],B:L1[2],
5
 
RETURN(E(L1,[])+        L*P(L1,[])+     L^2*M(L1,[])+   L^3*N(L1,[])+
6
 
L^4*W(L1,[])+   L^5*Z(L1,[])
7
 
)],
8
 
A:L2[1],B:L2[2],
9
 
E([],L2)-       L*P([],L2)+
10
 
L^2*(-M([],L2)+P([AL],[A])*P([],[AL,B]))+
11
 
L^3*(-N([],L2)+P([AL],[A])*M([],[AL,B])+M([AL],[A])*P([],[AL,B])
12
 
-P([AL],[A])*P([BL],[AL])*P([],[BL,B]))+
13
 
L^4*(
14
 
-W([],L2)+
15
 
M([AL],[A])*M([],[AL,B])+
16
 
P([AL],[A])*N([],[AL,B])+
17
 
N([AL],[A])*P([],[AL,B])+
18
 
P([AL],[A])*P([BL],[AL])*P([CL],[BL])*P([],[B,CL])-
19
 
P([AL],[BL])*P([],[A,AL])*M([BL],[B])-
20
 
P([AL],[BL])*M([],[A,AL])*P([BL],[B])-
21
 
M([AL],[BL])*P([],[A,AL])*P([BL],[B])
22
 
)
23
 
+L^5*(-P([],[DL,A])*P([AL],[CL])*P([BL],[B])*P([CL],[BL])*P([DL],[AL])+
24
 
M([],[BL,CL])*P([AL],[A])*P([BL],[AL])*P([CL],[B])+
25
 
P([],[AL,BL])*M([AL],[CL])*P([BL],[B])*P([CL],[A])+
26
 
P([],[BL,CL])*M([AL],[A])*P([BL],[B])*P([CL],[AL])+
27
 
P([],[CL,A])*P([AL],[BL])*M([BL],[B])*P([CL],[AL])-
28
 
N([],[BL,A])*P([AL],[B])*P([BL],[AL])-
29
 
P([],[BL,A])*N([AL],[B])*P([BL],[AL])-
30
 
M([],[BL,A])*M([AL],[B])*P([BL],[AL])-
31
 
P([],[BL,A])*P([AL],[B])*N([BL],[AL])-
32
 
M([],[BL,A])*P([AL],[B])*M([BL],[AL])-
33
 
P([],[BL,A])*M([AL],[B])*M([BL],[AL])+
34
 
P([],[AL,A])*W([AL],[B])+
35
 
W([],[AL,A])*P([AL],[B])+
36
 
M([],[AL,A])*N([AL],[B])+
37
 
N([],[AL,A])*M([AL],[B])-
38
 
Z([],[B,A])
39
 
))$
40
 
SHOW(G([I,J],[]))$
41
 
SHOW(G([],[I,J]))$
42
 
BATCH(EXPAND,10000,TENSOR)$
43
 
RATWTLVL:5$
44
 
SHOW(CONTRACT(RENAME(RATEXPAND(G([I,J],[])*G([],[J,K])))))$
45
 
S(X):=CONTRACT(RENAME(CONTRACT(RENAME(CONTRACT(RENAME(X))))))$