share/tensor/adsitter.dem

   1 /* Copyright (C) 2004 Viktor T. Toth <http://www.vttoth.com/>
   2  *
   3  * This program is free software; you can redistribute it and/or
   4  * modify it under the terms of the GNU General Public License as
   5  * published by the Free Software Foundation; either version 2 of
   6  * the License, or (at your option) any later version.
   7  *
   8  * This program is distributed in the hope that it will be
   9  * useful, but WITHOUT ANY WARRANTY; without even the implied
  10  * warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
  11  * PURPOSE.  See the GNU General Public License for more details.
  12  *
  13  * Classification of the anti-deSitter metric
  14  *
  15  */
  16
  17 if get('ctensor,'version)=false then load(ctensor);
  18 ("Computing the Petrov classification of the anti-deSitter metric")$
  19 ("We use frames")$
  20 cframe_flag:true;
  21 ("We set up the metric by hand")$
  22 ct_coords:[t,x,theta,phi];
  23 (lfg:ident(4),lfg[1,1]:-1,lfg)$
  24 fri:matrix([ 1/sqrt(2), -%i*a*cos(t/a)/sqrt(2),0,0],
  25            [-%i*1/sqrt(2),a*cos(t/a)/sqrt(2),0,0],
  26     [0,0,a*cos(t/a)*sinh(x)/sqrt(2), %i*a*cos(t/a)*sin(theta)*sinh(x)/sqrt(2)],
  27     [0,0,a*cos(t/a)*sinh(x)/sqrt(2),-%i*a*cos(t/a)*sin(theta)*sinh(x)/sqrt(2)]);
  28 cmetric(false);
  29 ug:invert(lg);
  30 weyl(false);
  31 nptetrad(false);
  32 ctrgsimp:true;
  33 psi(true);
  34 ("And obtain the Petrov class")$
  35 petrov();
  36
  37 /* End of demo -- comment line needed by MAXIMA to resume demo menu */