share/tensor/ctensor6.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  * Algebraic classification of metrics
  14  *
  15  */
  16
  17 if get('ctensor,'version)=false then load(ctensor);
  18 ("Computing the Petrov classification of the Schwarzschild metric")$
  19 ("We use frames, set up simplification flags")$
  20 (cframe_flag:true,gcd:spmod,ctrgsimp:true,ratwtlvl:false,ratfac:true);
  21 ("We use the stock Schwarzschild frame")$
  22 ct_coordsys(exteriorschwarzschild,all);
  23 ("The inverse metric tensor is needed by nptetrad()")$
  24 ug:invert(lg)$
  25 ("The basis for the calculations is the Weyl tensor")$
  26 weyl(false);
  27 ("We compute a Newman-Penrose null tetrad")$
  28 nptetrad(false);
  29 ("Now we're ready to compute the coefficients")$
  30 psi(true);
  31 ("And obtain the Petrov class")$
  32 petrov();
  33
  34 /* End of demo -- comment line needed by MAXIMA to resume demo menu */