share/tensor/tracefree-code/formQ0-examples-additional.mac

   1 load("formQ0.mac")$
   2 load(itensor)$
   3 dim:3$
   4 decsym(g,0,2,[],[sym(all)])$
   5 imetric(g)$
   6 flipflag:true$
   7 defcon(x,x,r2)$
   8 components(r2([],[]),r^2)$
   9 showtime: true $
  10 ishow(canform('kdels([a,b,c],[j,k,l])*Q0(x([],[a])*x([],[b])*x([],[c]))))$
  11 time3: time(%);
  12 ishow(canform('kdels([a,b,c,d],[j,k,l,m])*Q0(x([],[a])*x([],[b])*x([],[c])*x([],[d]))))$
  13 time4: time(%);
  14 ishow(canform('kdels([a,b,c,d,e],[j,k,l,m,n])*Q0(x([],[a])*x([],[b])*x([],[c])*x([],[d])*x([],[e]))))$
  15 time5: time(%);
  16 ishow(canform('kdels([a,b,c,d,e,f],[j,k,l,m,n,o])*Q0(x([],[a])*x([],[b])*x([],[c])*x([],[d])*x([],[e])*x([],[f]))))$
  17 time6: time(%);
  18 ishow(canform('kdels([a,b,c,d,e,f,h],[j,k,l,m,n,o,p])*Q0(x([],[a])*x([],[b])*x([],[c])*x([],[d])*x([],[e])*x([],[f])*x([],[h]))))$
  19 time7: time(%);
  20 [time3, time4, time5, time6, time7];