plots.nb

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

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

(* CreatedBy='Mathematica 9.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[     16420,        464]
NotebookOptionsPosition[     14514,        396]
NotebookOutlinePosition[     14857,        411]
CellTagsIndexPosition[     14814,        408]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Some basic plotting commands", "Section",
 CellChangeTimes->{{3.676379915049796*^9, 3.676379934912633*^9}, 
   3.6763842780818925`*^9, {3.678910954140424*^9, 3.67891095550366*^9}}],

Cell["\<\
Below are some quick code snippets to provide an idea of how to plot various \
quantities.\
\>", "Text",
 CellChangeTimes->{{3.6763842798468175`*^9, 3.676384308693612*^9}}],

Cell[CellGroupData[{

Cell["Find directories with data", "Subsection",
 CellChangeTimes->{{3.6763810374788733`*^9, 3.6763810552071257`*^9}}],

Cell["Work in the notebook directory", "Text",
 CellChangeTimes->{{3.676379906882577*^9, 3.676379941160733*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"SetDirectory", "[", 
   RowBox[{"NotebookDirectory", "[", "]"}], "]"}], ";"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.6109704409332256`*^9, 3.61097050007586*^9}, {
   3.627921033258557*^9, 3.627921048130514*^9}, 3.6279215246919727`*^9, {
   3.676379886482505*^9, 3.6763799037163377`*^9}, {3.676379944346958*^9, 
   3.6763799453136005`*^9}}],

Cell["\<\
Get a list of subdirectories with log.txt files, indicative of run data\
\>", "Text",
 CellChangeTimes->{{3.6763799463940277`*^9, 3.676379966325293*^9}, {
  3.676380000303673*^9, 3.676380021802626*^9}, {3.676380176053421*^9, 
  3.6763801770143204`*^9}, {3.6763803911193895`*^9, 3.6763804031049733`*^9}}],

Cell[BoxData[
 RowBox[{"directories", " ", "=", " ", 
  RowBox[{"DirectoryName", "/@", 
   RowBox[{"FileNames", "[", 
    RowBox[{"\"\<log.txt\>\"", ",", "\"\<\>\"", ",", "\[Infinity]"}], 
    "]"}]}]}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.67637996887175*^9, 3.676379997592099*^9}, {
  3.6763803845767183`*^9, 3.6763803879999356`*^9}, {3.6763810489849515`*^9, 
  3.676381050457323*^9}, {3.678903891624094*^9, 3.678903892765128*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Plot Reference FRW values", "Subsection",
 CellChangeTimes->{{3.676381081924012*^9, 3.676381090903887*^9}, {
  3.676381409867964*^9, 3.676381412692502*^9}}],

Cell["\<\
Output at https://github.com/jbcm627/cosmograph/blob/master/IO/io.cc#L193
If this shows up blank, the reference metric may not have been used.\
\>", "Text",
 CellChangeTimes->{{3.6763842498342404`*^9, 3.676384263817173*^9}, {
  3.6763844330598946`*^9, 3.6763844341803136`*^9}, {3.678910904025009*^9, 
  3.6789109224286876`*^9}}],

Cell[BoxData[{
 RowBox[{"frwfiles", " ", "=", " ", 
  RowBox[{"FileNames", "[", 
   RowBox[{"\"\<*.frwdat.*\>\"", ",", "directories", ",", "\[Infinity]"}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"frwdata", "=", 
   RowBox[{"Import", "/@", "frwfiles"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.676381108037866*^9, 3.676381166726838*^9}, {
  3.6763812705077963`*^9, 3.67638128609394*^9}, {3.6763840621515303`*^9, 
  3.676384062337332*^9}}],

Cell[BoxData[
 RowBox[{"ListPlot", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"frwdata", "\[LeftDoubleBracket]", 
    RowBox[{";;", ",", ";;", ",", 
     RowBox[{"{", 
      RowBox[{"1", ",", "2"}], "}"}]}], "\[RightDoubleBracket]"}], ",", 
   "\[IndentingNewLine]", 
   RowBox[{"Joined", "\[Rule]", "True"}], ",", 
   RowBox[{"Axes", "\[Rule]", "False"}], ",", 
   RowBox[{"Frame", "\[Rule]", "True"}], ",", "\[IndentingNewLine]", 
   RowBox[{"FrameLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
     "\"\<\!\(\*SubscriptBox[\(\[Phi]\), \(FRW\)]\)\>\"", ",", 
      "\"\<\!\(\*SubscriptBox[\(K\), \(FRW\)]\)\>\""}], "}"}]}], ",", 
   RowBox[{"PlotLegends", "\[Rule]", "frwfiles"}]}], "\[IndentingNewLine]", 
  "]"}]], "Input",
 CellChangeTimes->{{3.6763813103306446`*^9, 3.6763814229815764`*^9}, {
  3.6789108895318727`*^9, 3.6789108923805113`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Plot metric statistics", "Subsection",
 CellChangeTimes->{{3.6763840305524282`*^9, 3.6763840397323556`*^9}}],

Cell["\<\
Output at https://github.com/jbcm627/cosmograph/blob/master/IO/io.cc#L217\
\>", "Text",
 CellChangeTimes->{{3.6763842092922983`*^9, 3.6763842465045705`*^9}, 
   3.6763844388375034`*^9}],

Cell[BoxData[{
 RowBox[{"statfiles", " ", "=", " ", 
  RowBox[{"FileNames", "[", 
   RowBox[{
   "\"\<calculated.dat.gz\>\"", ",", "directories", ",", "\[Infinity]"}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"statdata", " ", "=", " ", 
   RowBox[{"Import", "/@", "statfiles"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.676384125616172*^9, 3.6763841539874916`*^9}}],

Cell[BoxData[
 RowBox[{"ListPlot", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"Range", "[", 
             RowBox[{"Length", "[", "#", "]"}], "]"}], "-", "1"}], ")"}], "/", 
          RowBox[{"Length", "[", "#", "]"}]}], ",", "#"}], "}"}], 
       "\[Transpose]"}], ")"}], "&"}], "/@", 
    RowBox[{"statdata", "\[LeftDoubleBracket]", 
     RowBox[{";;", ",", ";;", ",", "4"}], "\[RightDoubleBracket]"}]}], ",", 
   "\[IndentingNewLine]", 
   RowBox[{"Joined", "\[Rule]", "True"}], ",", 
   RowBox[{"Axes", "\[Rule]", "False"}], ",", 
   RowBox[{"Frame", "\[Rule]", "True"}], ",", "\[IndentingNewLine]", 
   RowBox[{"FrameLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
     "\"\<Simulation Step (% of total)\>\"", ",", 
      "\"\<\!\(\*SubscriptBox[\(\[Sigma]\), \(K\)]\)\>\""}], "}"}]}], ",", 
   RowBox[{"PlotLegends", "\[Rule]", "statfiles"}]}], "\[IndentingNewLine]", 
  "]"}]], "Input",
 CellChangeTimes->{{3.676384042545755*^9, 3.6763840425697737`*^9}, {
   3.676384161587376*^9, 3.6763841653335476`*^9}, {3.676384318373906*^9, 
   3.6763843344486113`*^9}, {3.6763843734063406`*^9, 3.676384418859645*^9}, 
   3.6789073190738974`*^9, {3.678907478907933*^9, 3.67890748219191*^9}, {
   3.6789097396165266`*^9, 3.6789097430025473`*^9}, {3.6789098737318993`*^9, 
   3.678909894955335*^9}, {3.6789108743150873`*^9, 3.6789108749596677`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Plot Constraint Violation", "Subsection",
 CellChangeTimes->{{3.6763814467527256`*^9, 3.676381450898116*^9}, 
   3.6763818124226866`*^9}],

Cell["\<\
Output at https://github.com/jbcm627/cosmograph/blob/master/IO/io.cc#L154\
\>", "Text",
 CellChangeTimes->{{3.6763841928401213`*^9, 3.676384201546498*^9}, 
   3.6763844455345354`*^9}],

Cell[BoxData[{
 RowBox[{"Hviolfiles", " ", "=", " ", 
  RowBox[{"FileNames", "[", 
   RowBox[{
   "\"\<H_violations.dat.gz\>\"", ",", "directories", ",", "\[Infinity]"}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Hviol", "=", 
   RowBox[{"Import", "/@", "Hviolfiles"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Mviolfiles", " ", "=", " ", 
  RowBox[{"FileNames", "[", 
   RowBox[{
   "\"\<M_violations.dat.gz\>\"", ",", "directories", ",", "\[Infinity]"}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Mviol", " ", "=", " ", 
   RowBox[{"Import", "/@", "Mviolfiles"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.676380663691797*^9, 3.6763806706422873`*^9}, {
  3.676381825482644*^9, 3.676381848690467*^9}, {3.67638262556114*^9, 
  3.676382643058236*^9}, {3.676384065489628*^9, 3.676384094571913*^9}}],

Cell[CellGroupData[{

Cell["Hamiltonian Constraint Violation", "Subsubsection",
 CellChangeTimes->{{3.6789107804486055`*^9, 3.6789107878863907`*^9}}],

Cell[BoxData[
 RowBox[{"ListPlot", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"Range", "[", 
             RowBox[{"Length", "[", "#", "]"}], "]"}], "-", "1"}], ")"}], "/", 
          RowBox[{"Length", "[", "#", "]"}]}], ",", "#"}], "}"}], 
       "\[Transpose]"}], ")"}], "&"}], "/@", 
    RowBox[{"Hviol", "\[LeftDoubleBracket]", 
     RowBox[{";;", ",", ";;", ",", "4"}], "\[RightDoubleBracket]"}]}], ",", 
   "\[IndentingNewLine]", 
   RowBox[{"Joined", "\[Rule]", "True"}], ",", 
   RowBox[{"Axes", "\[Rule]", "False"}], ",", 
   RowBox[{"Frame", "\[Rule]", "True"}], ",", "\[IndentingNewLine]", 
   RowBox[{"FrameLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<Simulation output step #\>\"", ",", "\"\<Max[H]\>\""}], 
     "}"}]}], ",", 
   RowBox[{"PlotLegends", "\[Rule]", "Hviolfiles"}]}], "\[IndentingNewLine]", 
  "]"}]], "Input",
 CellChangeTimes->{{3.6763825668337636`*^9, 3.6763826051340213`*^9}, {
   3.676383244021113*^9, 3.676383244172325*^9}, 3.6763843402998533`*^9, {
   3.67890396162869*^9, 3.6789040014638176`*^9}, {3.6789040529178777`*^9, 
   3.6789040802386346`*^9}, {3.6789041601736937`*^9, 
   3.6789041709095716`*^9}, {3.6789042327110972`*^9, 3.678904257518583*^9}, {
   3.6789073841103325`*^9, 3.678907433927961*^9}, {3.678907525432096*^9, 
   3.678907581840746*^9}, 3.678908014124953*^9, {3.678909763438777*^9, 
   3.678909763667985*^9}, {3.678910851463373*^9, 3.678910864749474*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Momentum Constraint Violation", "Subsubsection",
 CellChangeTimes->{{3.6789107804486055`*^9, 3.678910796241962*^9}}],

Cell[BoxData[
 RowBox[{"ListPlot", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"Range", "[", 
             RowBox[{"Length", "[", "#", "]"}], "]"}], "-", "1"}], ")"}], "/", 
          RowBox[{"Length", "[", "#", "]"}]}], ",", "#"}], "}"}], 
       "\[Transpose]"}], ")"}], "&"}], "/@", 
    RowBox[{"Mviol", "\[LeftDoubleBracket]", 
     RowBox[{";;", ",", ";;", ",", "4"}], "\[RightDoubleBracket]"}]}], ",", 
   "\[IndentingNewLine]", 
   RowBox[{"Joined", "\[Rule]", "True"}], ",", 
   RowBox[{"Axes", "\[Rule]", "False"}], ",", 
   RowBox[{"Frame", "\[Rule]", "True"}], ",", "\[IndentingNewLine]", 
   RowBox[{"FrameLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<Simulation output step #\>\"", ",", "\"\<Max[M]\>\""}], 
     "}"}]}], ",", 
   RowBox[{"PlotLegends", "\[Rule]", "Mviolfiles"}]}], "\[IndentingNewLine]", 
  "]"}]], "Input",
 CellChangeTimes->{{3.6763831286381717`*^9, 3.6763831289820766`*^9}, {
   3.6763832518387747`*^9, 3.676383255446208*^9}, 3.676384348346648*^9, {
   3.6789040211156263`*^9, 3.6789040243515606`*^9}, {3.678904216282689*^9, 
   3.6789042165439253`*^9}, {3.6789042499407096`*^9, 3.678904253600025*^9}, {
   3.678907355533436*^9, 3.678907368896544*^9}, {3.678907443023204*^9, 
   3.6789074434946327`*^9}, 3.6789076295802374`*^9, {3.6789108543830185`*^9, 
   3.678910870091254*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Plot 1-d, 2-d, 3-d initial data from \[CapitalDelta]\[Phi] field", \
"Subsection",
 CellChangeTimes->{{3.6763844518579907`*^9, 3.6763844773057413`*^9}, {
  3.6763847404258757`*^9, 3.6763847551883264`*^9}}],

Cell[CellGroupData[{

Cell["3 - D", "Subsubsection",
 CellChangeTimes->{{3.6763846570523834`*^9, 3.676384660706005*^9}}],

Cell[BoxData[{
 RowBox[{"\[CapitalDelta]\[Phi]files3D", "=", 
  RowBox[{"FileNames", "[", 
   RowBox[{
   "\"\<3D_DIFFphi.0.3d_grid.h5.gz\>\"", ",", "directories", ",", 
    "\[Infinity]"}], "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[CapitalDelta]\[Phi]3D", "=", 
   RowBox[{
    RowBox[{
     RowBox[{"Import", "[", 
      RowBox[{"#", ",", 
       RowBox[{"{", 
        RowBox[{"\"\<Datasets\>\"", ",", "\"\<DS1\>\""}], "}"}]}], "]"}], 
     "&"}], "/@", "\[CapitalDelta]\[Phi]files3D"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.676384545778428*^9, 3.6763845857854123`*^9}}],

Cell[BoxData[
 RowBox[{"GraphicsGrid", "[", 
  RowBox[{"{", "\[IndentingNewLine]", 
   RowBox[{"ListDensityPlot3D", "/@", "\[CapitalDelta]\[Phi]3D"}], 
   "\[IndentingNewLine]", "}"}], "]"}]], "Input",
 CellChangeTimes->{{3.6763846079839797`*^9, 3.676384647294714*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["2 - D", "Subsubsection",
 CellChangeTimes->{{3.676384665752746*^9, 3.676384667700614*^9}}],

Cell[BoxData[{
 RowBox[{"\[CapitalDelta]\[Phi]files2D", "=", 
  RowBox[{"FileNames", "[", 
   RowBox[{
   "\"\<2D_DIFFphi.0.2d_grid.h5.gz\>\"", ",", "directories", ",", 
    "\[Infinity]"}], "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[CapitalDelta]\[Phi]2D", "=", 
   RowBox[{
    RowBox[{
     RowBox[{"Import", "[", 
      RowBox[{"#", ",", 
       RowBox[{"{", 
        RowBox[{"\"\<Datasets\>\"", ",", "\"\<Dataset1\>\""}], "}"}]}], "]"}],
      "&"}], "/@", "\[CapitalDelta]\[Phi]files2D"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.676384700261*^9, 3.6763847168230896`*^9}, {
  3.676384808282786*^9, 3.6763848214613295`*^9}}],

Cell[BoxData[
 RowBox[{"GraphicsGrid", "[", 
  RowBox[{"{", "\[IndentingNewLine]", 
   RowBox[{"ListDensityPlot", "/@", "\[CapitalDelta]\[Phi]2D"}], 
   "\[IndentingNewLine]", "}"}], "]"}]], "Input",
 CellChangeTimes->{{3.6763847244644933`*^9, 3.67638472579235*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["1 - D", "Subsubsection",
 CellChangeTimes->{{3.676384665752746*^9, 3.6763846754262295`*^9}}],

Cell[BoxData[{
 RowBox[{"\[CapitalDelta]\[Phi]files1D", "=", 
  RowBox[{"FileNames", "[", 
   RowBox[{
   "\"\<1D_DIFFphi.strip.dat.gz\>\"", ",", "directories", ",", 
    "\[Infinity]"}], "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[CapitalDelta]\[Phi]1D", "=", 
   RowBox[{
    RowBox[{
     RowBox[{"Import", "[", "#", "]"}], "&"}], "/@", 
    "\[CapitalDelta]\[Phi]files1D"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.6763847759505186`*^9, 3.67638483907016*^9}}],

Cell[BoxData[
 RowBox[{"GraphicsGrid", "[", 
  RowBox[{
   RowBox[{"{", "\[IndentingNewLine]", 
    RowBox[{"ListPlot", "/@", 
     RowBox[{"\[CapitalDelta]\[Phi]1D", "\[LeftDoubleBracket]", 
      RowBox[{";;", ",", "1", ",", ";;"}], "\[RightDoubleBracket]"}]}], 
    "\[IndentingNewLine]", "}"}], ",", 
   RowBox[{"ImageSize", "\[Rule]", "500"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.676384843968258*^9, 3.676384901805817*^9}}]
}, Open  ]]
}, Open  ]]
}, Open  ]]
},
WindowSize->{1600, 827},
WindowMargins->{{-8, Automatic}, {Automatic, -8}},
FrontEndVersion->"10.2 for Microsoft Windows (64-bit) (July 29, 2015)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 186, 2, 70, "Section"],
Cell[768, 26, 182, 4, 30, "Text"],
Cell[CellGroupData[{
Cell[975, 34, 118, 1, 49, "Subsection"],
Cell[1096, 37, 112, 1, 30, "Text"],
Cell[1211, 40, 395, 8, 31, "Input",
 InitializationCell->True],
Cell[1609, 50, 313, 5, 30, "Text"],
Cell[1925, 57, 455, 9, 31, "Input",
 InitializationCell->True]
}, Open  ]],
Cell[CellGroupData[{
Cell[2417, 71, 162, 2, 49, "Subsection"],
Cell[2582, 75, 338, 6, 49, "Text"],
Cell[2923, 83, 457, 10, 52, "Input"],
Cell[3383, 95, 863, 19, 112, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[4283, 119, 114, 1, 49, "Subsection"],
Cell[4400, 122, 195, 4, 30, "Text"],
Cell[4598, 128, 382, 9, 52, "Input"],
Cell[4983, 139, 1521, 34, 112, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[6541, 178, 143, 2, 49, "Subsection"],
Cell[6687, 182, 193, 4, 30, "Text"],
Cell[6883, 188, 830, 19, 92, "Input"],
Cell[CellGroupData[{
Cell[7738, 211, 127, 1, 39, "Subsubsection"],
Cell[7868, 214, 1595, 35, 112, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[9500, 254, 122, 1, 39, "Subsubsection"],
Cell[9625, 257, 1500, 34, 112, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[11174, 297, 211, 3, 49, "Subsection"],
Cell[CellGroupData[{
Cell[11410, 304, 98, 1, 39, "Subsubsection"],
Cell[11511, 307, 591, 15, 52, "Input"],
Cell[12105, 324, 269, 5, 72, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[12411, 334, 96, 1, 39, "Subsubsection"],
Cell[12510, 337, 644, 16, 52, "Input"],
Cell[13157, 355, 266, 5, 72, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[13460, 365, 98, 1, 39, "Subsubsection"],
Cell[13561, 368, 477, 12, 52, "Input"],
Cell[14041, 382, 433, 9, 72, "Input"]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}
]
*)