tiny changes

1 files changed, 273 insertions, 171 deletions

diff --git a/tempering_analysis.nb b/tempering_analysis.nb
index 044b448..d5e07d8 100644
--- a/tempering_analysis.nb
+++ b/tempering_analysis.nb
@@ -10,61 +10,131 @@
 NotebookFileLineBreakTest
 NotebookFileLineBreakTest
 NotebookDataPosition[       158,          7]
-NotebookDataLength[     33282,        767]
-NotebookOptionsPosition[     31337,        730]
-NotebookOutlinePosition[     31675,        745]
-CellTagsIndexPosition[     31632,        742]
+NotebookDataLength[     37301,        869]
+NotebookOptionsPosition[     35385,        834]
+NotebookOutlinePosition[     35723,        849]
+CellTagsIndexPosition[     35680,        846]
 WindowFrame->Normal*)
 
 (* Beginning of Notebook Content *)
 Notebook[{
 Cell[BoxData[
  RowBox[{
+  RowBox[{"getPData", "[", "dir_", "]"}], ":=", 
+  RowBox[{
+   RowBox[{
+    RowBox[{"With", "[", 
+     RowBox[{
+      RowBox[{"{", 
+       RowBox[{
+        RowBox[{"data", "=", 
+         RowBox[{"Import", "[", "#", "]"}]}], ",", 
+        RowBox[{"metadata", "=", 
+         RowBox[{"ToExpression", "/@", 
+          RowBox[{"Rest", "[", 
+           RowBox[{"StringSplit", "[", 
+            RowBox[{
+             RowBox[{"FileBaseName", "[", "#", "]"}], ",", "\"\<_\>\""}], 
+            "]"}], "]"}]}]}]}], "}"}], ",", 
+      RowBox[{"<|", 
+       RowBox[{
+        RowBox[{"\"\<n\>\"", "\[Rule]", 
+         RowBox[{"metadata", "[", 
+          RowBox[{"[", "1", "]"}], "]"}]}], ",", 
+        RowBox[{"\"\<\[Beta]0\>\"", "\[Rule]", 
+         RowBox[{"metadata", "[", 
+          RowBox[{"[", "2", "]"}], "]"}]}], ",", 
+        RowBox[{"\"\<\[Beta]f\>\"", "\[Rule]", 
+         RowBox[{"metadata", "[", 
+          RowBox[{"[", "3", "]"}], "]"}]}], ",", 
+        RowBox[{"\"\<n\[Beta]\>\"", "\[Rule]", 
+         RowBox[{"metadata", "[", 
+          RowBox[{"[", "4", "]"}], "]"}]}], ",", 
+        RowBox[{"\"\<As\>\"", "\[Rule]", 
+         RowBox[{
+          RowBox[{"Rest", "[", "data", "]"}], "/", 
+          RowBox[{"First", "[", 
+           RowBox[{"First", "[", "data", "]"}], "]"}]}]}]}], "|>"}]}], "]"}], 
+    "&"}], "/@", 
+   RowBox[{"FileNames", "[", 
+    RowBox[{"FileNameJoin", "[", 
+     RowBox[{"{", 
+      RowBox[{"dir", ",", "\"\<probs_*.dat\>\""}], "}"}], "]"}], 
+    "]"}]}]}]], "Input",
+ CellChangeTimes->{{3.778948155736003*^9, 3.778948220094614*^9}, {
+   3.778948251774682*^9, 3.778948445776267*^9}, 3.7849915996103563`*^9, {
+   3.784991674882406*^9, 3.784991676833317*^9}},
+ CellLabel->"In[1]:=",ExpressionUUID->"6ca216e4-4a7b-4e5f-a718-f4e70c77d6a2"],
+
+Cell[BoxData[
+ RowBox[{
+  RowBox[{"getEData", "[", "dir_", "]"}], ":=", 
+  RowBox[{
+   RowBox[{
+    RowBox[{"With", "[", 
+     RowBox[{
+      RowBox[{"{", 
+       RowBox[{
+        RowBox[{"data", "=", 
+         RowBox[{"Import", "[", "#", "]"}]}], ",", 
+        RowBox[{"metadata", "=", 
+         RowBox[{"ToExpression", "/@", 
+          RowBox[{"Rest", "[", 
+           RowBox[{"StringSplit", "[", 
+            RowBox[{
+             RowBox[{"FileBaseName", "[", "#", "]"}], ",", "\"\<_\>\""}], 
+            "]"}], "]"}]}]}]}], "}"}], ",", 
+      RowBox[{"<|", 
+       RowBox[{
+        RowBox[{"\"\<n\>\"", "\[Rule]", 
+         RowBox[{"metadata", "[", 
+          RowBox[{"[", "1", "]"}], "]"}]}], ",", 
+        RowBox[{"\"\<\[Beta]0\>\"", "\[Rule]", 
+         RowBox[{"metadata", "[", 
+          RowBox[{"[", "2", "]"}], "]"}]}], ",", 
+        RowBox[{"\"\<\[Beta]f\>\"", "\[Rule]", 
+         RowBox[{"metadata", "[", 
+          RowBox[{"[", "3", "]"}], "]"}]}], ",", 
+        RowBox[{"\"\<n\[Beta]\>\"", "\[Rule]", 
+         RowBox[{"metadata", "[", 
+          RowBox[{"[", "4", "]"}], "]"}]}], ",", 
+        RowBox[{"\"\<Es\>\"", "\[Rule]", 
+         RowBox[{
+          RowBox[{"First", "@", 
+           RowBox[{"Rest", "[", "data", "]"}]}], "/", 
+          RowBox[{"First", "[", 
+           RowBox[{"First", "[", "data", "]"}], "]"}]}]}]}], "|>"}]}], "]"}], 
+    "&"}], "/@", 
+   RowBox[{"FileNames", "[", 
+    RowBox[{"FileNameJoin", "[", 
+     RowBox[{"{", 
+      RowBox[{"dir", ",", "\"\<energies_*.dat\>\""}], "}"}], "]"}], 
+    "]"}]}]}]], "Input",
+ CellChangeTimes->{{3.778948155736003*^9, 3.778948220094614*^9}, {
+  3.778948251774682*^9, 3.778948445776267*^9}, {3.7849915996103563`*^9, 
+  3.784991614545053*^9}, {3.784991671691276*^9, 3.784991674145091*^9}, {
+  3.784991874442992*^9, 3.784991874617395*^9}, {3.784992011338725*^9, 
+  3.784992012353281*^9}},
+ CellLabel->"In[14]:=",ExpressionUUID->"1ffc3385-0105-4d64-94b0-d5fc2c51fc67"],
+
+Cell[BoxData[
+ RowBox[{
   RowBox[{"getData", "[", "dir_", "]"}], ":=", 
   RowBox[{"Dataset", "[", 
    RowBox[{
-    RowBox[{
-     RowBox[{"With", "[", 
+    RowBox[{"Merge", "[", "First", "]"}], "/@", 
+    RowBox[{"GatherBy", "[", 
+     RowBox[{
+      RowBox[{"Join", "[", 
+       RowBox[{
+        RowBox[{"getPData", "[", "dir", "]"}], ",", 
+        RowBox[{"getEData", "[", "dir", "]"}]}], "]"}], ",", " ", 
       RowBox[{
        RowBox[{"{", 
-        RowBox[{
-         RowBox[{"data", "=", 
-          RowBox[{"Import", "[", "#", "]"}]}], ",", 
-         RowBox[{"metadata", "=", 
-          RowBox[{"ToExpression", "/@", 
-           RowBox[{"Rest", "[", 
-            RowBox[{"StringSplit", "[", 
-             RowBox[{
-              RowBox[{"FileBaseName", "[", "#", "]"}], ",", "\"\<_\>\""}], 
-             "]"}], "]"}]}]}]}], "}"}], ",", 
-       RowBox[{"<|", 
-        RowBox[{
-         RowBox[{"\"\<n\>\"", "\[Rule]", 
-          RowBox[{"metadata", "[", 
-           RowBox[{"[", "1", "]"}], "]"}]}], ",", 
-         RowBox[{"\"\<\[Beta]0\>\"", "\[Rule]", 
-          RowBox[{"metadata", "[", 
-           RowBox[{"[", "2", "]"}], "]"}]}], ",", 
-         RowBox[{"\"\<\[Beta]f\>\"", "\[Rule]", 
-          RowBox[{"metadata", "[", 
-           RowBox[{"[", "3", "]"}], "]"}]}], ",", 
-         RowBox[{"\"\<n\[Beta]\>\"", "\[Rule]", 
-          RowBox[{"metadata", "[", 
-           RowBox[{"[", "4", "]"}], "]"}]}], ",", 
-         RowBox[{"\"\<As\>\"", "\[Rule]", 
-          RowBox[{
-           RowBox[{"Rest", "[", "data", "]"}], "/", 
-           RowBox[{"First", "[", 
-            RowBox[{"First", "[", "data", "]"}], "]"}]}]}]}], "|>"}]}], "]"}],
-      "&"}], "/@", 
-    RowBox[{"FileNames", "[", 
-     RowBox[{"FileNameJoin", "[", 
-      RowBox[{"{", 
-       RowBox[{"dir", ",", "\"\<probs_*.dat\>\""}], "}"}], "]"}], "]"}]}], 
-   "]"}]}]], "Input",
- CellChangeTimes->{{3.778948155736003*^9, 3.778948220094614*^9}, {
-  3.778948251774682*^9, 3.778948445776267*^9}},
- CellLabel->"In[12]:=",ExpressionUUID->"6ca216e4-4a7b-4e5f-a718-f4e70c77d6a2"],
+        RowBox[{"#n", ",", "#\[Beta]0", ",", "#\[Beta]f", ",", "#n\[Beta]"}], 
+        "}"}], "&"}]}], "]"}]}], "]"}]}]], "Input",
+ CellChangeTimes->{{3.784991617461961*^9, 3.784991728472725*^9}},
+ CellLabel->"In[3]:=",ExpressionUUID->"10f87413-28fb-4b72-bd1b-eaeef3d02b15"],
 
 Cell[BoxData[
  RowBox[{
@@ -84,7 +154,7 @@ Cell[BoxData[
     RowBox[{"InverseErfc", "[", "A", "]"}], "2"]}]}]], "Input",
  CellChangeTimes->{{3.7789451133983593`*^9, 3.778945185413146*^9}, {
   3.778945446992125*^9, 3.778945447478743*^9}},
- CellLabel->"In[14]:=",ExpressionUUID->"f756d73a-4921-45d7-acd1-44fe21438a02"],
+ CellLabel->"In[4]:=",ExpressionUUID->"f756d73a-4921-45d7-acd1-44fe21438a02"],
 
 Cell[BoxData[
  RowBox[{
@@ -104,7 +174,7 @@ Cell[BoxData[
    "]"}]}]], "Input",
  CellChangeTimes->{{3.778948523412203*^9, 3.778948610089134*^9}, {
   3.7789486512266293`*^9, 3.778948654153482*^9}},
- CellLabel->"In[15]:=",ExpressionUUID->"c91e85b2-f90e-40d7-bf19-7a03fd9499bb"],
+ CellLabel->"In[5]:=",ExpressionUUID->"c91e85b2-f90e-40d7-bf19-7a03fd9499bb"],
 
 Cell[BoxData[
  RowBox[{
@@ -130,7 +200,7 @@ Cell[BoxData[
        RowBox[{"dat", "[", "\"\<As\>\"", "]"}]}]}], "}"}]}], "]"}]}]], "Input",\
 
  CellChangeTimes->{{3.7789537137281647`*^9, 3.7789537203664207`*^9}},
- CellLabel->"In[16]:=",ExpressionUUID->"1b068f0a-fa61-487d-893c-5b00f279e758"],
+ CellLabel->"In[6]:=",ExpressionUUID->"1b068f0a-fa61-487d-893c-5b00f279e758"],
 
 Cell[BoxData[
  RowBox[{
@@ -180,7 +250,32 @@ Cell[BoxData[
   3.778948688346896*^9}, {3.7789536224795513`*^9, 3.778953623110043*^9}, {
   3.778953673879937*^9, 3.778953685238799*^9}, {3.778953742582839*^9, 
   3.778953745606895*^9}},
- CellLabel->"In[17]:=",ExpressionUUID->"83b164e2-d1bf-4be5-bf92-d3645cb950ca"],
+ CellLabel->"In[7]:=",ExpressionUUID->"83b164e2-d1bf-4be5-bf92-d3645cb950ca"],
+
+Cell[BoxData[
+ RowBox[{
+  RowBox[{"Z0", "[", "n_", "]"}], ":=", 
+  RowBox[{
+   SuperscriptBox["2", 
+    RowBox[{
+     RowBox[{"(", 
+      RowBox[{"n", "-", "1"}], ")"}], 
+     RowBox[{"(", 
+      RowBox[{
+       RowBox[{"n", "/", "4"}], "+", "1"}], ")"}]}]], 
+   RowBox[{"Product", "[", 
+    RowBox[{
+     RowBox[{
+      SuperscriptBox["\[Pi]", 
+       RowBox[{"k", "/", "2"}]], 
+      RowBox[{"Gamma", "[", 
+       RowBox[{"k", "/", "2"}], "]"}]}], ",", 
+     RowBox[{"{", 
+      RowBox[{"k", ",", "2", ",", "n"}], "}"}]}], "]"}]}]}]], "Input",
+ CellChangeTimes->{{3.784991222937964*^9, 3.784991266472834*^9}, {
+  3.784991403570047*^9, 3.784991404177*^9}, {3.7849915386807756`*^9, 
+  3.784991539440803*^9}},
+ CellLabel->"In[8]:=",ExpressionUUID->"f81a2199-caac-407b-95bf-2211d7b265a5"],
 
 Cell[BoxData[
  RowBox[{
@@ -234,8 +329,10 @@ Cell[BoxData[
   3.779700575188902*^9}, {3.779700627710291*^9, 3.7797006535654173`*^9}, {
   3.779700716133662*^9, 3.7797008200306168`*^9}, {3.779700857824486*^9, 
   3.7797008584386597`*^9}, {3.7797008994884863`*^9, 3.779700901055118*^9}, {
-  3.779700939144706*^9, 3.779700939215414*^9}},
- CellLabel->"In[18]:=",ExpressionUUID->"cfeba3ce-dc48-40c8-b8b8-f34cf59a253d"],
+  3.779700939144706*^9, 3.779700939215414*^9}, {3.784991004834445*^9, 
+  3.7849910388807697`*^9}, {3.78499127505939*^9, 3.78499130611341*^9}, {
+  3.784991358538706*^9, 3.784991380665337*^9}},
+ CellLabel->"In[9]:=",ExpressionUUID->"cfeba3ce-dc48-40c8-b8b8-f34cf59a253d"],
 
 Cell[BoxData[
  RowBox[{
@@ -256,8 +353,14 @@ Cell[BoxData[
          RowBox[{"[", 
           RowBox[{"All", ",", "2"}], "]"}], "]"}], "]"}]}], "}"}], "]"}]}], 
    "]"}]}]], "Input",
- CellChangeTimes->{{3.779700960584425*^9, 3.7797010023114433`*^9}},
- CellLabel->"In[19]:=",ExpressionUUID->"a1830ef0-911d-442d-b5d2-b62c202fce89"],
+ CellChangeTimes->{{3.779700960584425*^9, 3.7797010023114433`*^9}, {
+  3.7849910515540743`*^9, 3.784991075832299*^9}, {3.784991369146078*^9, 
+  3.7849913764007483`*^9}, {3.784991547938244*^9, 3.784991556568716*^9}, {
+  3.7849917900181932`*^9, 3.784991803720964*^9}, {3.784991863633103*^9, 
+  3.784991884800845*^9}, {3.784992068826951*^9, 3.7849920762732773`*^9}, {
+  3.784992122226707*^9, 3.7849921237610617`*^9}, {3.7849921543708973`*^9, 
+  3.784992167089695*^9}, {3.784992239922765*^9, 3.784992242761147*^9}},
+ CellLabel->"In[26]:=",ExpressionUUID->"a1830ef0-911d-442d-b5d2-b62c202fce89"],
 
 Cell[BoxData[
  RowBox[{
@@ -266,7 +369,7 @@ Cell[BoxData[
   ";"}]], "Input",
  CellChangeTimes->{{3.778948447393362*^9, 3.778948476191824*^9}, {
   3.7795618755287123`*^9, 3.779561880399358*^9}},
- CellLabel->"In[37]:=",ExpressionUUID->"3b7ef50e-b362-4da3-b734-818bb17a7c2f"],
+ CellLabel->"In[34]:=",ExpressionUUID->"3b7ef50e-b362-4da3-b734-818bb17a7c2f"],
 
 Cell[CellGroupData[{
 
@@ -301,28 +404,28 @@ K9prVXjaPqLGeudBDSZNLYdYHNxHr9C4GjGrXHmBXWKyHA6pf26K6LcVLyG4
 jl2JYrLw36wHPi1Jdk+i5lhu6QdmVtnYpy04DEvtfMM7tag834ifO/rNuOGV
 uMKyvmNZfqKNMMvC23HG5Ge2kHpBEJKLo+uuhTjzWVmE/wGo7glY
   "],
- CellLabel->"In[46]:=",ExpressionUUID->"73943d9b-70f3-4562-ab84-3881106d7e0c"],
+ CellLabel->"In[35]:=",ExpressionUUID->"73943d9b-70f3-4562-ab84-3881106d7e0c"],
 
 Cell[BoxData[
  GraphicsBox[{{}, {{{}, {}, 
      {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.007333333333333334],
        AbsoluteThickness[1.6], LineBox[CompressedData["
-1:eJw1yX8s1HEcx3E/MstRuHK5k993x3EnSYzyfZuNobpiLW2t/Go3UUdIUoro
-bKoty2ZSGelO6+iIsNSZULj8TMf5I+t33fy4kB/tVPP+vLfX3nvs6RQvjjpp
-ZGBgQP3b/79+asrRbNTrVnkyte5PVFLkmkwRnYfWUmr5o9uvHtxA6yhxi+5g
-8OtS9CJ1Iir+fEH+XfQKtZaeYC6tqkbrqXLTwbn7wlq0IVR2/EhXJNWjjcH5
-+2xVv3Uj2gRCslVT1M9mtCmUPpb/pnW1oTdCS3e+3kb5Ak2DEu0H6zJJB9oC
-suRpQSPCTvRmyMoalhbmdq37iiW0sneZSYp60FagLSyhGex/g7YGpjGDX23a
-h6ZDgOWzXL/j/egtIAwyFNOnVOitoD9mmsGMGkDbQO/1yn2H/AbRDCiu2St8
-qSHeBvdUEZ3Lp4bQtqC0lbFaF4iZILMyn489M4xmgYwShl/QENtBzPNpUVrY
-CHo73DzM5bs8IbaH4KX8nW2sUbQDiJYiOWNiNDhC4Di7ljeEVjrCHx1LEid4
-h90J6hY3rGZmopVO4H92wk0xQLozmAyKk3WMMezOkBjXSStLRYMLtDdMeqe0
-k+4Cp70lFW/nSHeFA7HzfiF73mN3hRi7zyZhl9HAhoRJC+a3TtLZUJ9tZcTT
-k86BcBGPa+Kvxs4B+x2JffQcNHAhuO5i/3Qz6VyYKdJGSOdIdwN2T3gkJ3Ac
-uxuYa7p7Ui+hwR0s3UU+9DrS3SHiV3ac4AvpPGjx6RbFOkxg54HUpji09iga
-PMA3oFV34xrpHjDzsdqWaifdE85VqHYvaEn3hK8B8HDVVoOdD9GMxqu9YWgl
-H3I2USl3MkgXwBG1nD1cQ7oAQi2afPP6SPeCphXFcsMs6V4wVs4ueOoySf0F
-QLMyxA==
+1:eJw1yQ0s1HEcx/HrIswajjzHnftzd9y5ZdxM4/9ttpJs8jAqustDDRlTzUpT
+1JRqxYaURCs2auO6nplc6elwDqfcXSbjOlYpWm4YVs33990+++61Nyc9P/4w
+k8Fg0P/2/6+fjk6PLPAeKcyh122kF/MtqHZmKfoHHXpoSDNhuIr+TTePDozn
+TtegzXSgKtzxnU0Depk2DUtEOVFN6DVaNyE0ln1pRW+AEy+44XaadvRG6KqK
+KL92/hHaEvJKs5lBx56hraCAappfSOpE28CNhor2uuButC30ys0lIdWv0Juh
+pynmtiKtB20H383Hw54nvV13iT20uULI3cL3aAfwvzC7d2e0Cs2C2LoZfQez
+D+0IRpGsovpIP9oJtPcsrXavqNFbwM30cn9CqgbtDC6bWjzlEYNoFyia+1Z0
+aozYFazVU7RzxhDaDQo6Wm11X4ndgTU5LjGlDqM9QKNfGXY3EHtChnp1znmH
+Fr0V5hcGI60bib3ARctnem0aQXuDmTrgtJqCBjb0WgbXxnWhlWwIS+nsLbT7
+iJ0Dic0npZFJaCUHHtJVxWfuk+4DtRx3lmKRdB8wbPv55mnsJ+xc2GMuk4lv
+opVcUH/Qno2eIJ2CRn5xXJBgFDsFXfXsoXO5aPAFj2TVTFYb6b5wObNVZzVL
+uh/UcMUqIU+H3Q+uxPCk0jw08OBo8sqirZx0HjBTVl//MpHOB4nDxfpigR47
+HxIVS0+mM9AggIRdk3f+1JAugFssWYi+n3R/4DSrnaIZBuz+wChzFTiEoiEA
+MmuXri9lkR4AsrRMaVQ96UJwNIz6ZPWRLoQxZXmLYo10EVBTRpFc9Bm7CB7n
+HmQX7UNDIFRK+vXyS6QHQvzpBx2V3aSLgdo+YFw2kS6GbItp+3mPMfov2T42
+7A==
        "]]}, 
      {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.007333333333333334],
        AbsoluteThickness[1.6], LineBox[CompressedData["
@@ -454,8 +557,14 @@ ADfH
    3.7800874839642763`*^9}, {3.784987528733613*^9, 3.7849875536676188`*^9}, 
    3.7849876962409067`*^9, 3.7849880617038116`*^9, 3.784988269808079*^9, 
    3.7849886106727123`*^9, 3.784989205039175*^9, {3.784989237650683*^9, 
-   3.784989241626642*^9}, {3.7849892794456577`*^9, 3.784989301810355*^9}},
- CellLabel->"Out[46]=",ExpressionUUID->"520bf791-28f4-41af-9c2f-24068e5eb363"]
+   3.784989241626642*^9}, {3.7849892794456577`*^9, 3.784989301810355*^9}, 
+   3.784991010484539*^9, {3.7849910412866077`*^9, 3.7849910778183727`*^9}, 
+   3.784991309318139*^9, 3.784991386947817*^9, {3.784991549707535*^9, 
+   3.7849915579184523`*^9}, 3.784991736881611*^9, {3.784991879079419*^9, 
+   3.784991887673251*^9}, 3.7849919728655863`*^9, 3.784992017379757*^9, 
+   3.78499208028115*^9, 3.784992125434637*^9, {3.784992155783885*^9, 
+   3.784992168646801*^9}, 3.7849922441034927`*^9, 3.784992833706853*^9},
+ CellLabel->"Out[35]=",ExpressionUUID->"4e0edb5f-c69d-4119-94ee-c0c8715ba2f1"]
 }, Open  ]],
 
 Cell[CellGroupData[{
@@ -485,28 +594,28 @@ Cell[BoxData[
   3.780087398738736*^9}, {3.7800877026590443`*^9, 3.780087726448501*^9}, {
   3.78498770229777*^9, 3.7849877073353243`*^9}, {3.78498924891437*^9, 
   3.7849892518156357`*^9}},
- CellLabel->"In[44]:=",ExpressionUUID->"60949ea9-18bb-4b5b-868e-6a12e7d013f6"],
+ CellLabel->"In[36]:=",ExpressionUUID->"60949ea9-18bb-4b5b-868e-6a12e7d013f6"],
 
 Cell[BoxData[
  GraphicsBox[{{}, {{{}, {}, 
      {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.007333333333333334],
        AbsoluteThickness[1.6], LineBox[CompressedData["
-1:eJw10XtIU3EUB3DNTNRi5nO6zdfmY/O1RC2LugdF6WGIK4lgE7EWpmWmm6NE
-HWZJDhMU35kV7Z+StFkDK/PWfM18wzT/EPMZ2nC9SLOxHnh+Fy6XD9/fPef8
-fj+/jCsi6Q4rKyvq3/v/u/18oBrs5haNnTJq20vUPm1j5gOnCrSRijNms+cr
-a9HfqImY3pZsUSv6J5XyKa/tAKjRW1T8SuIxRlc72kKV88+I1UdfoK2hoe5r
-5tB6F9oGmFWpAZXmHrQttP2wdZAIe9F24BykL0rl9aPt4fxi59rO+UG0I5x4
-vTq46jaC3gMWo2bj7LQezYBiztOOldXxbSud4P6sOls1SbwX8mq+1yesj6Kd
-weDoYtSbx9Au8OtPx/4k5hDaFZo7vnTf3SLr3eBtjY3R3ULsDl4uNQkb1thf
-6QHpqunrkUKSM+H5k7IYYdZ7tCcsCOqavB8Re8GCSDPcs4X7U7LAcyEn5Wru
-AJoNPrrSCVsHknPgQrKjbNOX2Bsaew2J5oekvw9ctldYPmrwPMEXxt1EsXIt
-rqd9obQ700MsRYMfmCwjh28E6zD3g6j4aO3FXVgP/OFWT725uQnr0f4wXVe4
-eXOT5FzoV7TIx1g4L82FU8yWtTQBjTkPIg2qsOU7JOcBa2N336GXfZgHQKuu
-/FmhHu+bDoAu2cEZ6TDOA4Ega3q8HMEh/QNh8lxFetQU1oMgmJpclla9wXp0
-EFwLLP1cMETqB0N+9WmuzwyeNx0MCk0ab5ScF/AhfeBIgauQ7J8PlZKq26Zq
-nAcEMHMpsbZMTuYXwO+Z42y5mJxfCCi56mQJh/wfAvfqe+bGtGS+UFDkZ81l
-MMh8ofAqtMhpKZbcTxiUdBdHq0xkf2EwvMjwZBtIHg7mTIl6No70D4e6nIJO
-Nusd5hGQbCrRJPHRdARIc8XtupN66i9WdTxG
+1:eJw10ntI01EUB3BtyTKztJxuWm1zmzof06hpFvk7Ba3HVHQ5yWGmtZUoPaRM
+tHD7JywhfFQqJFoyJ6UWlhpp6EozH5mPVmpmuVpCIJZTtllmDzz3Bz9+fPje
+e+45lx/32BmZaoWdnR317/3/XX5GqXzzDQ/22Dlq2SZqIEopyVi8ip6mMrS6
+6ficErSZ+tbg5PFDWYm2UGLmlpjMJB36J6Xuqt+a5dWAXqIiGz1naiyNaHu4
+yWS10Hc8QdPgkGr+LUfwDO0AxqMpu99rO9F0kAqLKk8c70I7Qm7rn6q9wm60
+E3iabAFt9b1oZyiP/ehZwiNeBxnZRtqI/PWyNS6Q7+4mP5g3hHaFhV8K02Pr
+AHo9aOf5U+kG4g1g3nVFxvDqQ7tB3/a03uyIQTQD5l0KzBOF/Wh3GBrdrbtz
+sgftAeGrDkTbrK/QTAgJq3pk30lyFpilp8turcF+NZ4g32ZvnqnG+TRecC38
+bH7UFMk3gmln6H3vGbwPzSZoPXzPvUL3Ar0ZnB9+V4gKyX42TAh4OV+5mAMH
+Bhmy8MxmzPUcgMT0S0/VJOcCe2HlMI/RgTkXWgaSb6sSX2LuDbVpFxRNkXi+
+3htSOzjK5FycB3gwydCpJSyS80BsodmSZkl9PsTFlsYnuGI9PR+aqpdSWG9I
+LoBUo8O4yojWC8BQyL9M209yHxhPU1eMhOL/ofeB7ILQ+cpocr4vyOaUlOsH
+Ut8XQuTJlrwkXA9+UOMoOtUehuv1fjAQHTdCL8V+QQhl3c/rme1kfiH87hG3
+WRkk94c6pkIrZZH+/MG/WHpRHEPuJwDaDIZuyzDpLwBM16Oy/CiyPxDK99ma
+HfpJHgjn67LW5qwm+4OAVes1N5FN6geBdZbe/OkL9gMiWEw9Uj2xh8wnAspR
+4lT0mcwfDA/uWovfRZD+gyFhTDrJlvRSfwGOLjJQ
        "]]}, 
      {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.007333333333333334],
        AbsoluteThickness[1.6], LineBox[CompressedData["
@@ -626,64 +735,72 @@ oPtb9a38pvH77cxffAM2LA==
    3.780087378630178*^9, 3.780087399174458*^9}, {3.780087705097999*^9, 
    3.780087726696293*^9}, {3.784987697054266*^9, 3.7849877075856333`*^9}, 
    3.7849880628367*^9, 3.784988270798539*^9, 3.7849886153459663`*^9, 
-   3.78498920646163*^9, 3.784989252209024*^9},
- CellLabel->"Out[44]=",ExpressionUUID->"4b818f0b-341a-496c-8216-10c283d7a417"]
+   3.78498920646163*^9, 3.784989252209024*^9, 3.7849913833561907`*^9, 
+   3.784992834391604*^9},
+ CellLabel->"Out[36]=",ExpressionUUID->"8d2214e2-fecc-4e88-ad31-f8070f2e3de8"]
 }, Open  ]],
 
-Cell[BoxData[
- RowBox[{
-  RowBox[{"a", "=", 
-   RowBox[{
-   "Import", "[", 
-    "\"\<~/doc/research/hadamard/code/energies_32_0.250000_16.000000_64.dat\>\
-\"", "]"}]}], ";"}]], "Input",
- CellChangeTimes->{{3.779699855009756*^9, 3.779699866112481*^9}, {
-  3.784988180321373*^9, 3.784988193494982*^9}, {3.784988626865487*^9, 
-  3.784988627055429*^9}, {3.784989217849276*^9, 3.7849892179355593`*^9}},
- CellLabel->"In[40]:=",ExpressionUUID->"bb64ee8f-2ab2-44da-86d3-071e144df4ac"],
-
 Cell[CellGroupData[{
 
 Cell[BoxData[
- RowBox[{"ListPlot", "[", "a", "]"}]], "Input",
- CellChangeTimes->{{3.7796998674408073`*^9, 3.779699868848329*^9}},
- CellLabel->"In[41]:=",ExpressionUUID->"27200d6e-a8c7-4ded-bf66-0ba360a8d0c7"],
+ RowBox[{"ListPlot", "[", 
+  RowBox[{
+   RowBox[{"d", "[", 
+    RowBox[{
+     RowBox[{"Select", "[", 
+      RowBox[{
+       RowBox[{"KeyExistsQ", "[", 
+        RowBox[{"#", ",", "\"\<Es\>\""}], "]"}], "&"}], "]"}], ",", 
+     RowBox[{
+      RowBox[{"With", "[", 
+       RowBox[{
+        RowBox[{"{", 
+         RowBox[{"L", "=", 
+          SuperscriptBox["#n", "2"]}], "}"}], ",", 
+        RowBox[{"Thread", "[", 
+         RowBox[{"{", 
+          RowBox[{
+           RowBox[{"Ts", "[", "#", "]"}], ",", 
+           RowBox[{"#Es", "/", 
+            SuperscriptBox["#n", "2"]}]}], "}"}], "]"}]}], "]"}], "&"}]}], 
+    "]"}], ",", 
+   RowBox[{"Joined", "\[Rule]", "True"}]}], "]"}]], "Input",
+ CellChangeTimes->{{3.784992259322311*^9, 3.7849923458651857`*^9}, {
+  3.784992651707872*^9, 3.784992651985375*^9}, {3.784992846100277*^9, 
+  3.784992888545384*^9}},
+ CellLabel->"In[40]:=",ExpressionUUID->"ee8c7b42-c7b8-4a5f-a611-21e18e0fcbb2"],
 
 Cell[BoxData[
- GraphicsBox[{{}, {{
+ GraphicsBox[{{}, {{{}, {}, 
      {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
-      0.011000000000000001`], AbsoluteThickness[1.6], 
-      PointBox[{{1., 100000.}, {1., 100000.}}]}, 
-     {RGBColor[0.880722, 0.611041, 0.142051], PointSize[
-      0.011000000000000001`], AbsoluteThickness[1.6], PointBox[CompressedData["
-1:eJw1yQlIE3AUx/FZWuKBYmFGopk6p2jMI2vT6W9ublM3l1JDPOY8yprHjNSS
-lEaaOrAyDyhRIalwWi4rbR2KShBaZFiLYRkaYcdGrRBKoxP//wePx+d9/fI1
-6QfWMRgM/Nv/d21s8emhd548Dw7hkwe0xoHvU0F2E2t0hOekWmL0ZxO746K5
-pvDU71jizTjTaVDFVQiJvTCAsHclK0nE3sh8kCG95C0n3o7FGdvw8HIa8Q74
-m41WFnsfcQAUQtmSWqwgZkJwtNuDF5dBzEJv4ojGsiWTOASDy3OV7k+ziEMR
-r7Dp7pflEO+Eint98PU3JTEbC/UJnxChIg6Hoo6t9lmgjsC2OcYjZUMecSRa
-GseumHzziaNQ4KRobhui3gV290vmOU4BcTRGS5kd0R6FxLshsrCmb7gdJN6D
-pjCeqZWvJuZA52SSFl2j5kIbvP54rUsxcQyUHu9tj7OpY1Feaz7CvkrNwyuZ
-bdO9j9RxUEdxRFkBJcTxeCjwahMpqQF3e8lceCuxFvjD+1GhmKSdj6n+pkT7
-L7Tz8ZnfNN/jVUp6Am6ft2TM8om1CZDH1Lhai2gX4EJE+WG9jnYB2seXvq7o
-aRfC4MCr75qmXQiZxLJy4gPtibDrVbHq7MtIT4R09ecLqy8xQwRv20ozk0O7
-COZacZ4qjXYxcrP7q0yHaBcjzCFA2XOSdgmONbYX6ztol2BDV6DnZT3tSdg4
-dKtzYJT2JET2jQzrn9GeDLcaF+fKt7QnI1epbT6+THsKmL9GjX0OGtJTYJid
-XyzxJGZIUVl1M6svkHYpXLh+Y6ejaJfh7pvo0Gwh7TJUt/jozqbRngpDEDd1
-fy7tqXA05DjrSmmXw6KYGc+ppl0OnTMYWxto3wvrqkOAa6tm4i9L0ywS
-       "]]}}, {
+      0.011000000000000001`], AbsoluteThickness[1.6], LineBox[CompressedData["
+
+1:eJw1yX0s1AEcx/FDGnYpDiEPh7tzx7k7tZA//D7rLIfRhppal5iHabOOGmVy
+t86GEu0ILSG3rtWK1LhleUhbam0eWjJbs1aJJFft6KZQzX2/22ffvfYOyDqV
+kmPLYrHwb///5o0zlU87U+vOfR3c9Humi5PX7jFEnmc4i8YnhWvk78ydrvOG
+3ohFq82M+FWK7bYCsoU5lHRfa2kj/2FEHwslxyfILKyvCj4YN8i2iJs2coal
+36zegnInYYlCSd4KvW+IpvQi2QFMk6OqpofsBI9bBq1lhszGiWxfryr7Jaud
+8daUoDWJyDvwup6t9E4ku2DJkCyLzCe7YkGXxhdryRyoDENl8hayG2rb/A9o
+HpLdcaZ5Jt7vOdkDqfcCB8unyDvRP8Yz2M6RPTG7V72fMZO9wOZXOETZmKz2
+xs/RB7+Pscm7UN3smDjiQfZBuC4uy5lL9kWRpDsqQUj2Q/5S8WSYjOyPHqkj
+VxFJ5oKdXpVrF2O1hguFy6/0z3LqAXip4ipH46kHIFo3tYGD1AMh9xz4EpNG
+PRAd4yv6tXTqQRh7YbN8REk9CMUJmp6GTOo81KV0rj7Loc7DjdbLI1fzqfMR
+nBW9Ii2gzof5LKu8RkVdALnbmxjLaeoCNJq7J+1KqAejKFyd51pKPRi7f2xX
+vCujLsRsv/90n5q6EM31NfaxF6iLcDJqjtVfQV0EPWOXIaikHoLMxAhebzX1
+EKznBmUMXaIeihXZozH7WuqheLywJ9b9CnUxhtsHCkJ11MUQ6bL3meuph6HP
+5+hdy1XqYejglDY0NlGXwIcpVR++Rl2CQv3tT97XqUthKZ9wb2yhLkXSvJ8x
+qZW6DKKby8l9babBv1tsSh4=
+       "]]}}, 
+    {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.011000000000000001`],
+      AbsoluteThickness[1.6]}, {
      {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
-      0.011000000000000001`], AbsoluteThickness[1.6]}, 
-     {RGBColor[0.880722, 0.611041, 0.142051], PointSize[
-      0.011000000000000001`], AbsoluteThickness[1.6]}, {}, {}, {}}, {
+      0.011000000000000001`], AbsoluteThickness[1.6]}, {}}, {
      {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
-      0.011000000000000001`], AbsoluteThickness[1.6]}, 
-     {RGBColor[0.880722, 0.611041, 0.142051], PointSize[
-      0.011000000000000001`], AbsoluteThickness[1.6]}, {}, {}, {}}}, {{}, {}}},
+      0.011000000000000001`], AbsoluteThickness[1.6]}, {}}}, {{}, {}}},
   AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
   Axes->{True, True},
   AxesLabel->{None, None},
-  AxesOrigin->{0., 0},
+  AxesOrigin->{0, -0.8513855643374515},
   DisplayFunction->Identity,
   Frame->{{False, False}, {False, False}},
   FrameLabel->{{None, None}, {None, None}},
@@ -702,7 +819,7 @@ fy7tqXA05DjrSmmXw6KYGc+ppl0OnTMYWxto3wvrqkOAa6tm4i9L0ywS
          Part[#, 1]], 
         (Identity[#]& )[
          Part[#, 2]]}& )}},
-  PlotRange->{{0., 64.}, {-1.655204112573*^7, 857321.899835}},
+  PlotRange->{{0, 16.}, {-0.9882480184724375, -0.8585888513971874}},
   PlotRangeClipping->True,
   PlotRangePadding->{{
      Scaled[0.02], 
@@ -710,26 +827,13 @@ fy7tqXA05DjrSmmXw6KYGc+ppl0OnTMYWxto3wvrqkOAa6tm4i9L0ywS
      Scaled[0.05], 
      Scaled[0.05]}},
   Ticks->{Automatic, Automatic}]], "Output",
- CellChangeTimes->{3.7796998691208344`*^9, 3.784988194503076*^9, 
-  3.7849886277986927`*^9, 3.784989220788183*^9},
- CellLabel->"Out[41]=",ExpressionUUID->"96361a65-d311-41a5-b251-e9f9bd1db01c"]
-}, Open  ]],
-
-Cell[CellGroupData[{
-
-Cell[BoxData[
- RowBox[{
-  RowBox[{"16", "/", "64"}], "//", "N"}]], "Input",
- CellChangeTimes->{{3.7796999517774973`*^9, 3.779699965625023*^9}},
- CellLabel->"In[12]:=",ExpressionUUID->"85d1832b-b409-4be7-b8d2-26df59eff8fb"],
-
-Cell[BoxData["0.25`"], "Output",
- CellChangeTimes->{{3.779699954939776*^9, 3.779699966125688*^9}},
- CellLabel->"Out[12]=",ExpressionUUID->"5dbe4cc3-7c89-4c96-870e-625a9bf27876"]
+ CellChangeTimes->{{3.784992312221224*^9, 3.784992346103744*^9}, 
+   3.7849926522781267`*^9, {3.7849928355405827`*^9, 3.7849928888523283`*^9}},
+ CellLabel->"Out[40]=",ExpressionUUID->"0fd34c5b-9bf7-4c9c-b984-ed2f76da4c8d"]
 }, Open  ]]
 },
 WindowSize->{1916, 1196},
-WindowMargins->{{Automatic, -1916}, {Automatic, 0}},
+WindowMargins->{{Automatic, -1918}, {Automatic, 2}},
 FrontEndVersion->"12.0 for Linux x86 (64-bit) (April 8, 2019)",
 StyleDefinitions->"Default.nb"
 ]
@@ -744,30 +848,28 @@ CellTagsIndex->{}
 *)
 (*NotebookFileOutline
 Notebook[{
-Cell[558, 20, 1783, 46, 101, "Input",ExpressionUUID->"6ca216e4-4a7b-4e5f-a718-f4e70c77d6a2"],
-Cell[2344, 68, 669, 18, 59, "Input",ExpressionUUID->"f756d73a-4921-45d7-acd1-44fe21438a02"],
-Cell[3016, 88, 723, 18, 31, "Input",ExpressionUUID->"c91e85b2-f90e-40d7-bf19-7a03fd9499bb"],
-Cell[3742, 108, 805, 24, 31, "Input",ExpressionUUID->"1b068f0a-fa61-487d-893c-5b00f279e758"],
-Cell[4550, 134, 1754, 48, 31, "Input",ExpressionUUID->"83b164e2-d1bf-4be5-bf92-d3645cb950ca"],
-Cell[6307, 184, 1969, 53, 52, "Input",ExpressionUUID->"cfeba3ce-dc48-40c8-b8b8-f34cf59a253d"],
-Cell[8279, 239, 713, 20, 31, "Input",ExpressionUUID->"a1830ef0-911d-442d-b5d2-b62c202fce89"],
-Cell[8995, 261, 337, 7, 31, "Input",ExpressionUUID->"3b7ef50e-b362-4da3-b734-818bb17a7c2f"],
-Cell[CellGroupData[{
-Cell[9357, 272, 1299, 31, 39, "Input",ExpressionUUID->"73943d9b-70f3-4562-ab84-3881106d7e0c"],
-Cell[10659, 305, 8153, 152, 245, "Output",ExpressionUUID->"520bf791-28f4-41af-9c2f-24068e5eb363"]
-}, Open  ]],
+Cell[558, 20, 1786, 46, 78, "Input",ExpressionUUID->"6ca216e4-4a7b-4e5f-a718-f4e70c77d6a2"],
+Cell[2347, 68, 1946, 49, 78, "Input",ExpressionUUID->"1ffc3385-0105-4d64-94b0-d5fc2c51fc67"],
+Cell[4296, 119, 678, 17, 31, "Input",ExpressionUUID->"10f87413-28fb-4b72-bd1b-eaeef3d02b15"],
+Cell[4977, 138, 668, 18, 59, "Input",ExpressionUUID->"f756d73a-4921-45d7-acd1-44fe21438a02"],
+Cell[5648, 158, 722, 18, 31, "Input",ExpressionUUID->"c91e85b2-f90e-40d7-bf19-7a03fd9499bb"],
+Cell[6373, 178, 804, 24, 31, "Input",ExpressionUUID->"1b068f0a-fa61-487d-893c-5b00f279e758"],
+Cell[7180, 204, 1753, 48, 31, "Input",ExpressionUUID->"83b164e2-d1bf-4be5-bf92-d3645cb950ca"],
+Cell[8936, 254, 786, 23, 39, "Input",ExpressionUUID->"f81a2199-caac-407b-95bf-2211d7b265a5"],
+Cell[9725, 279, 2112, 55, 52, "Input",ExpressionUUID->"cfeba3ce-dc48-40c8-b8b8-f34cf59a253d"],
+Cell[11840, 336, 1157, 26, 31, "Input",ExpressionUUID->"a1830ef0-911d-442d-b5d2-b62c202fce89"],
+Cell[13000, 364, 337, 7, 31, "Input",ExpressionUUID->"3b7ef50e-b362-4da3-b734-818bb17a7c2f"],
 Cell[CellGroupData[{
-Cell[18849, 462, 1028, 25, 39, "Input",ExpressionUUID->"60949ea9-18bb-4b5b-868e-6a12e7d013f6"],
-Cell[19880, 489, 7361, 140, 245, "Output",ExpressionUUID->"4b818f0b-341a-496c-8216-10c283d7a417"]
+Cell[13362, 375, 1299, 31, 39, "Input",ExpressionUUID->"73943d9b-70f3-4562-ab84-3881106d7e0c"],
+Cell[14664, 408, 8586, 158, 245, "Output",ExpressionUUID->"4e0edb5f-c69d-4119-94ee-c0c8715ba2f1"]
 }, Open  ]],
-Cell[27256, 632, 476, 10, 31, "Input",ExpressionUUID->"bb64ee8f-2ab2-44da-86d3-071e144df4ac"],
 Cell[CellGroupData[{
-Cell[27757, 646, 208, 3, 31, "Input",ExpressionUUID->"27200d6e-a8c7-4ded-bf66-0ba360a8d0c7"],
-Cell[27968, 651, 2914, 63, 209, "Output",ExpressionUUID->"96361a65-d311-41a5-b251-e9f9bd1db01c"]
+Cell[23287, 571, 1028, 25, 39, "Input",ExpressionUUID->"60949ea9-18bb-4b5b-868e-6a12e7d013f6"],
+Cell[24318, 598, 7415, 141, 245, "Output",ExpressionUUID->"8d2214e2-fecc-4e88-ad31-f8070f2e3de8"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[30919, 719, 222, 4, 31, "Input",ExpressionUUID->"85d1832b-b409-4be7-b8d2-26df59eff8fb"],
-Cell[31144, 725, 177, 2, 35, "Output",ExpressionUUID->"5dbe4cc3-7c89-4c96-870e-625a9bf27876"]
+Cell[31770, 744, 947, 26, 39, "Input",ExpressionUUID->"ee8c7b42-c7b8-4a5f-a611-21e18e0fcbb2"],
+Cell[32720, 772, 2649, 59, 224, "Output",ExpressionUUID->"0fd34c5b-9bf7-4c9c-b984-ed2f76da4c8d"]
 }, Open  ]]
 }
 ]