add a mathematica notebook computing Christoffels.

1 files changed, 1153 insertions, 0 deletions

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]]
+}
+]
+*)
+