saddle-lib-test.nb

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

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

(* CreatedBy='Mathematica 11.2' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     37875,        820]
NotebookOptionsPosition[     35064,        765]
NotebookOutlinePosition[     35404,        780]
CellTagsIndexPosition[     35361,        777]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"SetDirectory", "[", 
   RowBox[{"NotebookDirectory", "[", "]"}], "]"}], ";"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.7673569687731447`*^9, 
  3.7673569687738857`*^9}},ExpressionUUID->"7e8350ee-f09c-4d78-ac3b-\
7f39904190c4"],

Cell["\<\
Disable current notebook FileOutlineCache and CellChangeTimes.  See: \
https://stackoverflow.com/questions/2816628/version-control-of-mathematica-\
notebooks \
\>", "Text",ExpressionUUID->"caed0e72-f004-4254-b57a-ccde69ef9676"],

Cell[CellGroupData[{

Cell["Library", "Title",ExpressionUUID->"77d08537-f7f3-4451-808d-7153544399d8"],

Cell[BoxData[
 RowBox[{
  RowBox[{"<<", "\"\<mma-lib/saddle-point.m\>\""}], ";"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.782562775018255*^9, 
  3.78256277677363*^9}},ExpressionUUID->"7484096c-cb41-4df2-b847-\
1c41437d49ea"],

Cell[BoxData[
 RowBox[{
  RowBox[{"removeTriangle", "[", "mat_", "]"}], ":=", 
  RowBox[{"Array", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{"#1", "\[GreaterEqual]", "#2"}], ",", 
       RowBox[{"mat", "[", 
        RowBox[{"[", "##", "]"}], "]"}], ",", "0"}], "]"}], "&"}], ",", 
    RowBox[{"Dimensions", "[", "mat", "]"}]}], "]"}]}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.782656605472694*^9, 3.782656611787201*^9}, {
  3.782656830385612*^9, 
  3.7826569154717007`*^9}},ExpressionUUID->"5789a40d-5130-461c-b144-\
b55de59e0464"],

Cell[CellGroupData[{

Cell["Create random saddle point at zero", "Section",
 CellChangeTimes->{{3.782561716415669*^9, 3.782561751325232*^9}, {
  3.782583609295702*^9, 
  3.782583627213561*^9}},ExpressionUUID->"588fbb2d-502a-4553-92c3-\
4c3a45b5df12"],

Cell["\<\
Independent test of the saddle-lib program using an explicitly constructed \
Hessian.  We use blockSize=1, since blockSize has been well-vetted.  Also, \
the hessian matrix is dense.\
\>", "Text",
 CellChangeTimes->{{3.782560899417079*^9, 3.7825610135091877`*^9}, {
  3.7825836769597282`*^9, 
  3.782583715750764*^9}},ExpressionUUID->"b49dc30d-f685-4917-888e-\
8925caca7789"],

Cell[BoxData[
 RowBox[{
  RowBox[{"makeZeroSaddle", "[", 
   RowBox[{"constraints_", ",", "largeShifts_", ",", "nTotal_"}], "]"}], ":=", 
  RowBox[{"Block", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"n1", "=", 
       RowBox[{"nTotal", "-", "largeShifts", "-", "constraints"}]}], ",", 
      RowBox[{"large", "=", "4"}]}], "}"}], ",", 
    RowBox[{
     RowBox[{"hess", "=", 
      RowBox[{"DiagonalMatrix", "[", 
       RowBox[{"Join", "[", 
        RowBox[{
         RowBox[{"Table", "[", 
          RowBox[{"0", ",", 
           RowBox[{"{", "constraints", "}"}]}], "]"}], ",", 
         RowBox[{"RandomReal", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"-", "1"}], ",", "1"}], "}"}], "/", "large"}], ",", 
           "largeShifts"}], "]"}], ",", 
         RowBox[{"RandomReal", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{"-", "1"}], ",", "1"}], "}"}], " ", ",", 
           RowBox[{"nTotal", "-", "constraints", "-", "largeShifts"}]}], 
          "]"}]}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"gauge", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{"i", "\[Equal]", "j"}], ",", "1", ",", "0"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "constraints"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"j", ",", "nTotal"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"x0", " ", "=", " ", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"i", ">", "constraints"}], "&&", 
           RowBox[{"i", "<=", 
            RowBox[{"constraints", "+", "largeShifts"}]}]}], ",", "large", 
          ",", "1"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "nTotal"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"grad", "=", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"x0", ".", "hess"}], "+", 
         RowBox[{"hess", ".", "x0"}]}], ")"}], "/", "2"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"nc", "=", 
      RowBox[{"Sqrt", "[", "2", "]"}]}], ";", 
     RowBox[{"latticeDimensions", "=", 
      RowBox[{"{", 
       RowBox[{"GCD", "[", 
        RowBox[{"constraints", ",", "nTotal"}], "]"}], "}"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"oo", "=", 
      RowBox[{"RandomVariate", "[", 
       RowBox[{"CircularRealMatrixDistribution", "[", "nTotal", "]"}], 
       "]"}]}], ";"}]}], "]"}]}]], "Input",
 CellChangeTimes->CompressedData["
1:eJwdzTlIwnEAxXERKQdTrCEwsTQiTHAwITLQGjo3KyFrKjuWDkiipqDBkEAb
cigJDKzAyArPkogIizBMyII0osMhJFGyUiiSfu8/PD7Lg69waKp7hE6j0fhk
UC+12xzSVMuDPG2HvixvB/bqo/vwjOVwQ+u42gOXQrNHkNV1HICeau0pDBU0
13BXc3UDGw8bXmDKm/iA3KjtE/b5z/9gjmks2iHeBsUc2BzOUMq4cR4UfFgq
IN3cXAWzKk4dfGUPNMHtxwUVzOeEanhQaO+B7MC7GZ7EVMuQdRHxQtPTsA+u
/pQp9oiVX1uUJ/lNnosYMRglUGZKdEBdbZpSyWwLu4n+1rcIFPwq3qGlJpOC
Uq5hxYOf3r0GJ0wMtR+d7xJKRfKccUdc7A8WQxF/uvyeGFgX8eD84IYYrvTM
SeBl0GmNEbX1SUqly6mLE93yzjFYemqbhM+jjBn4DwZi5rQ=
  
  "],ExpressionUUID->"8128455f-421f-4879-bdbd-fb8f6064fdbe"],

Cell[BoxData[
 RowBox[{"makeZeroSaddle", "[", 
  RowBox[{"12", ",", "4", ",", "60"}], "]"}]], "Input",
 CellChangeTimes->{{3.782562149584385*^9, 3.7825621689429417`*^9}, {
  3.7825705660447607`*^9, 3.782570571604999*^9}, {3.782583509032957*^9, 
  3.7825835142447844`*^9}, {3.782589110365363*^9, 3.7825891167148933`*^9}, {
  3.782589264405539*^9, 3.782589292943904*^9}, {3.782649817034506*^9, 
  3.782649823058633*^9}, {3.7826498806649837`*^9, 3.782649882956093*^9}, {
  3.782649941637237*^9, 3.782649945364552*^9}, {3.7826508673561697`*^9, 
  3.782650867596459*^9}, {3.7826563611718817`*^9, 3.782656362242104*^9}, {
  3.782656513024809*^9, 3.782656518228669*^9}, {3.782657064661723*^9, 
  3.782657068284443*^9}, {3.7826591322753487`*^9, 3.782659132737564*^9}, {
  3.78266030684516*^9, 3.782660307315263*^9}, {3.782661140419877*^9, 
  3.7826611408795357`*^9}},ExpressionUUID->"77587275-3b49-4031-a9b4-\
9c149aeaac7b"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Block", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"hess", "=", "hess"}], ",", 
     RowBox[{"grad", "=", "grad"}], ",", 
     RowBox[{"gauge", "=", "gauge"}]}], "}"}], ",", 
   RowBox[{"Chop", "[", 
    RowBox[{"findDelta", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"hess", ",", "grad", ",", "gauge"}], "}"}], ",", 
      RowBox[{"largeShiftCutoff", "\[Rule]", "2"}], ",", 
      RowBox[{"Method", "\[Rule]", "\"\<Dense\>\""}]}], "]"}], "]"}]}], 
  "]"}]], "Input",
 CellChangeTimes->{
  3.782562435513649*^9, {3.782562813524234*^9, 3.7825628354367933`*^9}, {
   3.782562909287883*^9, 3.782562910970261*^9}, {3.782563461629524*^9, 
   3.7825634651587067`*^9}, {3.7825704759391127`*^9, 
   3.7825705463209257`*^9}, {3.7825721665905724`*^9, 
   3.7825721676784353`*^9}, {3.782578548774455*^9, 3.782578553016066*^9}, {
   3.782581986219469*^9, 3.782582008100095*^9}, {3.7825825380329113`*^9, 
   3.782582569371564*^9}, {3.78258260902708*^9, 3.7825826345987864`*^9}, {
   3.7825827030536737`*^9, 3.782582718608828*^9}, {3.7825827585838337`*^9, 
   3.782582869875939*^9}, {3.782584755701001*^9, 3.7825847590952578`*^9}, {
   3.7825849558778267`*^9, 3.782584982346477*^9}, {3.782589074326005*^9, 
   3.782589075888097*^9}, {3.7826549083024473`*^9, 
   3.782654923258662*^9}},ExpressionUUID->"c8b6ea22-a1ba-4c97-aac5-\
8b2874ed978b"],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"applyCutoff3:  \"\>", "\[InvisibleSpace]", "4", 
   "\[InvisibleSpace]", "\<\" of \"\>", "\[InvisibleSpace]", "48", 
   "\[InvisibleSpace]", "\<\" zeros, between \"\>", "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{"36", ",", "46"}], "}"}]}],
  SequenceForm["applyCutoff3:  ", 4, " of ", 48, " zeros, between ", {36, 46}],
  Editable->False]], "Print",
 CellChangeTimes->{{3.7825851350658503`*^9, 3.7825851450076923`*^9}, {
   3.7825890764887857`*^9, 3.782589119715064*^9}, {3.7825892683322983`*^9, 
   3.782589297298403*^9}, 3.782649827132422*^9, 3.782649949468762*^9, 
   3.782649982765266*^9, {3.7826508832860613`*^9, 3.782650897302432*^9}, {
   3.782654902053808*^9, 3.782654924775929*^9}, 3.782655303553157*^9, {
   3.782656218589714*^9, 3.782656245846547*^9}, 3.782656365471498*^9, 
   3.7826565220969276`*^9, 3.782657021827548*^9, 3.782657073583559*^9, 
   3.7826582392108727`*^9, 3.782659137045724*^9, 3.782659295561084*^9, 
   3.7826593668561993`*^9, 3.782660314727646*^9, 3.782660932718391*^9, {
   3.782661124254121*^9, 
   3.7826611456520853`*^9}},ExpressionUUID->"142be610-21f4-445d-b0a6-\
ba30920fb274"],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"shifts norms and rescale:  \"\>", "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{"1.`", ",", "6.6332495807108`", ",", "1.`"}], "}"}]}],
  SequenceForm["shifts norms and rescale:  ", {1., 6.6332495807108, 1.}],
  Editable->False]], "Print",
 CellChangeTimes->{{3.7825851350658503`*^9, 3.7825851450076923`*^9}, {
   3.7825890764887857`*^9, 3.782589119715064*^9}, {3.7825892683322983`*^9, 
   3.782589297298403*^9}, 3.782649827132422*^9, 3.782649949468762*^9, 
   3.782649982765266*^9, {3.7826508832860613`*^9, 3.782650897302432*^9}, {
   3.782654902053808*^9, 3.782654924775929*^9}, 3.782655303553157*^9, {
   3.782656218589714*^9, 3.782656245846547*^9}, 3.782656365471498*^9, 
   3.7826565220969276`*^9, 3.782657021827548*^9, 3.782657073583559*^9, 
   3.7826582392108727`*^9, 3.782659137045724*^9, 3.782659295561084*^9, 
   3.7826593668561993`*^9, 3.782660314727646*^9, 3.782660932718391*^9, {
   3.782661124254121*^9, 
   3.782661145653695*^9}},ExpressionUUID->"8b1bca75-7853-476d-9cb8-\
9f14d120e0ba"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", 
   ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
   "0", ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "0.9999999999999999`"}], ",", 
   RowBox[{"-", "0.9999999999999999`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "0.9999999999999999`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "0.9999999999999999`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.`"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.782585116679306*^9, 3.782585145010666*^9}, {
   3.782589076491521*^9, 3.7825891197182093`*^9}, {3.782589268334833*^9, 
   3.782589297300885*^9}, 3.782649827140768*^9, 3.782649949471655*^9, 
   3.782649982771448*^9, {3.782650883290009*^9, 3.782650897305552*^9}, {
   3.7826549020555983`*^9, 3.7826549247793283`*^9}, 3.78265530355704*^9, {
   3.78265621859307*^9, 3.782656245849204*^9}, 3.782656365475668*^9, 
   3.782656522100181*^9, 3.782657021830366*^9, 3.782657073587343*^9, 
   3.782658239214233*^9, 3.782659137049383*^9, 3.782659295564555*^9, 
   3.782659366859021*^9, 3.782660314730866*^9, 3.782660932722849*^9, {
   3.7826611242582912`*^9, 
   3.782661145655909*^9}},ExpressionUUID->"ec83d663-d7e5-4054-9491-\
3da493055aa8"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Block", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"hess", "=", 
      RowBox[{"oo", ".", "hess", ".", 
       RowBox[{"Transpose", "[", "oo", "]"}]}]}], ",", 
     RowBox[{"grad", "=", 
      RowBox[{"oo", ".", "grad"}]}], ",", 
     RowBox[{"gauge", "=", 
      RowBox[{"SparseArray", "[", 
       RowBox[{"gauge", ".", 
        RowBox[{"Transpose", "[", "oo", "]"}]}], "]"}]}]}], "}"}], ",", 
   RowBox[{
    RowBox[{"hess", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"hess", "+", 
        RowBox[{"Transpose", "[", "hess", "]"}]}], ")"}], "/", "2"}]}], ";", 
    RowBox[{"Chop", "[", 
     RowBox[{
      RowBox[{"findDelta", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"hess", ",", "grad", ",", "gauge"}], "}"}], ",", 
        RowBox[{"largeShiftCutoff", "\[Rule]", "1"}], ",", 
        RowBox[{"rescaleCutoff", "\[Rule]", "2"}], ",", 
        RowBox[{"Method", "\[Rule]", "\"\<Dense\>\""}]}], "]"}], ".", "oo"}], 
     "]"}]}]}], "]"}]], "Input",
 CellChangeTimes->{
  3.782562435513649*^9, {3.782562813524234*^9, 3.7825628354367933`*^9}, {
   3.782562909287883*^9, 3.782562910970261*^9}, {3.782563461629524*^9, 
   3.7825634651587067`*^9}, {3.7825704759391127`*^9, 
   3.7825705463209257`*^9}, {3.7825721665905724`*^9, 
   3.7825721676784353`*^9}, {3.782578548774455*^9, 3.782578553016066*^9}, {
   3.782581986219469*^9, 3.782582008100095*^9}, {3.7825825380329113`*^9, 
   3.782582569371564*^9}, {3.78258260902708*^9, 3.7825826345987864`*^9}, {
   3.7825827030536737`*^9, 3.782582718608828*^9}, {3.7825827585838337`*^9, 
   3.782582869875939*^9}, {3.782584755701001*^9, 3.7825847590952578`*^9}, {
   3.7825850319203167`*^9, 3.7825850914003067`*^9}, {3.7825851798539143`*^9, 
   3.782585180033781*^9}, {3.782589085956164*^9, 3.782589087554303*^9}, {
   3.782650011940165*^9, 3.7826500434328327`*^9}, {3.78265076361721*^9, 
   3.7826507690566683`*^9}, {3.782650827695301*^9, 3.782650841951812*^9}, {
   3.782650917844472*^9, 3.782650977633123*^9}, {3.782654831001894*^9, 
   3.782654854259512*^9}, {3.7826549344939537`*^9, 3.782654935929555*^9}, {
   3.7826553145441*^9, 3.782655392396159*^9}, {3.7826554303379707`*^9, 
   3.782655445683961*^9}, {3.782655673150926*^9, 3.782655675358139*^9}, {
   3.782655716663412*^9, 3.7826557692651978`*^9}, {3.78265582772617*^9, 
   3.7826558335995197`*^9}, {3.782655960697959*^9, 3.782655965437128*^9}, {
   3.782656445857448*^9, 3.782656453279531*^9}, {3.7826564895231133`*^9, 
   3.7826564958916693`*^9}, {3.782656559463686*^9, 3.782656598674707*^9}, {
   3.782656713483946*^9, 3.782656733169656*^9}, {3.782656948066584*^9, 
   3.782656967966096*^9}, {3.782657146392611*^9, 3.782657148154387*^9}, {
   3.782657358071356*^9, 3.782657364805691*^9}, {3.782657397188362*^9, 
   3.78265742834508*^9}, {3.782657458582782*^9, 3.782657482025961*^9}, {
   3.7826575208110323`*^9, 3.78265755162647*^9}, {3.782657607146769*^9, 
   3.782657621423006*^9}, 3.782657665061451*^9, {3.7826598756548433`*^9, 
   3.7826598768646307`*^9}, {3.782661168323925*^9, 3.782661199028166*^9}, {
   3.7826613298636017`*^9, 3.782661330444346*^9}, {3.7826614154853487`*^9, 
   3.782661444314023*^9}},ExpressionUUID->"7eee7402-8768-40f3-bd8d-\
427eaaa99c77"],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"applyCutoff3:  \"\>", "\[InvisibleSpace]", "4", 
   "\[InvisibleSpace]", "\<\" of \"\>", "\[InvisibleSpace]", "48", 
   "\[InvisibleSpace]", "\<\" zeros, between \"\>", "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{"36", ",", "46"}], "}"}]}],
  SequenceForm["applyCutoff3:  ", 4, " of ", 48, " zeros, between ", {36, 46}],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.78265762290689*^9, 3.782657667594302*^9, 3.7826582490942993`*^9, 
   3.782659140934462*^9, 3.7826592978835382`*^9, 3.78265937470713*^9, {
   3.782659872104104*^9, 3.7826598774892683`*^9}, 3.7826603172228727`*^9, 
   3.782660937997993*^9, {3.782661129559134*^9, 3.782661200119714*^9}, {
   3.782661313859033*^9, 3.782661360862439*^9}, {3.782661417512732*^9, 
   3.782661444778261*^9}},ExpressionUUID->"9aeb4762-51ad-42ed-914c-\
50f5bd103df7"],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"shifts norms and rescale:  \"\>", "\[InvisibleSpace]", 
   RowBox[{"{", 
    RowBox[{
    "1.813615607700426`", ",", "6.63324958071084`", ",", 
     "1.0000000000000702`"}], "}"}]}],
  SequenceForm[
  "shifts norms and rescale:  ", {1.813615607700426, 6.63324958071084, 
   1.0000000000000702`}],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.78265762290689*^9, 3.782657667594302*^9, 3.7826582490942993`*^9, 
   3.782659140934462*^9, 3.7826592978835382`*^9, 3.78265937470713*^9, {
   3.782659872104104*^9, 3.7826598774892683`*^9}, 3.7826603172228727`*^9, 
   3.782660937997993*^9, {3.782661129559134*^9, 3.782661200119714*^9}, {
   3.782661313859033*^9, 3.782661360862439*^9}, {3.782661417512732*^9, 
   3.782661444780324*^9}},ExpressionUUID->"d9afd4f1-92ac-47d5-a513-\
85db1bbdd34f"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", 
   ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
   "0", ",", 
   RowBox[{"-", "1.0000000000000264`"}], ",", 
   RowBox[{"-", "0.9999999999999981`"}], ",", 
   RowBox[{"-", "1.0000000000000449`"}], ",", 
   RowBox[{"-", "1.0000000000000195`"}], ",", 
   RowBox[{"-", "1.000000000000015`"}], ",", 
   RowBox[{"-", "0.9999999999999872`"}], ",", 
   RowBox[{"-", "1.0000000000000249`"}], ",", 
   RowBox[{"-", "0.9999999999999889`"}], ",", 
   RowBox[{"-", "1.0000000000000002`"}], ",", 
   RowBox[{"-", "1.0000000000000024`"}], ",", 
   RowBox[{"-", "0.999999999999999`"}], ",", 
   RowBox[{"-", "1.0000000000000313`"}], ",", 
   RowBox[{"-", "0.9999999999999855`"}], ",", 
   RowBox[{"-", "0.9999999999999871`"}], ",", 
   RowBox[{"-", "1.000000000000027`"}], ",", 
   RowBox[{"-", "1.0000000000000049`"}], ",", 
   RowBox[{"-", "1.0000000000000093`"}], ",", 
   RowBox[{"-", "1.0000000000000129`"}], ",", 
   RowBox[{"-", "1.0000000000000022`"}], ",", 
   RowBox[{"-", "0.9999999999999932`"}], ",", 
   RowBox[{"-", "0.9999999999999908`"}], ",", 
   RowBox[{"-", "1.0000000000000102`"}], ",", 
   RowBox[{"-", "1.0000000000000042`"}], ",", 
   RowBox[{"-", "1.0000000000003477`"}], ",", 
   RowBox[{"-", "1.0000000000000084`"}], ",", 
   RowBox[{"-", "1.0000000000000215`"}], ",", 
   RowBox[{"-", "1.0000000000000124`"}], ",", 
   RowBox[{"-", "0.9999999999999982`"}], ",", 
   RowBox[{"-", "1.0000000000000773`"}], ",", 
   RowBox[{"-", "1.000000000000003`"}], ",", 
   RowBox[{"-", "0.9999999999999899`"}], ",", 
   RowBox[{"-", "1.0000000000000195`"}], ",", 
   RowBox[{"-", "0.9999999999999852`"}], ",", 
   RowBox[{"-", "1.0000000000000004`"}], ",", 
   RowBox[{"-", "0.9999999999998483`"}], ",", 
   RowBox[{"-", "0.9999999999999917`"}], ",", 
   RowBox[{"-", "1.0000000000000044`"}], ",", 
   RowBox[{"-", "1.000000000000011`"}], ",", 
   RowBox[{"-", "0.9999999999999855`"}], ",", 
   RowBox[{"-", "0.9999999999999594`"}], ",", 
   RowBox[{"-", "0.9999999999999853`"}], ",", 
   RowBox[{"-", "0.9999999999999981`"}], ",", 
   RowBox[{"-", "0.999999999999999`"}], ",", 
   RowBox[{"-", "0.9999999999998511`"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.7826575247351027`*^9, 3.782657552452633*^9}, 
   3.782657622909482*^9, 3.782657667597128*^9, 3.7826582490976257`*^9, 
   3.782659140937252*^9, 3.782659297886698*^9, 3.782659374709334*^9, {
   3.782659872108581*^9, 3.782659877493594*^9}, 3.7826603172303753`*^9, 
   3.7826609380021133`*^9, {3.782661129563051*^9, 3.782661200123313*^9}, {
   3.782661313862755*^9, 3.7826613608672323`*^9}, {3.782661417516595*^9, 
   3.782661444782874*^9}},ExpressionUUID->"1b55b683-3da5-4077-be59-\
86bc3955046a"]
}, Open  ]],

Cell[BoxData["StringQ"], "Input",
 CellChangeTimes->{{3.7826626914965353`*^9, 
  3.782662698503046*^9}},ExpressionUUID->"818e24ff-dc61-4d09-a373-\
e51326e36275"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Block", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"hess", "=", 
      RowBox[{"SparseArray", "[", 
       RowBox[{"removeTriangle", "[", 
        RowBox[{"oo", ".", "hess", ".", 
         RowBox[{"Transpose", "[", "oo", "]"}]}], "]"}], "]"}]}], ",", 
     RowBox[{"grad", "=", 
      RowBox[{"oo", ".", "grad"}]}], ",", 
     RowBox[{"gauge", "=", 
      RowBox[{"SparseArray", "[", 
       RowBox[{"gauge", ".", 
        RowBox[{"Transpose", "[", "oo", "]"}]}], "]"}]}]}], "}"}], ",", 
   RowBox[{"Chop", "[", 
    RowBox[{
     RowBox[{"findDelta", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"hess", ",", "grad", ",", "gauge"}], "}"}], ",", 
       RowBox[{"largeShiftOptions", "\[Rule]", 
        RowBox[{"{", 
         RowBox[{"eigenPairs", "\[Rule]", 
          RowBox[{"-", "20"}]}], "}"}]}], ",", 
       RowBox[{"largeShiftCutoff", "\[Rule]", "1"}], ",", 
       RowBox[{"rescaleCutoff", "\[Rule]", "2"}], ",", 
       RowBox[{"Method", "\[Rule]", "\"\<External\>\""}], ",", 
       RowBox[{"externalAction", "\[Rule]", "\"\<remote\>\""}], ",", 
       RowBox[{"processes", "\[Rule]", "2"}]}], "]"}], ".", "oo"}], "]"}]}], 
  "]"}]], "Input",
 CellChangeTimes->{
  3.782562435513649*^9, {3.782562813524234*^9, 3.7825628354367933`*^9}, {
   3.782562909287883*^9, 3.782562910970261*^9}, {3.782563461629524*^9, 
   3.7825634651587067`*^9}, {3.7825704759391127`*^9, 
   3.7825705463209257`*^9}, {3.7825721665905724`*^9, 
   3.7825721676784353`*^9}, {3.782578548774455*^9, 3.782578553016066*^9}, {
   3.782581986219469*^9, 3.782582008100095*^9}, {3.7825825380329113`*^9, 
   3.782582569371564*^9}, {3.78258260902708*^9, 3.7825826345987864`*^9}, {
   3.7825827030536737`*^9, 3.782582718608828*^9}, {3.7825827585838337`*^9, 
   3.782582869875939*^9}, {3.782584755701001*^9, 3.7825847590952578`*^9}, {
   3.7825850319203167`*^9, 3.7825850914003067`*^9}, {3.7825851798539143`*^9, 
   3.782585180033781*^9}, {3.782589085956164*^9, 3.782589087554303*^9}, {
   3.782650011940165*^9, 3.7826500434328327`*^9}, {3.78265076361721*^9, 
   3.7826507690566683`*^9}, {3.782650827695301*^9, 3.782650841951812*^9}, {
   3.782650917844472*^9, 3.782650977633123*^9}, {3.782654831001894*^9, 
   3.782654854259512*^9}, {3.7826549344939537`*^9, 3.782654935929555*^9}, {
   3.7826553145441*^9, 3.782655392396159*^9}, {3.7826554303379707`*^9, 
   3.782655445683961*^9}, {3.782655673150926*^9, 3.782655675358139*^9}, {
   3.782655716663412*^9, 3.7826557692651978`*^9}, {3.78265582772617*^9, 
   3.7826558335995197`*^9}, {3.782655960697959*^9, 3.782655965437128*^9}, {
   3.782656445857448*^9, 3.782656453279531*^9}, {3.7826564895231133`*^9, 
   3.7826564958916693`*^9}, {3.782656559463686*^9, 3.782656598674707*^9}, {
   3.782656713483946*^9, 3.782656733169656*^9}, {3.782656948066584*^9, 
   3.782656967966096*^9}, {3.782657146392611*^9, 3.782657148154387*^9}, {
   3.782657358071356*^9, 3.782657364805691*^9}, {3.782657397188362*^9, 
   3.78265742834508*^9}, {3.782657458582782*^9, 3.782657482025961*^9}, {
   3.782657922457035*^9, 3.782657954956361*^9}, {3.782659093281973*^9, 
   3.7826591127463503`*^9}, {3.782659180246702*^9, 3.782659232624591*^9}, {
   3.782659323532857*^9, 3.782659337721921*^9}, {3.782659391827957*^9, 
   3.782659440874511*^9}, {3.782659486272327*^9, 3.782659651967956*^9}, {
   3.782659692819592*^9, 3.7826597644620953`*^9}, {3.782659834464961*^9, 
   3.7826598403092203`*^9}, {3.782659884980027*^9, 3.782659960634584*^9}, {
   3.782660051477244*^9, 3.782660068384879*^9}, {3.782660123625054*^9, 
   3.782660125299101*^9}, {3.782660229562516*^9, 3.782660269933401*^9}, {
   3.7826604522757883`*^9, 3.782660472406783*^9}, {3.782660737551478*^9, 
   3.782660739292865*^9}, {3.782661207840393*^9, 3.782661237492841*^9}, {
   3.782661489346595*^9, 3.7826614896918*^9}, {3.7826616150306787`*^9, 
   3.782661635683399*^9}, {3.782661739201438*^9, 3.782661740793182*^9}, {
   3.7826618055586243`*^9, 
   3.782661810353057*^9}},ExpressionUUID->"cf2e1a15-8016-42d8-bcde-\
aa0802c3be56"],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Running \"\>", "\[InvisibleSpace]", "\<\"pshifts\"\>", 
   "\[InvisibleSpace]", "\<\" at \"\>", 
   "\[InvisibleSpace]", \
"\<\"samson:lattice/gauge-field-analyze/saddle-lib\"\>"}],
  SequenceForm[
  "Running ", "pshifts", " at ", 
   "samson:lattice/gauge-field-analyze/saddle-lib"],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.7826560440283613`*^9, 3.7826562554537477`*^9, 3.782656375649253*^9, 
   3.782656454286672*^9, {3.782656496678426*^9, 3.782656525608779*^9}, {
   3.782656567504281*^9, 3.782656578992988*^9}, 3.7826567355771923`*^9, 
   3.782656969604823*^9, 3.782657028663678*^9, 3.782657076452814*^9, 
   3.7826571491955557`*^9, 3.782657366127198*^9, {3.782657398692021*^9, 
   3.782657429916443*^9}, {3.782657472182067*^9, 3.782657482863453*^9}, 
   3.782657671650832*^9, 3.782657793845943*^9, {3.782657925486659*^9, 
   3.7826579560394163`*^9}, 3.782658260218503*^9, {3.782659095478405*^9, 
   3.782659113762329*^9}, {3.7826591747489443`*^9, 3.782659233488926*^9}, {
   3.782659318482155*^9, 3.782659338648018*^9}, {3.782659382842613*^9, 
   3.782659441993248*^9}, {3.7826594887197943`*^9, 3.7826596526932087`*^9}, {
   3.7826597186616583`*^9, 3.782659765455974*^9}, 3.782659841270282*^9, 
   3.782659886815648*^9, {3.782659938044525*^9, 3.7826599616450977`*^9}, {
   3.782660057531056*^9, 3.782660069388386*^9}, 3.782660126102231*^9, 
   3.782660239484991*^9, 3.7826602707747793`*^9, 3.78266033801071*^9, {
   3.782660454032401*^9, 3.782660473792946*^9}, 3.7826607405971813`*^9, 
   3.782660950963111*^9, {3.7826612200048847`*^9, 3.782661238544505*^9}, {
   3.782661478123623*^9, 3.7826614906276197`*^9}, 3.782661636905671*^9, 
   3.7826617413842583`*^9, 
   3.7826618111308813`*^9},ExpressionUUID->"ef1c6de6-db8b-4780-bfb7-\
65aff542d63d"],

Cell[BoxData[
 StyleBox["\<\"Opening file \
lattice/gauge-field-analyze/saddle-lib/hess-grad-gauge.json n=60, gauge=12 \
blockSize=1 partitions=12\\nSetting OpenMP & OpenBLAS to 2 threads\\nOpening \
file lattice/gauge-field-analyze/saddle-lib/hess.mtx with 1830 nonzero \
elements.\\nOpening file lattice/gauge-field-analyze/saddle-lib/grad.dat for \
a vector of length 60\\nOpening file \
lattice/gauge-field-analyze/saddle-lib/gauge.mtx with 720 nonzero \
elements.\\nOpening file lattice/gauge-field-analyze/saddle-lib/gauge.mtx \
with 720 nonzero elements.\\ncutoffNullSpace:  4 of 20 zeros between [2, \
12]\\nOpening output file \
lattice/gauge-field-analyze/saddle-lib/shifts0.dat\\nMatrix hessian ops=152, \
total process localTime=0.02 mpiTime=0.00\\n    global 0.1 GFLOP/s, \
stream=1.1 GB/s\\nMatrix gauge ops=1619, total process localTime=0.07 \
mpiTime=0.00\\n    global 0.1 GFLOP/s, stream=0.7 GB/s\\nMatrix gaugeT \
ops=1619, total process localTime=0.07 mpiTime=0.01\\n    global 0.1 GFLOP/s, \
stream=0.6 GB/s\\nshifts.c:  input/output in 0.07 sec (0 wall)\\nshifts.c:  \
overall time 0.30 sec (0 wall)\\n\"\>",
  StripOnInput->False,
  FontColor->RGBColor[0, 0, 1]]], "Print",
 CellChangeTimes->{
  3.7826560440283613`*^9, 3.7826562554537477`*^9, 3.782656375649253*^9, 
   3.782656454286672*^9, {3.782656496678426*^9, 3.782656525608779*^9}, {
   3.782656567504281*^9, 3.782656578992988*^9}, 3.7826567355771923`*^9, 
   3.782656969604823*^9, 3.782657028663678*^9, 3.782657076452814*^9, 
   3.7826571491955557`*^9, 3.782657366127198*^9, {3.782657398692021*^9, 
   3.782657429916443*^9}, {3.782657472182067*^9, 3.782657482863453*^9}, 
   3.782657671650832*^9, 3.782657793845943*^9, {3.782657925486659*^9, 
   3.7826579560394163`*^9}, 3.782658260218503*^9, {3.782659095478405*^9, 
   3.782659113762329*^9}, {3.7826591747489443`*^9, 3.782659233488926*^9}, {
   3.782659318482155*^9, 3.782659338648018*^9}, {3.782659382842613*^9, 
   3.782659441993248*^9}, {3.7826594887197943`*^9, 3.7826596526932087`*^9}, {
   3.7826597186616583`*^9, 3.782659765455974*^9}, 3.782659841270282*^9, 
   3.782659886815648*^9, {3.782659938044525*^9, 3.7826599616450977`*^9}, {
   3.782660057531056*^9, 3.782660069388386*^9}, 3.782660126102231*^9, 
   3.782660239484991*^9, 3.7826602707747793`*^9, 3.78266033801071*^9, {
   3.782660454032401*^9, 3.782660473792946*^9}, 3.7826607405971813`*^9, 
   3.782660950963111*^9, {3.7826612200048847`*^9, 3.782661238544505*^9}, {
   3.782661478123623*^9, 3.7826614906276197`*^9}, 3.782661636905671*^9, 
   3.7826617413842583`*^9, 
   3.782661812708956*^9},ExpressionUUID->"17e75830-6a9d-4e3c-9518-\
949ce59ac784"],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"findDelta shifts time (seconds):  total=\"\>", 
   "\[InvisibleSpace]", "1.869543`6.4222504563288965"}],
  SequenceForm[
  "findDelta shifts time (seconds):  total=", 1.869543`6.4222504563288965],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.7826560440283613`*^9, 3.7826562554537477`*^9, 3.782656375649253*^9, 
   3.782656454286672*^9, {3.782656496678426*^9, 3.782656525608779*^9}, {
   3.782656567504281*^9, 3.782656578992988*^9}, 3.7826567355771923`*^9, 
   3.782656969604823*^9, 3.782657028663678*^9, 3.782657076452814*^9, 
   3.7826571491955557`*^9, 3.782657366127198*^9, {3.782657398692021*^9, 
   3.782657429916443*^9}, {3.782657472182067*^9, 3.782657482863453*^9}, 
   3.782657671650832*^9, 3.782657793845943*^9, {3.782657925486659*^9, 
   3.7826579560394163`*^9}, 3.782658260218503*^9, {3.782659095478405*^9, 
   3.782659113762329*^9}, {3.7826591747489443`*^9, 3.782659233488926*^9}, {
   3.782659318482155*^9, 3.782659338648018*^9}, {3.782659382842613*^9, 
   3.782659441993248*^9}, {3.7826594887197943`*^9, 3.7826596526932087`*^9}, {
   3.7826597186616583`*^9, 3.782659765455974*^9}, 3.782659841270282*^9, 
   3.782659886815648*^9, {3.782659938044525*^9, 3.7826599616450977`*^9}, {
   3.782660057531056*^9, 3.782660069388386*^9}, 3.782660126102231*^9, 
   3.782660239484991*^9, 3.7826602707747793`*^9, 3.78266033801071*^9, {
   3.782660454032401*^9, 3.782660473792946*^9}, 3.7826607405971813`*^9, 
   3.782660950963111*^9, {3.7826612200048847`*^9, 3.782661238544505*^9}, {
   3.782661478123623*^9, 3.7826614906276197`*^9}, 3.782661636905671*^9, 
   3.7826617413842583`*^9, 
   3.7826618129785843`*^9},ExpressionUUID->"240ea087-643f-43f6-ae84-\
65167f3ae102"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", 
   ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
   "0", ",", 
   RowBox[{"-", "0.9999999999999999`"}], ",", 
   RowBox[{"-", "0.9999999999999525`"}], ",", 
   RowBox[{"-", "0.9999999999999997`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.000000000000064`"}], ",", 
   RowBox[{"-", "0.9999999999999986`"}], ",", 
   RowBox[{"-", "1.0000000000000484`"}], ",", 
   RowBox[{"-", "0.9999999999999993`"}], ",", 
   RowBox[{"-", "0.999999999999999`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "0.9999999999997277`"}], ",", 
   RowBox[{"-", "1.0000000000000475`"}], ",", 
   RowBox[{"-", "0.9999999999999704`"}], ",", 
   RowBox[{"-", "0.9999999999999885`"}], ",", 
   RowBox[{"-", "1.0000000000000067`"}], ",", 
   RowBox[{"-", "0.9999999999999998`"}], ",", 
   RowBox[{"-", "0.9999999999999988`"}], ",", 
   RowBox[{"-", "0.9999999999999994`"}], ",", 
   RowBox[{"-", "0.9999999999999999`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "0.9999999999999899`"}], ",", 
   RowBox[{"-", "0.9999999999999997`"}], ",", 
   RowBox[{"-", "0.9999999999999989`"}], ",", 
   RowBox[{"-", "1.0000000000009017`"}], ",", 
   RowBox[{"-", "0.9999999999999992`"}], ",", 
   RowBox[{"-", "0.9999999999999992`"}], ",", 
   RowBox[{"-", "0.9999999999999999`"}], ",", 
   RowBox[{"-", "1.`"}], ",", 
   RowBox[{"-", "1.0000000000000004`"}], ",", 
   RowBox[{"-", "0.9999999999999998`"}], ",", 
   RowBox[{"-", "0.9999999999999974`"}], ",", 
   RowBox[{"-", "1.0000000000000042`"}], ",", 
   RowBox[{"-", "1.0000000000000016`"}], ",", 
   RowBox[{"-", "0.9999999999999994`"}], ",", 
   RowBox[{"-", "0.9999999999999994`"}], ",", 
   RowBox[{"-", "1.0000000000000004`"}], ",", 
   RowBox[{"-", "0.9999999999999989`"}], ",", 
   RowBox[{"-", "1.0000000000000004`"}], ",", 
   RowBox[{"-", "0.9999999999999992`"}], ",", 
   RowBox[{"-", "1.0000000000000067`"}], ",", 
   RowBox[{"-", "0.9999999999999891`"}], ",", 
   RowBox[{"-", "0.9999999999999991`"}], ",", 
   RowBox[{"-", "1.0000000000000004`"}], ",", 
   RowBox[{"-", "0.9999999999999996`"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.782656044039515*^9, 3.782656255464892*^9, 3.7826563756602173`*^9, 
   3.7826564542981853`*^9, {3.78265649669221*^9, 3.7826565256252527`*^9}, {
   3.782656567517809*^9, 3.782656579007235*^9}, 3.782656735595517*^9, 
   3.782656969622855*^9, 3.782657028680867*^9, 3.782657076457417*^9, 
   3.782657149199078*^9, 3.7826573661308002`*^9, {3.782657398695157*^9, 
   3.782657429920418*^9}, {3.782657472196245*^9, 3.782657482867399*^9}, 
   3.782657671670463*^9, 3.7826577938491488`*^9, {3.7826579255013638`*^9, 
   3.782657956053795*^9}, 3.782658260233121*^9, {3.782659095493382*^9, 
   3.78265911376683*^9}, {3.7826591747518682`*^9, 3.78265923349337*^9}, {
   3.7826593184966497`*^9, 3.7826593386646967`*^9}, {3.782659382846299*^9, 
   3.782659442009965*^9}, {3.782659488736351*^9, 3.78265965269774*^9}, {
   3.78265971866533*^9, 3.78265976546028*^9}, 3.7826598412739067`*^9, 
   3.7826598868195763`*^9, {3.782659938048023*^9, 3.7826599616493607`*^9}, {
   3.7826600575356903`*^9, 3.782660069392701*^9}, 3.782660126105638*^9, 
   3.7826602394888372`*^9, 3.78266027077871*^9, 3.782660338015094*^9, {
   3.782660454036304*^9, 3.782660473797304*^9}, 3.782660740601551*^9, 
   3.782660950966507*^9, {3.782661220009088*^9, 3.782661238547531*^9}, {
   3.782661478126342*^9, 3.782661490631761*^9}, 3.782661663150096*^9, 
   3.7826617423809967`*^9, 
   3.782661812980236*^9},ExpressionUUID->"dc4ff695-c52e-4390-818d-\
d862c6a13932"]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.7826596828646803`*^9, 
  3.782659684840878*^9}},ExpressionUUID->"834e5c3b-ac9f-47a0-a6fa-\
2346d845975b"]
}, Open  ]]
}, Open  ]]
},
WindowSize->{1001, 642},
WindowMargins->{{41, Automatic}, {Automatic, 44}},
FrontEndVersion->"11.2 for Linux x86 (64-bit) (September 10, 2017)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 280, 7, 47, "Input",ExpressionUUID->"7e8350ee-f09c-4d78-ac3b-7f39904190c4",
 InitializationCell->True],
Cell[841, 29, 237, 4, 58, "Text",ExpressionUUID->"caed0e72-f004-4254-b57a-ccde69ef9676"],
Cell[CellGroupData[{
Cell[1103, 37, 79, 0, 98, "Title",ExpressionUUID->"77d08537-f7f3-4451-808d-7153544399d8"],
Cell[1185, 39, 244, 6, 47, "Input",ExpressionUUID->"7484096c-cb41-4df2-b847-1c41437d49ea",
 InitializationCell->True],
Cell[1432, 47, 593, 16, 47, "Input",ExpressionUUID->"5789a40d-5130-461c-b144-b55de59e0464",
 InitializationCell->True],
Cell[CellGroupData[{
Cell[2050, 67, 228, 4, 68, "Section",ExpressionUUID->"588fbb2d-502a-4553-92c3-4c3a45b5df12"],
Cell[2281, 73, 385, 8, 58, "Text",ExpressionUUID->"b49dc30d-f685-4917-888e-8925caca7789"],
Cell[2669, 83, 3223, 85, 216, "Input",ExpressionUUID->"8128455f-421f-4879-bdbd-fb8f6064fdbe"],
Cell[5895, 170, 916, 14, 31, "Input",ExpressionUUID->"77587275-3b49-4031-a9b4-9c149aeaac7b"],
Cell[CellGroupData[{
Cell[6836, 188, 1381, 29, 55, "Input",ExpressionUUID->"c8b6ea22-a1ba-4c97-aac5-8b2874ed978b"],
Cell[CellGroupData[{
Cell[8242, 221, 1171, 20, 25, "Print",ExpressionUUID->"142be610-21f4-445d-b0a6-ba30920fb274"],
Cell[9416, 243, 1055, 18, 25, "Print",ExpressionUUID->"8b1bca75-7853-476d-9cb8-9f14d120e0ba"]
}, Open  ]],
Cell[10486, 264, 2405, 61, 82, "Output",ExpressionUUID->"ec83d663-d7e5-4054-9491-3da493055aa8"]
}, Open  ]],
Cell[CellGroupData[{
Cell[12928, 330, 3241, 61, 78, "Input",ExpressionUUID->"7eee7402-8768-40f3-bd8d-427eaaa99c77"],
Cell[CellGroupData[{
Cell[16194, 395, 872, 16, 25, "Print",ExpressionUUID->"9aeb4762-51ad-42ed-914c-50f5bd103df7"],
Cell[17069, 413, 841, 18, 25, "Print",ExpressionUUID->"d9afd4f1-92ac-47d5-a513-85db1bbdd34f"]
}, Open  ]],
Cell[17925, 434, 2803, 57, 82, "Output",ExpressionUUID->"1b55b683-3da5-4077-be59-86bc3955046a"]
}, Open  ]],
Cell[20743, 494, 161, 3, 31, "Input",ExpressionUUID->"818e24ff-dc61-4d09-a373-e51326e36275"],
Cell[CellGroupData[{
Cell[20929, 501, 4039, 72, 101, "Input",ExpressionUUID->"cf2e1a15-8016-42d8-bcde-aa0802c3be56"],
Cell[CellGroupData[{
Cell[24993, 577, 1814, 31, 25, "Print",ExpressionUUID->"ef1c6de6-db8b-4780-bfb7-65aff542d63d"],
Cell[26810, 610, 2644, 41, 400, "Print",ExpressionUUID->"17e75830-6a9d-4e3c-9518-949ce59ac784"],
Cell[29457, 653, 1726, 28, 25, "Print",ExpressionUUID->"240ea087-643f-43f6-ae84-65167f3ae102"]
}, Open  ]],
Cell[31198, 684, 3669, 71, 82, "Output",ExpressionUUID->"dc4ff695-c52e-4390-818d-d862c6a13932"]
}, Open  ]],
Cell[34882, 758, 154, 3, 31, "Input",ExpressionUUID->"834e5c3b-ac9f-47a0-a6fa-2346d845975b"]
}, Open  ]]
}, Open  ]]
}
]
*)