share/tensor/tracefree-code/formQ0-examples.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 ishow(canform('kdels([a,b,c],[j,k,l])*Q0(x([],[a])*x([],[b])*x([],[c]))))$
  10 /* output shown in paper:
  11                                          %1  %2 %3      j k l     2
  12                                       3 x   g      kdels         r
  13            %1  %2  %3      j k l                        %1 %2 %3
  14 (%t9)     x   x   x   kdels         - -----------------------------
  15                            %1 %2 %3                 5
  16  */
  17 ishow(canform('kdels([a,b,c,d],[j,k,l,m])
  18                      *Q0(x([],[a])*x([],[b])*x([],[c])*x([],[d]))))$
  19 /* output shown in paper:
  20           %1 %2  %3 %4      j k l m      4
  21        3 g      g      kdels            r
  22                             %1 %2 %3 %4
  23 (%t10) -----------------------------------
  24                        35
  25            %1  %2  %3 %4      j k l m      2
  26         6 x   x   g      kdels            r
  27                               %1 %2 %3 %4       %1  %2  %3  %4      j k l m
  28       - ------------------------------------ + x   x   x   x   kdels
  29                          7                                          %1 %2 %3 %4
  30  */