share/tensor/rtest_atensor.mac

   1 (if get('atensor,'version)=false then load(atensor),done);
   2 done;
   3 init_atensor(clifford);
   4 done;
   5 atensimp(u.u);
   6 sf(u,u);
   7 atensimp(u.v-v.u);
   8 2*(u.v)-2*sf(u,v);
   9 init_atensor(symplectic);
  10 done;
  11 atensimp(v.u+u.v);
  12 2*(u.v)-2*af(u,v);
  13 atensimp(w.(u.v));
  14 -2*u*af(v,w)-2*v*af(u,w)+u.v.w;
  15 init_atensor(clifford,3);
  16 done;
  17 asymbol;
  18 v;
  19 atensimp(v[1].v[1]);
  20 1;
  21 init_atensor(quaternion);
  22 done;
  23 adim;
  24 2;
  25 aform;
  26 matrix([-1,0],[0,-1]);
  27 atensimp(v[1].v[1]);
  28 -1;
  29 atensimp(v[2].v[2]);
  30 -1;
  31 atensimp((v[1].v[2]).(v[1].v[2]));
  32 -1;
  33 init_atensor(lie_envelop);
  34 done;
  35 lbr(u,v):=u.v-v.u;
  36 lbr(u,v):=u.v-v.u;
  37 atensimp(lbr(u,lbr(v,w)))+atensimp(lbr(v,lbr(w,u)))+atensimp(lbr(w,lbr(u,v)));
  38 0;
  39
  40 (reset(dotscrules,dotexptsimp,dotdistrib),0);
  41 0$