From 3ba26f3f986478b616160d06e679cf38ad015a13 Mon Sep 17 00:00:00 2001 From: David Hilditch Date: Fri, 24 Feb 2023 15:51:03 +0000 Subject: add a mathematica notebook computing Christoffels. --- math_notebooks/Christoffels.nb | 1153 ++++++++++++++++++++++++++++++++++++++++ 1 file changed, 1153 insertions(+) create mode 100644 math_notebooks/Christoffels.nb diff --git a/math_notebooks/Christoffels.nb b/math_notebooks/Christoffels.nb new file mode 100644 index 0000000..be80b41 --- /dev/null +++ b/math_notebooks/Christoffels.nb @@ -0,0 +1,1153 @@ +(* Content-type: application/vnd.wolfram.mathematica *) + +(*** Wolfram Notebook File ***) +(* http://www.wolfram.com/nb *) + +(* CreatedBy='Mathematica 13.1' *) + +(*CacheID: 234*) +(* Internal cache information: +NotebookFileLineBreakTest +NotebookFileLineBreakTest +NotebookDataPosition[ 158, 7] +NotebookDataLength[ 42377, 1145] +NotebookOptionsPosition[ 37437, 1058] +NotebookOutlinePosition[ 37878, 1076] +CellTagsIndexPosition[ 37835, 1073] +WindowFrame->Normal*) + +(* Beginning of Notebook Content *) +Notebook[{ +Cell["\<\ +Christoffel Computation +David Hilditch & Anton Khirnov, Feb. 2023.\ +\>", "Subsubsection", + CellChangeTimes->{{3.886238482476939*^9, 3.8862385051531477`*^9}, { + 3.886238556279541*^9, + 3.886238570895651*^9}},ExpressionUUID->"373d7a92-84df-4cc8-94b3-\ +ee4da185a64c"], + +Cell[CellGroupData[{ + +Cell["Define coordinates, metric and inverse metric", "Text", + CellChangeTimes->{{3.88623852761099*^9, + 3.886238544018025*^9}},ExpressionUUID->"7b3be704-efcd-4b3c-a35c-\ +5fafcf85c005"], + +Cell[BoxData[ + RowBox[{ + RowBox[{"X", "=", + RowBox[{"{", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "}"}]}], ";"}]], "Input", + CellChangeTimes->{{3.886171931510335*^9, 3.886171933623374*^9}, { + 3.886171972930284*^9, 3.886171975545806*^9}, {3.88617265745574*^9, + 3.886172663671729*^9}, {3.886173418537775*^9, 3.8861734198173847`*^9}}, + CellLabel->"In[1]:=",ExpressionUUID->"33ae6975-ec1d-4e40-8c00-336b4482e361"], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{ + RowBox[{"g", "=", + RowBox[{"{", + RowBox[{ + RowBox[{"{", + RowBox[{ + RowBox[{"gtt", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", + RowBox[{"gtx", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", + RowBox[{"gty", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", + RowBox[{"gtz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"gtx", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", + RowBox[{"gxx", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", + RowBox[{"gxy", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", + RowBox[{"gxz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"gty", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", + RowBox[{"gxy", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", + RowBox[{"gyy", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", + RowBox[{"gyz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"gtz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", + RowBox[{"gxz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", + RowBox[{"gyz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", + RowBox[{"gzz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}], "}"}]}], + "}"}]}], ";", + RowBox[{"g", "//", "MatrixForm"}]}]], "Input", + CellChangeTimes->{{3.8861719282262793`*^9, 3.886171929989098*^9}, { + 3.886171987992716*^9, 3.88617200715347*^9}, {3.886172316594161*^9, + 3.886172333075539*^9}, {3.886172647434791*^9, 3.886172651808694*^9}, { + 3.886172682276737*^9, 3.886172909660861*^9}, {3.8861729438025417`*^9, + 3.886172948517481*^9}, {3.8861733110110073`*^9, 3.8861734651163197`*^9}, { + 3.8861744456232347`*^9, 3.886174447108243*^9}, 3.886174614151684*^9}, + CellLabel->"In[2]:=",ExpressionUUID->"bfbd2c04-efbd-49b8-940f-e8e2a48cee58"], + +Cell[BoxData[ + TagBox[ + RowBox[{"(", "\[NoBreak]", GridBox[{ + { + RowBox[{"gtt", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], + RowBox[{"gtx", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], + RowBox[{"gty", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], + RowBox[{"gtz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}, + { + RowBox[{"gtx", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], + RowBox[{"gxx", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], + RowBox[{"gxy", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], + RowBox[{"gxz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}, + { + RowBox[{"gty", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], + RowBox[{"gxy", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], + RowBox[{"gyy", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], + RowBox[{"gyz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}, + { + RowBox[{"gtz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], + RowBox[{"gxz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], + RowBox[{"gyz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], + RowBox[{"gzz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]} + }, + GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Baseline}}}, + GridBoxSpacings->{"Columns" -> { + Offset[0.27999999999999997`], { + Offset[0.7]}, + Offset[0.27999999999999997`]}, "Rows" -> { + Offset[0.2], { + Offset[0.4]}, + Offset[0.2]}}], "\[NoBreak]", ")"}], + Function[BoxForm`e$, + MatrixForm[BoxForm`e$]]]], "Output", + CellChangeTimes->{{3.886172749895281*^9, 3.8861727752028303`*^9}, { + 3.886172893537417*^9, 3.886172912569895*^9}, 3.886172949054069*^9, + 3.886173465656674*^9, 3.886174204828216*^9, 3.886174447409831*^9, + 3.8861746149973183`*^9, 3.886238766426577*^9}, + CellLabel-> + "Out[2]//MatrixForm=",ExpressionUUID->"afb62946-6d22-40db-87c7-\ +06d8aaa5eafb"] +}, Open ]], + +Cell[BoxData[ + RowBox[{ + RowBox[{"ginv", "=", + RowBox[{ + RowBox[{"g", "//", "Inverse"}], "//", "Simplify"}]}], ";"}]], "Input", + CellChangeTimes->{{3.886173014702764*^9, 3.88617303331397*^9}}, + CellLabel->"In[3]:=",ExpressionUUID->"ff8d5500-b0da-4dea-aed3-005d3f38ce38"] +}, Closed]], + +Cell[CellGroupData[{ + +Cell["Compute Christoffels (ddd & udd)", "Text", + CellChangeTimes->{{3.8862386105171843`*^9, + 3.8862386296497107`*^9}},ExpressionUUID->"0e39b00c-5443-4cc9-9420-\ +976408ee456d"], + +Cell[BoxData[ + RowBox[{ + RowBox[{"Christoffelsddd", "=", + RowBox[{ + RowBox[{"Table", "[", + RowBox[{ + RowBox[{ + RowBox[{"1", "/", "2"}], + RowBox[{"(", + RowBox[{ + RowBox[{ + SubscriptBox["\[PartialD]", + RowBox[{"X", "[", + RowBox[{"[", "\[Mu]", "]"}], "]"}]], + RowBox[{"g", "[", + RowBox[{"[", + RowBox[{"\[Nu]", ",", "\[Delta]"}], "]"}], "]"}]}], "+", + RowBox[{ + SubscriptBox["\[PartialD]", + RowBox[{"X", "[", + RowBox[{"[", "\[Nu]", "]"}], "]"}]], + RowBox[{"g", "[", + RowBox[{"[", + RowBox[{"\[Mu]", ",", "\[Delta]"}], "]"}], "]"}]}], "-", + RowBox[{ + SubscriptBox["\[PartialD]", + RowBox[{"X", "[", + RowBox[{"[", "\[Delta]", "]"}], "]"}]], + RowBox[{"g", "[", + RowBox[{"[", + RowBox[{"\[Mu]", ",", "\[Nu]"}], "]"}], "]"}]}]}], ")"}]}], ",", + RowBox[{"{", + RowBox[{"\[Delta]", ",", "4"}], "}"}], ",", + RowBox[{"{", + RowBox[{"\[Mu]", ",", "4"}], "}"}], ",", + RowBox[{"{", + RowBox[{"\[Nu]", ",", "4"}], "}"}]}], "]"}], "//", "Simplify"}]}], + ";"}]], "Input", + CellChangeTimes->{{3.8861743287685547`*^9, 3.886174365030443*^9}}, + CellLabel->"In[4]:=",ExpressionUUID->"ac869ef8-5c04-4ff7-b6ff-063dc324d8f4"], + +Cell[BoxData[ + RowBox[{ + RowBox[{"Christoffelsudd", "=", + RowBox[{ + RowBox[{"Table", "[", + RowBox[{ + RowBox[{ + RowBox[{"1", "/", "2"}], + RowBox[{"Sum", "[", + RowBox[{ + RowBox[{ + RowBox[{"ginv", "[", + RowBox[{"[", + RowBox[{"\[Alpha]", ",", "\[Delta]"}], "]"}], "]"}], + RowBox[{"(", + RowBox[{ + RowBox[{ + SubscriptBox["\[PartialD]", + RowBox[{"X", "[", + RowBox[{"[", "\[Mu]", "]"}], "]"}]], + RowBox[{"g", "[", + RowBox[{"[", + RowBox[{"\[Nu]", ",", "\[Delta]"}], "]"}], "]"}]}], "+", + RowBox[{ + SubscriptBox["\[PartialD]", + RowBox[{"X", "[", + RowBox[{"[", "\[Nu]", "]"}], "]"}]], + RowBox[{"g", "[", + RowBox[{"[", + RowBox[{"\[Mu]", ",", "\[Delta]"}], "]"}], "]"}]}], "-", + RowBox[{ + SubscriptBox["\[PartialD]", + RowBox[{"X", "[", + RowBox[{"[", "\[Delta]", "]"}], "]"}]], + RowBox[{"g", "[", + RowBox[{"[", + RowBox[{"\[Mu]", ",", "\[Nu]"}], "]"}], "]"}]}]}], ")"}]}], + ",", + RowBox[{"{", + RowBox[{"\[Delta]", ",", "4"}], "}"}]}], "]"}]}], ",", + RowBox[{"{", + RowBox[{"\[Alpha]", ",", "4"}], "}"}], ",", + RowBox[{"{", + RowBox[{"\[Mu]", ",", "4"}], "}"}], ",", + RowBox[{"{", + RowBox[{"\[Nu]", ",", "4"}], "}"}]}], "]"}], "//", "Simplify"}]}], + ";"}]], "Input", + CellChangeTimes->{{3.8861729804852133`*^9, 3.886173003793852*^9}, { + 3.886173056335821*^9, 3.886173105859511*^9}, {3.886173175235386*^9, + 3.88617327140637*^9}, 3.886173478420867*^9, {3.886173835805682*^9, + 3.8861738550901203`*^9}, {3.886174098114808*^9, 3.886174113642325*^9}, { + 3.886174235751086*^9, 3.8861742386193647`*^9}}, + CellLabel->"In[5]:=",ExpressionUUID->"e7e016e2-c00b-42d2-91be-a38628e97843"] +}, Closed]], + +Cell[CellGroupData[{ + +Cell["\<\ +Rules for twist-free axisymmetry with reflection symmetry, on the symmetry \ +axis and equatorial plane\ +\>", "Text", + CellChangeTimes->{{3.8862386444103394`*^9, + 3.88623870518185*^9}},ExpressionUUID->"c11437e5-f609-4db8-bb2d-\ +e4ba921912d9"], + +Cell[BoxData[ + RowBox[{ + RowBox[{"KillMeOnAxisRules", "=", + RowBox[{"{", + RowBox[{ + RowBox[{"gtx", "->", + RowBox[{"Function", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "}"}], ",", "0"}], + "]"}]}], ",", + RowBox[{"gty", "->", + RowBox[{"Function", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "}"}], ",", "0"}], + "]"}]}], ",", + RowBox[{"gxy", "->", + RowBox[{"Function", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "}"}], ",", "0"}], + "]"}]}], ",", + RowBox[{"gyz", "->", + RowBox[{"Function", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "}"}], ",", "0"}], + "]"}]}], ",", + RowBox[{"gxz", "\[Rule]", + RowBox[{"Function", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "}"}], ",", "0"}], + "]"}]}]}], "}"}]}], ";"}]], "Input", + CellChangeTimes->{{3.88617348439834*^9, 3.886173509870784*^9}, { + 3.886173553075624*^9, 3.8861735746880627`*^9}, {3.886173613403329*^9, + 3.886173656011073*^9}, {3.886173711669825*^9, 3.886173730830943*^9}, { + 3.886174771825536*^9, 3.886174773975914*^9}}, + CellLabel->"In[6]:=",ExpressionUUID->"1f3a5dca-b140-4e51-93fd-93ff73818e14"], + +Cell[BoxData[ + RowBox[{ + RowBox[{"KillMeEquatorRules", "=", + RowBox[{"{", + RowBox[{ + RowBox[{"gtz", "->", + RowBox[{"Function", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "}"}], ",", "0"}], + "]"}]}], ",", + RowBox[{"gty", "->", + RowBox[{"Function", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "}"}], ",", "0"}], + "]"}]}], ",", + RowBox[{"gxy", "->", + RowBox[{"Function", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "}"}], ",", "0"}], + "]"}]}], ",", + RowBox[{"gyz", "->", + RowBox[{"Function", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "}"}], ",", "0"}], + "]"}]}], ",", + RowBox[{ + RowBox[{"gxz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], "->", "0"}]}], + "}"}]}], ";"}]], "Input", + CellChangeTimes->{{3.886174790055666*^9, 3.886174820279228*^9}, { + 3.886174952284816*^9, 3.886174960675373*^9}}, + CellLabel->"In[7]:=",ExpressionUUID->"9461b4d0-c922-4076-a54e-0d59230005dd"] +}, Closed]], + +Cell[CellGroupData[{ + +Cell["Examine metric on axis and in equator :", "Text", + CellChangeTimes->{{3.886238741422525*^9, + 3.8862387597451887`*^9}},ExpressionUUID->"9b144fb2-3580-49bd-9bde-\ +8f9e9d66dc8a"], + +Cell[CellGroupData[{ + +Cell[BoxData[{ + RowBox[{ + RowBox[{"g", "/.", "KillMeEquatorRules"}], "//", + "MatrixForm"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"g", "/.", "KillMeOnAxisRules"}], "//", "MatrixForm"}]}], "Input", + CellChangeTimes->{{3.886239099478282*^9, 3.8862391391676826`*^9}}, + CellLabel->"In[12]:=",ExpressionUUID->"bd17db4b-99a8-4a6c-9eb7-49accaf01cdd"], + +Cell[BoxData[ + TagBox[ + RowBox[{"(", "\[NoBreak]", GridBox[{ + { + RowBox[{"gtt", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], + RowBox[{"gtx", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], "0", "0"}, + { + RowBox[{"gtx", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], + RowBox[{"gxx", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], "0", "0"}, + {"0", "0", + RowBox[{"gyy", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], "0"}, + {"0", "0", "0", + RowBox[{"gzz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]} + }, + GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Baseline}}}, + GridBoxSpacings->{"Columns" -> { + Offset[0.27999999999999997`], { + Offset[0.7]}, + Offset[0.27999999999999997`]}, "Rows" -> { + Offset[0.2], { + Offset[0.4]}, + Offset[0.2]}}], "\[NoBreak]", ")"}], + Function[BoxForm`e$, + MatrixForm[BoxForm`e$]]]], "Output", + CellChangeTimes->{{3.8862391317376204`*^9, 3.8862391395839252`*^9}}, + CellLabel-> + "Out[12]//MatrixForm=",ExpressionUUID->"bb34dbee-88b7-4abe-9f81-\ +02f49b097c80"], + +Cell[BoxData[ + TagBox[ + RowBox[{"(", "\[NoBreak]", GridBox[{ + { + RowBox[{"gtt", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], "0", "0", + RowBox[{"gtz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}, + {"0", + RowBox[{"gxx", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], "0", "0"}, + {"0", "0", + RowBox[{"gyy", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], "0"}, + { + RowBox[{"gtz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], "0", "0", + RowBox[{"gzz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]} + }, + GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Baseline}}}, + GridBoxSpacings->{"Columns" -> { + Offset[0.27999999999999997`], { + Offset[0.7]}, + Offset[0.27999999999999997`]}, "Rows" -> { + Offset[0.2], { + Offset[0.4]}, + Offset[0.2]}}], "\[NoBreak]", ")"}], + Function[BoxForm`e$, + MatrixForm[BoxForm`e$]]]], "Output", + CellChangeTimes->{{3.8862391317376204`*^9, 3.8862391395856543`*^9}}, + CellLabel-> + "Out[13]//MatrixForm=",ExpressionUUID->"d86ab8bb-5d04-472e-a292-\ +412d5e32b5d3"] +}, Open ]] +}, Closed]], + +Cell[CellGroupData[{ + +Cell["Collect the relevant components on Axis:", "Text", + CellChangeTimes->{{3.886238862895907*^9, + 3.886238894746866*^9}},ExpressionUUID->"ac3861a3-e496-4de1-ad12-\ +e033917d7b88"], + +Cell[BoxData[ + RowBox[{ + RowBox[{ + RowBox[{ + RowBox[{"{", + RowBox[{ + RowBox[{"Christoffelsddd", "[", + RowBox[{"[", + RowBox[{"1", ",", "1", ",", "1"}], "]"}], "]"}], ",", + RowBox[{"Christoffelsddd", "[", + RowBox[{"[", + RowBox[{"4", ",", "1", ",", "1"}], "]"}], "]"}], ",", + RowBox[{"Christoffelsddd", "[", + RowBox[{"[", + RowBox[{"1", ",", "4", ",", "4"}], "]"}], "]"}], ",", + RowBox[{"Christoffelsddd", "[", + RowBox[{"[", + RowBox[{"4", ",", "4", ",", "4"}], "]"}], "]"}], ",", + RowBox[{"Christoffelsddd", "[", + RowBox[{"[", + RowBox[{"1", ",", "1", ",", "4"}], "]"}], "]"}], ",", + RowBox[{"Christoffelsddd", "[", + RowBox[{"[", + RowBox[{"4", ",", "1", ",", "4"}], "]"}], "]"}]}], "}"}], "/.", + "KillMeOnAxisRules"}], "//", "Simplify"}], ";"}]], "Input", + CellChangeTimes->{{3.886173999514325*^9, 3.8861740245071373`*^9}, + 3.886174252714652*^9, 3.886174383803364*^9, {3.8861749956068773`*^9, + 3.886174997474677*^9}, {3.886175279412593*^9, 3.886175386312916*^9}, { + 3.886175446572671*^9, 3.886175449889043*^9}, 3.886239866253716*^9}, + CellLabel->"In[62]:=",ExpressionUUID->"49dcf80e-4d26-47f8-aad5-d314e4df0e69"], + +Cell["\<\ +Double check that we can get the udd part of the metric by multiplying with \ +the inverse of just one block of the inverse:\ +\>", "Text", + CellChangeTimes->{{3.886238889485704*^9, 3.886238900333441*^9}, { + 3.886239153187326*^9, + 3.8862391961989613`*^9}},ExpressionUUID->"92d83e45-e318-4d91-be89-\ +af13b6883cab"], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"g", "/.", "KillMeOnAxisRules"}]], "Input", + CellChangeTimes->{{3.8862392389049063`*^9, 3.886239242557054*^9}}, + CellLabel->"In[14]:=",ExpressionUUID->"de8bd87c-f1a9-413f-9eb9-806bac2b9685"], + +Cell[BoxData[ + RowBox[{"{", + RowBox[{ + RowBox[{"{", + RowBox[{ + RowBox[{"gtt", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", "0", ",", "0", + ",", + RowBox[{"gtz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"0", ",", + RowBox[{"gxx", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", "0", ",", + "0"}], "}"}], ",", + RowBox[{"{", + RowBox[{"0", ",", "0", ",", + RowBox[{"gyy", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", "0"}], "}"}], + ",", + RowBox[{"{", + RowBox[{ + RowBox[{"gtz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", "0", ",", "0", + ",", + RowBox[{"gzz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}], "}"}]}], + "}"}]], "Output", + CellChangeTimes->{3.886239242979628*^9}, + CellLabel->"Out[14]=",ExpressionUUID->"76943539-2335-4cc0-9e7e-b0a0226fa3ad"] +}, Open ]], + +Cell[BoxData[ + RowBox[{ + RowBox[{"XAxi", "=", + RowBox[{"{", + RowBox[{"t", ",", "z"}], "}"}]}], ";"}]], "Input", + CellChangeTimes->{{3.886239392262786*^9, 3.886239399279251*^9}}, + CellLabel->"In[21]:=",ExpressionUUID->"01ae6f5e-d085-4a3b-928b-925b844bd0ac"], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{ + RowBox[{"gOnAxisBlock", "=", + RowBox[{"{", + RowBox[{ + RowBox[{"{", + RowBox[{ + RowBox[{"gtt", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", + RowBox[{"gtz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"gtz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", + RowBox[{"gzz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}], "}"}]}], + "}"}]}], ";", + RowBox[{"gOnAxisBlock", "//", "MatrixForm"}]}]], "Input", + CellChangeTimes->{{3.886239203828929*^9, 3.886239300860549*^9}}, + CellLabel->"In[19]:=",ExpressionUUID->"8d6186c2-bad4-4e3e-9b47-61d5bbf88a87"], + +Cell[BoxData[ + TagBox[ + RowBox[{"(", "\[NoBreak]", GridBox[{ + { + RowBox[{"gtt", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], + RowBox[{"gtz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}, + { + RowBox[{"gtz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], + RowBox[{"gzz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]} + }, + GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Baseline}}}, + GridBoxSpacings->{"Columns" -> { + Offset[0.27999999999999997`], { + Offset[0.7]}, + Offset[0.27999999999999997`]}, "Rows" -> { + Offset[0.2], { + Offset[0.4]}, + Offset[0.2]}}], "\[NoBreak]", ")"}], + Function[BoxForm`e$, + MatrixForm[BoxForm`e$]]]], "Output", + CellChangeTimes->{3.886239301125842*^9}, + CellLabel-> + "Out[19]//MatrixForm=",ExpressionUUID->"2977af70-299b-46e5-9314-\ +010cb4214309"] +}, Open ]], + +Cell[BoxData[ + RowBox[{ + RowBox[{"gOnAxisBlockinv", "=", + RowBox[{ + RowBox[{"gOnAxisBlock", "//", "Inverse"}], "//", "Simplify"}]}], + ";"}]], "Input", + CellChangeTimes->{{3.886239279053372*^9, 3.8862393025269327`*^9}}, + CellLabel->"In[20]:=",ExpressionUUID->"26c014b5-335e-41f3-a343-c5c438670da3"], + +Cell[BoxData[ + RowBox[{ + RowBox[{"ChristoffelsAxisBlockddd", "=", + RowBox[{ + RowBox[{"Table", "[", + RowBox[{ + RowBox[{ + RowBox[{"1", "/", "2"}], + RowBox[{"(", + RowBox[{ + RowBox[{ + SubscriptBox["\[PartialD]", + RowBox[{"XAxi", "[", + RowBox[{"[", "\[Mu]", "]"}], "]"}]], + RowBox[{"gOnAxisBlock", "[", + RowBox[{"[", + RowBox[{"\[Nu]", ",", "\[Delta]"}], "]"}], "]"}]}], "+", + RowBox[{ + SubscriptBox["\[PartialD]", + RowBox[{"XAxi", "[", + RowBox[{"[", "\[Nu]", "]"}], "]"}]], + RowBox[{"gOnAxisBlock", "[", + RowBox[{"[", + RowBox[{"\[Mu]", ",", "\[Delta]"}], "]"}], "]"}]}], "-", + RowBox[{ + SubscriptBox["\[PartialD]", + RowBox[{"XAxi", "[", + RowBox[{"[", "\[Delta]", "]"}], "]"}]], + RowBox[{"gOnAxisBlock", "[", + RowBox[{"[", + RowBox[{"\[Mu]", ",", "\[Nu]"}], "]"}], "]"}]}]}], ")"}]}], ",", + RowBox[{"{", + RowBox[{"\[Delta]", ",", "2"}], "}"}], ",", + RowBox[{"{", + RowBox[{"\[Mu]", ",", "2"}], "}"}], ",", + RowBox[{"{", + RowBox[{"\[Nu]", ",", "2"}], "}"}]}], "]"}], "//", "Simplify"}]}], + ";"}]], "Input", + CellChangeTimes->{{3.8862393767955112`*^9, 3.886239447571889*^9}}, + CellLabel->"In[25]:=",ExpressionUUID->"de4e5987-1579-4a27-b9fc-70eb429269ae"], + +Cell[BoxData[ + RowBox[{ + RowBox[{"ChristoffelsAxisBlockudd", "=", + RowBox[{"Table", "[", + RowBox[{ + RowBox[{"Sum", "[", + RowBox[{ + RowBox[{ + RowBox[{"gOnAxisBlockinv", "[", + RowBox[{"[", + RowBox[{"\[Alpha]", ",", "\[Delta]"}], "]"}], "]"}], + RowBox[{"ChristoffelsAxisBlockddd", "[", + RowBox[{"[", + RowBox[{"\[Delta]", ",", "\[Mu]", ",", "\[Nu]"}], "]"}], "]"}]}], + ",", + RowBox[{"{", + RowBox[{"\[Delta]", ",", "2"}], "}"}]}], "]"}], ",", + RowBox[{"{", + RowBox[{"\[Alpha]", ",", "2"}], "}"}], ",", + RowBox[{"{", + RowBox[{"\[Mu]", ",", "2"}], "}"}], ",", + RowBox[{"{", + RowBox[{"\[Nu]", ",", "2"}], "}"}]}], "]"}]}], ";"}]], "Input", + CellChangeTimes->{{3.8862394559284563`*^9, 3.886239526238406*^9}, { + 3.886239680673327*^9, 3.886239685179586*^9}}, + CellLabel->"In[36]:=",ExpressionUUID->"640e7de3-3183-400a-8794-02a2c9d32492"], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{ + RowBox[{ + RowBox[{"{", + RowBox[{ + RowBox[{ + RowBox[{"ChristoffelsAxisBlockudd", "[", + RowBox[{"[", + RowBox[{"1", ",", "1", ",", "1"}], "]"}], "]"}], "-", + RowBox[{"Christoffelsudd", "[", + RowBox[{"[", + RowBox[{"1", ",", "1", ",", "1"}], "]"}], "]"}]}], ",", + RowBox[{ + RowBox[{"ChristoffelsAxisBlockudd", "[", + RowBox[{"[", + RowBox[{"1", ",", "1", ",", "2"}], "]"}], "]"}], "-", + RowBox[{"Christoffelsudd", "[", + RowBox[{"[", + RowBox[{"1", ",", "1", ",", "4"}], "]"}], "]"}]}], ",", + RowBox[{ + RowBox[{"ChristoffelsAxisBlockudd", "[", + RowBox[{"[", + RowBox[{"1", ",", "2", ",", "2"}], "]"}], "]"}], "-", + RowBox[{"Christoffelsudd", "[", + RowBox[{"[", + RowBox[{"1", ",", "4", ",", "4"}], "]"}], "]"}]}], ",", + RowBox[{ + RowBox[{"ChristoffelsAxisBlockudd", "[", + RowBox[{"[", + RowBox[{"2", ",", "1", ",", "1"}], "]"}], "]"}], "-", + RowBox[{"Christoffelsudd", "[", + RowBox[{"[", + RowBox[{"4", ",", "1", ",", "1"}], "]"}], "]"}]}], ",", + RowBox[{ + RowBox[{"ChristoffelsAxisBlockudd", "[", + RowBox[{"[", + RowBox[{"2", ",", "1", ",", "2"}], "]"}], "]"}], "-", + RowBox[{"Christoffelsudd", "[", + RowBox[{"[", + RowBox[{"4", ",", "1", ",", "4"}], "]"}], "]"}]}], ",", + RowBox[{ + RowBox[{"ChristoffelsAxisBlockudd", "[", + RowBox[{"[", + RowBox[{"2", ",", "2", ",", "2"}], "]"}], "]"}], "-", + RowBox[{"Christoffelsudd", "[", + RowBox[{"[", + RowBox[{"4", ",", "4", ",", "4"}], "]"}], "]"}]}]}], "}"}], "/.", + "KillMeOnAxisRules"}], "//", "Simplify"}]], "Input", + CellChangeTimes->{{3.886239550225062*^9, 3.88623958206875*^9}, { + 3.886239692524212*^9, 3.8862398032153063`*^9}, {3.886240229048933*^9, + 3.886240259620287*^9}}, + CellLabel->"In[74]:=",ExpressionUUID->"06a60ae0-778e-4912-85dc-d32667a00362"], + +Cell[BoxData[ + RowBox[{"{", + RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], + "}"}]], "Output", + CellChangeTimes->{{3.886239569587137*^9, 3.8862395822670918`*^9}, { + 3.8862396864278193`*^9, 3.886239734674528*^9}, {3.8862397688325167`*^9, + 3.8862398035105886`*^9}, {3.886240252178359*^9, 3.886240264435549*^9}}, + CellLabel->"Out[74]=",ExpressionUUID->"9359abe8-5143-4c48-81db-e7d805acb053"] +}, Open ]], + +Cell["All is well!", "Text", + CellChangeTimes->{{3.886239880987774*^9, + 3.886239885684226*^9}},ExpressionUUID->"221e1d75-6a4b-4bc8-9548-\ +1d2a9bac055f"] +}, Closed]], + +Cell[CellGroupData[{ + +Cell["Collect the relevant components at the equator:", "Text", + CellChangeTimes->{{3.886238862895907*^9, 3.886238894746866*^9}, { + 3.8862399139244413`*^9, + 3.886239918455309*^9}},ExpressionUUID->"0c813465-1111-450c-abee-\ +81cee2d5e354"], + +Cell[BoxData[ + RowBox[{ + RowBox[{ + RowBox[{ + RowBox[{"{", + RowBox[{ + RowBox[{"Christoffelsddd", "[", + RowBox[{"[", + RowBox[{"1", ",", "1", ",", "1"}], "]"}], "]"}], ",", + RowBox[{"Christoffelsddd", "[", + RowBox[{"[", + RowBox[{"2", ",", "1", ",", "1"}], "]"}], "]"}], ",", + RowBox[{"Christoffelsddd", "[", + RowBox[{"[", + RowBox[{"1", ",", "2", ",", "2"}], "]"}], "]"}], ",", + RowBox[{"Christoffelsddd", "[", + RowBox[{"[", + RowBox[{"2", ",", "2", ",", "2"}], "]"}], "]"}], ",", + RowBox[{"Christoffelsddd", "[", + RowBox[{"[", + RowBox[{"1", ",", "1", ",", "2"}], "]"}], "]"}], ",", + RowBox[{"Christoffelsddd", "[", + RowBox[{"[", + RowBox[{"2", ",", "1", ",", "2"}], "]"}], "]"}]}], "}"}], "/.", + "KillMeEquatorRules"}], "//", "Simplify"}], ";"}]], "Input", + CellChangeTimes->{{3.886173999514325*^9, 3.8861740245071373`*^9}, + 3.886174252714652*^9, 3.886174383803364*^9, {3.8861749956068773`*^9, + 3.886174997474677*^9}, {3.886175279412593*^9, 3.886175386312916*^9}, { + 3.886175446572671*^9, 3.886175449889043*^9}, 3.886239866253716*^9, { + 3.8862399323925657`*^9, 3.886239950169238*^9}}, + CellLabel->"In[63]:=",ExpressionUUID->"d230e578-d7a3-4696-81a1-55a1e97d698f"], + +Cell["\<\ +Double check that we can get the udd part of the metric by multiplying with \ +the inverse of just one block of the inverse:\ +\>", "Text", + CellChangeTimes->{{3.886238889485704*^9, 3.886238900333441*^9}, { + 3.886239153187326*^9, + 3.8862391961989613`*^9}},ExpressionUUID->"1b418052-d781-4b3a-886c-\ +17105256f20d"], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"g", "/.", "KillMeEquatorRules"}]], "Input", + CellChangeTimes->{{3.8862392389049063`*^9, 3.886239242557054*^9}, { + 3.886239955096569*^9, 3.8862399559882517`*^9}}, + CellLabel->"In[75]:=",ExpressionUUID->"49bce66e-c2ac-4737-a18b-83b9af6efec3"], + +Cell[BoxData[ + RowBox[{"{", + RowBox[{ + RowBox[{"{", + RowBox[{ + RowBox[{"gtt", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", + RowBox[{"gtx", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", "0", ",", + "0"}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"gtx", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", + RowBox[{"gxx", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", "0", ",", + "0"}], "}"}], ",", + RowBox[{"{", + RowBox[{"0", ",", "0", ",", + RowBox[{"gyy", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", "0"}], "}"}], + ",", + RowBox[{"{", + RowBox[{"0", ",", "0", ",", "0", ",", + RowBox[{"gzz", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}], "}"}]}], + "}"}]], "Output", + CellChangeTimes->{3.886239242979628*^9, 3.886239956472385*^9, + 3.886240293779134*^9}, + CellLabel->"Out[75]=",ExpressionUUID->"dcbac780-95b9-4925-9a81-911c83eabaee"] +}, Open ]], + +Cell[BoxData[ + RowBox[{ + RowBox[{"XEq", "=", + RowBox[{"{", + RowBox[{"t", ",", "x"}], "}"}]}], ";"}]], "Input", + CellChangeTimes->{{3.886239392262786*^9, 3.886239399279251*^9}, { + 3.8862399619870443`*^9, 3.886239964013056*^9}}, + CellLabel->"In[76]:=",ExpressionUUID->"2f1d3424-d70f-481b-a66e-b1449c0fcc65"], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{ + RowBox[{"gEquatorBlock", "=", + RowBox[{"{", + RowBox[{ + RowBox[{"{", + RowBox[{ + RowBox[{"gtt", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", + RowBox[{"gtx", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"gtx", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", + RowBox[{"gxx", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}], "}"}]}], + "}"}]}], ";", + RowBox[{"gEquatorBlock", "//", "MatrixForm"}]}]], "Input", + CellChangeTimes->{{3.886239203828929*^9, 3.886239300860549*^9}, { + 3.8862399680307703`*^9, 3.886239991444262*^9}}, + CellLabel->"In[77]:=",ExpressionUUID->"4c2fad4b-2113-406e-9140-500b5acb0e14"], + +Cell[BoxData[ + TagBox[ + RowBox[{"(", "\[NoBreak]", GridBox[{ + { + RowBox[{"gtt", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], + RowBox[{"gtx", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}, + { + RowBox[{"gtx", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], + RowBox[{"gxx", "[", + RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]} + }, + GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Baseline}}}, + GridBoxSpacings->{"Columns" -> { + Offset[0.27999999999999997`], { + Offset[0.7]}, + Offset[0.27999999999999997`]}, "Rows" -> { + Offset[0.2], { + Offset[0.4]}, + Offset[0.2]}}], "\[NoBreak]", ")"}], + Function[BoxForm`e$, + MatrixForm[BoxForm`e$]]]], "Output", + CellChangeTimes->{3.886239301125842*^9, 3.886239992311648*^9, + 3.8862402203987093`*^9, 3.88624029644741*^9}, + CellLabel-> + "Out[77]//MatrixForm=",ExpressionUUID->"4debc2a1-1480-4b9c-9db2-\ +294463da8105"] +}, Open ]], + +Cell[BoxData[ + RowBox[{ + RowBox[{"gEquatorBlockinv", "=", + RowBox[{ + RowBox[{"gEquatorBlock", "//", "Inverse"}], "//", "Simplify"}]}], + ";"}]], "Input", + CellChangeTimes->{{3.886239279053372*^9, 3.8862393025269327`*^9}, { + 3.8862400002414837`*^9, 3.886240007499395*^9}}, + CellLabel->"In[78]:=",ExpressionUUID->"9f98fd8a-6fa5-4aaa-b85b-5ca08fae24a6"], + +Cell[BoxData[ + RowBox[{ + RowBox[{"ChristoffelsEquatorBlockddd", "=", + RowBox[{ + RowBox[{"Table", "[", + RowBox[{ + RowBox[{ + RowBox[{"1", "/", "2"}], + RowBox[{"(", + RowBox[{ + RowBox[{ + SubscriptBox["\[PartialD]", + RowBox[{"XEq", "[", + RowBox[{"[", "\[Mu]", "]"}], "]"}]], + RowBox[{"gEquatorBlock", "[", + RowBox[{"[", + RowBox[{"\[Nu]", ",", "\[Delta]"}], "]"}], "]"}]}], "+", + RowBox[{ + SubscriptBox["\[PartialD]", + RowBox[{"XEq", "[", + RowBox[{"[", "\[Nu]", "]"}], "]"}]], + RowBox[{"gEquatorBlock", "[", + RowBox[{"[", + RowBox[{"\[Mu]", ",", "\[Delta]"}], "]"}], "]"}]}], "-", + RowBox[{ + SubscriptBox["\[PartialD]", + RowBox[{"XEq", "[", + RowBox[{"[", "\[Delta]", "]"}], "]"}]], + RowBox[{"gEquatorBlock", "[", + RowBox[{"[", + RowBox[{"\[Mu]", ",", "\[Nu]"}], "]"}], "]"}]}]}], ")"}]}], ",", + RowBox[{"{", + RowBox[{"\[Delta]", ",", "2"}], "}"}], ",", + RowBox[{"{", + RowBox[{"\[Mu]", ",", "2"}], "}"}], ",", + RowBox[{"{", + RowBox[{"\[Nu]", ",", "2"}], "}"}]}], "]"}], "//", "Simplify"}]}], + ";"}]], "Input", + CellChangeTimes->{{3.8862393767955112`*^9, 3.886239447571889*^9}, { + 3.8862400174952917`*^9, 3.886240060859981*^9}}, + CellLabel->"In[80]:=",ExpressionUUID->"f8b887e2-1f20-49ec-8a4e-6bdd7b719e87"], + +Cell[BoxData[ + RowBox[{ + RowBox[{"ChristoffelsEquatorBlockudd", "=", + RowBox[{"Table", "[", + RowBox[{ + RowBox[{"Sum", "[", + RowBox[{ + RowBox[{ + RowBox[{"gEquatorBlockinv", "[", + RowBox[{"[", + RowBox[{"\[Alpha]", ",", "\[Delta]"}], "]"}], "]"}], + RowBox[{"ChristoffelsEquatorBlockddd", "[", + RowBox[{"[", + RowBox[{"\[Delta]", ",", "\[Mu]", ",", "\[Nu]"}], "]"}], "]"}]}], + ",", + RowBox[{"{", + RowBox[{"\[Delta]", ",", "2"}], "}"}]}], "]"}], ",", + RowBox[{"{", + RowBox[{"\[Alpha]", ",", "2"}], "}"}], ",", + RowBox[{"{", + RowBox[{"\[Mu]", ",", "2"}], "}"}], ",", + RowBox[{"{", + RowBox[{"\[Nu]", ",", "2"}], "}"}]}], "]"}]}], ";"}]], "Input", + CellChangeTimes->{{3.8862394559284563`*^9, 3.886239526238406*^9}, { + 3.886239680673327*^9, 3.886239685179586*^9}, {3.8862400721140003`*^9, + 3.88624009482312*^9}}, + CellLabel->"In[81]:=",ExpressionUUID->"c5674cf3-b1f7-4213-ab57-9659993e8b80"], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{ + RowBox[{ + RowBox[{"{", + RowBox[{ + RowBox[{ + RowBox[{"ChristoffelsEquatorBlockudd", "[", + RowBox[{"[", + RowBox[{"1", ",", "1", ",", "1"}], "]"}], "]"}], "-", + RowBox[{"Christoffelsudd", "[", + RowBox[{"[", + RowBox[{"1", ",", "1", ",", "1"}], "]"}], "]"}]}], ",", + RowBox[{ + RowBox[{"ChristoffelsEquatorBlockudd", "[", + RowBox[{"[", + RowBox[{"1", ",", "1", ",", "2"}], "]"}], "]"}], "-", + RowBox[{"Christoffelsudd", "[", + RowBox[{"[", + RowBox[{"1", ",", "1", ",", "2"}], "]"}], "]"}]}], ",", + "\[IndentingNewLine]", + RowBox[{ + RowBox[{"ChristoffelsEquatorBlockudd", "[", + RowBox[{"[", + RowBox[{"1", ",", "2", ",", "2"}], "]"}], "]"}], "-", + RowBox[{"Christoffelsudd", "[", + RowBox[{"[", + RowBox[{"1", ",", "2", ",", "2"}], "]"}], "]"}]}], ",", + "\[IndentingNewLine]", + RowBox[{ + RowBox[{"ChristoffelsEquatorBlockudd", "[", + RowBox[{"[", + RowBox[{"2", ",", "1", ",", "1"}], "]"}], "]"}], "-", + RowBox[{"Christoffelsudd", "[", + RowBox[{"[", + RowBox[{"2", ",", "1", ",", "1"}], "]"}], "]"}]}], ",", + "\[IndentingNewLine]", + RowBox[{ + RowBox[{"ChristoffelsEquatorBlockudd", "[", + RowBox[{"[", + RowBox[{"2", ",", "1", ",", "2"}], "]"}], "]"}], "-", + RowBox[{"Christoffelsudd", "[", + RowBox[{"[", + RowBox[{"2", ",", "1", ",", "2"}], "]"}], "]"}]}], ",", + "\[IndentingNewLine]", + RowBox[{ + RowBox[{"ChristoffelsEquatorBlockudd", "[", + RowBox[{"[", + RowBox[{"2", ",", "2", ",", "2"}], "]"}], "]"}], "-", + RowBox[{"Christoffelsudd", "[", + RowBox[{"[", + RowBox[{"2", ",", "2", ",", "2"}], "]"}], "]"}]}]}], "}"}], "/.", + "KillMeEquatorRules"}], "//", "Simplify"}]], "Input", + CellChangeTimes->{{3.886239550225062*^9, 3.88623958206875*^9}, { + 3.886239692524212*^9, 3.8862398032153063`*^9}, {3.886240109340292*^9, + 3.886240207302966*^9}}, + CellLabel->"In[82]:=",ExpressionUUID->"cd24a2c2-cec1-4e13-8ee2-25af439fceaa"], + +Cell[BoxData[ + RowBox[{"{", + RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], + "}"}]], "Output", + CellChangeTimes->{{3.886239569587137*^9, 3.8862395822670918`*^9}, { + 3.8862396864278193`*^9, 3.886239734674528*^9}, {3.8862397688325167`*^9, + 3.8862398035105886`*^9}, 3.886240207733821*^9, 3.886240302195738*^9}, + CellLabel->"Out[82]=",ExpressionUUID->"33e9b9d5-5aeb-406f-b752-17986db4beb9"] +}, Open ]], + +Cell["All is well!", "Text", + CellChangeTimes->{{3.886239880987774*^9, + 3.886239885684226*^9}},ExpressionUUID->"b4df87a0-82a6-4fbd-a126-\ +8cedb791ad67"] +}, Closed]] +}, +CellGrouping->Manual, +WindowSize->{1440, 852}, +WindowMargins->{{0, Automatic}, {0, Automatic}}, +Magnification:>1.5 Inherited, +FrontEndVersion->"13.1 for Linux x86 (64-bit) (June 16, 2022)", +StyleDefinitions->"Default.nb", +ExpressionUUID->"4f361d4c-d8b0-44e2-b48c-8cc05b0570b8" +] +(* End of Notebook Content *) + +(* Internal cache information *) +(*CellTagsOutline +CellTagsIndex->{} +*) +(*CellTagsIndex +CellTagsIndex->{} +*) +(*NotebookFileOutline +Notebook[{ +Cell[558, 20, 276, 7, 107, "Subsubsection",ExpressionUUID->"373d7a92-84df-4cc8-94b3-ee4da185a64c"], +Cell[CellGroupData[{ +Cell[859, 31, 186, 3, 44, "Text",ExpressionUUID->"7b3be704-efcd-4b3c-a35c-5fafcf85c005"], +Cell[1048, 36, 425, 8, 45, "Input",ExpressionUUID->"33ae6975-ec1d-4e40-8c00-336b4482e361"], +Cell[CellGroupData[{ +Cell[1498, 48, 2332, 53, 173, "Input",ExpressionUUID->"bfbd2c04-efbd-49b8-940f-e8e2a48cee58"], +Cell[3833, 103, 2223, 56, 147, "Output",ExpressionUUID->"afb62946-6d22-40db-87c7-06d8aaa5eafb"] +}, Open ]], +Cell[6071, 162, 277, 6, 45, "Input",ExpressionUUID->"ff8d5500-b0da-4dea-aed3-005d3f38ce38"] +}, Closed]], +Cell[CellGroupData[{ +Cell[6385, 173, 178, 3, 44, "Text",ExpressionUUID->"0e39b00c-5443-4cc9-9420-976408ee456d"], +Cell[6566, 178, 1392, 39, 49, "Input",ExpressionUUID->"ac869ef8-5c04-4ff7-b6ff-063dc324d8f4"], +Cell[7961, 219, 2017, 52, 80, "Input",ExpressionUUID->"e7e016e2-c00b-42d2-91be-a38628e97843"] +}, Closed]], +Cell[CellGroupData[{ +Cell[10015, 276, 253, 6, 44, "Text",ExpressionUUID->"c11437e5-f609-4db8-bb2d-e4ba921912d9"], +Cell[10271, 284, 1420, 39, 77, "Input",ExpressionUUID->"1f3a5dca-b140-4e51-93fd-93ff73818e14"], +Cell[11694, 325, 1202, 35, 77, "Input",ExpressionUUID->"9461b4d0-c922-4076-a54e-0d59230005dd"] +}, Closed]], +Cell[CellGroupData[{ +Cell[12933, 365, 183, 3, 44, "Text",ExpressionUUID->"9b144fb2-3580-49bd-9bde-8f9e9d66dc8a"], +Cell[CellGroupData[{ +Cell[13141, 372, 352, 7, 77, "Input",ExpressionUUID->"bd17db4b-99a8-4a6c-9eb7-49accaf01cdd"], +Cell[13496, 381, 1214, 33, 147, "Output",ExpressionUUID->"bb34dbee-88b7-4abe-9f81-02f49b097c80"], +Cell[14713, 416, 1214, 33, 147, "Output",ExpressionUUID->"d86ab8bb-5d04-472e-a292-412d5e32b5d3"] +}, Open ]] +}, Closed]], +Cell[CellGroupData[{ +Cell[15976, 455, 182, 3, 44, "Text",ExpressionUUID->"ac3861a3-e496-4de1-ad12-e033917d7b88"], +Cell[16161, 460, 1255, 29, 77, "Input",ExpressionUUID->"49dcf80e-4d26-47f8-aad5-d314e4df0e69"], +Cell[17419, 491, 324, 7, 52, "Text",ExpressionUUID->"92d83e45-e318-4d91-be89-af13b6883cab"], +Cell[CellGroupData[{ +Cell[17768, 502, 213, 3, 45, "Input",ExpressionUUID->"de8bd87c-f1a9-413f-9eb9-806bac2b9685"], +Cell[17984, 507, 995, 29, 83, "Output",ExpressionUUID->"76943539-2335-4cc0-9e7e-b0a0226fa3ad"] +}, Open ]], +Cell[18994, 539, 264, 6, 45, "Input",ExpressionUUID->"01ae6f5e-d085-4a3b-928b-925b844bd0ac"], +Cell[CellGroupData[{ +Cell[19283, 549, 768, 20, 45, "Input",ExpressionUUID->"8d6186c2-bad4-4e3e-9b47-61d5bbf88a87"], +Cell[20054, 571, 948, 27, 94, "Output",ExpressionUUID->"2977af70-299b-46e5-9314-010cb4214309"] +}, Open ]], +Cell[21017, 601, 306, 7, 45, "Input",ExpressionUUID->"26c014b5-335e-41f3-a343-c5c438670da3"], +Cell[21326, 610, 1444, 39, 115, "Input",ExpressionUUID->"de4e5987-1579-4a27-b9fc-70eb429269ae"], +Cell[22773, 651, 959, 25, 77, "Input",ExpressionUUID->"640e7de3-3183-400a-8794-02a2c9d32492"], +Cell[CellGroupData[{ +Cell[23757, 680, 2000, 51, 205, "Input",ExpressionUUID->"06a60ae0-778e-4912-85dc-d32667a00362"], +Cell[25760, 733, 413, 7, 51, "Output",ExpressionUUID->"9359abe8-5143-4c48-81db-e7d805acb053"] +}, Open ]], +Cell[26188, 743, 154, 3, 52, "Text",ExpressionUUID->"221e1d75-6a4b-4bc8-9548-1d2a9bac055f"] +}, Closed]], +Cell[CellGroupData[{ +Cell[26379, 751, 240, 4, 44, "Text",ExpressionUUID->"0c813465-1111-450c-abee-81cee2d5e354"], +Cell[26622, 757, 1308, 30, 77, "Input",ExpressionUUID->"d230e578-d7a3-4696-81a1-55a1e97d698f"], +Cell[27933, 789, 324, 7, 52, "Text",ExpressionUUID->"1b418052-d781-4b3a-886c-17105256f20d"], +Cell[CellGroupData[{ +Cell[28282, 800, 265, 4, 45, "Input",ExpressionUUID->"49bce66e-c2ac-4737-a18b-83b9af6efec3"], +Cell[28550, 806, 1036, 29, 83, "Output",ExpressionUUID->"dcbac780-95b9-4925-9a81-911c83eabaee"] +}, Open ]], +Cell[29601, 838, 314, 7, 45, "Input",ExpressionUUID->"2f1d3424-d70f-481b-a66e-b1449c0fcc65"], +Cell[CellGroupData[{ +Cell[29940, 849, 821, 21, 45, "Input",ExpressionUUID->"4c2fad4b-2113-406e-9140-500b5acb0e14"], +Cell[30764, 872, 1018, 28, 94, "Output",ExpressionUUID->"4debc2a1-1480-4b9c-9db2-294463da8105"] +}, Open ]], +Cell[31797, 903, 359, 8, 45, "Input",ExpressionUUID->"9f98fd8a-6fa5-4aaa-b85b-5ca08fae24a6"], +Cell[32159, 913, 1498, 40, 115, "Input",ExpressionUUID->"f8b887e2-1f20-49ec-8a4e-6bdd7b719e87"], +Cell[33660, 955, 1016, 26, 77, "Input",ExpressionUUID->"c5674cf3-b1f7-4213-ab57-9659993e8b80"], +Cell[CellGroupData[{ +Cell[34701, 985, 2135, 55, 237, "Input",ExpressionUUID->"cd24a2c2-cec1-4e13-8ee2-25af439fceaa"], +Cell[36839, 1042, 413, 7, 51, "Output",ExpressionUUID->"33e9b9d5-5aeb-406f-b752-17986db4beb9"] +}, Open ]], +Cell[37267, 1052, 154, 3, 52, "Text",ExpressionUUID->"b4df87a0-82a6-4fbd-a126-8cedb791ad67"] +}, Closed]] +} +] +*) + -- cgit v1.2.3