Axisymmetric evolution code for the Einstein equations on hyperboloidal slices by Rinne.nb

(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 12.1' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     31953,        698]
NotebookOptionsPosition[     28124,        628]
NotebookOutlinePosition[     28590,        645]
CellTagsIndexPosition[     28547,        642]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["\<\
Notebook for :  
Axisymmetric evolution code for the Einstein equations on hyperboloidal \
slices by Rinne\
\>", "Title",
 CellGroupingRules->{"GroupTogetherGrouping", 10000.},
 CellChangeTimes->{{3.8217294718661222`*^9, 3.821729502043838*^9}, 
   3.8217295683553457`*^9, 3.821732221863006*^9, {3.821811203295624*^9, 
   3.8218112036548367`*^9}, {3.822041768762134*^9, 3.822041781290854*^9}, {
   3.822048009752109*^9, 3.8220480569607773`*^9}, {3.822504877727007*^9, 
   3.8225048874039307`*^9}, {3.822565447028533*^9, 3.822565457908236*^9}, {
   3.822822107930128*^9, 3.822822113128089*^9}, {3.823176459993638*^9, 
   3.823176470345335*^9}, {3.823359338529508*^9, 3.823359361520629*^9}, {
   3.823422938182802*^9, 3.823422938584877*^9}, {3.823431880221382*^9, 
   3.823431897644202*^9}, {3.82343570279051*^9, 3.823435714197551*^9}, {
   3.8279532410864677`*^9, 3.8279532445579853`*^9}, {3.827991128705386*^9, 
   3.827991132320956*^9}, {3.827992663416297*^9, 3.8279926679446497`*^9}, 
   3.828002571005631*^9, {3.828435487218956*^9, 3.828435494634301*^9}, {
   3.828455314498684*^9, 3.828455317626096*^9}, {3.82845929819172*^9, 
   3.828459314167211*^9}, {3.82846914237554*^9, 3.828469146910945*^9}, {
   3.828471572443717*^9, 3.8284715890511923`*^9}, {3.828529438373527*^9, 
   3.8285294469068527`*^9}, {3.828533473276997*^9, 3.8285334845643597`*^9}, {
   3.828605529086955*^9, 3.828605544687042*^9}, {3.8286904376033077`*^9, 
   3.828690443890785*^9}, {3.8287894528609467`*^9, 3.828789480877036*^9}, {
   3.828799096740452*^9, 3.828799108533012*^9}, {3.828814117044565*^9, 
   3.8288141187243*^9}, {3.828858846895076*^9, 3.82885884811847*^9}, {
   3.829979043091065*^9, 3.829979061081265*^9}, {3.830076319791185*^9, 
   3.830076322806036*^9}, {3.830081835983293*^9, 3.830081845630418*^9}, {
   3.8300821479108353`*^9, 3.830082150190091*^9}, {3.830083730839076*^9, 
   3.830083739182637*^9}, 3.8554992941324177`*^9, 3.8773025369493303`*^9},
 TextAlignment->Center,ExpressionUUID->"3d6666ab-568a-4851-baef-ae695b51f0cd"],

Cell["\<\
Geoff Cope
University of Utah
                                                                              \
                               \|01040f\|01042d\|01044c\|010432\|010442\
\|010432\|010449\|010445\|01042e\|01043b\|010428 \|010432\|010442 \|01040f\
\|01042d\|01043b\|01042b
May 15, 2021\
\>", "Author",
 CellGroupingRules->{"GroupTogetherGrouping", 10000.},
 CellChangeTimes->{{3.821729507652328*^9, 3.821729529657844*^9}, 
   3.8220480097522783`*^9, {3.822048070235491*^9, 3.822048084689715*^9}, {
   3.8225654346519403`*^9, 3.8225654361641483`*^9}, {3.8231764728567343`*^9, 
   3.823176482313031*^9}, {3.8233640113399878`*^9, 3.823364023347041*^9}, {
   3.823422996286858*^9, 3.823422999926001*^9}, {3.823602893791738*^9, 
   3.823602896188092*^9}, {3.8236833287515574`*^9, 3.823683328958127*^9}, {
   3.823966903508623*^9, 3.823966903859309*^9}, {3.823966992570949*^9, 
   3.823966994785722*^9}, {3.825859229609981*^9, 3.825859230168909*^9}, {
   3.850737092718917*^9, 3.850737113741796*^9}, 3.8521317157482023`*^9, {
   3.852133617349394*^9, 3.852133621456737*^9}, 3.852135095983374*^9, {
   3.855275349439538*^9, 3.855275352632153*^9}, 3.855276459218237*^9, {
   3.87730250808038*^9, 3.87730253694945*^9}},
 TextAlignment->Center,ExpressionUUID->"c94396a1-e700-429a-ac7f-26930c3cd713"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Hyperlink To Documentation", "Section",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.},
 CellChangeTimes->{{3.820944675605734*^9, 3.8209446931410646`*^9}, 
   3.877302532218264*^9},ExpressionUUID->"25f71703-2bc9-42df-9f5d-\
1f6e38728f77"],

Cell[BoxData[
 RowBox[{"Hyperlink", "[", 
  RowBox[{
  "\"\<GeneralRelativityTensors Documentation and Download\>\"", ",", 
   "\[IndentingNewLine]", 
   "\"\<https://github.com/BlackHolePerturbationToolkit\>\""}], 
  "]"}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.},
 CellChangeTimes->{{3.820944619768764*^9, 3.820944681268063*^9}, 
   3.877302532218608*^9},
 CellLabel->"In[42]:=",ExpressionUUID->"60620864-fa1e-46bd-98a1-e52d473c5e77"],

Cell[BoxData[
 TemplateBox[{
  "\"GeneralRelativityTensors Documentation and Download\"", 
   "https://github.com/BlackHolePerturbationToolkit"},
  "HyperlinkURL"]], "Output",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.},
 CellChangeTimes->{{3.820944653562152*^9, 3.820944681268403*^9}, 
   3.821729905862431*^9, 3.821732641380642*^9, 3.8218069931429567`*^9, 
   3.821811425912986*^9, 3.822042652654992*^9, {3.822048793766124*^9, 
   3.8220488009116488`*^9}, 3.82250502887521*^9, 3.822505578829629*^9, 
   3.822565490613895*^9, 3.8225662097813597`*^9, 3.822606717894772*^9, 
   3.822739433713368*^9, 3.822822210464344*^9, 3.8233594284652433`*^9, 
   3.8233607803938913`*^9, 3.823360880234103*^9, 3.823361224481102*^9, 
   3.82336209870968*^9, 3.823364030542309*^9, 3.823423099366008*^9, 
   3.823432318595525*^9, 3.823435951894755*^9, 3.8234364866582403`*^9, 
   3.8279535191386757`*^9, 3.827991252970625*^9, 3.827992769438249*^9, 
   3.8284359064920053`*^9, 3.828455531917368*^9, 3.828459339591674*^9, 
   3.828469241787932*^9, 3.828471664239974*^9, 3.8284721155283823`*^9, 
   3.828530571609256*^9, 3.828533466311865*^9, 3.828605666606776*^9, 
   3.828611773790597*^9, 3.82861654226705*^9, 3.828690472516346*^9, 
   3.828693683718005*^9, 3.828696145626916*^9, 3.828697685976324*^9, 
   3.828789577672448*^9, 3.828799204299177*^9, 3.828803846851424*^9, 
   3.828806278746703*^9, 3.828814185433984*^9, 3.8288338117285624`*^9, 
   3.8288341327645693`*^9, 3.828858949744699*^9, 3.828859274913085*^9, 
   3.8288629564377623`*^9, 3.8288631705507107`*^9, 3.828864762563582*^9, 
   3.8288658476756067`*^9, 3.828867772240056*^9, 3.829979175106444*^9, 
   3.830076411662816*^9, 3.830082279965604*^9, 3.830082624032351*^9, 
   3.877302532218936*^9},
 CellLabel->"Out[42]=",ExpressionUUID->"8db0f510-4cee-4b37-8b38-2fb042b2b514"]
}, Closed]],

Cell[CellGroupData[{

Cell["Hyperlink  To Article", "Section",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.},
 CellChangeTimes->{
  3.820937794883616*^9, {3.820938266275806*^9, 3.82093829464357*^9}, 
   3.820943883376583*^9, {3.8210324492405148`*^9, 3.821032453951701*^9}, {
   3.8217297279218817`*^9, 3.8217297313782063`*^9}, {3.822593613363134*^9, 
   3.822593616181035*^9}, {3.822736144359599*^9, 3.8227361450547543`*^9}, {
   3.823359417697976*^9, 3.823359420672374*^9}, {3.828691206350932*^9, 
   3.828691209871169*^9}, {3.828703099384598*^9, 3.828703101328167*^9}, {
   3.855499374752931*^9, 
   3.855499376597884*^9}},ExpressionUUID->"cbf1b83d-a087-4194-ae25-\
c5034b331158"],

Cell[BoxData[
 RowBox[{"Hyperlink", "[", 
  RowBox[{
  "\"\<Axisymmetric evolution code for the Einstein equations on \
hyperboloidal slices by Rinne\>\"", ",", "\[IndentingNewLine]", 
   "\"\<https://arxiv.org/abs/0910.0139\>\""}], "]"}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.},
 CellChangeTimes->{{3.8554993277202263`*^9, 3.855499344853116*^9}},
 CellLabel->"In[2]:=",ExpressionUUID->"ecb12663-7450-4eb3-b230-be611c0044ba"],

Cell[BoxData[
 TemplateBox[{
  "\"Axisymmetric evolution code for the Einstein equations on hyperboloidal \
slices by Rinne\"", "https://arxiv.org/abs/0910.0139"},
  "HyperlinkURL"]], "Output",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.},
 CellChangeTimes->{3.855499346685338*^9},
 CellLabel->"Out[2]=",ExpressionUUID->"4b68460b-c448-4278-a631-6961ebbed2f9"]
}, Closed]],

Cell[CellGroupData[{

Cell["Utilities and Package Load", "Section",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.},
 CellChangeTimes->{{3.820937646765551*^9, 3.8209376725449123`*^9}, 
   3.828806689104821*^9},ExpressionUUID->"0e4cc30c-24f6-4803-9550-\
70d5277cad1d"],

Cell[BoxData[{
 RowBox[{"ClearAll", "[", "\"\<Global`*\>\"", "]"}], "\[IndentingNewLine]", 
 RowBox[{"<<", "Utilities`CleanSlate`"}], "\[IndentingNewLine]", 
 RowBox[{"CleanSlate", "[", "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ClearInOut", "[", "]"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"pdConv", "[", "f_", "]"}], ":=", 
  RowBox[{"TraditionalForm", "[", 
   RowBox[{"f", "/.", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"Derivative", "[", "inds__", "]"}], "[", "g_", "]"}], "[", 
      "vars__", "]"}], "\[RuleDelayed]", 
     RowBox[{"Apply", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"Defer", "[", 
         RowBox[{"D", "[", 
          RowBox[{
           RowBox[{"g", "[", "vars", "]"}], ",", "##"}], "]"}], "]"}], "&"}], 
       ",", 
       RowBox[{
        RowBox[{"Transpose", "[", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"{", "vars", "}"}], ",", 
           RowBox[{"{", "inds", "}"}]}], "}"}], "]"}], "/.", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{
           RowBox[{"{", 
            RowBox[{"var_", ",", "0"}], "}"}], "\[RuleDelayed]", 
           RowBox[{"Sequence", "[", "]"}]}], ",", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"var_", ",", "1"}], "}"}], "\[RuleDelayed]", 
           RowBox[{"{", "var", "}"}]}]}], "}"}]}]}], "]"}]}]}], 
   "]"}]}]}], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.},
 CellChangeTimes->{
  3.8209376549030037`*^9, {3.821032517164062*^9, 3.821032518258361*^9}, 
   3.828806689105043*^9},
 CellLabel->"In[43]:=",ExpressionUUID->"2149b6f5-1146-424e-a6fd-d224d23ae853"],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"  (CleanSlate) Contexts purged: \"\>", "\[InvisibleSpace]", 
   RowBox[{"{", "\<\"Global`\"\>", "}"}]}],
  SequenceForm["  (CleanSlate) Contexts purged: ", {"Global`"}],
  Editable->False]], "Print",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.},
 CellChangeTimes->{
  3.820932754769992*^9, 3.820933427630714*^9, 3.820937654903143*^9, 
   3.821729906190301*^9, 3.821732641685704*^9, 3.8218069933778133`*^9, 
   3.82181142618342*^9, 3.8220426528487463`*^9, {3.822048793940605*^9, 
   3.822048801089246*^9}, 3.822505029191699*^9, 3.822505579049889*^9, 
   3.822565490838564*^9, 3.822566210008327*^9, 3.822606718180009*^9, 
   3.8226371008270607`*^9, 3.822739433887854*^9, 3.822761800089258*^9, 
   3.8228172997642097`*^9, 3.822822210725956*^9, 3.822988687208267*^9, 
   3.823359428732101*^9, 3.823360780582996*^9, 3.823360880477087*^9, 
   3.823361224742215*^9, 3.8233620990886593`*^9, 3.8233640308219624`*^9, 
   3.823423099601778*^9, 3.8234323187773438`*^9, 3.823435952157167*^9, 
   3.823436486849058*^9, 3.827953519486154*^9, 3.82799125321449*^9, 
   3.82799276999289*^9, 3.828435906665306*^9, 3.828455532084074*^9, 
   3.828459339797955*^9, 3.828469241959221*^9, 3.8284716644571867`*^9, 
   3.828472115722972*^9, 3.828530571883527*^9, 3.828533467717411*^9, 
   3.828605666780765*^9, 3.828611773924569*^9, 3.82861654251404*^9, 
   3.828690472852253*^9, 3.828693683937971*^9, 3.8286961458514643`*^9, 
   3.828697686995554*^9, 3.828789578115527*^9, 3.828799204762722*^9, 
   3.8288038473503647`*^9, 3.828806280080041*^9, 3.828806689105138*^9, 
   3.8288141858769903`*^9, 3.828833812193956*^9, 3.828834133179943*^9, 
   3.828858950064399*^9, 3.828859275539399*^9, 3.828862957119445*^9, 
   3.828863171652452*^9, 3.828864762870928*^9, 3.828865847990973*^9, 
   3.828867772550511*^9, 3.8299791755203648`*^9, 3.830076412132945*^9, 
   3.830082280334453*^9, 3.830082632338097*^9},
 CellLabel->
  "During evaluation of \
In[43]:=",ExpressionUUID->"8f445517-9389-48de-9d7d-a0f0f3d2b9ba"],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"  (CleanSlate) Approximate kernel memory recovered: \"\>", 
   "\[InvisibleSpace]", "\<\"136 Kb\"\>"}],
  SequenceForm[
  "  (CleanSlate) Approximate kernel memory recovered: ", "136 Kb"],
  Editable->False]], "Print",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.},
 CellChangeTimes->{
  3.820932754769992*^9, 3.820933427630714*^9, 3.820937654903143*^9, 
   3.821729906190301*^9, 3.821732641685704*^9, 3.8218069933778133`*^9, 
   3.82181142618342*^9, 3.8220426528487463`*^9, {3.822048793940605*^9, 
   3.822048801089246*^9}, 3.822505029191699*^9, 3.822505579049889*^9, 
   3.822565490838564*^9, 3.822566210008327*^9, 3.822606718180009*^9, 
   3.8226371008270607`*^9, 3.822739433887854*^9, 3.822761800089258*^9, 
   3.8228172997642097`*^9, 3.822822210725956*^9, 3.822988687208267*^9, 
   3.823359428732101*^9, 3.823360780582996*^9, 3.823360880477087*^9, 
   3.823361224742215*^9, 3.8233620990886593`*^9, 3.8233640308219624`*^9, 
   3.823423099601778*^9, 3.8234323187773438`*^9, 3.823435952157167*^9, 
   3.823436486849058*^9, 3.827953519486154*^9, 3.82799125321449*^9, 
   3.82799276999289*^9, 3.828435906665306*^9, 3.828455532084074*^9, 
   3.828459339797955*^9, 3.828469241959221*^9, 3.8284716644571867`*^9, 
   3.828472115722972*^9, 3.828530571883527*^9, 3.828533467717411*^9, 
   3.828605666780765*^9, 3.828611773924569*^9, 3.82861654251404*^9, 
   3.828690472852253*^9, 3.828693683937971*^9, 3.8286961458514643`*^9, 
   3.828697686995554*^9, 3.828789578115527*^9, 3.828799204762722*^9, 
   3.8288038473503647`*^9, 3.828806280080041*^9, 3.828806689105138*^9, 
   3.8288141858769903`*^9, 3.828833812193956*^9, 3.828834133179943*^9, 
   3.828858950064399*^9, 3.828859275539399*^9, 3.828862957119445*^9, 
   3.828863171652452*^9, 3.828864762870928*^9, 3.828865847990973*^9, 
   3.828867772550511*^9, 3.8299791755203648`*^9, 3.830076412132945*^9, 
   3.830082280334453*^9, 3.8300826323414707`*^9},
 CellLabel->
  "During evaluation of \
In[43]:=",ExpressionUUID->"8dbf193b-88cf-408c-a6c3-e7ecdafcf426"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"\<\"Utilities`CleanSlate`\"\>", ",", "\<\"VariationalMethods`\"\>",
    ",", "\<\"GeneralRelativityTensors`\"\>", 
   ",", "\<\"GeneralRelativityTensors`CommonTensors`\"\>", 
   ",", "\<\"GeneralRelativityTensors`TensorDerivatives`\"\>", 
   ",", "\<\"GeneralRelativityTensors`TensorManipulation`\"\>", 
   ",", "\<\"GeneralRelativityTensors`TensorDefinitions`\"\>", 
   ",", "\<\"GeneralRelativityTensors`Utils`\"\>", 
   ",", "\<\"DocumentationSearch`\"\>", ",", "\<\"ResourceLocator`\"\>", 
   ",", "\<\"System`\"\>", ",", "\<\"Global`\"\>"}], "}"}]], "Output",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.},
 CellChangeTimes->{
  3.820932754773623*^9, 3.820933427635436*^9, 3.8209376549033947`*^9, 
   3.821729906193787*^9, 3.821732641691049*^9, 3.821806993386058*^9, 
   3.821811426187071*^9, 3.822042652852829*^9, {3.822048793944077*^9, 
   3.822048801092909*^9}, 3.822505029195156*^9, 3.8225055790570173`*^9, 
   3.8225654908422747`*^9, 3.8225662100125237`*^9, 3.8226067181902227`*^9, 
   3.822637100839163*^9, 3.8227394339100733`*^9, 3.822761800103311*^9, 
   3.822817299776285*^9, 3.822822210739685*^9, 3.822988687217742*^9, 
   3.823359428736412*^9, 3.823360780587545*^9, 3.8233608804810534`*^9, 
   3.823361224746982*^9, 3.823362099093968*^9, 3.823364030838743*^9, 
   3.823423099606819*^9, 3.8234323187826138`*^9, 3.8234359521606693`*^9, 
   3.8234364868537683`*^9, 3.8279535194907703`*^9, 3.827991253218878*^9, 
   3.827992769998845*^9, 3.8284359066702137`*^9, 3.828455532087866*^9, 
   3.82845933980285*^9, 3.828469241964273*^9, 3.828471664461882*^9, 
   3.828472115727516*^9, 3.8285305718886538`*^9, 3.82853346772556*^9, 
   3.8286056667871323`*^9, 3.828611773929594*^9, 3.828616542519825*^9, 
   3.828690472858699*^9, 3.828693683944024*^9, 3.82869614585636*^9, 
   3.828697687001523*^9, 3.8287895781234426`*^9, 3.828799204770609*^9, 
   3.828803847356839*^9, 3.828806280086404*^9, 3.828806689105328*^9, 
   3.828814185882605*^9, 3.828833812202352*^9, 3.828834133188683*^9, 
   3.828858950068405*^9, 3.828859275549589*^9, 3.828862957127707*^9, 
   3.828863171659013*^9, 3.8288647628760147`*^9, 3.82886584799824*^9, 
   3.8288677725571413`*^9, 3.8299791755258827`*^9, 3.830076412139241*^9, 
   3.830082280340658*^9, 3.830082632344269*^9},
 CellLabel->"Out[0]=",ExpressionUUID->"6958ad4d-bb31-4ea3-aeab-09df344dbb41"],

Cell[BoxData[
 RowBox[{"<<", "GeneralRelativityTensors`"}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.},
 CellChangeTimes->{3.8209376549035053`*^9, 3.828806689105419*^9},
 CellLabel->"In[2]:=",ExpressionUUID->"00e10672-0c1f-492b-9d5d-a8531d6214e2"],

Cell[BoxData[
 RowBox[{"Needs", "[", "\"\<VariationalMethods`\>\"", "]"}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.},
 CellChangeTimes->{3.8287882026748323`*^9},
 CellLabel->"In[3]:=",ExpressionUUID->"006f065b-9f32-4f3c-9323-3f51d1b0b8ef"]
}, Closed]],

Cell[CellGroupData[{

Cell["Custom Notation", "Section",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.},
 CellChangeTimes->{{3.8225114042235403`*^9, 3.822511435665856*^9}, 
   3.823176577809782*^9},ExpressionUUID->"7b375b49-ee04-4b8f-9690-\
3260f9cb9414"],

Cell[BoxData[{
 RowBox[{"Clear", "[", "dtReplace", "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"dtReplace", " ", "=", " ", 
   RowBox[{"{", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"Dt", "[", "r", "]"}], "\[Rule]", " ", "dr"}], " ", ",", " ", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Dt", "[", "t", "]"}], "\[Rule]", " ", "dt"}], " ", ",", " ", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Dt", "[", "u", "]"}], "\[Rule]", " ", "du"}], " ", ",", " ", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Dt", "[", "v", "]"}], "\[Rule]", " ", "dv"}], " ", ",", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Dt", "[", "x", "]"}], "\[Rule]", " ", "dx"}], " ", ",", " ", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Dt", "[", "y", "]"}], "\[Rule]", " ", "dy"}], " ", ",", " ", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Dt", "[", "z", "]"}], "\[Rule]", " ", "dz"}], " ", ",", "  ", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Dt", "[", "\[Tau]", "]"}], "\[Rule]", " ", "d\[Tau]"}], " ", 
     ",", " ", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Dt", "[", "\[Theta]", "]"}], " ", "\[Rule]", "d\[Theta]"}], 
     ",", " ", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Dt", "[", "\[Phi]", "]"}], " ", "\[Rule]", "d\[Phi]"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Dt", "[", "\[Eta]", "]"}], " ", "\[Rule]", "d\[Eta]"}], ",", 
     " ", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Dt", "[", "\[Chi]", "]"}], "\[Rule]", " ", "d\[Chi]"}], " ", 
     ",", " ", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Dt", "[", "\[Rho]", "]"}], " ", "\[Rule]", " ", "d\[Rho]"}], 
     " ", ",", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Dt", "[", "\[CurlyPhi]", "]"}], " ", "\[Rule]", " ", 
      "d\[CurlyPhi]"}], " ", ",", " ", "\[IndentingNewLine]", " ", 
     RowBox[{
      RowBox[{"Dt", "[", "\[ScriptT]", "]"}], " ", "\[Rule]", " ", 
      "d\[ScriptT]"}], " ", ",", " ", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Dt", "[", "\[ScriptZ]", "]"}], " ", "\[Rule]", " ", 
      "d\[ScriptZ]"}], ",", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Dt", "[", "T", "]"}], "\[Rule]", " ", "dT"}], " ", ",", " ", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Dt", "[", "U", "]"}], "\[Rule]", " ", "dU"}], " ", ",", " ", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Dt", "[", "V", "]"}], "\[Rule]", " ", "dV"}], " ", ",", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Dt", "[", "\[ScriptR]", "]"}], "\[Rule]", " ", "d\[ScriptR]"}],
      "  ", ",", " ", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Dt", "[", "\[CapitalPhi]", "]"}], "\[Rule]", " ", 
      "d\[CapitalPhi]"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"dtReplace", " ", "//", " ", "TableForm"}]}], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.},
 CellChangeTimes->{{3.82317656063256*^9, 3.823176577809904*^9}, {
  3.828606009159594*^9, 3.828606024422215*^9}, {3.828701721392323*^9, 
  3.828701725375505*^9}, {3.829979800080254*^9, 3.8299798262997217`*^9}},
 CellLabel->"In[4]:=",ExpressionUUID->"82edc7de-deef-4b9a-90f8-5d6133ff2655"],

Cell[BoxData[
 TagBox[
  TagBox[GridBox[{
     {
      RowBox[{
       RowBox[{"Dt", "[", "r", "]"}], "\[Rule]", "dr"}]},
     {
      RowBox[{
       RowBox[{"Dt", "[", "t", "]"}], "\[Rule]", "dt"}]},
     {
      RowBox[{
       RowBox[{"Dt", "[", "u", "]"}], "\[Rule]", "du"}]},
     {
      RowBox[{
       RowBox[{"Dt", "[", "v", "]"}], "\[Rule]", "dv"}]},
     {
      RowBox[{
       RowBox[{"Dt", "[", "x", "]"}], "\[Rule]", "dx"}]},
     {
      RowBox[{
       RowBox[{"Dt", "[", "y", "]"}], "\[Rule]", "dy"}]},
     {
      RowBox[{
       RowBox[{"Dt", "[", "z", "]"}], "\[Rule]", "dz"}]},
     {
      RowBox[{
       RowBox[{"Dt", "[", "\[Tau]", "]"}], "\[Rule]", "d\[Tau]"}]},
     {
      RowBox[{
       RowBox[{"Dt", "[", "\[Theta]", "]"}], "\[Rule]", "d\[Theta]"}]},
     {
      RowBox[{
       RowBox[{"Dt", "[", "\[Phi]", "]"}], "\[Rule]", "d\[Phi]"}]},
     {
      RowBox[{
       RowBox[{"Dt", "[", "\[Eta]", "]"}], "\[Rule]", "d\[Eta]"}]},
     {
      RowBox[{
       RowBox[{"Dt", "[", "\[Chi]", "]"}], "\[Rule]", "d\[Chi]"}]},
     {
      RowBox[{
       RowBox[{"Dt", "[", "\[Rho]", "]"}], "\[Rule]", "d\[Rho]"}]},
     {
      RowBox[{
       RowBox[{"Dt", "[", "\[CurlyPhi]", "]"}], "\[Rule]", "d\[CurlyPhi]"}]},
     {
      RowBox[{
       RowBox[{"Dt", "[", "\[ScriptT]", "]"}], "\[Rule]", "d\[ScriptT]"}]},
     {
      RowBox[{
       RowBox[{"Dt", "[", "\[ScriptZ]", "]"}], "\[Rule]", "d\[ScriptZ]"}]},
     {
      RowBox[{
       RowBox[{"Dt", "[", "T", "]"}], "\[Rule]", "dT"}]},
     {
      RowBox[{
       RowBox[{"Dt", "[", "U", "]"}], "\[Rule]", "dU"}]},
     {
      RowBox[{
       RowBox[{"Dt", "[", "V", "]"}], "\[Rule]", "dV"}]},
     {
      RowBox[{
       RowBox[{"Dt", "[", "\[ScriptR]", "]"}], "\[Rule]", "d\[ScriptR]"}]},
     {
      RowBox[{
       RowBox[{"Dt", "[", "\[CapitalPhi]", "]"}], "\[Rule]", 
       "d\[CapitalPhi]"}]}
    },
    GridBoxAlignment->{"Columns" -> {{Left}}, "Rows" -> {{Baseline}}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.5599999999999999]}, 
        Offset[0.27999999999999997`]}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}}],
   Column],
  Function[BoxForm`e$, 
   TableForm[BoxForm`e$]]]], "Output",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.},
 CellChangeTimes->{{3.8231765606326637`*^9, 3.82317657781*^9}, 
   3.8233594287980833`*^9, 3.823360780614534*^9, 3.8233608805437098`*^9, 
   3.823361224810658*^9, 3.823362099159925*^9, 3.8233640309198627`*^9, 
   3.8234230996755753`*^9, 3.823432318844563*^9, 3.823435952223415*^9, 
   3.8234364868866034`*^9, 3.827953519527835*^9, 3.8279912532863617`*^9, 
   3.827992770095962*^9, 3.828435906737782*^9, 3.828455532145804*^9, 
   3.82845933986725*^9, 3.8284692420170603`*^9, 3.828471664529696*^9, 
   3.82847211575672*^9, 3.828530571960513*^9, 3.828533469228689*^9, 
   3.828605666816091*^9, 3.8286060250588503`*^9, 3.8286117739948387`*^9, 
   3.828616542589233*^9, 3.828690472957294*^9, 3.828693684019373*^9, 
   3.828694593299365*^9, 3.828696145930181*^9, 3.8286976871734247`*^9, 
   3.828701726000832*^9, 3.8287895781937857`*^9, 3.828799204839447*^9, 
   3.828803847429542*^9, 3.828806280308897*^9, 3.828814185992796*^9, 
   3.8288338122924232`*^9, 3.8288341333989553`*^9, 3.828858950097949*^9, 
   3.828859275654873*^9, 3.8288629572392063`*^9, 3.828863171871059*^9, 
   3.828864762961931*^9, 3.8288658488547573`*^9, 3.8288677726704407`*^9, 
   3.8299791756359386`*^9, {3.829979804930752*^9, 3.829979814038682*^9}, 
   3.8300764121955757`*^9, 3.830082280437125*^9},
 CellLabel->
  "Out[6]//TableForm=",ExpressionUUID->"e25565a5-2371-48a9-ae1e-fe77124e28b6"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"M", " ", "/:", " ", 
    RowBox[{"Dt", "[", "M", "]"}], " ", "=", " ", "0"}], " ", ";"}], "  ", 
  RowBox[{"(*", " ", 
   RowBox[{
   "Mass", " ", "of", " ", "Schwarzschild", " ", "is", " ", "constant"}], " ",
    "*)"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{
   "We", " ", "use", " ", "M", " ", "for", " ", "mass", " ", "and", " ", "m", 
    " ", "for", " ", "null", " ", "vector"}], " ", "*)"}], " "}]], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.},
 CellChangeTimes->{{3.82869475118476*^9, 3.828694761903348*^9}, 
   3.828695041652215*^9, {3.828697403799129*^9, 3.8286974141262493`*^9}},
 CellLabel->"In[7]:=",ExpressionUUID->"5a996066-ba9d-4c2b-8955-1aab763ce070"]
}, Closed]],

Cell[CellGroupData[{

Cell["Line Element and Metric Functions", "Section",
 CellGroupingRules->{"GroupTogetherGrouping", 10003.},
 CellChangeTimes->{{3.820937548919169*^9, 3.8209375728447866`*^9}, 
   3.8209383129401073`*^9, 3.8220481344968433`*^9, {3.8225115057038183`*^9, 
   3.822511507783984*^9}, 
   3.822763734568725*^9},ExpressionUUID->"69166c17-09ce-4a98-905f-\
cf54b53a37ff"],

Cell[BoxData[{
 RowBox[{"Clear", "[", "lineToMetric", "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"lineToMetric", "[", 
    RowBox[{"lineelement_", " ", ",", " ", "differentials_"}], "]"}], ":=", 
   " ", "\[IndentingNewLine]", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{"i", "\[Equal]", "j"}], " ", ",", 
       RowBox[{"Coefficient", "[", 
        RowBox[{"lineelement", ",", " ", 
         RowBox[{
          RowBox[{"differentials", "[", 
           RowBox[{"[", "i", "]"}], "]"}], " ", 
          RowBox[{"differentials", "[", 
           RowBox[{"[", "j", "]"}], "]"}]}]}], "]"}], ",", 
       "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"1", "/", "2"}], ")"}], "*", 
        RowBox[{"Coefficient", "[", 
         RowBox[{"lineelement", ",", " ", 
          RowBox[{
           RowBox[{"differentials", "[", 
            RowBox[{"[", "i", "]"}], "]"}], " ", 
           RowBox[{"differentials", "[", 
            RowBox[{"[", "j", "]"}], "]"}]}]}], "]"}]}]}], "]"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", 
       RowBox[{"Length", "[", "differentials", "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "1", ",", 
       RowBox[{"Length", "[", "differentials", "]"}]}], "}"}]}], "]"}]}], 
  "  ", ";"}]}], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10003.},
 CellChangeTimes->{3.821567412203915*^9, 3.821641945174967*^9, 
  3.821641981470141*^9, 3.8220481344970007`*^9, 3.822763734569355*^9},
 CellLabel->"In[4]:=",ExpressionUUID->"8c11b0fe-dbf1-4361-a426-0cfcb695cba7"],

Cell[BoxData[{
 RowBox[{"Clear", "[", "metricToLine", "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"metricToLine", "[", 
   RowBox[{"metric_", ",", " ", "differentials_"}], "]"}], ":=", " ", 
  "\[IndentingNewLine]", 
  RowBox[{"Sum", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"metric", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", 
     RowBox[{"differentials", "[", 
      RowBox[{"[", "i", "]"}], "]"}], " ", 
     RowBox[{"differentials", "[", 
      RowBox[{"[", "j", "]"}], "]"}]}], " ", ",", " ", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", 
      RowBox[{"Length", "[", "differentials", "]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"j", ",", "1", ",", 
      RowBox[{"Length", "[", "differentials", "]"}]}], "}"}]}], 
   "]"}]}]}], "Input",
 CellGroupingRules->{"GroupTogetherGrouping", 10003.},
 CellChangeTimes->{{3.822762557092309*^9, 3.8227626322008677`*^9}, 
   3.8227637345698223`*^9, 3.822988709740047*^9},
 CellLabel->"In[6]:=",ExpressionUUID->"74c3693d-1f50-406e-b772-cf55d33e19d4"]
}, Closed]]
},
WindowSize->{1278, 675},
WindowMargins->{{Automatic, 0}, {Automatic, 0}},
FrontEndVersion->"13.0 for Mac OS X x86 (64-bit) (February 4, 2022)",
StyleDefinitions->FrontEnd`FileName[{"Report"}, "StandardReport.nb", 
  CharacterEncoding -> "UTF-8"],
ExpressionUUID->"520c3445-0d73-40a0-a6d4-55443903a185"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 2031, 30, 221, "Title",ExpressionUUID->"3d6666ab-568a-4851-baef-ae695b51f0cd",
 CellGroupingRules->{"GroupTogetherGrouping", 10000.}],
Cell[2614, 54, 1310, 22, 94, "Author",ExpressionUUID->"c94396a1-e700-429a-ac7f-26930c3cd713",
 CellGroupingRules->{"GroupTogetherGrouping", 10000.}]
}, Open  ]],
Cell[CellGroupData[{
Cell[3961, 81, 251, 4, 69, "Section",ExpressionUUID->"25f71703-2bc9-42df-9f5d-1f6e38728f77",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.}],
Cell[4215, 87, 460, 10, 63, "Input",ExpressionUUID->"60620864-fa1e-46bd-98a1-e52d473c5e77",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.}],
Cell[4678, 99, 1828, 28, 43, "Output",ExpressionUUID->"8db0f510-4cee-4b37-8b38-2fb042b2b514",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.}]
}, Closed]],
Cell[CellGroupData[{
Cell[6543, 132, 668, 11, 55, "Section",ExpressionUUID->"cbf1b83d-a087-4194-ae25-c5034b331158",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.}],
Cell[7214, 145, 450, 8, 63, "Input",ExpressionUUID->"ecb12663-7450-4eb3-b230-be611c0044ba",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.}],
Cell[7667, 155, 368, 7, 43, "Output",ExpressionUUID->"4b68460b-c448-4278-a631-6961ebbed2f9",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.}]
}, Closed]],
Cell[CellGroupData[{
Cell[8072, 167, 251, 4, 55, "Section",ExpressionUUID->"0e4cc30c-24f6-4803-9550-70d5277cad1d",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.}],
Cell[8326, 173, 1667, 45, 189, "Input",ExpressionUUID->"2149b6f5-1146-424e-a6fd-d224d23ae853",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.}],
Cell[9996, 220, 2037, 33, 24, "Print",ExpressionUUID->"8f445517-9389-48de-9d7d-a0f0f3d2b9ba",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.}],
Cell[12036, 255, 2058, 34, 24, "Print",ExpressionUUID->"8dbf193b-88cf-408c-a6c3-e7ecdafcf426",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.}],
Cell[14097, 291, 2378, 36, 84, "Output",ExpressionUUID->"6958ad4d-bb31-4ea3-aeab-09df344dbb41",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.}],
Cell[16478, 329, 268, 4, 41, "Input",ExpressionUUID->"00e10672-0c1f-492b-9d5d-a8531d6214e2",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.}],
Cell[16749, 335, 261, 4, 41, "Input",ExpressionUUID->"006f065b-9f32-4f3c-9323-3f51d1b0b8ef",
 CellGroupingRules->{"GroupTogetherGrouping", 10002.}]
}, Closed]],
Cell[CellGroupData[{
Cell[17047, 344, 240, 4, 55, "Section",ExpressionUUID->"7b375b49-ee04-4b8f-9690-3260f9cb9414",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.}],
Cell[17290, 350, 3262, 74, 525, "Input",ExpressionUUID->"82edc7de-deef-4b9a-90f8-5d6133ff2655",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.}],
Cell[20555, 426, 3696, 99, 411, "Output",ExpressionUUID->"e25565a5-2371-48a9-ae1e-fe77124e28b6",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.}],
Cell[24254, 527, 735, 16, 41, "Input",ExpressionUUID->"5a996066-ba9d-4c2b-8955-1aab763ce070",
 CellGroupingRules->{"GroupTogetherGrouping", 10001.}]
}, Closed]],
Cell[CellGroupData[{
Cell[25026, 548, 362, 6, 55, "Section",ExpressionUUID->"69166c17-09ce-4a98-905f-cf54b53a37ff",
 CellGroupingRules->{"GroupTogetherGrouping", 10003.}],
Cell[25391, 556, 1659, 41, 126, "Input",ExpressionUUID->"8c11b0fe-dbf1-4361-a426-0cfcb695cba7",
 CellGroupingRules->{"GroupTogetherGrouping", 10003.}],
Cell[27053, 599, 1055, 26, 84, "Input",ExpressionUUID->"74c3693d-1f50-406e-b772-cf55d33e19d4",
 CellGroupingRules->{"GroupTogetherGrouping", 10003.}]
}, Closed]]
}
]
*)