MathTensor Information (Long)

*To*: mathgroup at yoda.physics.unc.edu*Subject*: MathTensor Information (Long)*From*: Steven M. Christensen <stevec>*Date*: Mon, 28 Jun 93 02:43:17 EDT*Organization*: MathSolutions, Inc.

[Due to a mailer glitch, this message was accidentally appended to a previous message and did not get to the entire list. Sorry if you get two copies.] Following is my periodic posting of information about MathTensor. Thanks go to members of MathGroup who have passed this information onto interested colleagues. Steve Christensen ******************************************************* MathTensor(tm) Information Sheet (Retail Pricing) ******************************************************* Prices Valid After June 15, 1993. MathSolutions offers discounts to academic institutions inside and outside the USA, dealers, and students. Contact MathSolutions for rates. See the end of this email for functionality new to version 2.1.5. *********************************************************** Hardware Requirements: MathTensor requires approximately 2 megabytes of disk space. It is recommended that a workstation have at least 8 megabytes of RAM memory. The Macintosh version runs on machines with 8 megabytes or more. DOS versions require at least 5 megabytes of memory, but more is highly recommended. The Windows 3.x version works with 8 megabytes of RAM or more. Software Requirements: MathTensor requires Mathematica 1.2, 2.0, or greater. Contact: MathSolutions, Inc. P.O. Box 16175 Chapel Hill, NC 27516 Telephone/Answering Machine/FAX: 919-967-9853 Email: mathtensor at wri.com * MathTensor Product Price List: (Prices include media/manual/shipping via first class mail, Upgrade contracts provide upgrades to MathTensor for one year - a minimum of two upgrades - from date of shipment.) Macintosh/DOS/Windows: No Upgrade $550.00 With Upgrade $700.00 Single Processor Workstations: No Upgrade $800.00 With Upgrade $1000.00 Convex/Multi-Processor Servers No Upgrade $3000.00 With Upgrade $4000.00 Extra Manuals: $30.00 Federal Express Delivery add: $15.00 (in US) $50.00 (outside US) Upgrades and updates are only sent to users with Upgrade Contracts. Volume discounts available. Consulting contracts for MathTensor and Mathematica users are also available. ** North Carolina residents add 6% sales tax ** Purchase orders are accepted. ** Payment terms - NET 15 ** Prices, availability and details of functionality of any products or services mentioned on this page are subject to change without notice. ** No credit card or letter of credit order Availability: MathTensor is available in 1/4" tape (QIC-24) tar format, MS-DOS 3.5" floppies, Macintosh 800K floppies, and Sun 3.5" bar format floppies. Other formats will be made available as needed but may require more time for shipment. Copy Protection: MathTensor, like Mathematica, requires a password on UNIX machines. This password must be obtained from MathSolutions, Inc. Authors: Leonard Parker, Ph.D. and Steven M. Christensen, Ph.D., founders of MathSolutions, Inc. -------------------- Order Form Cut Here --------------------------------- MathTensor 2.1.5 Ordering Information (June 15, 1993) Price Number Total o Personal Computers [ ] Macintosh (without upgrade, 800K 3.5" floppies) $550 _____ _______ [ ] Macintosh (with upgrade, 800K 3.5" floppies) $700 _____ _______ [ ] DOS/Windows (without upgrade, 1.44 Meg 3.5" floppies) $550 _____ _______ [ ] DOS/Windows (with upgrade, 1.44 Meg 3.5" floppies) $700 _____ _______ o Single Processor Workstations** [ ] NeXT Workstation (without upgrade, 1.44 Meg 3.5" DOS floppies) $800 _____ _______ [ ] NeXT Workstation (with upgrade, 1.44 Meg 3.5" DOS floppies) $1000 _____ _______ [ ] Sun Sparcstation (without upgrade, 3.5" Sun bar floppies) $800 _____ _______ [ ] Sun Sparcstation (with upgrade, 3.5" Sun bar floppies) $1000 _____ _______ [ ] UNIX Workstation (without upgrade, 1/4" QIC-24 tar tape) $800 _____ _______ [ ] UNIX Workstation (with upgrade, 1/4" QIC-24 tar tape) $1000 _____ _______ [ ] UNIX Workstation (without upgrade, 1/4" QIC-150 tar tape) $800 _____ _______ [ ] UNIX Workstation (with upgrade, 1/4" QIC-150 tar tape) $1000 _____ _______ Specify UNIX Workstation Type _______________________________ [ ] VMS Workstation (without upgrade, 1/4" QIC-24 tar tape) $800 _____ _______ [ ] VMS Workstation (with upgrade, 1/4" QIC-24 tar tape) $1000 _____ _______ Give the Mathematica $MachineID by typing $MachineID while in Mathematica __________________________ If your workstation for MathTensor does not have a $MachineID (usually if you have a network license), type hostid at the UNIX or VMS prompt and enter it here ______________________ o Convex/Multiprocessor Servers** [ ] Convex (without upgrade, 1/4" QIC-24 tar tape) $3000 _____ _______ [ ] Convex (with upgrade, 1/4" QIC-24 tar tape) $4000 _____ _______ [ ] Multiprocessor Server (without upgrade, 1/4" QIC-24 tar tape) $3000 _____ _______ [ ] Multiprocessor Server (with upgrade, 1/4" QIC-24 tar tape) $4000 _____ _______ Specify Machine Type _______________________________ Give the Mathematica $MachineID by typing $MachineID while in Mathematica __________________________ If your workstation for MathTensor does not have a $MachineID (usually if you have a network license), type hostid at the UNIX prompt and enter it here ______________________ [ ] Extra Manuals $30 _____ _______ [ ] Federal Express Delivery (give shipping address and phone number) $15 (US) _______ $50 (Outside US) _______ Subtotal _______ North Carolina residents add 6% sales tax _______ Order Total (Send PO or check in US Dollars) allow 4 weeks for delivery, no credit card orders, terms NET 15) _______ All prices include shipping and one manual. Prices are subject to change without notice. Upgrade gives the purchaser one year of MathTensor updates (minimum of 2). ** If your computer is on the Internet, it is also possible for us to do a remote installation. Contact MSI. o Shipping Information: Name: ________________________________________________ Address: __________________________________________________________ __________________________________________________________ City: ___________________________ State: _________________________ Postal/Zip Code: __________________ Country: _____________________ PO Number __________________________ Telephone: ________________________________ Fax: ________________________________ Email: ________________________________________ --------------------- End of Order Form ---------------------------------- Some of the users of MathTensor can be found at: Wake Forest University The University of Maryland Insituto de Fisica Fundamental - Madrid, Spain City College of New York The University of Winnipeg -Canada Stanford University Caltech The University of California - Santa Barbara Texas A&M University The University of North Carolina - Chapel Hill Universitat Konstanz - Germany Schlumberger KK - Japan Sumitomo Corporation - Japan Oakland University Louisiana Tech University Polaroid Lawrence Livermore Labs Los Alamos Labs The University of New Brunswick -Canada United Technologies Research Center Universitat Hannover - Germany Reed College Yale University The University of Bergen - Norway Cotton, Inc. NASA - Langley NASA - Goddard NASA - JPL National Center for Supercomputing Applications Wolfram Research The University of Wisconsin - Milwaukee The Hebrew University of Jerusalem -Israel Martin Marietta Cornell University Central Connecticut State University Universite de Liege - Belgium Pune University - India Cal State University - Fullerton Queen's University - Canada Utah State University Southwestern University - Texas University of South Florida Deutsches Klimarenzentrum - Hamburg, Germany University of Oklahoma - Norman Montana State University University of Chicago University of Milan - Italy Wellesley College University of San Francisco University of Cologne - Germany University of California - San Diego University of Berne - Switzerland University d. Bundeswehr - Hamburg, Germany University of Thessaloniki - Greece University of California - San Diego Royal Insitute of Technology (KTH) - Sweden Univeritat d. Bundeswehr - Hamburg University of Cologne - German Hamamatsu Photonics - Japan Linkoeping University - Sweden University of Washington - Seattle University of Cincinnati W.R. Grace, Co. Aberdeen Proving Ground Nissin-High Voltage Co, Ltd. - Japan Ibaraki University - Japan University of California - Berkeley CNR-FISBAT - Italy Austin College Mercyhurst College Universidad De Puerto Rico Unviversity of Alberta University of Oregon NeXT Computer Corporation Queen's University - Canada Tokyo Institute of Technology - Japan Univ. Autonoma Metropolitana - Mexico S3IS - France PICA Software - Australia RITME Informatique - France Mount Sinai Hospital - New York IVIC - Venezuela Carnegie Mellon University Choong Puk National University - Korea Microsoft, Inc. University of Texas - Austin General Tire and many others ....... MathTensor has been mentioned in articles in MacWorld, the Mathematica Journal, Computers In Physics, Science, and PC Magazine. "The idea of packages and special tools has been taken to a kind of modern limit in the form of MathTensor, a tensor- analysis package built on Mathematica." - Richard Crandall, Howard Vollum Professor of Science, Reed College and Chief Scientist, NeXT Computer, Inc., from Computers in Physics, Nov/Dec 1991. "I would like to take this opportunity to tell you how pleased I have been with MathTensor. I have been using it to study higher-dimensional cosmological solutions of Einstein's Equations. It has been an indispensable tool for me. You have done a great job." - Dwight Vincent, University of Winnipeg - Canada, Feb 1992. "Regarding my comments on MathTensor, you can be assured that they were coming out of my heart, and that they are TRUE. Yes, it is a very nice and useful tool. Many people, like myself, who are getting older and loosing the patience (and in some cases the ability to do long tedious computations) to deal with the "nitty gritty" of things can find here REAL HELP. Also, for those of us who are daring and curious it opens up a door to tackle problems that we otherwise wouldn't." Juan Perez Mercader, Insituto de Fisica Fundamental - Madrid, Spain, Feb 1992. "MathTensor is an outstanding package for doing tensor calculus. Best of all is the support. Questions receive an immediate reply by email from both Leonard and Steve, including advice which frequently goes far beyond the original point. Part of the MathTensor package is the continuing insight of two researchers who use this mathematics as their tool." - George Ruppeiner, Associate Professor of Physics, University of South Florida, Sept 1992. -----------------------General Information-------------------------- MathTensor Tensor analysis is extensively used in applications in physics, mathematics, engineering, and many other areas of scientific research. Problems involving tensors often are extraordinarily large and can be some of the most difficult computations in all of science. Equations with thousands of terms are common and can only be manipulated by computer mathematics systems like Mathematica. MathTensor is the largest Mathematica package yet developed outside of Wolfram Research. It adds over 250 new functions and objects to Mathematica to give the user both elementary and advanced tensor analysis functionality. MathTensor is a general tool for handling both indicial and concrete tensor indices. Standard objects like Riemann tensor, Ricci tensor, metric and others are built into the system along with common functions like the covariant derivative, index commutation, raising and lowering of indices, and various differential forms operations. MathTensor has been under development by Leonard Parker and Steven M. Christensen since the first alpha test release of Mathematica. It contains over 25,000 lines of Mathematica code contained in nearly 100 files totalling approximately 2.0 Megabytes of disk space. MathTensor will run on any machine that runs Mathematica and has sufficient RAM memory (generally 8 Megabytes or more) and disk space for file storage and swap. MathTensor runs under versions 1.2 and 2.X of Mathematica. MathSolutions, Inc. was formed by Leonard Parker and Steven M. Christensen. MathSolutions has its office in Chapel Hill, North Carolina. Both Parker and Christensen are theoretical physicists, specializing in research in Einstein's Special and General Theory of Relativity, quantum field theory, black hole theory, and cosmology. Christensen is a Contributing Editor to the Mathematica Journal and is the founder and Moderator of the MathGroup, Mathematica mailing list. Parker obtained his Ph.D. in Physics from Harvard University and Christensen his Ph.D. in Physics from the University of Texas at Austin. MathTensor is a trademark of MathSolutions, Inc. Mathematica is a registered trademark of Wolfram Research, Inc. ----------------------MathTensor Examples ---------------------------- MathTensor provides commands for simplifying and manipulating tensor expressions, as well as a knowledge base of transformation rules and definitions required for dealing with some of the more important tensors. MathTensor is designed to work along with the functions of Mathematica to provide users with the functions and objects they need to devise their own custom tensor analysis programs. MathTensor provides most of the basic structures needed for doing tensor computations and for programming new functions. As MathTensor is used and special applications are developed, they will be added to future versions of MathTensor. The following pages give a few very elementary examples of how MathTensor works. (* First, after starting Mathematica 2.X, we load a file which in turn loads many other files containing the MathTensor function and object definitions. *) In[1]:= <<MathTensor.m ====================================================== MathTensor (TM) 2.1.5 (UNIX (R)) (January 1, 1993) by Leonard Parker and Steven M. Christensen Copyright (c) 1991-1993 MathSolutions, Inc. Runs with Mathematica (R) Versions 1.2, 2.0, or 2.1 Licensed to machine sunny. ====================================================== No unit system is chosen. If you want one, you must edit the file called Conventions.m, or enter a command to interactively set units. Units: {} Sign conventions: Rmsign = 1 Rcsign = 1 MetricgSign = 1 DetgSign = -1 TensorForm turned on, ShowTime turned off, MetricgFlag = True. ========================================= (* MathTensor has a number of tensors already defined. Lower indices are entered with an "l" in front, while upper indices are entered with a "u" in front. In output lines the indices are placed where they should be on the tensor. *) In[2]:= RiemannR[la,lb,lc,ld] Out[2]= R abcd In[3]:= RiemannR[la,ub,lc,ud] b d Out[3]= R a c (* Operations like covariant differentiation are done simply with the CD function, with the derivative index placed after the tensor. *) In[4]:= CD[RicciR[la,lb],uc] c Out[4]= R ab; (* The standard summation convention is recognized. Here indices are summed using the standard Mathematica substitution operation which in the case below simply renames the upper index. *) In[5]:= % /. uc->ub b Out[5]= R ab; (* The DefineTensor function permits you to define your own tensors. The input name of the object defined below is "tensor", and its print name, which appears in output lines, is "t". The last argument indicates that it will have two indices, which upon interchange result in multiplication of the object by a weight factor of 1 -- that is, it is a symmetric tensor. *) In[6]:= DefineTensor[tensor,"t",{{2,1},1}] PermWeight::sym: Symmetries of t assigned PermWeight::def: Object t defined In[7]:= tensor[la,lb] Out[7]= t ab (* MathTensor now automatically reorders symmetric indices into lexical order. *) In[8]:= tensor[lb,la] Out[8]= t ab (* The standard symmetries of the Riemann tensor are built into its definition. *) In[9]:= RiemannR[lb,la,lc,ld] Out[9]= -R abcd In[10]:= RiemannR[lc,ld,la,lb] Out[10]= R abcd (* MathTensor knows that the appropriate sum of indices on the Riemann tensor gives the Ricci tensor. *) In[11]:= RiemannR[la,lb,ua,lc] Out[11]= R bc (* The sum of the indices on the Ricci tensor gives the Riemann Scalar. *) In[12]:= RicciR[la,ua] Out[12]= R (* Now we define a new tensor with four indices and no symmetries. *) In[13]:= DefineTensor[T,"T",{{1,2,3,4},1}] PermWeight::sym: Symmetries of T assigned PermWeight::def: Object T defined (* We produce a complicated product of seven of these tensor with multiple summations of indices and then add it to another similar object. *) In[14]:= SevenTensorTest := T[la,lb,uc,ud] T[lc,ld,ue,uf] T[le,lf,ug,uh] T[lg,li,ui,uj] * T[lh,lj,uk,ul] T[lk,ll,um,un] T[lm,ln,ua,ub] - T[la,lb,uc,ud] T[lc,le,ue,uf] * T[ld,lf,ug,uh] T[lg,lh,ui,uj] T[li,lj,uk,ul] T[lk,ll,um,un] T[lm,ln,ua,ub] (* Trying to find some simplification of SevenTensorTest is not easy by hand. *) In[15]:= SevenTensorTest cd ef gh ij kl mn ab Out[15]= T T T T T T T - ab cd ef gi hj kl mn cd ef gh ij kl mn ab > T T T T T T T ab ce df gh ij kl mn (* But the MathTensor command Tsimplify rapidly finds that the two terms are equal. *) In[16]:= Tsimplify[%] Out[16]= 0 (* Some terms differ only by the renaming of summation indices. *) In[17]:= T[la,lb,lc,ld] RiemannR[ua,ub,uc,ud] + T[le,lf,lg,lh] RiemannR[ue,uf,ug,uh] abcd efgh Out[17]= R T + R T abcd efgh (* MathTensor's canonicalization functions can rename indices and combine terms. *) In[18]:= Canonicalize[%] pqrs Out[18]= 2 R T pqrs (* MathTensor can symmetrize or antisymmetrize pairs of indices. *) In[19]:= Symmetrize[T[la,lb,lc,ld], {la,lb}] T + T abcd bacd Out[19]= ------------- 2 In[20]:= Expand[%] T T abcd bacd Out[20]= ----- + ----- 2 2 In[21]:= Antisymmetrize[%,{lc,ld}] T T T T abcd abdc bacd badc ----- - ----- + ----- - ----- 2 2 2 2 Out[21]= ----------------------------- 2 In[22]:= Expand[%] T T T T abcd abdc bacd badc Out[22]= ----- - ----- + ----- - ----- 4 4 4 4 (* MathTensor understands how to convert covariant to ordinary partial derivatives with affine connection terms added. *) In[23]:= CD[RicciR[la,lb],lc] Out[23]= R ab;c In[24]:= CDtoOD[%] p p Out[24]= R - G R - G R ab,c bc pa ac pb (* Using positive and negative index values, MathTensor can deal with concrete contravariant or covariant indices. *) In[25]:= RiemannR[1,2,3,4] 1234 Out[25]= R In[26]:= RiemannR[-1,2,-3,4] 4 2 Out[26]= R 3 1 (* We can set the dimension of the spacetime to some value, like 4. *) In[27]:= Dimension = 4 Out[27]= 4 (* Then using the MakeSum function, we can explicitly write out sums in terms of concrete indices. *) In[28]:= MakeSum[RicciR[la,lb] RicciR[lc,ub]] 1 2 3 4 Out[28]= R R + R R + R R + R R 1a c 2a c 3a c 4a c (* Now we define tensor T with two indices that are symmetric. *) In[29]:= DefineTensor[T,"T",{{2,1},1}] PermWeight::sym: Symmetries of T assigned PermWeight::def: Object T defined (* Then the lower components of T can be defined in terms of the components of other tensors like the Ricci tensor. *) In[30]:= SetComponents[T[la,lb],RicciR[la,lc] RicciR[lb,uc]] Components assigned to T (* We can ask for a specific covariant component of T. *) In[31]:= T[-1,-1] 1 2 3 4 Out[31]= R R + R R + R R + R R 11 1 21 1 31 1 41 1 (* One example application built into MathTensor does variations with respect to the metric tensor of structures that are functions of the metric tensor, Metricg. The variation of the square root of the determinant of the metric times the Riemann scalar gives terms in the variation of the metric, called h. *) In[32]:= Sqrt[Detg] ScalarR Out[32]= Sqrt[g] R In[33]:= Variation[%,Metricg] MetricgFlag::off: MetricgFlag is turned off by this operation pq p q Out33]= Sqrt[g] h - Sqrt[g] h + pq; p ;q pq Sqrt[g] R g h pq pq ----------------- - Sqrt[g] R h 2 pq (* A more complicated variation is just as easily found. *) In[34]:= Sqrt[Detg] RicciR[la,lb] RicciR[ua,ub] ab Out[34]= Sqrt[g] R R ab In[35]:= Variation[%,Metricg] pq r -(Sqrt[g] h R ) ;r pq p rq Out[35]= --------------------- + Sqrt[g] h R - 2 q; pr p qr r pq Sqrt[g] h R Sqrt[g] h R p ; qr pq;r ------------------ - ------------------ + 2 2 p qr Sqrt[g] h R q pr p ;qr Sqrt[g] h R - ------------------ + pq;r 2 pq rs Sqrt[g] g R R h rs pq qr p ----------------------- - Sqrt[g] R R h - 2 pq r pr q Sqrt[g] R R h pq r (* MathTensor's ApplyRules function permits the user to define large sets of rules that can be applied to expressions to simplify them. One set of rules that is provided as an example are the RiemannRules. Familiar rules are applied in the next two examples. MathTensor includes several functions DefUnique and RuleUnique that help the user devise their own rules and save them for later use. *) In[36]:= CD[RicciR[la,lb],ub] b Out[36]= R ab; In[37]:= ApplyRules[%,RiemannRules] R ;a Out[37]= --- 2 In[38]:= RiemannR[la,lb,lc,ld] RiemannR[ua,uc,ub,ud] acbd Out[38]= R R abcd In[39]:= ApplyRules[%,RiemannRules] pqrs R R pqrs Out[39]= ----------- 2 (* Suppose we define a tensor with four indices. *) In[40]:= DefineTensor[tensor,"t",{{1,2,3,4},1}] PermWeight::sym: Symmetries of t assigned PermWeight::def: Object t defined (* MathTensor's multiple index facility permits the user to add indices with one, two or three primes to be used as extra non-spacetime indices. *) In[41]:= AddIndexTypes In[42]:= tensor[ala,blb,clc,ld] Out[42]= t a'b''c'''d (* Tools for building rules involving these extra indices are provided. Future releases of MathTensor will extend this functionality. *) (* Shipped with MathTensor is the file Components.m which may be run separately from MathTensor. Components.m takes a file like SchwarzschildIn.m, listed next, and computes the components of the affine connection, the Riemann tensor, Ricci tensor, Riemann scalar, Weyl tensor, Einstein tensor and several other objects. SchwarzschildIn.m is the input file containing information about the famous Schwarzschild metric. This metric represents the curved space of a gravitating spherically symmetric object with mass M in otherwise empty space. *) ---------------------- SchwarzschildIn.m file listing -------------------- (* Copyright (c) 1992 MathSolutions, Inc.*) (* SchwarzschildIn.m *) Dimension = 4 x/: x[1] = r x/: x[2] = theta x/: x[3] = phi x/: x[4] = t Metricg/: Metricg[-1, -1] = (1 - (2*G*M)/r)^(-1) Metricg/: Metricg[-2, -1] = 0 Metricg/: Metricg[-3, -1] = 0 Metricg/: Metricg[-4, -1] = 0 Metricg/: Metricg[-2, -2] = r^2 Metricg/: Metricg[-3, -2] = 0 Metricg/: Metricg[-4, -2] = 0 Metricg/: Metricg[-3, -3] = r^2*Sin[theta]^2 Metricg/: Metricg[-4, -3] = 0 Metricg/: Metricg[-4, -4] = -(1 - (2*G*M)/r) Rmsign = 1 Rcsign = 1 CalcEinstein = 1 CalcRiemann = 1 CalcWeyl = 1 SetOptions[Expand, Trig->True] SetOptions[Together, Trig->True] CompSimp[a_] := Expand[Together[a/.CompSimpRules[1] ] ] CompSimpRules[1] = {} (* End of file SchwarzschildIn.m *) --------------------------------------------------------------------- (* If we load Components.m into Mathematica we can run the following command to produce from SchwarzschildIn.m two new files, SchwarzschildOut.m, which contains results that can be used immediately in MathTensor, and SchwarzschildOut.out which can be printed. The computation below takes from just a few seconds to less than one minute on typical workstations. More complex metrics can take a bit longer. *) In[1]:= <<Components.m ==================================================== MathTensor (TM) 2.1.5 (UNIX (R)) (January 1, 1993) Components Package by Leonard Parker and Steven M. Christensen Copyright (c) 1991-1993 MathSolutions, Inc. Runs with Mathematica (R) Versions 1.2, 2.0, and 2.1. Licensed to machine sunny. ==================================================== In[2]:= Components["SchwarzschildIn.m","SchwarzschildOut.m", "SchwarzschildOut.out"] The following tensors have been calculated and stored in the file SchwarzschildOut.m in InputForm, and in the file SchwarzschildOut.out in OutputForm: Metricg MatrixMetricgLower MatrixMetricgUpper Detg AffineG[ua,lb,lc] RicciR[la,lb] ScalarR EinsteinG[la,lb,lc,ld] RiemannR[la,lb,lc,ld] WeylC[la,lb,lc,ld] You can edit SchwarzschildOut.out to print a record of the results. (* The listing of the output file is too long to give here, but it gives the calculated values of the components of the Riemann and related tensors in the Schwarzschild spacetime. *) The names of the various commands and other objects available in MathTensor is given below: A List of MathTensor Functions and Objects (Pre-2.1.5, See 2.1.5 Announcement for details of new functions.) Absorb Epsilon InvertFast Absorbg EpsilonProductTensor Kdelta AbsorbKdelta EpsilonProductTensorRule KdeltaRule AbsorbRule EpsilonSign Kill AddIndexTypes EpsilonToEpsDownRule Lap AffineG EpsilonToEpsUpRule LieD AffineToMetric EpsUp LieDtoCD AffineToMetricRule EpsUpToEpsDownRule LieDtoCDrule AllSymmetries EpsUpToEpsilonRule Lightc Antisymmetrize EqApart LorentzGaugeRule ApplyRules EqApply Lower ApplyRulesRepeated EqCancel Lowera Arglist EqCollect LowerAllPairs ArglistAllTypes EqDivide LowerAllTypes AskSignsFlag EqExpand Lowerb AskSignsProcedure EqExpandAll Lowerc BianchiFirstPairRule EqFactor LowerIndexAllTypesQ BianchiSecondPairRule EqFactorTerms LowerIndexaQ CanAll EqMinus LowerIndexbQ CanApplyRules EqPlus LowerIndexcQ CanApplyRulesFast EqPower LowerIndexQ CanApplyRulesRepeated EqReverse MakeAllSymmetries CanDum EqSimplify MakePermWeightGroup Cannn EqSolve MakeSum CannnDum EqSubtract Matchlist CanNonInvert EqTimes Matchlista Canonicalize EqTogether MatchlistAllTypes CanSame EqTwoDivide Matchlistb CanSuperApplyRules EqTwoPlus Matchlistc CD EqTwoSubtract MatchlistOrd CDtoOD EqTwoTimes $MathTensorVersionNumber ClearComponents esuUnits MaxwellA ClearUnits EvaluateODFlag MaxwellB CommuteCD Evenlist MaxwellCyclicEquation Components Explode MaxwellCyclicRule CountNewDums FirstCubicRiemannRule MaxwellDivergenceEquation CoXD FirstQuadraticRiemannRule MaxwellDivergenceRule DefineForm FreeList MaxwellE DefineTensor FtoC MaxwellF DefUnique FtoCrule MaxwellJ Detg GaussianUnits Maxwellk1 DetgSign GenLap Maxwellk3 Dimension GravitationalUnits Maxwellrho Downdummylist hbar MaxwellT Downlist HeavisideLorentzUnits MaxwellTexpression Downuserlist HodgeStar MaxwellTtoFrule DualStar Implode MaxwellVectorPotentialRule Dum IndexAllTypesQ Metricg DumAllTypes IndexaQ MetricgFlag EinsteinG IndexbQ MetricgFlagOff EinsteinToRicciRule IndexcQ MetricgFlagOn emuUnits IndexQ MetricgSign Eps0 IndexTypes Mu0 EpsDown IndicesAndNotOrderedQ NaturalUnits EpsDownToEpsilonRule InListQ NegIntegerQ EpsDownToEpsUpRule Invert NewtonG NonTensorPart Symmetries OD SymmetriesOfSymbol Oddlist Symmetrize OrderedArgsQ SyntaxCheck Pair SyntaxCheckOff PairAllTypes SyntaxCheckOn PairAllTypesQ TensorForm PairaQ TensorPart PairbQ TensorPartSameQ PaircQ TensorQ Pairdum TensorSimp PairQ TensorSimpAfter PermWeight TraceFreeRicciR PIntegrate TraceFreeRicciToRicciRule PosIntegerQ Tsimplify PrettyOff TsimplifyAfter PrettyOn Units Raise Unlist Raisea Updowndummylist RaiseAllTypes Updummylist Raiseb Uplist Raisec UpLo RankForm UpLoa RationalizedGaussianUnit UpLob RationalizedMKSUnits UpLoc Rcsign UpperIndexAllTypesQ RicciR UpperIndexaQ RicciSquared UpperIndexbQ RicciToAffine UpperIndexcQ RicciToAffineRule UpperIndexQ RicciToTraceFreeRicciRule Upuserlist RiemannCyclicFirstThreeRule Var RiemannCyclicSecondThreeRule Varg RiemannR Variation RiemannRules VariationalDerivative RiemannSquared VectorA RiemannToAffine VectorAFlag RiemannToAffineRule WeylC RiemannToWeylRule WeylToRiemannRule Rmsign XD Rulelists XDtoCDflag RuleUnique XP RuleUniqueAllTypes ZeroFormQ ScalarR ScalarRtoAffine ScalarRtoAffineRule SecondCubicRiemannRule SecondQuadraticRiemannRule SetAntisymmetric SetComponents SetSymmetric ShowNumbers ShowTime SIUnits SuperApplyRules SwapDum o New in version 2.1.5: **************** * Ttransform * **************** Ttransform is an easy-to-use function for performing general coordinate transformations on tensorial objects. Here is an example from gravitational physics: Example: General Relativity --- Kruskal coordinates Transform the Schwarzschild metric of a black hole from Schwarzschild to Kruskal coordinates. This transformation does not involve the angular coordinates, so it is sufficient to work in a two dimensional spacetime. In Schwarzschild coordinates r, t, the Schwarzschild line element is ds^2 = (1-(2 M)/r)^(-1) dr^2 - (1-(2 M)/r)) dt^2 . Kruskal gave the transformation X = (r/(2 M)-1)^(1/2) exp[r/(4 M)] cosh[t/(4 M)] T = (r/(2 M)-1)^(1/2) exp[r/ 4 M)] sinh[t/(4 M)] to Kruskal coordinates X, T. Use Ttransform to calculate the form of the metric tensor in Kruskal coordinates. (* Clear any previous assignments. *) In[2]:= Clear[ga,gb] (* Define the dimension. *) In[3]:= Dimension = 2 Out[3]= 2 (* Define the coordinates. *) In[4]:= coords = {r,t} Out[4]= {r, t} (* Define the transformations giving the kruskal coordinates {X,T} in terms of {r,t}. *) In[5]:= trans = {(r/(2 M) - 1)^(1/2) Exp[r/(4 M)] Cosh[t/(4 M)], (r/(2 M) - 1)^(1/2) Exp[r/(4 M)] Sinh[t/(4 M)]} r/(4 M) r t r/(4 M) r t Out[5]= {E Sqrt[-1 + ---] Cosh[---], E Sqrt[-1 + ---] Sinh[---]} 2 M 4 M 2 M 4 M (* Let ga[la,lb] be the Schwarzschild metric for r,t. The independent covariant components of the metric are assigned below. *) In[6]:= DefineTensor[ga,{{2,1},1}] (* Set -2, -2 component. *) In[7]:= ga[-2,-2] = -(1 - (2 M)/r) 2 M Out[7]= -1 + --- r (* Set -2, -1 component. *) In[8]:= ga[-2,-1] = 0 Out[8]= 0 (* Set -1, -1 component. *) In[9]:= ga[-1,-1] = (1 - (2 M)/r)^(-1) 1 Out[9]= ------- 2 M 1 - --- r (* gb will denote the transformed metric. *) In[10]:= DefineTensor[gb,{{2,1},1}] Calculate gb from ga. The transformation of coordinates and metric are in the same direction, from Schwarzschild to Kruskal, so the last argument of Ttransform is 1. The metric components of gb[la,lb] refer to the Kruskal coordinates X, T, but they will appear as functions of r, t, since ga[la,lb] and the right-hand-sides of the transformation equations are given as functions of r, t. In fact, that is the usual way that the Kruskal metric is written, since r and t cannot be expressed in a simple algebraic form in terms of X, T. In[11]:= Ttransform[gb,ga[la,lb],coords,trans,1] Components assigned to gb (* Display the result. *) In[12]:= Table[gb[-i,-j],{i,2},{j,2}] 3 3 32 M -32 M Out[12]= {{----------, 0}, {0, ----------}} r/(2 M) r/(2 M) E r E r This result means that the line element in Kruskal coordinates has the well-known form, ds^2 = 32 M^3 r^(-1) Exp[-r/(2 M)] (dX^2 - dT^2). ********************** * OrderCD * ********************** In previous versions of MathTensor there has been only one function that could be used to reorder covariant derivatives on a tensor object, CommuteCD. A number of users of MathTensor have asked for a function that will reorder all the indices into alphabetical order at once. The new OrderCD function does this. (* Define a two index tensor. *) In[2]:= DefineTensor[A,"A",{{1,2},1}] (* Take several covariant derivatives. In this case they are not in alphabetical order. *) In[3]:= CD[A[la,lb],lf,le,ld,lc] Out[3]= A ab;fedc (* The OrderCD function reorders the covariant derivatives and generates the appropriate Riemann tensor terms. This expression is still equal to A , but now contains A . All derivatives ab;fedc ab;cdef on all terms are put in alphabetical order as well. *) In[4]:= OrderCD[%] p p p p Out[4]= A R + A R + A R + A R + pb;f a cd;e pb;e a cd;f pb;d a ce;f pb;f a de;c p p p p > A R + A R + A R + A R + pb;c a de;f pb;e a df;c pb;d a ef;c pb;c a ef;d p p p p > A R + A R + A R + A R + ap;f b cd;e ap;e b cd;f ap;d b ce;f ap;f b de;c p p p p > A R + A R + A R + A R + ap;c b de;f ap;e b df;c ap;d b ef;c ap;c b ef;d p p p p > A R + A R + A R + A R + ab;p ced ;f ab;p dfe ;c pb a cd;ef pb a de;cf p p p p > A R + A R + A R + A R + A + pb a ef;cd ap b cd;ef ap b de;cf ap b ef;cd ab;cdef p p p p > A R + A R + A R + A R + pb;ef a cd pb;df a ce pb;de a cf pb;cf a de p p p p > A R + A R + A R + A R + pb;ce a df pb;cd a ef ap;ef b cd ap;df b ce p p p p > A R + A R + A R + A R + ap;de b cf ap;cf b de ap;ce b df ap;cd b ef p q p q p q p q > A R R + A R R + A R R + A R R + pq a ef b cd pq a df b ce pq a ce b df pq a cd b ef p q p q p q p q > A R R + A R R - A R R - A R R + pb qfa cde ap qfb cde pb qea cdf ap qeb cdf p pq pq p > A R + A R R + A R R + A R + ab;pf ced pb qadf ce ap qbdf ce ab;pe cfd p p q p q p > A R + A R R + A R R + A R + ab;pd cfe pb qda cfe ap qdb cfe ab;pc dfe pq pq > A R R + A R R pb qacd ef ap qbcd ef ********************** * Components * ********************** The Components function is used to compute the components of the Riemann tensor, Ricci tensor, Einstein tensor, affine connection, and Weyl tensor given the spacetime dimension and metric tensor. The latest version of Components permits you to set up the simplification rules you want to use to put the results of Components into the desired form. You can also treat metrics having a small parameter perturbatively by series expansion. As a simple example, we consider perturbations of flat spacetime using the series expansion feature of Components to perturb about a given metric. For a perfectly general perturbation it is probably simpler to work with symbolic indices within a MathTensor session. But for more special perturbations, it is convenient to use Components. In this example, the most general purely time-dependent perturbation about flat spacetime is considered in the linear approximation. The coordinates are called x1, x2, x3, and t. The ten components of the perturbation of the symmetric metric tensor are called h11[t], h21[t], h31[t], h41[t], h22[t], h32[t], h42[t], h33[t], h43[t], h44[t]. They are taken to depend only on the time coordinate t. The definition of CompSimp, which determines the simplification and other operations done by the Components function, expands to first order in a small parameter q which is introduced as a factor in each perturbation term of the metric. Here is the input file, which is called exIn.m. (* exIn.m *) Dimension = 4 x/: x[1] = x1 x/: x[2] = x2 x/: x[3] = x3 x/: x[4] = t Metricg/: Metricg[-1, -1] = 1 + q h11[t] Metricg/: Metricg[-2, -1] = q h21[t] Metricg/: Metricg[-3, -1] = q h31[t] Metricg/: Metricg[-4, -1] = q h41[t] Metricg/: Metricg[-2, -2] = 1 + q h22[t] Metricg/: Metricg[-3, -2] = q h32[t] Metricg/: Metricg[-4, -2] = q h41[t] Metricg/: Metricg[-3, -3] = 1 + q h33[t] Metricg/: Metricg[-4, -3] = q h43[t] Metricg/: Metricg[-4, -4] = -1 + q h44[t] Rmsign = 1 Rcsign = 1 CalcEinstein = 1 CalcRiemann = 1 CalcWeyl = 1 CompSimpRules[1] = {} SetOptions[Expand, Trig->True] SetOptions[Together, Trig->True] CompSimp[a_] := Together[Series[Expand[a/.CompSimpRules[1]], {q,0,1}]] (* end of file *) Suppose that in a Mathematica session you have loaded Components.m and issued the command Components["exIn.m","exOut.m","exOut.out"] to produce the file exOut.m suitable for input to a later MathTensor session, and the file exOut.out with output in a more human readable form. Here is an example MathTensor session which displays and manipulates the calculated tensor components. (* Load the output of the Components computation above. *) In[1]:= << exOut.m MetricgFlag has been turned off. (* Display the independent Ricci tensor components. *) In[2]:= Simplify[Table[RicciR[-i,-j],{i,4},{j,i}] ] h11''[t] q 2 Out[2]= {{---------- + O[q] }, 2 h21''[t] q 2 h22''[t] q 2 {---------- + O[q] , ---------- + O[q] }, 2 2 h31''[t] q 2 h32''[t] q 2 {---------- + O[q] , ---------- + O[q] , 2 2 h33''[t] q 2 ---------- + O[q] }, 2 2 2 2 {O[q] , O[q] , O[q] , (-h11''[t] - h22''[t] - h33''[t]) q 2 ----------------------------------- + O[q] }} 2 (* Here is the scalar curvature invariant. *) In[3]:= ScalarR 2 Out[3]= (h11''[t] + h22''[t] + h33''[t]) q + O[q] (* Here are three components of the Riemann tensor, namely, R_{1212}, R_{2323}, and R_{3434}. *) In[4]:= Simplify[Table[RiemannR[-i,-i-1,-i,-i-1],{i,1,3}] ] 2 2 -(h33''[t] q) 2 Out[4]= {O[q] , O[q] , ------------- + O[q] } 2 (* Here is the invariant square of the Ricci tensor. *) In[5]:= MakeSum[RicciR[la,lb] RicciR[ua,ub] ] 2 2 2 2 h11''[t] h21''[t] h22''[t] h31''[t] Out[5]= (--------- + --------- + --------- + --------- + 4 2 4 2 2 h32''[t] -h11''[t] h22''[t] h33''[t] 2 --------- + (--------- - -------- - --------) + 2 2 2 2 2 h33''[t] 2 3 ---------) q + O[q] 4 >From the displayed components of the Ricci tensor, it is clear that in vacuum, where all the components of the Ricci tensor must vanish, there is no purely time-dependent perturbation of Minkowski space which is zero at early times and non-zero at later times. On the other hand, there are non-zero vacuum perturbations which are linear in t, but these are not physical. ********************* * Faster Loading * ********************* On some computers with low memory or slow disks, MathTensor can be slow to load. In version 2.1.5, a special loading file called MathTensorFast.m has been added. This loads a smaller version of MathTensor which initially with contains fewer functions. If a function is needed that is not loaded by MathTensorFast.m, it is loaded when used the first time. This system has shown a speedup in loading of about 40% on most computers. ***************************************** * PairSymmetrize and PairAntisymmetrize * ***************************************** The standard MathTensor Symmetrize and Antisymmetrize functions permit symmetrization over a list of indices as in the following examples. (* Define some arbitrary four index tensor. *) In[2]:= DefineTensor[A,"A",{{1,2,3,4},1}] (* With the usual Symmetrize or Antisymmetrize functions, you can symmetrize over two indices. *) In[3]:= Symmetrize[A[la,lb,lc,ld],{la,lc}] A + A abcd cbad Out[3]= ------------- 2 (* You can also symmetrize or antisymmetrize over a list of indices. *) In[4]:= Symmetrize[A[la,lb,lc,ld],{la,lb,lc,ld}] A A A A A A A A abcd abdc acbd acdb adbc adcb bacd badc Out[4]= ----- + ----- + ----- + ----- + ----- + ----- + ----- + ----- + 24 24 24 24 24 24 24 24 A A A A A A A A A bcad bcda bdac bdca cabd cadb cbad cbda cdab > ----- + ----- + ----- + ----- + ----- + ----- + ----- + ----- + ----- + 24 24 24 24 24 24 24 24 24 A A A A A A A cdba dabc dacb dbac dbca dcab dcba > ----- + ----- + ----- + ----- + ----- + ----- + ----- 24 24 24 24 24 24 24 (* With PairSymmetrize you can manipulate pairs of indices. The resulting expression is symmetric under interchange of the index pair {a,b} with the pair {c,d}. *) In[5]:= PairSymmetrize[A[la,lb,lc,ld],{{la,lb},{lc,ld}}] A A abcd cdab Out[5]= ----- + ----- 2 2 (* PairAntisymmetrize does the expected thing. *) In[6]:= PairAntisymmetrize[A[la,lb,lc,ld],{{la,lb},{lc,ld}}] A A abcd cdab Out[6]= ----- - ----- 2 2 (* Define a new six index tensor. *) In[7]:= DefineTensor[B,"B",{{1,2,3,4,5,6},1}] (* You can manipulate more than one pair of indices. The resulting expression is antisymmetric under interchange of any of the indicated pairs of indices with each other, such as exchange of the pair {a,b} with either of the other two pairs, {c,d} or {e,f}. *) In[8]:= PairAntisymmetrize[B[la,lb,lc,ld,le,lf],{{la,lb},{lc,ld},{le,lf}}] B B B B B B abcdef abefcd cdabef cdefab efabcd efcdab Out[8]= ------- - ------- - ------- + ------- + ------- - ------- 6 6 6 6 6 6 (* These new functions can also work on triples and other sets of indices. *) In[9]:= PairSymmetrize[B[la,lb,lc,ld,le,lf],{{la,lb,lc},{ld,le,lf}}] B B abcdef defabc Out[9]= ------- + ------- 2 2 *********************** * Bug Fixes * *********************** All known MathTensor problems have been fixed in 2.1.5.

**update of Classes.m**

**Teaching numerical analysis**

**MathTensor Information (Long)**

**Axes ticks in Plot**