From 046b292abd7ce75e8491f33fd5ab8d949c202885 Mon Sep 17 00:00:00 2001 From: Jaron Kent-Dobias Date: Mon, 14 Oct 2019 20:27:43 -0400 Subject: initialized repo --- tempering_analysis.nb | 706 ++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 706 insertions(+) create mode 100644 tempering_analysis.nb (limited to 'tempering_analysis.nb') diff --git a/tempering_analysis.nb b/tempering_analysis.nb new file mode 100644 index 0000000..1b6c9aa --- /dev/null +++ b/tempering_analysis.nb @@ -0,0 +1,706 @@ +(* Content-type: application/vnd.wolfram.mathematica *) + +(*** Wolfram Notebook File ***) +(* http://www.wolfram.com/nb *) + +(* CreatedBy='Mathematica 12.0' *) + +(*CacheID: 234*) +(* Internal cache information: +NotebookFileLineBreakTest +NotebookFileLineBreakTest +NotebookDataPosition[ 158, 7] +NotebookDataLength[ 29108, 698] +NotebookOptionsPosition[ 27171, 661] +NotebookOutlinePosition[ 27505, 676] +CellTagsIndexPosition[ 27462, 673] +WindowFrame->Normal*) + +(* Beginning of Notebook Content *) +Notebook[{ +Cell[BoxData[ + RowBox[{ + RowBox[{"getData", "[", "dir_", "]"}], ":=", + RowBox[{"Dataset", "[", + 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[{"\"\\"", "\[Rule]", + RowBox[{"metadata", "[", + RowBox[{"[", "1", "]"}], "]"}]}], ",", + RowBox[{"\"\<\[Beta]0\>\"", "\[Rule]", + RowBox[{"metadata", "[", + RowBox[{"[", "2", "]"}], "]"}]}], ",", + RowBox[{"\"\<\[Beta]f\>\"", "\[Rule]", + RowBox[{"metadata", "[", + RowBox[{"[", "3", "]"}], "]"}]}], ",", + RowBox[{"\"\\"", "\[Rule]", + RowBox[{"metadata", "[", + RowBox[{"[", "4", "]"}], "]"}]}], ",", + RowBox[{"\"\\"", "\[Rule]", + RowBox[{ + RowBox[{"Rest", "[", "data", "]"}], "/", + RowBox[{"First", "[", + RowBox[{"First", "[", "data", "]"}], "]"}]}]}]}], "|>"}]}], "]"}], + "&"}], "/@", + RowBox[{"FileNames", "[", + RowBox[{"FileNameJoin", "[", + RowBox[{"{", + RowBox[{"dir", ",", "\"\\""}], "}"}], "]"}], "]"}]}], + "]"}]}]], "Input", + CellChangeTimes->{{3.778948155736003*^9, 3.778948220094614*^9}, { + 3.778948251774682*^9, 3.778948445776267*^9}}, + CellLabel->"In[1]:=",ExpressionUUID->"6ca216e4-4a7b-4e5f-a718-f4e70c77d6a2"], + +Cell[BoxData[ + RowBox[{ + RowBox[{"d", "=", + RowBox[{"getData", "[", "\"\<~/doc/research/hadamard/code/\>\"", "]"}]}], + ";"}]], "Input", + CellChangeTimes->{{3.778948447393362*^9, 3.778948476191824*^9}, { + 3.7795618755287123`*^9, 3.779561880399358*^9}}, + CellLabel->"In[15]:=",ExpressionUUID->"3b7ef50e-b362-4da3-b734-818bb17a7c2f"], + +Cell[BoxData[ + RowBox[{ + RowBox[{ + RowBox[{"c", "[", + RowBox[{"\[Beta]i_", ",", "\[Beta]j_"}], "]"}], "[", "A_", "]"}], ":=", + RowBox[{"2", + FractionBox[ + RowBox[{"(", + RowBox[{ + SuperscriptBox["\[Beta]i", "2"], "+", + SuperscriptBox["\[Beta]j", "2"]}], ")"}], + SuperscriptBox[ + RowBox[{"(", + RowBox[{"\[Beta]i", "-", "\[Beta]j"}], ")"}], "2"]], + SuperscriptBox[ + RowBox[{"InverseErfc", "[", "A", "]"}], "2"]}]}]], "Input", + CellChangeTimes->{{3.7789451133983593`*^9, 3.778945185413146*^9}, { + 3.778945446992125*^9, 3.778945447478743*^9}}, + CellLabel->"In[3]:=",ExpressionUUID->"f756d73a-4921-45d7-acd1-44fe21438a02"], + +Cell[BoxData[ + RowBox[{ + RowBox[{"Ts", "[", "dat_", "]"}], ":=", + RowBox[{"Range", "[", + RowBox[{ + RowBox[{"dat", "[", "\"\<\[Beta]0\>\"", "]"}], ",", + RowBox[{"dat", "[", "\"\<\[Beta]f\>\"", "]"}], ",", + RowBox[{ + RowBox[{"(", + RowBox[{ + RowBox[{"dat", "[", "\"\<\[Beta]f\>\"", "]"}], "-", + RowBox[{"dat", "[", "\"\<\[Beta]0\>\"", "]"}]}], ")"}], "/", + RowBox[{"(", + RowBox[{ + RowBox[{"dat", "[", "\"\\"", "]"}], "-", "1"}], ")"}]}]}], + "]"}]}]], "Input", + CellChangeTimes->{{3.778948523412203*^9, 3.778948610089134*^9}, { + 3.7789486512266293`*^9, 3.778948654153482*^9}}, + CellLabel->"In[4]:=",ExpressionUUID->"c91e85b2-f90e-40d7-bf19-7a03fd9499bb"], + +Cell[BoxData[ + RowBox[{ + RowBox[{"Cstmp", "[", "dat_", "]"}], ":=", + RowBox[{"With", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"ts", "=", + RowBox[{"Ts", "[", "dat", "]"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{ + RowBox[{ + RowBox[{"{", + RowBox[{"#", ",", + RowBox[{"#", "+", + RowBox[{"ts", "[", + RowBox[{"[", "2", "]"}], "]"}], "-", + RowBox[{"ts", "[", + RowBox[{"[", "1", "]"}], "]"}]}]}], "}"}], "&"}], "/@", "ts"}], + ",", + RowBox[{"Last", "/@", + RowBox[{"dat", "[", "\"\\"", "]"}]}]}], "}"}]}], "]"}]}]], "Input",\ + + CellChangeTimes->{{3.7789537137281647`*^9, 3.7789537203664207`*^9}}, + CellLabel->"In[5]:=",ExpressionUUID->"1b068f0a-fa61-487d-893c-5b00f279e758"], + +Cell[BoxData[ + RowBox[{ + RowBox[{"Cs", "[", "dat_", "]"}], ":=", + RowBox[{"With", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"ts", "=", + RowBox[{"Ts", "[", "dat", "]"}]}], "}"}], ",", + RowBox[{"MapThread", "[", + RowBox[{ + RowBox[{ + RowBox[{"{", + RowBox[{ + RowBox[{"Mean", "[", + RowBox[{"{", + RowBox[{ + RowBox[{"#1", "[", + RowBox[{"[", "1", "]"}], "]"}], ",", + RowBox[{"#1", "[", + RowBox[{"[", "2", "]"}], "]"}]}], "}"}], "]"}], ",", + RowBox[{ + RowBox[{"c", "[", + RowBox[{ + RowBox[{"#1", "[", + RowBox[{"[", "1", "]"}], "]"}], ",", + RowBox[{"#1", "[", + RowBox[{"[", "2", "]"}], "]"}]}], "]"}], "[", "#2", "]"}]}], + "}"}], "&"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{ + RowBox[{ + RowBox[{"{", + RowBox[{"#", ",", + RowBox[{"#", "+", + RowBox[{"ts", "[", + RowBox[{"[", "2", "]"}], "]"}], "-", + RowBox[{"ts", "[", + RowBox[{"[", "1", "]"}], "]"}]}]}], "}"}], "&"}], "/@", + RowBox[{"Most", "[", "ts", "]"}]}], ",", + RowBox[{"Last", "/@", + RowBox[{"dat", "[", "\"\\"", "]"}]}]}], "}"}]}], "]"}]}], + "]"}]}]], "Input", + CellChangeTimes->{{3.7789479000625353`*^9, 3.7789480157738333`*^9}, { + 3.778948503834503*^9, 3.778948518745183*^9}, {3.778948666921554*^9, + 3.778948688346896*^9}, {3.7789536224795513`*^9, 3.778953623110043*^9}, { + 3.778953673879937*^9, 3.778953685238799*^9}, {3.778953742582839*^9, + 3.778953745606895*^9}}, + CellLabel->"In[6]:=",ExpressionUUID->"83b164e2-d1bf-4be5-bf92-d3645cb950ca"], + +Cell[BoxData[ + RowBox[{ + RowBox[{"\[CapitalDelta]Ss", "[", "dat_", "]"}], ":=", + RowBox[{"With", "[", + RowBox[{ + RowBox[{"{", + RowBox[{ + RowBox[{"ts", "=", + RowBox[{"1", "/", + RowBox[{"Ts", "[", "dat", "]"}]}]}], ",", + RowBox[{"cs", "=", + RowBox[{"Cs", "[", "dat", "]"}]}]}], "}"}], ",", + RowBox[{ + RowBox[{ + RowBox[{"{", + RowBox[{ + RowBox[{"Mean", "[", + RowBox[{"1", "/", "#1"}], "]"}], ",", + FractionBox[ + RowBox[{"#2", + RowBox[{"(", + RowBox[{ + RowBox[{"#1", "[", + RowBox[{"[", "1", "]"}], "]"}], "-", + RowBox[{"#1", "[", + RowBox[{"[", "2", "]"}], "]"}]}], ")"}]}], + RowBox[{"Mean", "[", "#1", "]"}]]}], "}"}], "&"}], "@@@", + RowBox[{"Thread", "[", + RowBox[{"{", + RowBox[{ + RowBox[{"Rest", "@", + RowBox[{"Rest", "[", + RowBox[{"FoldList", "[", + RowBox[{ + RowBox[{ + RowBox[{"{", + RowBox[{ + RowBox[{"#1", "[", + RowBox[{"[", "2", "]"}], "]"}], ",", "#2"}], "}"}], "&"}], + ",", + RowBox[{"{", + RowBox[{"0", ",", "0"}], "}"}], ",", "ts"}], "]"}], "]"}]}], ",", + RowBox[{ + RowBox[{"Cs", "[", "dat", "]"}], "[", + RowBox[{"[", + RowBox[{"All", ",", "2"}], "]"}], "]"}]}], "}"}], "]"}]}]}], + "]"}]}]], "Input", + CellChangeTimes->{{3.779561994673626*^9, 3.7795620342479887`*^9}, { + 3.779700498628436*^9, 3.779700513292344*^9}, {3.7797005747090063`*^9, + 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[7]:=",ExpressionUUID->"cfeba3ce-dc48-40c8-b8b8-f34cf59a253d"], + +Cell[BoxData[ + RowBox[{ + RowBox[{"Ss", "[", "dat_", "]"}], ":=", + RowBox[{"With", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"ss", "=", + RowBox[{"\[CapitalDelta]Ss", "[", "dat", "]"}]}], "}"}], ",", + RowBox[{"Thread", "[", + RowBox[{"{", + RowBox[{ + RowBox[{"ss", "[", + RowBox[{"[", + RowBox[{"All", ",", "1"}], "]"}], "]"}], ",", + RowBox[{"Accumulate", "[", + RowBox[{"ss", "[", + RowBox[{"[", + RowBox[{"All", ",", "2"}], "]"}], "]"}], "]"}]}], "}"}], "]"}]}], + "]"}]}]], "Input", + CellChangeTimes->{{3.779700960584425*^9, 3.7797010023114433`*^9}}, + CellLabel->"In[8]:=",ExpressionUUID->"a1830ef0-911d-442d-b5d2-b62c202fce89"], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"ListPlot", "[", + RowBox[{ + RowBox[{"d", "[", + RowBox[{ + RowBox[{"Select", "[", + RowBox[{ + RowBox[{"#n\[Beta]", "\[Equal]", "64"}], "&"}], "]"}], ",", + RowBox[{ + RowBox[{"With", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"L", "=", + SuperscriptBox["#n", "2"]}], "}"}], ",", + RowBox[{ + RowBox[{ + RowBox[{"{", + RowBox[{"#1", ",", + RowBox[{"#2", "/", "L"}]}], "}"}], "&"}], "@@@", + RowBox[{"Ss", "[", "#", "]"}]}]}], "]"}], "&"}]}], "]"}], ",", + RowBox[{"Joined", "\[Rule]", "True"}], ",", + RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]], "Input", + CellChangeTimes->CompressedData[" +1:eJwdxVkogwEAB/A5krlTJtLYis2WI02KUZMkVxkxtQcWGqYpV59cL2SON0tR +yssaNiXKPvHg2DLTjBqtkCixpmwefA9zfP/v4dePp9LKO0NZLBaXhsuF9tZj +gU9mSxpXYm/tHYEPAk2jeNMh1uErBzWPnUob54TedMRPxoS6vxhvf7OleKzD +dY83iKcHrJ7de8avTuoFawQHXiy1L3xiSdhY6im9ecbCw1X75izMefoWYLZ8 +kcCVkt5RrP/9FZ7R83fdItyToijFZHtDGW6g8mvw2k4V89e5dh1HiTRGXNKY +t4sVg6o9vBzQWbAjfJrElLvPitXcCWY9GenCCS2GaxyRUxFipf/x1DMf6t/T +cVuoPQMXFHUP46Xg8ghutsQ5cJrRw3wz5y7xC32yy7rYMqx5nKLis32yZ9IQ +xLlbYnYC/luNwReNRdGF9PyXe+apFRMPJw5wM7E2IHH66eUt1S78NjR5i7s+ +ZHf4H8z/9Bs= + "], + CellLabel->"In[24]:=",ExpressionUUID->"73943d9b-70f3-4562-ab84-3881106d7e0c"], + +Cell[BoxData[ + GraphicsBox[{{}, {{{}, {}, + {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.009166666666666668], + AbsoluteThickness[1.6], LineBox[CompressedData[" +1:eJw1yX8s1HEcx/ETF83k5yHhHOe4H+gPrl/q+/4uZiyUa/0Y0oysbjL6pd3W +roZOIY21uq06TtouKrKS/Pj6kUlEWBx1lx8RaZOZTItq9/68t9fee+zJS86I +S93AYrGof/v/TTdClbRaC7yVmZTJU9Rg1vHrlpdz0POURB4uOaa5hV6kVG8i +KtxoNXqZesxyKSgylqFXqXfVmTHnhI/Qa1TxXHxCdF8V2gwK6dhDP5Nr0OZQ +Vq/cN1RUh2ZDktz31HDLS7QlOLqGhzR4vEZvgnbNh7Rr+c1oa+i5Q28LW2LQ +NtDaOLdy8mkb2hbo9uXyqNMdJivtYKYpcl21uxNtD7/P5s5eEHWhHeCPQTLt +F9GNdoQqqcWL2aUetBM8VOmlZ/L60Rzg0DH8zsRBtDNM8DM50gMf0S6QU3E+ +sU8+gnYFRa3YXtChR28Bg75rIbphFO0GiqXy9kLFGHor1MasqJ24n9Du0P1A +6pxTS+wBhqwWdv6ez2hPqDdeWZe1EXOhuJFJ2AwGk8ELDs6nWgW/QjNeoGZu +msvFRuw8GDfqZ7bfRzM8CGFEFiqLL9i94erC14aWXWjGG3R1ig10Ouk+sChW +2rHvke4Dzb7OmsRe0vlA9w5pmlZJ54M297uLg3Acuy8sJzlcCj6MZnxBln6E +7laSLoA43f6dsTrSBXDjRBSbNUC6H3i8l7LdV0n3g5WxJ2mlvAns/iDLrnHa +EYlm/CF7I2exPYN0IXRpW6m2UtKFMBqleu5YT7oI8nUjau0Y6SII8wxvMqyR +LoZqZuDXXt4kdjGU3BWo+2k0SODHRHK8OoV0CVSNJnwbzCU9AGJSspLNKkkP +AKvbjdyCDtID4dmkTYB2kvRAuCgdDo03m8IeBJV5FUm2XDQTBMKjb2XToVPU +X9S5NVs= + "]]}, + {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.009166666666666668], + AbsoluteThickness[1.6], LineBox[CompressedData[" +1:eJw1yX8s1HEcx/E7ThzLb+dHPzgO5/w4LX8ws+/bRmGtH7ejov7x41r0Q4XR +Vl0WhkK1hNZ2t0w/Fl2IkR9Hd3RuhUl+XXLnKpviosyPdap5f97ba+899mQn +nxekmdBoNOrf/v+tG6fo5pPC5qx0ast6KvGx4dC4/Dr6O8Uar32jDa9AL1Hi +KpnkGa0avUI9yRBtOAul6HVqW/thw6OROrSRajAYLMskz9F0aL6x4mR58iXa +FLjuitqiy81oM4ADa1Ghila0OSQsaPg29q/RTLhyocCHUdmFtoKNuYzfnfQe +9HZ44CK1og/0om0gpkiiXitRbFlsC9Epqam3BX1oO5hoiko+5/4WbQ/3wnLd +GnUqtAP8yI9sMWtQox0hMzp40arkHdoJVJv5cYOSQTQLQjezhJw7w2hnaLk5 +15M4Mop2gUaFXYcjfRLtCnuYX/7U359Cu0HwMWWhNkOD3gFdZ4fmlhw+oXeC +UWcbr6on3gXiyL44Zdg0ejdUi/PUAjmxO4xqRBHxUZ+3DB6g7ui0iulGyz3A +OLauYrjNYGeDLqTZGCJAy9lQmnv0OKOEdE84YlNe2N9NuifkD49kM5dJ9wJZ +39WP/d5a7F7gnGJ9dyEeDRxYXSr+KS0gnQMsGHI0viLdG6bFeYOpetK9oWIs +27zKWofdB5Zlvx6+D0fLfcC78dtUkoh0X1hripU1lZPuC6y8+banbaRzgXlt +XhqhJZ0LTgf9fFctZrH7QULvAHuKj5b7wdcTrq1VQtJ5kCXcn3Qml3QeSLiJ +iy9qSPeHW5UXeZou0v0h7UNpjlZHegBYGNodTBh67AFwqUbPm+GgIRCKmcy9 +pvtID4QVs/lTp0WkB8Fsf2yJsYD0ILDJdy0qqyOdD8pMw2akknQ+NEzkUOk6 +PfUX0Ek+tw== + "]]}, + {RGBColor[0.560181, 0.691569, 0.194885], PointSize[0.009166666666666668], + AbsoluteThickness[1.6], LineBox[CompressedData[" +1:eJw1yQ0s1HEcx3EPM0+JRLg4d8cd5zhjlZrW/6th5rHSNbPGJuqOPGx1pZKx +47pkhdtCSY2tVBhNHnZbd5eHQiQaIeT5LJlMxES1+/6+22ffvfZmxqWeStDT +0dGh/u3/194XaiMlf7utLpXSeoYak0ufUM0S9CLl2Tg/dCGsEL1C3Widmki0 +e4Beo4K8k9tFehXoTWoipHjysegZepvKEl3RlfJq0Loguns2vFdRj9YH312q +pNdLDWgDsG2hC/MCm9GGID9XwB1oUKCNIfbyoKzlgBJtCox3Iysp0Wq0GdCD +FyzFjFa0OchurnrZL7RpnWUBPuvCTUZzB3oPHHsY5/xU+h5tCYKmyNC1E13o +vVBZTS1W0j6grcDvuubtn/EetDW4Lxey419+RO+DUSo5LehVH9oGVjkRvqFF +n9C2wORH9tSI+9F2cC+qPZ3bP4CmgXxOZptTOozeD7k53aUxJmNoe2i0zj9t +rCF2gNHOECo9exxNB/XvXFq8xQTaEURWkuKEYjQwYOmHe/kRs29aqxggvL9x +cD4EDUzYybhlEyAlnQlb/ApOrJJ0Fuhe2vn+aJ10FkQF1E378SaxO4GA1tlf +H4NWOYHXNcWWUQHpziCYO/rZQE26MxxvelEj/kk6GzJX7TVBjlPY2RBs2jHo +FYYGDmzOVym2rpLOgYzooXxxJekuQMllcdJe0l1ALTkclblOuiuIL1bF8enT +2F1BKEmsO+OPBi5EaLKECUmkc4FeIpgWF5HuBoll4fpFTaS7QZl3WrD5V9J5 +IIpfZi3skM4DfolDNZc1g90daOzZ0Tf+aJU7ONf6dJsISfeA8pHsX3m3SfeA +Q8MnF1afk86H0NrdYNRJOh/ilX315zWke0JwQJGyy3AWuyeUKwLvCDiz1F84 +ATis + "]]}}, { + {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.009166666666666668], + AbsoluteThickness[1.6]}, + {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.009166666666666668], + AbsoluteThickness[1.6]}, + {RGBColor[0.560181, 0.691569, 0.194885], PointSize[0.009166666666666668], + AbsoluteThickness[1.6]}}, { + {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.009166666666666668], + AbsoluteThickness[1.6]}, + {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.009166666666666668], + AbsoluteThickness[1.6]}, + {RGBColor[0.560181, 0.691569, 0.194885], PointSize[0.009166666666666668], + AbsoluteThickness[1.6]}, {}, {}, {}, {}, {}, {}, {}}, { + {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.009166666666666668], + AbsoluteThickness[1.6]}, + {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.009166666666666668], + AbsoluteThickness[1.6]}, + {RGBColor[0.560181, 0.691569, 0.194885], PointSize[0.009166666666666668], + AbsoluteThickness[1.6]}, {}, {}, {}, {}, {}, {}, {}}}, {{}, {}}}, + AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], + Axes->{True, True}, + AxesLabel->{None, None}, + AxesOrigin->{0.052083333333333336`, 0}, + DisplayFunction->Identity, + Frame->{{False, False}, {False, False}}, + FrameLabel->{{None, None}, {None, None}}, + FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}}, + GridLines->{None, None}, + GridLinesStyle->Directive[ + GrayLevel[0.5, 0.4]], + ImageSize->{983., Automatic}, + Method->{ + "OptimizePlotMarkers" -> True, + "CoordinatesToolOptions" -> {"DisplayFunction" -> ({ + (Identity[#]& )[ + Part[#, 1]], + (Identity[#]& )[ + Part[#, 2]]}& ), "CopiedValueFunction" -> ({ + (Identity[#]& )[ + Part[#, 1]], + (Identity[#]& )[ + Part[#, 2]]}& )}}, + PlotRange->{{0.052083333333333336`, 15.875}, {0, 0.6609294624557724}}, + PlotRangeClipping->True, + PlotRangePadding->{{ + Scaled[0.02], + Scaled[0.02]}, { + Scaled[0.02], + Scaled[0.05]}}, + Ticks->{Automatic, Automatic}]], "Output", + CellChangeTimes->{ + 3.778946374557211*^9, 3.7789464676677723`*^9, 3.778946508497999*^9, + 3.778946604369417*^9, {3.7789466690966673`*^9, 3.778946701500992*^9}, { + 3.778946731994258*^9, 3.778946749349575*^9}, 3.7789467962636337`*^9, + 3.778946833931172*^9, 3.778947082771483*^9, {3.778947327915792*^9, + 3.778947345892826*^9}, 3.778948029994212*^9, 3.7789481053151817`*^9, { + 3.778948692166095*^9, 3.7789487038141737`*^9}, 3.778953585064103*^9, + 3.778953624385591*^9, 3.778953682031535*^9, 3.778953746826598*^9, { + 3.77956194067076*^9, 3.779561968768979*^9}, 3.779624858106914*^9, { + 3.779624924084155*^9, 3.779624933643355*^9}, {3.779625037660088*^9, + 3.779625061656746*^9}, 3.779699807199407*^9, {3.779700827690935*^9, + 3.7797008596365232`*^9}, 3.7797009048423357`*^9, 3.779700941258234*^9, { + 3.779701005894755*^9, 3.779701033694045*^9}, 3.78008725076254*^9, + 3.780087374491346*^9, 3.780087426664627*^9, {3.780087476814011*^9, + 3.7800874839642763`*^9}}, + CellLabel->"Out[24]=",ExpressionUUID->"2a29681a-d1fe-4e5f-9729-68c4489a550c"] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"ListPlot", "[", + RowBox[{ + RowBox[{"d", "[", + RowBox[{ + RowBox[{"Select", "[", + RowBox[{ + RowBox[{"#n", "\[Equal]", "20"}], "&"}], "]"}], ",", + RowBox[{ + RowBox[{"With", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"L", "=", + SuperscriptBox["#n", "2"]}], "}"}], ",", + RowBox[{ + RowBox[{ + RowBox[{"{", + RowBox[{"#1", ",", + RowBox[{"#2", "/", "L"}]}], "}"}], "&"}], "@@@", + RowBox[{"Cs", "[", "#", "]"}]}]}], "]"}], "&"}]}], "]"}], ",", + RowBox[{"Joined", "\[Rule]", "True"}], ",", + RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]], "Input", + CellChangeTimes->{{3.779700477940461*^9, 3.77970048742832*^9}, { + 3.780087268580113*^9, 3.7800872686230392`*^9}, {3.780087382286021*^9, + 3.780087398738736*^9}, {3.7800877026590443`*^9, 3.780087726448501*^9}}, + CellLabel->"In[28]:=",ExpressionUUID->"60949ea9-18bb-4b5b-868e-6a12e7d013f6"], + +Cell[BoxData[ + GraphicsBox[{{}, {{{}, {}, + {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.009166666666666668], + AbsoluteThickness[1.6], LineBox[CompressedData[" +1:eJw11HtM01cUB/AKoqwBZBQYrVQohUJ59gFlZo7fUYnoJjiRanAoOhyyhMIS +QAjhoWxmWwfTITBlEswMDt2gzOoEgQoVBLX0QUEQHQpThqNmILKB1Top557k +l/5xcz/3e869KeeTzPhPbWg0Wsqbb/E3XLxYRmrHBrOJ/TqRolnrIVX3zPPk +qLaQOpi6WBPU924fSywuxymtZrFMlPQOa2Pt2Wqq+tRizVBHxnNuv+9zluqz +rj+nitn0vdKsemqZ1Zun5k5XbVYuNFJL55mpJ2OsrdemLlJWPtVCKSvWRuW+ +9Ttl3a6hQaR+/fhxRQuu20DbF3aDkqNtmG857FugMWNzO/B8O7hqu10656pG +fyVEmW7uSl64jnntIX0y6VLSoW706FD93+6xsNwe9BxgfkX9DVlrL/Wj1XOE +87kZ3NJVt9FbBRaTeF3NBg3udwb5ld0Zt+734frb4Js+ZPT5RYueC9ipx1Yb +n+rwfBfgVDozW90NmJcBpo/2CEpPGKg0q+cKrlmFDrNDBvTcYLpt7z/7OvvR +c4db11eVdMYY0XMH12Mqp5F5I3rvQH5t2abtowPoeYAisKZr2+MB9Jgwu98t +49nwEHosiAxzNoZLJtFjgcuOgpvFUjtY8lbDRtP0+pcN0+h5gjlvm0Cd/gd6 +bHBhtDzu3km8NXCeFSPPvDKE978G/m07U/76ZT/O0wvqvrKzGCJIHm+IkSb4 +nmwYxPzesGWuxWO2fRA9DgSaeflp0UacNwfkTeWyzZUGfB8cKHIc8PHar0fP +B4o7sisf9OrR84HZle5/v/DWo8eFuOTBjj8L9Ohx4dEphX1eoR7754LMeCzW +MZx4vjDN/YudeFeHni9wH11+MCzVoecHBXJGmWqZDj0/oD+tqihqJ/ftB6LU +eedxiw49HmiVEtvmFSQfD6ZeZF+eSiL5/CF4rmyhJ4J4/nCRMh4MjSOeP+Qz +M5jL24kXAKzvqIhmPskXABZVQlZGMfH4wP+6oU7yinh8aM1RpuqOEI8PtvZO +VRx70m8gePEaNeGJxAuELrPMYeEQ6TcIytP5kY1FxAsCuuFVMBuIFwSKuAMX +xHuIFwzBuskv451Iv8HQntmlKmWTfCHgbz9SMexP7iMEZtTP6aY8ch8hcMY5 +8aeUKdJvKJhlW7wzfyD5QuFcKUMvDCJeGNAmaPqfGcQLA2/FZ7Uf2OrwPYZB +zFWTo1uydmk/UwCWUV7WvXHd0nvaKoAOz197RHLs97AAHrrHMtR38HylAGZC +5Mk2Wux/QgC/NQ2pMtdhHqYQvrnQJTvRqUVPCO9F8drOleL/wWEhxERvam1a +q0VPCC0BfIGvhnhCUOQ3869JcF5MEbSUvHt/wEA8EeRkpyRdKiGe6M18dp5W +iYkngnrVSO/nN7ToiUCebLznMUjyiSF+jpPs8SHxxBByoFutohNPDEcjxAl3 +a4knhl1p6W7f9vehJwZD9BOlvkZD/Q9wPq4r + "]]}, + {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.009166666666666668], + AbsoluteThickness[1.6], LineBox[CompressedData[" +1:eJw90Q0slHEcB/ATksmSq+Q4584d55yTVpHTnl8X0yytNOta4rALa05UitmN +arI1sZizvLVpLiNmctZVOk4vwnHI2fRmnWRRWNx2i1Xz+/dsz5599t3v+395 +2MmZsfJNNBqN+vv++248k5T41Fx3ed9VasNmKroldNlHWYKep2byXjDyt95D +L1PzH34ZNw80oFcpi5U/kRjVjLZSwrHiqczSdvQ65d120+Kl6kTbACNvii7T +a9G2EBcQrCpP6UbbQ7HTktCrpAftAOEVMVKOVo92BLu5K2OJn16ineDhj+ns +lfuv0c6QfmF0RNrVj94GCzWBSV6NAxsucIEhU+F3g2kIvR2SNGUJjLphtCuY +juRkyUtH0HR4owo3xSuM6B1wVyyZrPIcRe+Exo7aB8fixtC7oP/4YlhICrEb +9Gmesc6UjaN3g7ktOXei9x3aHSw1P/v5Z6fRDJC9CgmO59rChj3g81pT6Ee1 +BXNPsPhffFv4aBLNhEiJm/PRVrI/L0ioVNS5OxCzgDnLnW6twfOAN4hiIzLE +DWidNxjcyrSKVJKzITdGWVmdg/ehY8NprUF2uQUNHLgjr4r66kjmOcCnqtl0 +iQFzH7CyFsfjltA6H1DervwycInkXLC6zCbSaKSfC+G5+9eK0kg/D+qb3g9H +KMg8D4aUCwdi95HcF1qzn/PoPfj/dL6wnmrfZZ4l/X6wfEgvEUeTfj8YLKiQ +bkkjOR9spIZvqSfJPB+KdUyu8AYa/EFf23y41kjW94fsazyTJJOsLwAoigxz +zSK5AIwl9JwnBaQ/ADzSBeeu3yJ5AOyZ2buy2kH6hdCZKhA9/p8LQczyszuo +If2B8FQ9qPatJ3kgMNoTmgs4pF8E1YrC9SgZ2b8I5Bn5J3q4JA8CVe9vRaea +zAeB+bwzUxNjoP4Ab0s0ug== + "]]}}, { + {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.009166666666666668], + AbsoluteThickness[1.6]}, + {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.009166666666666668], + AbsoluteThickness[1.6]}}, { + {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.009166666666666668], + AbsoluteThickness[1.6]}, + {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.009166666666666668], + AbsoluteThickness[1.6]}, {}, {}, {}}, { + {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.009166666666666668], + AbsoluteThickness[1.6]}, + {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.009166666666666668], + AbsoluteThickness[1.6]}, {}, {}, {}}}, {{}, {}}}, + AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], + Axes->{True, True}, + AxesLabel->{None, None}, + AxesOrigin->{0, 0}, + DisplayFunction->Identity, + Frame->{{False, False}, {False, False}}, + FrameLabel->{{None, None}, {None, None}}, + FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}}, + GridLines->{None, None}, + GridLinesStyle->Directive[ + GrayLevel[0.5, 0.4]], + Method->{ + "OptimizePlotMarkers" -> True, + "CoordinatesToolOptions" -> {"DisplayFunction" -> ({ + (Identity[#]& )[ + Part[#, 1]], + (Identity[#]& )[ + Part[#, 2]]}& ), "CopiedValueFunction" -> ({ + (Identity[#]& )[ + Part[#, 1]], + (Identity[#]& )[ + Part[#, 2]]}& )}}, + PlotRange->{{0, 19.9}, {0, 2.666747640711386}}, + PlotRangeClipping->True, + PlotRangePadding->{{ + Scaled[0.02], + Scaled[0.02]}, { + Scaled[0.02], + Scaled[0.05]}}, + Ticks->{Automatic, Automatic}]], "Output", + CellChangeTimes->{ + 3.7797004888467417`*^9, {3.780087263365243*^9, 3.780087269048139*^9}, { + 3.780087378630178*^9, 3.780087399174458*^9}, {3.780087705097999*^9, + 3.780087726696293*^9}}, + CellLabel->"Out[28]=",ExpressionUUID->"1c27332b-7aa0-4f6c-9d7f-d912fe119acd"] +}, Open ]], + +Cell[BoxData[ + RowBox[{ + RowBox[{"a", "=", + RowBox[{ + "Import", "[", + "\"\<~/doc/research/hadamard/code/energies_24_0.200000_20.000000_100.dat\>\ +\"", "]"}]}], ";"}]], "Input", + CellChangeTimes->{{3.779699855009756*^9, 3.779699866112481*^9}}, + CellLabel->"In[8]:=",ExpressionUUID->"bb64ee8f-2ab2-44da-86d3-071e144df4ac"], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"ListPlot", "[", "a", "]"}]], "Input", + CellChangeTimes->{{3.7796998674408073`*^9, 3.779699868848329*^9}}, + CellLabel->"In[9]:=",ExpressionUUID->"27200d6e-a8c7-4ded-bf66-0ba360a8d0c7"], + +Cell[BoxData[ + GraphicsBox[{{}, {{ + {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.009166666666666668], + AbsoluteThickness[1.6], PointBox[{{1., 1.*^7}, {1., 1.*^7}}]}, + {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.009166666666666668], + AbsoluteThickness[1.6], PointBox[CompressedData[" +1:eJw1y38w13cAx/FvrUkWKm7JaX5m5Ed+fP1qystX+NIX3y9fvnw1oqVrd8vO +OhX+4MZNxbo1scba6vyaHcPXTCzlaCXbMGFHJu3K3He5r/zKj7Tt+r4+d+/7 +3OP9fL8sk1Mjj64XCATv/Xf+/7/6NH7LF6c8/bNmOrQXyLPOzrOvpnWRteZe +vdpHb0FPxz1JywJtjIyKzvhM02dam8Bo9ESxzJc2w6GiqfsD8bQFLjZfM/n8 +DG0F/+ynf5ddom2wDgOJdvW0LcSNu2rRTdvhXHHITp0JejfSriSPFi3Sjrg+ +8qZaT39Wa2cMRo0PJ1nSLrD3FfVeF9KuSIrw7rIMpt1gcNc1qyaWdkdR2Bej +8cdpIZRiqYX4DO0BYanM9lQ+7Yl01QfdM8W0F1Tdwbd/Lqe9EXgzrGOtgfZB +f07V45J2ei8EzY2/Vd6j30G/2UqN1xDti/LWo0uBE/Q+vLyprhtU0/uR8yDZ +ZmWe9oOr0+Hlupc0IO/TMdbdNPfK2UDq93vSNVu1FvjjxU7zG4Wm7P5IsGrV +77JiF8Fju6Spdje7CGPJ5Qv73NgDkHIio6rQhz0AJQ0t1YVgP4C1HoUkIpj9 +AAy8TBTDYeyBcLyiFLnK2QMxNtOqc0PJHgQ9Pc2RqiPsQfjk0ccO8gz2YHxU +XzGe+w17MAIK7psKa9nFONvn4y1qYRfjsm7c9P529hCkmST13OlkD0HFbNcx +g7vsoXBwNvcY6WEPhUmTuPRQL/tBxIy9tun47+wHEffHG5WSQXYJCttz654P +s0vg7OM2lj/CHgbVZO/A4gP2MCz/MlQmHmcPhzD8VMfZCfZwfFnfmFn9F3sE +4nozXOsfs0fgV72Sa5cn2aX4c4Pw8MkprSHFjoW8HR5qvpei6XxJqfofrW9J +YWPkNH5hmnsZ3t7Yfdtaw70Mcw+XhitnuJchMdFo3GqWexkePTUwuDDHfSRW +f5r2npnnPhI1q51Pgha5j8SnitO9Rc+5j4TRpe7coSXuozDwoeHKlhXuo/BZ +ivvV8FXuo2C7N7u84AX3UXDRKNf617iX47TI/Lsgwbx2L4eNsi7CYL3W2XJM +nj9XnLpB61tyqFzc27bpai2IRqt7p1ypz3003r96rPmkMffR2BzTZvaVKffR +SE9osP3anPsYbN+cUBJozX0M7iSFfavZxX0M6oZKml634z4GoVWqgR/tuVfg +h4K0t/QduFcgty1z0c2RewWeqcvG7Zy4V8BmY9I6Q2fuYzFnq9j6hEYsIJhP +aNnDfSwe5kgN8124j4VFyjbVu67cxyHZzzNL4jbf8S8dmScE + "]]}}, { + {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.009166666666666668], + AbsoluteThickness[1.6]}, + {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.009166666666666668], + AbsoluteThickness[1.6]}, {}, {}, {}}, { + {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.009166666666666668], + AbsoluteThickness[1.6]}, + {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.009166666666666668], + AbsoluteThickness[1.6]}, {}, {}, {}}}, {{}, {}}}, + AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], + Axes->{True, True}, + AxesLabel->{None, None}, + AxesOrigin->{0., -4.400895677366526*^9}, + DisplayFunction->Identity, + Frame->{{False, False}, {False, False}}, + FrameLabel->{{None, None}, {None, None}}, + FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}}, + GridLines->{None, None}, + GridLinesStyle->Directive[ + GrayLevel[0.5, 0.4]], + Method->{ + "OptimizePlotMarkers" -> True, + "CoordinatesToolOptions" -> {"DisplayFunction" -> ({ + (Identity[#]& )[ + Part[#, 1]], + (Identity[#]& )[ + Part[#, 2]]}& ), "CopiedValueFunction" -> ({ + (Identity[#]& )[ + Part[#, 1]], + (Identity[#]& )[ + Part[#, 2]]}& )}}, + PlotRange->{{0., 100.}, {-5.689886659390474*^9, -4.468737307999365*^9}}, + PlotRangeClipping->True, + PlotRangePadding->{{ + Scaled[0.02], + Scaled[0.02]}, { + Scaled[0.05], + Scaled[0.05]}}, + Ticks->{Automatic, Automatic}]], "Output", + CellChangeTimes->{3.7796998691208344`*^9}, + CellLabel->"Out[9]=",ExpressionUUID->"161e209e-ab11-4d64-99a1-d2cc4b5ffda5"] +}, 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"] +}, Open ]] +}, +WindowSize->{1277, 1076}, +WindowMargins->{{Automatic, 2}, {Automatic, 2}}, +FrontEndVersion->"12.0 for Linux x86 (64-bit) (April 8, 2019)", +StyleDefinitions->"Default.nb" +] +(* End of Notebook Content *) + +(* Internal cache information *) +(*CellTagsOutline +CellTagsIndex->{} +*) +(*CellTagsIndex +CellTagsIndex->{} +*) +(*NotebookFileOutline +Notebook[{ +Cell[558, 20, 1782, 46, 124, "Input",ExpressionUUID->"6ca216e4-4a7b-4e5f-a718-f4e70c77d6a2"], +Cell[2343, 68, 337, 7, 31, "Input",ExpressionUUID->"3b7ef50e-b362-4da3-b734-818bb17a7c2f"], +Cell[2683, 77, 668, 18, 59, "Input",ExpressionUUID->"f756d73a-4921-45d7-acd1-44fe21438a02"], +Cell[3354, 97, 722, 18, 31, "Input",ExpressionUUID->"c91e85b2-f90e-40d7-bf19-7a03fd9499bb"], +Cell[4079, 117, 804, 24, 31, "Input",ExpressionUUID->"1b068f0a-fa61-487d-893c-5b00f279e758"], +Cell[4886, 143, 1753, 48, 55, "Input",ExpressionUUID->"83b164e2-d1bf-4be5-bf92-d3645cb950ca"], +Cell[6642, 193, 1968, 53, 83, "Input",ExpressionUUID->"cfeba3ce-dc48-40c8-b8b8-f34cf59a253d"], +Cell[8613, 248, 712, 20, 31, "Input",ExpressionUUID->"a1830ef0-911d-442d-b5d2-b62c202fce89"], +Cell[CellGroupData[{ +Cell[9350, 272, 1317, 33, 39, "Input",ExpressionUUID->"73943d9b-70f3-4562-ab84-3881106d7e0c"], +Cell[10670, 307, 6446, 123, 630, "Output",ExpressionUUID->"2a29681a-d1fe-4e5f-9729-68c4489a550c"] +}, Open ]], +Cell[CellGroupData[{ +Cell[17153, 435, 979, 25, 39, "Input",ExpressionUUID->"60949ea9-18bb-4b5b-868e-6a12e7d013f6"], +Cell[18135, 462, 4725, 95, 244, "Output",ExpressionUUID->"1c27332b-7aa0-4f6c-9d7f-d912fe119acd"] +}, Open ]], +Cell[22875, 560, 330, 8, 31, "Input",ExpressionUUID->"bb64ee8f-2ab2-44da-86d3-071e144df4ac"], +Cell[CellGroupData[{ +Cell[23230, 572, 207, 3, 31, "Input",ExpressionUUID->"27200d6e-a8c7-4ded-bf66-0ba360a8d0c7"], +Cell[23440, 577, 3276, 68, 214, "Output",ExpressionUUID->"161e209e-ab11-4d64-99a1-d2cc4b5ffda5"] +}, Open ]], +Cell[CellGroupData[{ +Cell[26753, 650, 222, 4, 31, "Input",ExpressionUUID->"85d1832b-b409-4be7-b8d2-26df59eff8fb"], +Cell[26978, 656, 177, 2, 35, "Output",ExpressionUUID->"5dbe4cc3-7c89-4c96-870e-625a9bf27876"] +}, Open ]] +} +] +*) + -- cgit v1.2.3-54-g00ecf