MVlookup.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[     76683,       1756]
NotebookOptionsPosition[     73146,       1637]
NotebookOutlinePosition[     73485,       1652]
CellTagsIndexPosition[     73442,       1649]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"cost", " ", "of", " ", "MV", " ", "sumcheck"}], ",", " ", 
    RowBox[{
     RowBox[{"relative", " ", "to", " ", "the", " ", "size"}], "|", "H", "|", 
     " ", 
     RowBox[{
      RowBox[{"of", " ", "the", " ", "n"}], "-", 
      RowBox[{"dimensional", " ", "hypercube"}]}]}]}], "\n", "*)"}], "\n", 
  RowBox[{
   RowBox[{
    RowBox[{"ClearAll", "@", "sumcheck"}], ";"}], "\n", 
   RowBox[{
    RowBox[{"sumcheck", "[", 
     RowBox[{"n_", ",", "d_", ",", "numVars_", ",", 
      RowBox[{"{", 
       RowBox[{"Qm_", ",", "Qs_", ",", "Qa_"}], "}"}]}], "]"}], ":=", "\n", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"1", " ", "-", " ", 
       RowBox[{"1", "/", 
        RowBox[{"2", "^", 
         RowBox[{"(", "n", ")"}]}]}]}], ")"}], " ", "*", " ", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{
        RowBox[{"d", "*", "numVars"}], " ", "+", " ", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"d", "+", "1"}], ")"}], " ", "*", " ", "Qm"}]}], ",", " ", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"numVars", "+", " ", "d", " ", "+", " ", "1"}], ")"}], " ", 
        "*", " ", "Qs"}], ",", " ", 
       RowBox[{
        RowBox[{"d", "*", "numVars"}], " ", "+", " ", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"d", "+", "1"}], ")"}], "*", 
         RowBox[{"(", 
          RowBox[{"Qa", " ", "+", " ", "1"}], ")"}]}]}]}], "}"}], 
     " "}]}]}]}]], "Code",
 CellChangeTimes->{{3.8750892597671213`*^9, 3.8750893986050997`*^9}, {
   3.8750894332159452`*^9, 3.875089434032494*^9}, {3.875089522618288*^9, 
   3.8750895264495363`*^9}, {3.875090719855582*^9, 3.875090921681652*^9}, 
   3.875091020506954*^9, 3.875095525984243*^9, {3.875464062927775*^9, 
   3.875464063514955*^9}, {3.8754642664194927`*^9, 3.875464274987112*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"From", " ", "Berstein", " ", "and", " ", "Lange"}], ",", " ", 
    "2007", ",", " ", 
    RowBox[{
     RowBox[{
      RowBox[{
      "Analysis", " ", "of", " ", "EC", " ", "single", " ", "scalar", " ", 
       "multiplication", "\[IndentingNewLine]", "costPerScalarBit"}], " ", 
      "=", " ", "8"}], ";", " ", 
     RowBox[{"(*", " ", 
      RowBox[{"InvEdwards", " ", "coordinates"}], " ", "*)"}], "\n", 
     RowBox[{"costPerScalarMul", " ", "=", " ", 
      RowBox[{
       RowBox[{"8", "*", "255", "\n", "costPippengerPerScalar"}], " ", "=", 
       " ", 
       RowBox[{"costPerScalarMul", "/", "12"}]}]}]}]}], " ", 
   RowBox[{"(*", " ", 
    RowBox[{
     RowBox[{
      RowBox[{"1", " ", "/", " ", "log"}], "|", "H", "|"}], ",", " ", 
     RowBox[{
      RowBox[{"with", " ", "|", "H", "|"}], "=", 
      RowBox[{"2", "^", "12"}]}]}], " ", "*)"}], " ", "\n", "*)"}], "\n", 
  "\n", 
  RowBox[{"(*", " ", "\n", 
   RowBox[{
    RowBox[{
     RowBox[{
     "Cost", " ", "of", " ", "Pippenger", " ", "in", " ", "terms", " ", "of", 
      " ", "field", " ", 
      RowBox[{"multiplications", ".", "\n", "Based"}], " ", "on", " ", 
      "Marcin", " ", "single"}], "-", 
     RowBox[{"threaded", " ", "benchmarks", " ", "of", " ", "Pippenger", "\n", 
      RowBox[{"https", ":"}]}]}], "//", 
    RowBox[{
     RowBox[{
      RowBox[{"github", ".", "com"}], "/", "Orbis"}], "-", 
     RowBox[{
      RowBox[{
       RowBox[{"Tertius", "/", "Orbis"}], "/", "issues"}], "/", "97"}]}]}], 
   "\n", "*)"}], "\n", "\n", "\n", 
  RowBox[{
   RowBox[{
    RowBox[{"instanceSizes", " ", "=", " ", 
     RowBox[{"Range", "[", 
      RowBox[{"12", ",", "18"}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"benchesPippenger", "=", " ", 
     RowBox[{"{", 
      RowBox[{
      "46.010", ",", " ", "82.382", ",", " ", "153.67", ",", " ", "279.35", 
       ",", " ", "522.13", ",", " ", "1003.4", ",", " ", "1869.5"}], "}"}]}], 
    ";"}], " ", 
   RowBox[{"(*", " ", 
    RowBox[{
     RowBox[{
      RowBox[{"for", " ", "N"}], "=", "12"}], ",", " ", "13", ",", "...", ",",
      " ", "18"}], "*)"}], "\n", 
   RowBox[{
    RowBox[{"benchesFieldMul", " ", "=", " ", 
     RowBox[{"{", 
      RowBox[{
      "21.110", ",", " ", "42.315", ",", " ", "84.777", ",", " ", "169.70", 
       ",", " ", "339.25", ",", " ", "679.27", ",", " ", "1360.6"}], "}"}]}], 
    ";"}], " ", 
   RowBox[{"(*", " ", 
    RowBox[{"for", " ", 
     RowBox[{"2", "^", "8"}], " ", "*", " ", "N"}], " ", "*)"}], "\n", 
   RowBox[{"equivalentMuls", " ", "=", " ", 
    RowBox[{
     RowBox[{"benchesPippenger", "/", "benchesFieldMul"}], " ", "*", " ", 
     RowBox[{"2", "^", "8"}]}]}]}]}]], "Code",
 CellChangeTimes->{
  3.876385782402144*^9, {3.876546678241295*^9, 3.876546833037793*^9}, {
   3.8765470262945023`*^9, 3.876547038754353*^9}, {3.87655572188908*^9, 
   3.876556138055821*^9}, {3.8765561847248993`*^9, 3.876556222062521*^9}, 
   3.876556258761971*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "557.9611558503079`", ",", "498.39990547087325`", ",", 
   "464.03529259115084`", ",", "421.41190335886864`", ",", 
   "394.00229918938834`", ",", "378.1565504144155`", ",", 
   "351.75069822137294`"}], "}"}]], "Output",
 CellChangeTimes->{
  3.876556259442092*^9, 3.876556717781522*^9, 3.8765601122094593`*^9, 
   3.8765601718109207`*^9, 3.876560948777672*^9, {3.876565089141786*^9, 
   3.876565141901601*^9}, 3.876572676005343*^9, 3.8769972609926777`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"maxM", " ", "=", " ", "200"}], ";"}], "\n", 
 RowBox[{
  RowBox[{"domainSize", " ", "=", " ", "12"}], ";"}], "\n", 
 RowBox[{"costPippengerPerScalar", " ", "=", " ", 
  RowBox[{"equivalentMuls", "[", 
   RowBox[{"[", 
    RowBox[{
     RowBox[{
      RowBox[{"Position", "[", 
       RowBox[{"instanceSizes", ",", " ", "domainSize"}], "]"}], "//", 
      "Flatten"}], "//", "First"}], "]"}], "]"}]}]}], "Code",
 CellChangeTimes->{{3.876386163553298*^9, 3.8763861664544897`*^9}, {
  3.876546999625012*^9, 3.876546999759676*^9}, {3.876556226246236*^9, 
  3.876556226492296*^9}, {3.8765567156743507`*^9, 3.876556716055079*^9}, {
  3.8765601100124598`*^9, 3.876560110598651*^9}, {3.876560170382987*^9, 
  3.876560170719102*^9}, {3.876560946908489*^9, 3.876560947221896*^9}, {
  3.876565087382625*^9, 3.876565139674706*^9}, {3.876572674254801*^9, 
  3.876572674499219*^9}}],

Cell[BoxData["557.9611558503079`"], "Output",
 CellChangeTimes->{{3.876556228208119*^9, 3.8765562564138317`*^9}, {
  3.8765567178043118`*^9, 3.87655671821242*^9}, {3.8765601122130632`*^9, 
  3.876560112567276*^9}, {3.876560171830325*^9, 3.876560172126446*^9}, {
  3.876560948797903*^9, 3.8765609491039953`*^9}, {3.8765650891547327`*^9, 
  3.876565142228589*^9}, {3.876572676018879*^9, 3.876572676427266*^9}, {
  3.876997261037283*^9, 3.8769972627925158`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Based on the logarithmic derivative", "Section",
 CellChangeTimes->{{3.876385540206491*^9, 3.8763855496770477`*^9}}],

Cell[CellGroupData[{

Cell["small number of columns", "Subsection",
 CellChangeTimes->{{3.87657263997535*^9, 3.876572644391529*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"arithComp", " ", "=", " ", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{"4", "M"}], " ", "+", " ", "1"}], ",", "1", ",", 
     RowBox[{"M", "+", "1"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"costsumcheck", " ", "=", " ", 
  RowBox[{
   RowBox[{
    RowBox[{"sumcheck", "[", 
     RowBox[{"n", ",", " ", 
      RowBox[{"M", "+", "3"}], ",", 
      RowBox[{"M", "+", "4"}], ",", " ", "arithComp"}], " ", "]"}], "/", 
    RowBox[{"(", 
     RowBox[{"1", " ", "-", " ", 
      RowBox[{"1", "/", 
       RowBox[{"2", "^", "n"}]}]}], ")"}]}], "//", 
   "Simplify"}]}], "\[IndentingNewLine]", 
 RowBox[{"cost", " ", "=", " ", 
  RowBox[{
   RowBox[{"costsumcheck", "+", " ", 
    RowBox[{"{", 
     RowBox[{"7", ",", " ", "1", ",", " ", 
      RowBox[{
       RowBox[{"2", "M"}], "+", "1"}]}], "}"}]}], "//", 
   "Simplify"}]}]}], "Code",
 CellChangeTimes->{{3.8750909259659147`*^9, 3.8750909755423727`*^9}, {
   3.875091011817265*^9, 3.875091011971953*^9}, {3.875091761347869*^9, 
   3.8750917685033712`*^9}, {3.875095922978674*^9, 3.875095924237027*^9}, {
   3.875102731687755*^9, 3.875102732249934*^9}, {3.8751027766648827`*^9, 
   3.8751027771588306`*^9}, {3.8751035450282717`*^9, 3.875103545106289*^9}, {
   3.875103575541709*^9, 3.875103575695881*^9}, {3.87544903570872*^9, 
   3.875449051990718*^9}, {3.875459791604109*^9, 3.875459793447543*^9}, {
   3.875460098494342*^9, 3.875460146550344*^9}, 3.8754622749212112`*^9, {
   3.8754623550682497`*^9, 3.875462355922514*^9}, {3.8754628105165787`*^9, 
   3.8754628283143*^9}, {3.875462860582975*^9, 3.875462874171104*^9}, {
   3.8754629739336987`*^9, 3.8754629980591993`*^9}, {3.875464304445507*^9, 
   3.8754643097529507`*^9}, {3.875464374646179*^9, 3.875464379549264*^9}, {
   3.875464667348877*^9, 3.875464667406427*^9}, {3.875464789979189*^9, 
   3.87546480291568*^9}, {3.875464900566935*^9, 3.875464925993936*^9}, {
   3.8754665630601683`*^9, 3.875466563364303*^9}, {3.875603537872497*^9, 
   3.875603540199489*^9}, {3.87560358270327*^9, 3.875603659183642*^9}, {
   3.875604134438981*^9, 3.875604134819813*^9}, 3.875604513449641*^9, {
   3.8757017826018*^9, 3.87570182248929*^9}, {3.876295114711844*^9, 
   3.876295115023265*^9}, {3.876385582696113*^9, 3.876385587334529*^9}, {
   3.8763857218627453`*^9, 3.8763857221737823`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"16", "+", 
    RowBox[{"24", " ", "M"}], "+", 
    RowBox[{"5", " ", 
     SuperscriptBox["M", "2"]}]}], ",", 
   RowBox[{"2", " ", 
    RowBox[{"(", 
     RowBox[{"4", "+", "M"}], ")"}]}], ",", 
   RowBox[{"20", "+", 
    RowBox[{"13", " ", "M"}], "+", 
    RowBox[{"2", " ", 
     SuperscriptBox["M", "2"]}]}]}], "}"}]], "Output",
 CellChangeTimes->{{3.876385715009602*^9, 3.876385722596904*^9}, 
   3.8763867511917953`*^9, 3.876388874745801*^9, 3.876545897342609*^9, 
   3.876546877593207*^9, 3.876547012008567*^9, 3.8765562564425907`*^9, 
   3.8765567179227877`*^9, 3.876560112246744*^9, 3.876560171852413*^9, 
   3.8765609488236647`*^9, {3.876565089189497*^9, 3.87656514196437*^9}, 
   3.876572676140774*^9, 3.876997261144877*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"23", "+", 
    RowBox[{"24", " ", "M"}], "+", 
    RowBox[{"5", " ", 
     SuperscriptBox["M", "2"]}]}], ",", 
   RowBox[{"9", "+", 
    RowBox[{"2", " ", "M"}]}], ",", 
   RowBox[{"21", "+", 
    RowBox[{"15", " ", "M"}], "+", 
    RowBox[{"2", " ", 
     SuperscriptBox["M", "2"]}]}]}], "}"}]], "Output",
 CellChangeTimes->{{3.876385715009602*^9, 3.876385722596904*^9}, 
   3.8763867511917953`*^9, 3.876388874745801*^9, 3.876545897342609*^9, 
   3.876546877593207*^9, 3.876547012008567*^9, 3.8765562564425907`*^9, 
   3.8765567179227877`*^9, 3.876560112246744*^9, 3.876560171852413*^9, 
   3.8765609488236647`*^9, {3.876565089189497*^9, 3.87656514196437*^9}, 
   3.876572676140774*^9, 3.876997261145756*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["large number of columns", "Subsection",
 CellChangeTimes->{{3.876385600508483*^9, 3.876385604712728*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"arithComp", " ", "=", " ", 
   RowBox[{"{", 
    RowBox[{"4", ",", "1", ",", "2"}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"costsumcheckLargeM", " ", "=", " ", 
  RowBox[{
   RowBox[{"M", " ", "*", " ", 
    RowBox[{
     RowBox[{"sumcheck", "[", 
      RowBox[{
       RowBox[{"n", "+", "m"}], ",", " ", "4", ",", " ", "5", ",", " ", 
       "arithComp"}], "]"}], "/", 
     RowBox[{"(", 
      RowBox[{"1", " ", "-", " ", 
       RowBox[{"1", "/", 
        RowBox[{"2", "^", 
         RowBox[{"(", 
          RowBox[{"n", "+", "m"}], ")"}]}]}]}], ")"}]}]}], "//", 
   "Simplify"}]}], "\[IndentingNewLine]", 
 RowBox[{"costLargeM", " ", "=", " ", 
  RowBox[{
   RowBox[{"costsumcheckLargeM", " ", "+", " ", 
    RowBox[{"M", " ", "*", " ", 
     RowBox[{"{", 
      RowBox[{"4", ",", " ", "0", ",", " ", "2"}], "}"}]}]}], "//", 
   "Simplify"}]}]}], "Code",
 CellChangeTimes->{{3.8750909259659147`*^9, 3.8750909755423727`*^9}, {
   3.875091011817265*^9, 3.875091011971953*^9}, {3.875091761347869*^9, 
   3.8750917685033712`*^9}, {3.875095922978674*^9, 3.875095924237027*^9}, {
   3.875102731687755*^9, 3.875102732249934*^9}, {3.8751027766648827`*^9, 
   3.8751027771588306`*^9}, {3.8751035450282717`*^9, 3.875103545106289*^9}, {
   3.875103575541709*^9, 3.875103575695881*^9}, {3.87544903570872*^9, 
   3.875449051990718*^9}, {3.875459791604109*^9, 3.875459793447543*^9}, {
   3.875460098494342*^9, 3.875460146550344*^9}, 3.8754622749212112`*^9, {
   3.8754623550682497`*^9, 3.875462355922514*^9}, {3.8754628105165787`*^9, 
   3.8754628283143*^9}, {3.875462860582975*^9, 3.875462874171104*^9}, {
   3.8754629739336987`*^9, 3.8754629980591993`*^9}, {3.875464304445507*^9, 
   3.8754643097529507`*^9}, {3.875464374646179*^9, 3.875464379549264*^9}, {
   3.875464667348877*^9, 3.875464667406427*^9}, {3.875464789979189*^9, 
   3.87546480291568*^9}, {3.875464900566935*^9, 3.875464925993936*^9}, {
   3.8754665630601683`*^9, 3.875466563364303*^9}, {3.875603537872497*^9, 
   3.875603540199489*^9}, {3.87560358270327*^9, 3.875603659183642*^9}, {
   3.875604134438981*^9, 3.875604134819813*^9}, 3.875604513449641*^9, {
   3.8757017826018*^9, 3.87570182248929*^9}, {3.876294935922762*^9, 
   3.876294976871594*^9}, {3.8762950096121683`*^9, 3.87629504286515*^9}, {
   3.876295088768346*^9, 3.876295091983234*^9}, {3.876295140653759*^9, 
   3.8762952201241207`*^9}, {3.8762954318300037`*^9, 3.876295438845056*^9}, {
   3.876368111846792*^9, 3.876368114205063*^9}, {3.876385621768506*^9, 
   3.876385712184557*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"40", " ", "M"}], ",", 
   RowBox[{"10", " ", "M"}], ",", 
   RowBox[{"35", " ", "M"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.876385715044609*^9, 3.876385715142146*^9}, 
   3.87638675119805*^9, 3.8763888747507277`*^9, 3.8765458973800983`*^9, 
   3.876546877627763*^9, 3.876547012024343*^9, 3.876556256464425*^9, 
   3.8765567179640636`*^9, 3.876560112250761*^9, 3.8765601718693943`*^9, 
   3.8765609488288097`*^9, {3.876565089207079*^9, 3.876565141969204*^9}, 
   3.876572676165784*^9, 3.876997261153358*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"44", " ", "M"}], ",", 
   RowBox[{"10", " ", "M"}], ",", 
   RowBox[{"37", " ", "M"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.876385715044609*^9, 3.876385715142146*^9}, 
   3.87638675119805*^9, 3.8763888747507277`*^9, 3.8765458973800983`*^9, 
   3.876546877627763*^9, 3.876547012024343*^9, 3.876556256464425*^9, 
   3.8765567179640636`*^9, 3.876560112250761*^9, 3.8765601718693943`*^9, 
   3.8765609488288097`*^9, {3.876565089207079*^9, 3.876565141969204*^9}, 
   3.876572676165784*^9, 3.8769972611576033`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["break even point", "Subsection",
 CellChangeTimes->{{3.876385748469411*^9, 3.876385752247199*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"costWithComm", " ", "=", " ", 
  RowBox[{"cost", " ", "+", " ", 
   RowBox[{"costPippengerPerScalar", "*", "2"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"costWithCommLarge", " ", "=", " ", 
   RowBox[{"costLargeM", " ", "+", " ", 
    RowBox[{"costPippengerPerScalar", "*", 
     RowBox[{"(", 
      RowBox[{"M", "+", "1"}], ")"}]}]}]}], "\n", "\n", 
  RowBox[{"(*", " ", 
   RowBox[{
   "break", " ", "even", " ", "point", " ", "with", " ", "linear", " ", 
    "strategy"}], " ", "*)"}]}], "\n", 
 RowBox[{"sol", " ", "=", " ", 
  RowBox[{"NSolve", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"costWithComm", "[", 
      RowBox[{"[", "1", "]"}], "]"}], " ", "\[Equal]", " ", 
     RowBox[{"costWithCommLarge", "[", 
      RowBox[{"[", "1", "]"}], "]"}]}], ",", "M"}], "]"}]}], "\n", 
 RowBox[{
  RowBox[{"breakEven", " ", "=", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"M", " ", "/.", " ", "sol"}], " ", "//", " ", "Last"}], "//", 
    "Floor"}]}], "\n", "\n", 
  RowBox[{"(*", " ", 
   RowBox[{"maximum", " ", "speedup"}], " ", "*)"}]}], "\n", 
 RowBox[{
  RowBox[{"ratio", " ", "=", " ", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"costWithCommLarge", "[", 
        RowBox[{"[", "1", "]"}], "]"}], "/", 
       RowBox[{"costWithComm", "[", 
        RowBox[{"[", "1", "]"}], "]"}]}], "/.", 
      RowBox[{"{", 
       RowBox[{"M", "\[Rule]", "x"}], "}"}]}], ",", "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"x", ",", "1", ",", " ", "breakEven"}], "}"}]}], "]"}]}], 
  ";"}], "\n", 
 RowBox[{
  RowBox[{"max", " ", "=", " ", 
   RowBox[{"Max", "@", "ratio"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"max", "//", "N"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Position", "[", 
    RowBox[{"ratio", ",", "max"}], "]"}], "//", "Flatten"}], "\n", "\n", 
  RowBox[{"(*", " ", 
   RowBox[{
   "range", " ", "for", " ", "0.8", " ", "of", " ", "the", " ", "maximum", 
    " ", "speedup"}], " ", "*)"}]}], "\n", 
 RowBox[{"threshold", " ", "=", " ", "1.5"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Position", "[", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"#", "\[GreaterEqual]", "threshold"}], ")"}], "&"}], "/@", 
      "ratio"}], ",", "True"}], "]"}], "//", "Flatten"}], "\n"}], "\n", 
 RowBox[{"ListPlot", "[", 
  RowBox[{"ratio", ",", 
   RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]}], "Code",
 CellChangeTimes->{{3.8763858044317017`*^9, 3.876385921151272*^9}, {
  3.876386743411701*^9, 3.876386838144216*^9}, {3.8763889126917133`*^9, 
  3.8763889257377234`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"1138.9223117006159`", "\[VeryThinSpace]", "+", 
    RowBox[{"24", " ", "M"}], "+", 
    RowBox[{"5", " ", 
     SuperscriptBox["M", "2"]}]}], ",", 
   RowBox[{"1124.9223117006159`", "\[VeryThinSpace]", "+", 
    RowBox[{"2", " ", "M"}]}], ",", 
   RowBox[{"1136.9223117006159`", "\[VeryThinSpace]", "+", 
    RowBox[{"15", " ", "M"}], "+", 
    RowBox[{"2", " ", 
     SuperscriptBox["M", "2"]}]}]}], "}"}]], "Output",
 CellChangeTimes->{{3.876385910047246*^9, 3.876385923439322*^9}, {
   3.8763867512289867`*^9, 3.8763868387229633`*^9}, {3.876388874777326*^9, 
   3.876388875009164*^9}, 3.876388926834021*^9, 3.876545897410365*^9, 
   3.876546877634027*^9, 3.876547012070548*^9, 3.876556256487472*^9, 
   3.876556717970819*^9, 3.876560112282793*^9, 3.8765601718768682`*^9, 
   3.8765609488515263`*^9, {3.876565089213463*^9, 3.8765651419895573`*^9}, 
   3.876572676171049*^9, 3.876997261206531*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{
    RowBox[{"44", " ", "M"}], "+", 
    RowBox[{"557.9611558503079`", " ", 
     RowBox[{"(", 
      RowBox[{"1", "+", "M"}], ")"}]}]}], ",", 
   RowBox[{
    RowBox[{"10", " ", "M"}], "+", 
    RowBox[{"557.9611558503079`", " ", 
     RowBox[{"(", 
      RowBox[{"1", "+", "M"}], ")"}]}]}], ",", 
   RowBox[{
    RowBox[{"37", " ", "M"}], "+", 
    RowBox[{"557.9611558503079`", " ", 
     RowBox[{"(", 
      RowBox[{"1", "+", "M"}], ")"}]}]}]}], "}"}]], "Output",
 CellChangeTimes->{{3.876385910047246*^9, 3.876385923439322*^9}, {
   3.8763867512289867`*^9, 3.8763868387229633`*^9}, {3.876388874777326*^9, 
   3.876388875009164*^9}, 3.876388926834021*^9, 3.876545897410365*^9, 
   3.876546877634027*^9, 3.876547012070548*^9, 3.876556256487472*^9, 
   3.876556717970819*^9, 3.876560112282793*^9, 3.8765601718768682`*^9, 
   3.8765609488515263`*^9, {3.876565089213463*^9, 3.8765651419895573`*^9}, 
   3.876572676171049*^9, 3.876997261212948*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"M", "\[Rule]", "1.0140872173493711`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"M", "\[Rule]", "114.57814395271221`"}], "}"}]}], "}"}]], "Output",\

 CellChangeTimes->{{3.876385910047246*^9, 3.876385923439322*^9}, {
   3.8763867512289867`*^9, 3.8763868387229633`*^9}, {3.876388874777326*^9, 
   3.876388875009164*^9}, 3.876388926834021*^9, 3.876545897410365*^9, 
   3.876546877634027*^9, 3.876547012070548*^9, 3.876556256487472*^9, 
   3.876556717970819*^9, 3.876560112282793*^9, 3.8765601718768682`*^9, 
   3.8765609488515263`*^9, {3.876565089213463*^9, 3.8765651419895573`*^9}, 
   3.876572676171049*^9, 3.876997261254324*^9}],

Cell[BoxData["114"], "Output",
 CellChangeTimes->{{3.876385910047246*^9, 3.876385923439322*^9}, {
   3.8763867512289867`*^9, 3.8763868387229633`*^9}, {3.876388874777326*^9, 
   3.876388875009164*^9}, 3.876388926834021*^9, 3.876545897410365*^9, 
   3.876546877634027*^9, 3.876547012070548*^9, 3.876556256487472*^9, 
   3.876556717970819*^9, 3.876560112282793*^9, 3.8765601718768682`*^9, 
   3.8765609488515263`*^9, {3.876565089213463*^9, 3.8765651419895573`*^9}, 
   3.876572676171049*^9, 3.876997261255728*^9}],

Cell[BoxData["3.6601636210353585`"], "Output",
 CellChangeTimes->{{3.876385910047246*^9, 3.876385923439322*^9}, {
   3.8763867512289867`*^9, 3.8763868387229633`*^9}, {3.876388874777326*^9, 
   3.876388875009164*^9}, 3.876388926834021*^9, 3.876545897410365*^9, 
   3.876546877634027*^9, 3.876547012070548*^9, 3.876556256487472*^9, 
   3.876556717970819*^9, 3.876560112282793*^9, 3.8765601718768682`*^9, 
   3.8765609488515263`*^9, {3.876565089213463*^9, 3.8765651419895573`*^9}, 
   3.876572676171049*^9, 3.8769972612643023`*^9}],

Cell[BoxData[
 RowBox[{"{", "14", "}"}]], "Output",
 CellChangeTimes->{{3.876385910047246*^9, 3.876385923439322*^9}, {
   3.8763867512289867`*^9, 3.8763868387229633`*^9}, {3.876388874777326*^9, 
   3.876388875009164*^9}, 3.876388926834021*^9, 3.876545897410365*^9, 
   3.876546877634027*^9, 3.876547012070548*^9, 3.876556256487472*^9, 
   3.876556717970819*^9, 3.876560112282793*^9, 3.8765601718768682`*^9, 
   3.8765609488515263`*^9, {3.876565089213463*^9, 3.8765651419895573`*^9}, 
   3.876572676171049*^9, 3.87699726126502*^9}],

Cell[BoxData["1.5`"], "Output",
 CellChangeTimes->{{3.876385910047246*^9, 3.876385923439322*^9}, {
   3.8763867512289867`*^9, 3.8763868387229633`*^9}, {3.876388874777326*^9, 
   3.876388875009164*^9}, 3.876388926834021*^9, 3.876545897410365*^9, 
   3.876546877634027*^9, 3.876547012070548*^9, 3.876556256487472*^9, 
   3.876556717970819*^9, 3.876560112282793*^9, 3.8765601718768682`*^9, 
   3.8765609488515263`*^9, {3.876565089213463*^9, 3.8765651419895573`*^9}, 
   3.876572676171049*^9, 3.876997261265718*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "3", ",", "4", ",", "5", ",", "6", ",", "7", ",", "8", ",", "9", ",", "10", 
   ",", "11", ",", "12", ",", "13", ",", "14", ",", "15", ",", "16", ",", 
   "17", ",", "18", ",", "19", ",", "20", ",", "21", ",", "22", ",", "23", 
   ",", "24", ",", "25", ",", "26", ",", "27", ",", "28", ",", "29", ",", 
   "30", ",", "31", ",", "32", ",", "33", ",", "34", ",", "35", ",", "36", 
   ",", "37", ",", "38", ",", "39", ",", "40", ",", "41", ",", "42", ",", 
   "43", ",", "44", ",", "45", ",", "46", ",", "47", ",", "48", ",", "49", 
   ",", "50", ",", "51", ",", "52", ",", "53", ",", "54", ",", "55", ",", 
   "56", ",", "57", ",", "58", ",", "59", ",", "60", ",", "61", ",", "62", 
   ",", "63", ",", "64", ",", "65", ",", "66", ",", "67", ",", "68", ",", 
   "69", ",", "70", ",", "71", ",", "72", ",", "73"}], "}"}]], "Output",
 CellChangeTimes->{{3.876385910047246*^9, 3.876385923439322*^9}, {
   3.8763867512289867`*^9, 3.8763868387229633`*^9}, {3.876388874777326*^9, 
   3.876388875009164*^9}, 3.876388926834021*^9, 3.876545897410365*^9, 
   3.876546877634027*^9, 3.876547012070548*^9, 3.876556256487472*^9, 
   3.876556717970819*^9, 3.876560112282793*^9, 3.8765601718768682`*^9, 
   3.8765609488515263`*^9, {3.876565089213463*^9, 3.8765651419895573`*^9}, 
   3.876572676171049*^9, 3.876997261266651*^9}],

Cell[BoxData[
 GraphicsBox[{{}, 
   {RGBColor[0.24720000000000014`, 0.24, 0.6], PointBox[CompressedData["
1:eJw11AlMFFccBvCBlcohLFlbpcihlHpziMghKN8Cyy7Xwi7sArsUK9Ei0Qoo
HpGqE8GgDTEIseJ6ISlqKYpApWK7dKG2KngWMKKirBZEDWFFRUTB2vh/k0wm
v/neN5mZl/dmpGcpV5pzHLfuw/n/9eNhCmm2Xzb28OJgCN2AfFK9a4b2NdkS
Db13rZTO42R7bMPAu1Abc3z0p3B7nr34wtEJZAe49555o736CdkJezsbvUfL
LMnT8avnZTdna2uyG0RaG1PTHBuyO7zLhn02mU0iz8S1Ydlwp455NiwyBv9+
Mc48F4eW9XZ5+NmS56PmtrJbIGP2ROuhKWdtg5i9If99nDP7jHkBsvv65xju
sOf5YHmkIHttMfNCmOqf/LszkNkXwWmqpIHb7H0XIb9QLbDNYfZDUKUxwyhg
9keO7fh7m1L2vQGwrM6+bXRlDkSVr7nz5p+tyItR2bu79bQfcxBOy/nzgS3s
/wUjs0W49KSceQkOBDp8t/reRPJSbF/t1Xx3NXMIOlMznwWPsfkAOuTqKS3F
ZB7wtcs/XjyT5WLcbXS83tJkQbkYz2a4ncxMIXOh2Oc0u3/PMM03H4qpg/MV
fj+w+Q9D62Zjd1oAy8OgWTPrvKBbQHk4Vpb0CJ0LyHw4Sov21JR7sFwCx+uj
uh1d5pRLEP140qazhWQuAhs3iJ7MDWB5BIYquEftT80ol+KrmoPG8qNkXoqs
+yaXHWqWy2D9yv23DCHLZSioE3eLWznKIxHisOIqV0jmI+Hhn2O1T8LyKHy9
XqIatGB5FBbeUuwxK35P6yMalgPpFRumkvlouBwMLnapYOsnBl+Uht3o8SLz
Maj1e1lzzDBGeSyWvlWqpIlkPhZVp2YYzj19R7kczV2yy935ZF6OrU03px91
ZXkcpFv1UX36t5THYYvoje5wGpmLR1xthfcZMzLioX/dZrI6MUrj46Fsu8fv
l5MN8YiN+wURI2+or4DEvjx8wo9kKJAoXBJnUJB5BfzKlW2ZZmSDAiWv6v8y
1Y1QX4nuYz6PJd+QoUSxR9VazTQyr8QCx3nzpv9D+49BiQTPqI27v2f7UQIy
qvd77g4nIwGNInGdiCPzCfjS40HX5/ph6ifg3N6cFyV5ZC4RCtEu/bogMhIh
DNSln373ivqJaCgsyg36g2z4MP5xer8wn8ypEL7FVeQmI0OF0Y4r5avsWF8F
6bY63Z3Ol9RXIXdbXUrWETKnRm7ttUKXVWSoEVLQMa3Hh8yrYdfkoa4ef0F9
NfqHOvR5bWQuCaFL1jhLD5CRBK8HufzEVWQ+CZXZx70b/Fn/w3jpivoYS9ZP
xo22uIsXuoaonwyrvOHlDtVkPhmGZN0p8XayIRncickTgxLIXAok+ta347NY
PwU3+y5Jd409p34KLBc/aL/ZTjak4PDKh/W3qsicBk55o+4lO8jQ4N4li6L3
GtbXQNBjm+voy/oa2AxcLbpvy/pahJVG/iTtN1FfC/HO9U2KP8m8Fv7Oj4wj
R8gGLSLEZbLAPDKXCrv+KxJhMuunIkDXHrlpEeunIlRlbvx2sinkP+75BZs=

     "]]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesLabel->{None, None},
  AxesOrigin->{1., 0.9931502293261688},
  Method->{},
  PlotRangeClipping->True]], "Output",
 CellChangeTimes->{{3.876385910047246*^9, 3.876385923439322*^9}, {
   3.8763867512289867`*^9, 3.8763868387229633`*^9}, {3.876388874777326*^9, 
   3.876388875009164*^9}, 3.876388926834021*^9, 3.876545897410365*^9, 
   3.876546877634027*^9, 3.876547012070548*^9, 3.876556256487472*^9, 
   3.876556717970819*^9, 3.876560112282793*^9, 3.8765601718768682`*^9, 
   3.8765609488515263`*^9, {3.876565089213463*^9, 3.8765651419895573`*^9}, 
   3.876572676171049*^9, 3.876997261345639*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Best strategy ", "Subsection",
 CellChangeTimes->{{3.876386294229916*^9, 3.876386298685363*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"bestCostWithComm", " ", "=", " ", 
   RowBox[{"Join", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"costWithComm", "[", 
         RowBox[{"[", "1", "]"}], "]"}], "/.", 
        RowBox[{"{", 
         RowBox[{"M", "\[Rule]", "x"}], "}"}]}], ",", 
       RowBox[{"{", 
        RowBox[{"x", ",", "1", ",", "breakEven"}], "}"}]}], "]"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"costWithCommLarge", "[", 
         RowBox[{"[", "1", "]"}], "]"}], "/.", 
        RowBox[{"{", 
         RowBox[{"M", "\[Rule]", "x"}], "}"}]}], ",", 
       RowBox[{"{", 
        RowBox[{"x", ",", 
         RowBox[{"breakEven", "+", "1"}], ",", "maxM"}], "}"}]}], "]"}]}], 
    "\[IndentingNewLine]", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"ListPlot", "[", "bestCostWithComm", "]"}]}], "Code",
 CellChangeTimes->{{3.8763863084943447`*^9, 3.8763864302263527`*^9}}],

Cell[BoxData[
 GraphicsBox[{{}, 
   {RGBColor[0.24720000000000014`, 0.24, 0.6], PointBox[CompressedData["
1:eJw11XlYzAkcx/GfY7fWscQuoYdZRzm2hNzSdwkRm0pW5/zmHsKu1u3Bz1rn
PitJjtUxhDZhY1VuI4lcWTlShilC0uGaNmNqre+nnqen59W7z3xn5o/mG+WP
AZqmgiBc/vj9/0/+qvGav2xpVLrXFsIv6JOLG23PToqF27J/2wp/hR4HO2K/
DXZie+2AZeyrO+Hu7JW74J7s0ATYmS1PgnvTSU/P4YUXDXBfdvVu+Fu2azLs
xl65F3Znl+2DB7DlKfBAduWf8CB29AHYgz3uIDyY3eYwPAT7v+ChZEpMdHpZ
nQ4PY78+Ag9nNxyFR7A7HINHsodmwKPYykzYk709Cx7Nvnsc9mJ3OwkTO+oU
W4LzT6N/xx5yFh1OOYc+hub6+X25MMaIzrYEnkcfy71LNjrb8gwWvLmfuIDO
tkTnoI/jHnkRnW3xzUUfz33AJXS2pctl9AncW+ahsy0NsODDve4KOttSexV9
InfrNXS2pdkN9EncHfLRJ1FWioud1eEmuu8nn0qAJba13z/ok7mfhiW21e8W
+hTuZbDEti4vQP+eu+NtdLY1Axb8uAfdQWdba2FhKvf4u2xiW8few9+jV8JG
9J2F2Ptz97mPPdtaB0voh4qwR1cVYx/A3ekB9mzrPVhCjzNhjx70EPtAvP5H
2LOtJlgKpAcFsatNAWbs2bYsWJj2yeFdS7Bn29bAEttUBRvRZ5RiH8T7HJjY
NvfH2LNNibARvfUT7KfzfgVMbFs1LLFNyjLs0e/Cwg+8n/wUe7YtG5bYphHP
sEf/GxZm8N7tOfZsWyossU3O5dij74WFYN73fIF9MDosBWNfgT16KiyE4P5L
7EPw/GApBM+/Env0bFgIxeuvwj4U7w8sheL9q8YevRoWwqiqPO+8m1cN9uz6
tbDEnpsPG8OopOHss0VfvMI+nArfJvt07A9TOOWXr0/NnAZL4ZT7cE6L6Uth
YzidKQiIfJfUuI+gY5eHXtt6sXEfQWlnnFw9Khr3EbT7aJNNBW1fYx9BO1Ke
VkUNgQU5Rcdf9WsXBpOc1sakpx9ZBUtyWr42zsE/pXEv//T5WXOtcS9S5E/y
gujXsEwkpcbbo7/jGzyeSMEhfeJueMKiSFM//vubo4IlkSZ4v5neegNsEOnj
p13WwcOwUSQPtzOOk2/DZpH69dizpKIOFhTU3XFd0cZub3FfQZ1azx7ZdxxM
Cmrb1D8+bxYsKsiudrBNvxmWFFRf0TnCPhM2KOidueFsSjFsVNDLO0+6TWjy
DvcVVHolT3rqDAtKun/ucMmaybBMSTePxY7pFQWTki6lLk7O2Q6LSjqbGN5c
fQaWlJQRO0bT7DFsUNLB9S65e+wtuK+k5OWtXMa4wWYl/RH1al1JICyoKEZ3
97m0BJapaH3YqYmyJJhUtMLfcOBcDiyqaMH4NS3lL2BJRbNHzppd36YW91Wk
cve7njAYNqoopJeHm2cobFaRf+dO0Q8kWFCTT5v66mX7YZmavJo/ntrlGkxq
GlJ36cjJV7CoJteqg+1COv6L+2rq+Tjm57pRsEFNnQsX3t6hhI1qcrgeOnjY
etisJvts2nbvECxoqCGzV+3CAlimIUtaixkd6mDSUKWh+nhG1zrc19CTuNud
grxhSUPFG08sfTsTNmjo1srE4tho2KihvPmrRw3KgM0aMs7UJ9wqggUtZUZM
qZ8nvMd9LR0KHCh3cIZJS3t9OhrTfWFRS7s8P8imzoMlLW0ZWLKqehts0NIG
l9zSTadho5Ykp7SxbqWwWUuLHDbvvW5nxX0dzf18wWdzXGGZjjTWYG2rQJh0
FFoz+lLaYljUUUBZj96+ibCko4lF9hteXIANOqL8yvIN5bBRR0Nzbk3q0+YD
7uvI7URW2mUPWNBTr8PxrfQhsExPTsmr5thJMOmp/Q7tjf37YFFPLX737T/+
Kizpqckv7pvLamCDnmoXfv3q1w423NdTVeR7/56jYLOeysRHRy8oYGEmPQjK
aa9aZ6P/AEgDEbA=
     "]]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  Method->{},
  PlotRange->{{0, 200.}, {0, 120950.19232591189`}},
  PlotRangeClipping->True,
  PlotRangePadding->{{4., 4.}, {2419.003846518238, 
   2419.003846518238}}]], "Output",
 CellChangeTimes->{
  3.8763864308388166`*^9, 3.876386751300392*^9, 3.8763888748586206`*^9, 
   3.876545897521863*^9, 3.876546877708362*^9, 3.876547012149145*^9, 
   3.876556256565036*^9, 3.876556718045046*^9, 3.876560112371064*^9, 
   3.876560171945764*^9, 3.876560948926437*^9, {3.876565089283977*^9, 
   3.876565142052981*^9}, 3.876572676253577*^9, 3.8769972614114237`*^9}]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.87638637542977*^9, 3.8763863858549757`*^9}}]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{

Cell["Based on the univariate strategy", "Section",
 CellChangeTimes->{{3.876385540206491*^9, 3.8763855496770477`*^9}, {
  3.8763859504295197`*^9, 3.876385954050096*^9}}],

Cell[CellGroupData[{

Cell["small number of columns", "Subsection",
 CellChangeTimes->{{3.87657262920255*^9, 3.8765726324865093`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"arithComp", " ", "=", " ", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{"2", 
       RowBox[{"(", 
        RowBox[{"M", "+", "1"}], ")"}]}], "+", "3"}], ",", "2", ",", "1"}], 
    "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"costsumcheckUV", " ", "=", " ", 
  RowBox[{
   RowBox[{
    RowBox[{"sumcheck", "[", 
     RowBox[{"n", ",", " ", 
      RowBox[{"M", "+", "3"}], ",", " ", 
      RowBox[{
       RowBox[{"2", "M"}], "+", "6"}], ",", " ", "arithComp"}], "]"}], "/", 
    RowBox[{"(", 
     RowBox[{"1", " ", "-", " ", 
      RowBox[{"1", "/", 
       RowBox[{"2", "^", "n"}]}]}], ")"}]}], "//", 
   "Simplify"}]}], "\[IndentingNewLine]", 
 RowBox[{"costUV", " ", "=", " ", 
  RowBox[{
   RowBox[{"costsumcheckUV", "+", " ", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{
       RowBox[{"2", 
        RowBox[{"(", 
         RowBox[{"M", "+", "1"}], ")"}]}], "+", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"2", "M"}], "+", "4"}], ")"}], "+", "2"}], ",", " ", "0", 
      ",", " ", 
      RowBox[{
       RowBox[{"3", "M"}], "+", "5"}]}], "}"}]}], "//", 
   "Simplify"}]}]}], "Code",
 CellChangeTimes->{{3.8750909259659147`*^9, 3.8750909755423727`*^9}, {
   3.875091011817265*^9, 3.875091011971953*^9}, {3.875091761347869*^9, 
   3.8750917685033712`*^9}, {3.875095922978674*^9, 3.875095924237027*^9}, {
   3.875102731687755*^9, 3.875102732249934*^9}, {3.8751027766648827`*^9, 
   3.8751027771588306`*^9}, {3.8751035450282717`*^9, 3.875103545106289*^9}, {
   3.875103575541709*^9, 3.875103575695881*^9}, {3.87544903570872*^9, 
   3.875449051990718*^9}, {3.875459791604109*^9, 3.875459793447543*^9}, {
   3.875460098494342*^9, 3.875460146550344*^9}, 3.8754622749212112`*^9, {
   3.8754623550682497`*^9, 3.875462355922514*^9}, {3.8754628105165787`*^9, 
   3.8754628283143*^9}, {3.875462860582975*^9, 3.875462874171104*^9}, {
   3.8754629739336987`*^9, 3.8754629980591993`*^9}, {3.875464304445507*^9, 
   3.8754643097529507`*^9}, {3.875464374646179*^9, 3.875464379549264*^9}, {
   3.875464667348877*^9, 3.875464667406427*^9}, {3.875464789979189*^9, 
   3.87546480291568*^9}, {3.875464900566935*^9, 3.875464925993936*^9}, {
   3.8754665630601683`*^9, 3.875466563364303*^9}, {3.875603537872497*^9, 
   3.875603540199489*^9}, {3.87560358270327*^9, 3.875603659183642*^9}, {
   3.875604134438981*^9, 3.875604134819813*^9}, 3.875604513449641*^9, {
   3.8757017826018*^9, 3.87570182248929*^9}, {3.876295114711844*^9, 
   3.876295115023265*^9}, {3.876385582696113*^9, 3.876385587334529*^9}, {
   3.8763857218627453`*^9, 3.8763857221737823`*^9}, {3.8763859659634027`*^9, 
   3.8763860077477903`*^9}, {3.876572531010734*^9, 3.876572583060375*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"38", "+", 
    RowBox[{"25", " ", "M"}], "+", 
    RowBox[{"4", " ", 
     SuperscriptBox["M", "2"]}]}], ",", 
   RowBox[{"20", "+", 
    RowBox[{"6", " ", "M"}]}], ",", 
   RowBox[{"2", " ", 
    RowBox[{"(", 
     RowBox[{"13", "+", 
      RowBox[{"7", " ", "M"}], "+", 
      SuperscriptBox["M", "2"]}], ")"}]}]}], "}"}]], "Output",
 CellChangeTimes->{{3.876385715009602*^9, 3.876385722596904*^9}, 
   3.876386008897937*^9, 3.876386751315452*^9, 3.8763888748686457`*^9, 
   3.876545897538424*^9, 3.876546877718638*^9, 3.876547012161355*^9, 
   3.876556256580185*^9, 3.876556718056282*^9, 3.876560112393311*^9, 
   3.8765601719581823`*^9, 3.8765609489446907`*^9, {3.876565089295424*^9, 
   3.876565142063918*^9}, 3.876572584959175*^9, 3.876572676264659*^9, 
   3.876572738584547*^9, 3.8769972614221373`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"46", "+", 
    RowBox[{"29", " ", "M"}], "+", 
    RowBox[{"4", " ", 
     SuperscriptBox["M", "2"]}]}], ",", 
   RowBox[{"20", "+", 
    RowBox[{"6", " ", "M"}]}], ",", 
   RowBox[{"31", "+", 
    RowBox[{"17", " ", "M"}], "+", 
    RowBox[{"2", " ", 
     SuperscriptBox["M", "2"]}]}]}], "}"}]], "Output",
 CellChangeTimes->{{3.876385715009602*^9, 3.876385722596904*^9}, 
   3.876386008897937*^9, 3.876386751315452*^9, 3.8763888748686457`*^9, 
   3.876545897538424*^9, 3.876546877718638*^9, 3.876547012161355*^9, 
   3.876556256580185*^9, 3.876556718056282*^9, 3.876560112393311*^9, 
   3.8765601719581823`*^9, 3.8765609489446907`*^9, {3.876565089295424*^9, 
   3.876565142063918*^9}, 3.876572584959175*^9, 3.876572676264659*^9, 
   3.876572738584547*^9, 3.8769972614236393`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["large number of columns", "Subsection",
 CellChangeTimes->{{3.876385600508483*^9, 3.876385604712728*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"arithComp", " ", "=", " ", 
   RowBox[{"{", 
    RowBox[{"5", ",", "2", ",", "1"}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"costsumcheckUVlarge", " ", "=", " ", 
  RowBox[{
   RowBox[{
    RowBox[{"(", 
     RowBox[{"M", "+", "1"}], ")"}], " ", "*", " ", 
    RowBox[{
     RowBox[{"sumcheck", "[", 
      RowBox[{
       RowBox[{"n", "+", "m"}], ",", " ", "3", ",", " ", "6", ",", " ", 
       "arithComp"}], "]"}], "/", 
     RowBox[{"(", 
      RowBox[{"1", " ", "-", " ", 
       RowBox[{"1", "/", 
        RowBox[{"2", "^", 
         RowBox[{"(", 
          RowBox[{"n", "+", "m"}], ")"}]}]}]}], ")"}]}]}], "//", 
   "Simplify"}]}], "\[IndentingNewLine]", 
 RowBox[{"costUVlarge", " ", "=", " ", 
  RowBox[{
   RowBox[{"costsumcheckUVlarge", "+", " ", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"M", "+", "1"}], ")"}], 
     RowBox[{"{", 
      RowBox[{"7", ",", " ", "0", ",", " ", "4"}], "}"}]}]}], "//", 
   "Simplify"}]}]}], "Code",
 CellChangeTimes->{{3.8750909259659147`*^9, 3.8750909755423727`*^9}, {
   3.875091011817265*^9, 3.875091011971953*^9}, {3.875091761347869*^9, 
   3.8750917685033712`*^9}, {3.875095922978674*^9, 3.875095924237027*^9}, {
   3.875102731687755*^9, 3.875102732249934*^9}, {3.8751027766648827`*^9, 
   3.8751027771588306`*^9}, {3.8751035450282717`*^9, 3.875103545106289*^9}, {
   3.875103575541709*^9, 3.875103575695881*^9}, {3.87544903570872*^9, 
   3.875449051990718*^9}, {3.875459791604109*^9, 3.875459793447543*^9}, {
   3.875460098494342*^9, 3.875460146550344*^9}, 3.8754622749212112`*^9, {
   3.8754623550682497`*^9, 3.875462355922514*^9}, {3.8754628105165787`*^9, 
   3.8754628283143*^9}, {3.875462860582975*^9, 3.875462874171104*^9}, {
   3.8754629739336987`*^9, 3.8754629980591993`*^9}, {3.875464304445507*^9, 
   3.8754643097529507`*^9}, {3.875464374646179*^9, 3.875464379549264*^9}, {
   3.875464667348877*^9, 3.875464667406427*^9}, {3.875464789979189*^9, 
   3.87546480291568*^9}, {3.875464900566935*^9, 3.875464925993936*^9}, {
   3.8754665630601683`*^9, 3.875466563364303*^9}, {3.875603537872497*^9, 
   3.875603540199489*^9}, {3.87560358270327*^9, 3.875603659183642*^9}, {
   3.875604134438981*^9, 3.875604134819813*^9}, 3.875604513449641*^9, {
   3.8757017826018*^9, 3.87570182248929*^9}, {3.876294935922762*^9, 
   3.876294976871594*^9}, {3.8762950096121683`*^9, 3.87629504286515*^9}, {
   3.876295088768346*^9, 3.876295091983234*^9}, {3.876295140653759*^9, 
   3.8762952201241207`*^9}, {3.8762954318300037`*^9, 3.876295438845056*^9}, {
   3.876368111846792*^9, 3.876368114205063*^9}, {3.876385621768506*^9, 
   3.876385712184557*^9}, {3.8763859817375793`*^9, 3.87638599991722*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"38", " ", 
    RowBox[{"(", 
     RowBox[{"1", "+", "M"}], ")"}]}], ",", 
   RowBox[{"20", " ", 
    RowBox[{"(", 
     RowBox[{"1", "+", "M"}], ")"}]}], ",", 
   RowBox[{"26", " ", 
    RowBox[{"(", 
     RowBox[{"1", "+", "M"}], ")"}]}]}], "}"}]], "Output",
 CellChangeTimes->{{3.876385715044609*^9, 3.876385715142146*^9}, 
   3.876386016882557*^9, 3.876386751338686*^9, 3.876388874890017*^9, 
   3.876545897562159*^9, 3.876546877744257*^9, 3.876547012184167*^9, 
   3.87655625660221*^9, 3.876556718078918*^9, 3.8765601124171124`*^9, 
   3.876560171975246*^9, 3.8765609489653463`*^9, {3.876565089312255*^9, 
   3.876565142080762*^9}, 3.87657267628866*^9, 3.87699726147366*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"45", " ", 
    RowBox[{"(", 
     RowBox[{"1", "+", "M"}], ")"}]}], ",", 
   RowBox[{"20", " ", 
    RowBox[{"(", 
     RowBox[{"1", "+", "M"}], ")"}]}], ",", 
   RowBox[{"30", " ", 
    RowBox[{"(", 
     RowBox[{"1", "+", "M"}], ")"}]}]}], "}"}]], "Output",
 CellChangeTimes->{{3.876385715044609*^9, 3.876385715142146*^9}, 
   3.876386016882557*^9, 3.876386751338686*^9, 3.876388874890017*^9, 
   3.876545897562159*^9, 3.876546877744257*^9, 3.876547012184167*^9, 
   3.87655625660221*^9, 3.876556718078918*^9, 3.8765601124171124`*^9, 
   3.876560171975246*^9, 3.8765609489653463`*^9, {3.876565089312255*^9, 
   3.876565142080762*^9}, 3.87657267628866*^9, 3.8769972614750357`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["break even point", "Subsection",
 CellChangeTimes->{{3.876385748469411*^9, 3.876385752247199*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"costWithCommUV", " ", "=", " ", 
  RowBox[{"costUV", " ", "+", " ", 
   RowBox[{"costPippengerPerScalar", " ", "*", " ", 
    RowBox[{"(", 
     RowBox[{"M", "+", "2"}], ")"}]}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"costWithCommUVlarge", " ", "=", " ", 
   RowBox[{"costUVlarge", " ", "+", " ", 
    RowBox[{"costPippengerPerScalar", " ", "*", " ", "2", "*", 
     RowBox[{"(", 
      RowBox[{"M", "+", "1"}], ")"}]}]}]}], "\n", "\n", 
  RowBox[{"(*", " ", 
   RowBox[{
   "break", " ", "even", " ", "point", " ", "with", " ", "linear", " ", 
    "strategy"}], " ", "*)"}]}], "\n", 
 RowBox[{"sol", " ", "=", " ", 
  RowBox[{"NSolve", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"costWithCommUV", "[", 
      RowBox[{"[", "1", "]"}], "]"}], " ", "\[Equal]", " ", 
     RowBox[{"costWithCommUVlarge", "[", 
      RowBox[{"[", "1", "]"}], "]"}]}], ",", "M"}], "]"}]}], "\n", 
 RowBox[{
  RowBox[{"breakEvenUV", " ", "=", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"M", " ", "/.", " ", "sol"}], " ", "//", " ", "Last"}], "//", 
    "Floor"}]}], "\n", "\n", 
  RowBox[{"(*", " ", 
   RowBox[{"maximum", " ", "speedup"}], " ", "*)"}]}], "\n", 
 RowBox[{
  RowBox[{"ratio", " ", "=", " ", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"costWithCommUVlarge", "[", 
        RowBox[{"[", "1", "]"}], "]"}], "/", 
       RowBox[{"costWithCommUV", "[", 
        RowBox[{"[", "1", "]"}], "]"}]}], "/.", 
      RowBox[{"{", 
       RowBox[{"M", "\[Rule]", "x"}], "}"}]}], ",", "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"x", ",", "1", ",", " ", "breakEvenUV"}], "}"}]}], "]"}]}], 
  ";"}], "\n", 
 RowBox[{
  RowBox[{"max", " ", "=", " ", 
   RowBox[{"Max", "@", "ratio"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"max", "//", "N"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Position", "[", 
    RowBox[{"ratio", ",", "max"}], "]"}], "//", "Flatten"}], "\n", "\n", 
  RowBox[{"(*", " ", 
   RowBox[{
   "range", " ", "for", " ", "0.8", " ", "of", " ", "the", " ", "maximum", 
    " ", "speedup"}], " ", "*)"}]}], "\n", 
 RowBox[{"max", " ", "*", " ", "0.80"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Position", "[", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"#", "\[GreaterEqual]", 
         RowBox[{"max", "*", "0.80"}]}], ")"}], "&"}], "/@", "ratio"}], ",", 
     "True"}], "]"}], "//", "Flatten"}], "\n"}], "\n", 
 RowBox[{"ListPlot", "[", 
  RowBox[{"ratio", ",", 
   RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]}], "Code",
 CellChangeTimes->{{3.8763858044317017`*^9, 3.876385921151272*^9}, {
  3.876386031692253*^9, 3.8763861394701977`*^9}, {3.876386871077985*^9, 
  3.8763868856803837`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"46", "+", 
    RowBox[{"29", " ", "M"}], "+", 
    RowBox[{"4", " ", 
     SuperscriptBox["M", "2"]}], "+", 
    RowBox[{"557.9611558503079`", " ", 
     RowBox[{"(", 
      RowBox[{"2", "+", "M"}], ")"}]}]}], ",", 
   RowBox[{"20", "+", 
    RowBox[{"6", " ", "M"}], "+", 
    RowBox[{"557.9611558503079`", " ", 
     RowBox[{"(", 
      RowBox[{"2", "+", "M"}], ")"}]}]}], ",", 
   RowBox[{"31", "+", 
    RowBox[{"17", " ", "M"}], "+", 
    RowBox[{"2", " ", 
     SuperscriptBox["M", "2"]}], "+", 
    RowBox[{"557.9611558503079`", " ", 
     RowBox[{"(", 
      RowBox[{"2", "+", "M"}], ")"}]}]}]}], "}"}]], "Output",
 CellChangeTimes->{{3.876385910047246*^9, 3.876385923439322*^9}, 
   3.8763861470907803`*^9, 3.876386751344759*^9, 3.876386886276246*^9, 
   3.8763888748981743`*^9, 3.876545897568553*^9, 3.876546877751233*^9, 
   3.876547012193437*^9, 3.876556256611607*^9, 3.876556718102483*^9, 
   3.876560112421912*^9, 3.8765601719939823`*^9, 3.876560948972204*^9, {
   3.876565089318849*^9, 3.876565142087329*^9}, 3.8765726762957087`*^9, 
   3.876997261525393*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"1160.9223117006159`", " ", 
    RowBox[{"(", 
     RowBox[{"1", "+", "M"}], ")"}]}], ",", 
   RowBox[{"1135.9223117006159`", " ", 
    RowBox[{"(", 
     RowBox[{"1", "+", "M"}], ")"}]}], ",", 
   RowBox[{"1145.9223117006159`", " ", 
    RowBox[{"(", 
     RowBox[{"1", "+", "M"}], ")"}]}]}], "}"}]], "Output",
 CellChangeTimes->{{3.876385910047246*^9, 3.876385923439322*^9}, 
   3.8763861470907803`*^9, 3.876386751344759*^9, 3.876386886276246*^9, 
   3.8763888748981743`*^9, 3.876545897568553*^9, 3.876546877751233*^9, 
   3.876547012193437*^9, 3.876556256611607*^9, 3.876556718102483*^9, 
   3.876560112421912*^9, 3.8765601719939823`*^9, 3.876560948972204*^9, {
   3.876565089318849*^9, 3.876565142087329*^9}, 3.8765726762957087`*^9, 
   3.876997261526895*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"M", "\[Rule]", "0.0017422993389621553`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"M", "\[Rule]", "143.48854666323803`"}], "}"}]}], "}"}]], "Output",\

 CellChangeTimes->{{3.876385910047246*^9, 3.876385923439322*^9}, 
   3.8763861470907803`*^9, 3.876386751344759*^9, 3.876386886276246*^9, 
   3.8763888748981743`*^9, 3.876545897568553*^9, 3.876546877751233*^9, 
   3.876547012193437*^9, 3.876556256611607*^9, 3.876556718102483*^9, 
   3.876560112421912*^9, 3.8765601719939823`*^9, 3.876560948972204*^9, {
   3.876565089318849*^9, 3.876565142087329*^9}, 3.8765726762957087`*^9, 
   3.876997261528151*^9}],

Cell[BoxData["143"], "Output",
 CellChangeTimes->{{3.876385910047246*^9, 3.876385923439322*^9}, 
   3.8763861470907803`*^9, 3.876386751344759*^9, 3.876386886276246*^9, 
   3.8763888748981743`*^9, 3.876545897568553*^9, 3.876546877751233*^9, 
   3.876547012193437*^9, 3.876556256611607*^9, 3.876556718102483*^9, 
   3.876560112421912*^9, 3.8765601719939823`*^9, 3.876560948972204*^9, {
   3.876565089318849*^9, 3.876565142087329*^9}, 3.8765726762957087`*^9, 
   3.876997261529421*^9}],

Cell[BoxData["1.7193552968081192`"], "Output",
 CellChangeTimes->{{3.876385910047246*^9, 3.876385923439322*^9}, 
   3.8763861470907803`*^9, 3.876386751344759*^9, 3.876386886276246*^9, 
   3.8763888748981743`*^9, 3.876545897568553*^9, 3.876546877751233*^9, 
   3.876547012193437*^9, 3.876556256611607*^9, 3.876556718102483*^9, 
   3.876560112421912*^9, 3.8765601719939823`*^9, 3.876560948972204*^9, {
   3.876565089318849*^9, 3.876565142087329*^9}, 3.8765726762957087`*^9, 
   3.876997261530648*^9}],

Cell[BoxData[
 RowBox[{"{", "11", "}"}]], "Output",
 CellChangeTimes->{{3.876385910047246*^9, 3.876385923439322*^9}, 
   3.8763861470907803`*^9, 3.876386751344759*^9, 3.876386886276246*^9, 
   3.8763888748981743`*^9, 3.876545897568553*^9, 3.876546877751233*^9, 
   3.876547012193437*^9, 3.876556256611607*^9, 3.876556718102483*^9, 
   3.876560112421912*^9, 3.8765601719939823`*^9, 3.876560948972204*^9, {
   3.876565089318849*^9, 3.876565142087329*^9}, 3.8765726762957087`*^9, 
   3.876997261531743*^9}],

Cell[BoxData["1.3754842374464955`"], "Output",
 CellChangeTimes->{{3.876385910047246*^9, 3.876385923439322*^9}, 
   3.8763861470907803`*^9, 3.876386751344759*^9, 3.876386886276246*^9, 
   3.8763888748981743`*^9, 3.876545897568553*^9, 3.876546877751233*^9, 
   3.876547012193437*^9, 3.876556256611607*^9, 3.876556718102483*^9, 
   3.876560112421912*^9, 3.8765601719939823`*^9, 3.876560948972204*^9, {
   3.876565089318849*^9, 3.876565142087329*^9}, 3.8765726762957087`*^9, 
   3.876997261532824*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "2", ",", "3", ",", "4", ",", "5", ",", "6", ",", "7", ",", "8", ",", "9", 
   ",", "10", ",", "11", ",", "12", ",", "13", ",", "14", ",", "15", ",", 
   "16", ",", "17", ",", "18", ",", "19", ",", "20", ",", "21", ",", "22", 
   ",", "23", ",", "24", ",", "25", ",", "26", ",", "27", ",", "28", ",", 
   "29", ",", "30", ",", "31", ",", "32", ",", "33", ",", "34", ",", "35", 
   ",", "36", ",", "37", ",", "38", ",", "39", ",", "40", ",", "41", ",", 
   "42", ",", "43", ",", "44", ",", "45", ",", "46", ",", "47", ",", "48", 
   ",", "49", ",", "50", ",", "51", ",", "52", ",", "53", ",", "54", ",", 
   "55", ",", "56", ",", "57", ",", "58", ",", "59", ",", "60", ",", "61", 
   ",", "62", ",", "63"}], "}"}]], "Output",
 CellChangeTimes->{{3.876385910047246*^9, 3.876385923439322*^9}, 
   3.8763861470907803`*^9, 3.876386751344759*^9, 3.876386886276246*^9, 
   3.8763888748981743`*^9, 3.876545897568553*^9, 3.876546877751233*^9, 
   3.876547012193437*^9, 3.876556256611607*^9, 3.876556718102483*^9, 
   3.876560112421912*^9, 3.8765601719939823`*^9, 3.876560948972204*^9, {
   3.876565089318849*^9, 3.876565142087329*^9}, 3.8765726762957087`*^9, 
   3.876997261534243*^9}],

Cell[BoxData[
 GraphicsBox[{{}, 
   {RGBColor[0.24720000000000014`, 0.24, 0.6], PointBox[CompressedData["
1:eJw11A0wFHgcxvG96pSOThciTlLnvOW4VN6yD7FSaK21rMXuapKXnLPJpLvY
E1LX9XZN0QtbKl16D9PpKpGcLdxsqQuVHJ2XXkTyVrm7uf/PzI757DPf+c3u
7O6c6CTB6gkcDmfnv4///v//18cNq1Y8VzoOctkT0C30OKq6MMQ8Ba3N1Xkp
piPMemjUzubdf0fWR6w6u0C8Y5TZCM4n3dLqX5JN8diV51dlP8ZsjotHnGys
wsgWmHFo1aagZPI89LQdTk5JJ1si5/nsSxszyVYYMPpm3oYMsg1MElNm96eS
7WDc9oGjE0+2h6S6412RmOwAXnavKIZHdkTTllnbmh3IX6Py8afpa43JC1Bk
3TlrHYfshDKtlukLn9HrW4g6w+yY52ryItQWqQzqzpIXg/NJv9HAbrIzUnMn
v/kxhewCH7cbXedEZFcYbz7dnepMdkOkyjZ80JjsDlGY+Un3MXr/l2C4SF0r
aSV7oNlHlie9SubikipgUXABGXC9+MwDGcxKwMhEMeUrGe2eyInU/GHpSbsn
vHT2dsy3oN0Liq7a634TaffCiO8t64yOYbYvRdXUSn9NDbNyKXr+vFfHK6bd
Gx+il6qf5NLuDevzdtNU8bT7oOuHgJGt/rT7oHTx+OFCe9p5OKiovtWuRzsP
Rt57PggH6PPqC5fVxnlDTcxKX3iHzyi9fZn2ZYDdA57mIO3LUAj3z3UyaPeD
Ir+7LS2adj/c9ZpnasqjfTkSq+cvGbCmfTlyXiN/XJf2FRhLd0ny7H/L9hXY
punrvnyfmeOPCYFii5grtPuD85d0jp+K9gDoqg88lWfTHoACv9v88/G0B2Ko
bHGQE5/2QFjZtFe8WEj7StzoLoxoNaF9JR490n8//hHtfNx02icJ62Lfd/DR
Uv/WobOBWclHh5rrcKKMuZKPTHHo5IOH6PchCAu37ZdWb6Y+CIubRI/NEqgP
Quy1VV+cEVAfhDttKXcTXakXQDdnR5bMgnoBHjZEabKmUi+ASchuYXP/G9YL
kG1kckTWwswJhta3Gf4zbzIjGMeKp/W8P82sDIb1UzMd/X3UB6P4GtdbnEG9
EG38uLn1sdQL0V9u6ZokoF6IvZrfjLju1AthHR+w0t2S+hCc0B5KjNGjPgRp
va8nXR0dYH0ImmxPHeN2MleGID/zwaMXjcwcEWQVezxrK5ghwgHDOPnt49SL
8HNTYuPILupFUA/5mQm+pz4Uz5KXP727hvpQpPeuPZ0eTH0oXlkn6Au41Ifi
4NKaswG21IfhfJl6evJM6sNglXL03pWJ1Ifhb23Bzvl9/awPw9XW6uaaVmaO
GOeK7/Cz6pghhou7w/bYcmalGF35x60URdSLIYxvrlLtoj4c216uHx7YRH04
fokziU1IoD4c5TPWXtASUx+OkoIp62t41EswKtojLHaiXgLDrpL+EgvqJai8
M9Kv0aNegjiFxajJ+GvWR8DRsOFMzgtmRKBWfqL0s1ZmZQRU9SHl19XMlRHo
qc+R/PQr9ZHoamgP+O4k9ZFw/Hi689b91Eei2aCspjyH+kisO5BYNDGV+iiY
Z2nSk1ZTHwUH+TqtUSH1UdDKbeo+6k19FNxsFA/jnKiXwldQmhM4j3opnBve
FvH1qZdis3ZmR9Ik6qUYfLD1XcmbPtbLsH+B7ZFJnczmMti/istNa2KGDINR
9nFTbzHLZBgb9horL2dWypD0MmgovZhZJYPJ9XbD6DzmShk8fg/Vl21lbpNB
1dh7asNGui+HwfZz+0oS6L4cunZb9g5H0H05FrV8GS0NpPtyNPoaXHvCpfty
hFsVrklzpPtySKsem9nNpfty1CpKjw3p0305ljnXVLRo9XH/AR792YE=
     "]]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesLabel->{None, None},
  AxesOrigin->{1., 1.001674394199215},
  Method->{},
  PlotRangeClipping->True]], "Output",
 CellChangeTimes->{{3.876385910047246*^9, 3.876385923439322*^9}, 
   3.8763861470907803`*^9, 3.876386751344759*^9, 3.876386886276246*^9, 
   3.8763888748981743`*^9, 3.876545897568553*^9, 3.876546877751233*^9, 
   3.876547012193437*^9, 3.876556256611607*^9, 3.876556718102483*^9, 
   3.876560112421912*^9, 3.8765601719939823`*^9, 3.876560948972204*^9, {
   3.876565089318849*^9, 3.876565142087329*^9}, 3.8765726762957087`*^9, 
   3.876997261575238*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["best strategy", "Subsection",
 CellChangeTimes->{{3.8763865017050533`*^9, 3.876386504201005*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"bestCostWithCommUV", " ", "=", " ", 
   RowBox[{"Join", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"costWithCommUV", "[", 
         RowBox[{"[", "1", "]"}], "]"}], "/.", 
        RowBox[{"{", 
         RowBox[{"M", "\[Rule]", "x"}], "}"}]}], ",", 
       RowBox[{"{", 
        RowBox[{"x", ",", "1", ",", "breakEvenUV"}], "}"}]}], "]"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"costWithCommUVlarge", "[", 
         RowBox[{"[", "1", "]"}], "]"}], "/.", 
        RowBox[{"{", 
         RowBox[{"M", "\[Rule]", "x"}], "}"}]}], ",", 
       RowBox[{"{", 
        RowBox[{"x", ",", 
         RowBox[{"breakEvenUV", "+", "1"}], ",", "maxM"}], "}"}]}], "]"}]}], 
    "\[IndentingNewLine]", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"ListPlot", "[", "bestCostWithCommUV", "]"}]}], "Code",
 CellChangeTimes->{{3.876386459981937*^9, 3.876386492348509*^9}}],

Cell[BoxData[
 GraphicsBox[{{}, 
   {RGBColor[0.24720000000000014`, 0.24, 0.6], PointBox[CompressedData["
1:eJw11gtYTekaB/BlJE1yaQzJ7WwjMjTGEZoo/VOSS+l+v+z7DhNHGMZxshg6
McZlhjEux9njmkmkRy6Z2CNCckvIZdgujUIkEoqO862/nqen59e7//tda33f
u9bqpZ0WZvhIkqRT73///1f5qfUZeOhAdqZuA/gPzPzn3LTc5M20HUK7Oq/Y
2GMH3QHi84+y6U9h6zZAX6DPpbso+Rt5dHfcLTktb4/Lp1XK9905QH8Gy6SU
//z4jwLaRfn+1oV0X2y0a12Qvv0o3U/p1/cPuj8KL4VNGZp+jHZT+l8rogdi
QO/N3w7wPEEPUo5nYzH9d6xPe/bvXjan6MHK8U07TbvD7hh+drpdQg9Rjjes
lB6K2Y4rt7YtOUsPU45/9HnaA2m5064t0F6gv0LwopaJPi4XaU/0j1lrbXrw
wcOV88spo0fgTvORB7NnXaK9lPP1Kae98cuOv57U2V+mRyrnX/HBPgiZ2O5l
atYVGsr1mHtVsQy0bhj6tiqognVf5fr0vsa6L45sSrTRN9LSKOV6Xb7O+ijM
Cljc5nbuDdb9lOu3/Cbrfvjiya5P4qb+ybq/cj1DbrHuj/tryp0vu99mfTSs
Q/YFn3awsj4a53qmd8iSaSkAv9sFlmXUf6gHYGed42rDlDusj8Hamzci/e/S
8hgsLt7m1DvuLuuBYj1aXKLlQKjXe26wTrjH+lixPkdP0vJYeE0923OT333W
x4n1mmeh5XHoMkqzOX5kJevjxfoNL6Tl8XjRqb6Ps/dfrE8Q69lQSMsTcL46
c+cVnwesB4n1zT9Gy0HILuz+xeqAKtaDxXqnnaHlYGSsys0NDa1mfaJY/0EV
tDwRWoP/kPbqh6yHiP1QU0UjBN6eFQdK0x7x8yFif2S/pS0hcG779YilSx4z
Hyr2S0rnGuZDUW9tPjJmKy2Hiv3T1/0J86G4sO+nUa2KaClM7Kf74U+ZD8Ou
TNfiY5W0HIagQ6+KQrxqmQ+D/+4/XZBBS+EYvuXY4j0XaIRj8PsLsqnbM+bD
0e+HZWPamWhLOFQLp2e9zaOlCHSeHfVxajONCLw/vckJE+qYj4CNRnXmzDra
EoHGyFZuhx/QUiTqxj1c1n/Yc+YjUe1zvsZ5MS1Hiv27pJy2ROLK5+v2zHF5
wXyU2M9VM2lE4XhH3fSrJ2g5SuzvAKd65qOQ99bNfdgkWooW+33bYRrRMD94
+WJ1u5fMR4v9b6OlLdFYftGy/3k+LcWIedDZNzAfg3mHl84JS6blGDEff+yj
LTGYvC1i+F77V8zHinlRaWjEImZFz6b2B2k5VszP/PavmY/F6G+rCqeaaClO
zNOtozTi4K7Lm3+2yxvm48R8eaXRljj0Cprn61ZKS/Fi3jb0bWQ+Hu08Alp+
v4CW48X8vblJW+LRpOpwovqrJuYTxDzGrqGRgIf21zMC62g5QczngYlvmU9A
xYstgTtyaClRzGtnh3fMJ6L4Vqq97RRaThTzO6uEtiRi3ymPUn3/ZuaTxDyX
L6WRhF/zWiwvekTLSSgNnzF7XV/JV8knIfBcpeuqaFpKRlFgTMWpTBrJGFlU
kokCWk5Ggbe356tHH/LJGHZwT/XTHi2YV2Pv4M/Wu4bQKjXcclaP++9CGmpk
ubZujMqn1Wq4vL8dBFfRshqbuj+Oz+j2kWKzGl3XJjm8C6YtavzsePH3/IW0
VQ3HZX6pOftpSYNltvt7VD6kVRrYLeh3LuFvLdlfg0WN69M/jaDVGjTPavul
wxJa1mBu7fzbfkdoswYvJ9etOPyctmgwvVKPGZ/bsL8Gj5Ov1k5OpiUtrh3d
fSd1Ga3SivtTfi4NLU7unLOlqZxWa8X9yv81LWuR/9Mow/c9WrG/Vty/ynxp
ixZb/uXg2tVIW7XiealZSks6rDJdqcraTat0yEw4PLa2jIYO6aHm3zwaaLVO
PF/nd7Nlfx3e366+LvahzTroBk0821ZPW3SI6zNkYGQmbdUpz+NdtKRHYPt3
T+9doFV6+NjcCxlQT0OPYa9P7k1zbs3+evH8LvCmZT1c7q2a0UJLm/XoWvFN
eWAGbdHD8Wz80JW/0Va9eN5fPUdLBjTv79PQ8zmtMuBltn2M0cmO/Q2oMT89
mDOCVhvE+0F9Mi0bcGPpoblei2izAWXzN934Lou2GHB65ndeZ0ppq0G8r33y
jJaM2J8U9C6208fsb0RO+ODkXz1pGLE10MlSlUirjdjg3aQatJCWjfhx8J0F
s7fTZiOWuBbfPVJCW4yQu2f72T6lrUbxfhjU0Z79TZhqO6vVGg9aZYKhMdZ4
M56GCfG1I0/2lmm1CWGVvftN2UrLJoy9brck7xRtNgHna6pfP6YtJngcLxvn
69iG/U3K++hQWkpBn90bHS7E0qoUdN+yINUpnUYKOv5iPJe0mVanwP6H8V9u
L6blFLRYOGhlzUPanIKGbzo9G9Legf1T8GTKm9B57rQ1BZXq23lF0bQ0CTcj
j3dsM8/B93+y061V
     "]]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  Method->{},
  PlotRange->{{0, 200.}, {0, 233345.38465182378`}},
  PlotRangeClipping->True,
  PlotRangePadding->{{4., 4.}, {4666.907693036475, 
   4666.907693036475}}]], "Output",
 CellChangeTimes->{{3.876386483448996*^9, 3.876386492962757*^9}, 
   3.876386751374456*^9, 3.87638887493642*^9, 3.876545897625967*^9, 
   3.876546877811241*^9, 3.876547012247064*^9, 3.876556256663608*^9, 
   3.876556718134581*^9, 3.8765601124712*^9, 3.876560172013389*^9, 
   3.876560949011223*^9, {3.876565089351614*^9, 3.876565142117313*^9}, 
   3.8765726763488083`*^9, 3.876997261633602*^9}]
}, Open  ]]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Comparison", "Section",
 CellChangeTimes->{{3.87638651335948*^9, 3.876386514714394*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"Show", "[", 
  RowBox[{"{", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ListPlot", "[", "bestCostWithCommUV", "]"}], ",", 
    "\[IndentingNewLine]", 
    RowBox[{"ListPlot", "[", "bestCostWithComm", "]"}]}], 
   "\[IndentingNewLine]", "}"}], "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ratioWithComm", "=", " ", 
   RowBox[{"bestCostWithCommUV", "/", "bestCostWithComm"}]}], ";"}], "\n", 
 RowBox[{"ListPlot", "[", "\[IndentingNewLine]", "\t", "ratioWithComm", 
  "\[IndentingNewLine]", "]"}], "\n", 
 RowBox[{
  RowBox[{"minRatio", " ", "=", " ", 
   RowBox[{"ratioWithComm", "//", "Min"}]}], ";"}], "\n", 
 RowBox[{"minRatio", "//", "N"}], "\n", 
 RowBox[{"Position", "[", 
  RowBox[{"ratioWithComm", ",", " ", "minRatio"}], "]"}], "\n", 
 RowBox[{
  RowBox[{"maxRatio", " ", "=", " ", 
   RowBox[{"ratioWithComm", "//", "Max"}]}], ";"}], "\n", 
 RowBox[{"maxRatio", "//", "N"}], "\n", 
 RowBox[{"Position", "[", 
  RowBox[{"ratioWithComm", ",", " ", "maxRatio"}], "]"}], "\n", 
 RowBox[{"ratioWithComm", "//", "Last"}], "\n", 
 RowBox[{
  RowBox[{"Position", "[", 
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{"#", " ", "\[GreaterEqual]", " ", "2"}], ")"}], "&"}], "/@", 
     "ratioWithComm"}], ",", " ", "True"}], "]"}], "//", "Flatten"}], "\n", 
 RowBox[{
  RowBox[{"Position", "[", 
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{"#", " ", "\[GreaterEqual]", " ", "3"}], ")"}], "&"}], "/@", 
     "ratioWithComm"}], ",", " ", "True"}], "]"}], "//", "Flatten"}]}], "Code",\

 CellChangeTimes->{{3.876386528957674*^9, 3.876386614328298*^9}, {
  3.8763927745290747`*^9, 3.8763928021296797`*^9}, {3.876547109386219*^9, 
  3.876547280034383*^9}, {3.8765607172293777`*^9, 3.876560792368044*^9}, {
  3.876564969012892*^9, 3.876564977905294*^9}}],

Cell[BoxData[
 GraphicsBox[{{{}, 
    {RGBColor[0.24720000000000014`, 0.24, 0.6], PointBox[CompressedData["
1:eJw11gtYTekaB/BlJE1yaQzJ7WwjMjTGEZoo/VOSS+l+v+z7DhNHGMZxshg6
McZlhjEux9njmkmkRy6Z2CNCckvIZdgujUIkEoqO862/nqen59e7//tda33f
u9bqpZ0WZvhIkqRT73///1f5qfUZeOhAdqZuA/gPzPzn3LTc5M20HUK7Oq/Y
2GMH3QHi84+y6U9h6zZAX6DPpbso+Rt5dHfcLTktb4/Lp1XK9905QH8Gy6SU
//z4jwLaRfn+1oV0X2y0a12Qvv0o3U/p1/cPuj8KL4VNGZp+jHZT+l8rogdi
QO/N3w7wPEEPUo5nYzH9d6xPe/bvXjan6MHK8U07TbvD7hh+drpdQg9Rjjes
lB6K2Y4rt7YtOUsPU45/9HnaA2m5064t0F6gv0LwopaJPi4XaU/0j1lrbXrw
wcOV88spo0fgTvORB7NnXaK9lPP1Kae98cuOv57U2V+mRyrnX/HBPgiZ2O5l
atYVGsr1mHtVsQy0bhj6tiqognVf5fr0vsa6L45sSrTRN9LSKOV6Xb7O+ijM
Cljc5nbuDdb9lOu3/Cbrfvjiya5P4qb+ybq/cj1DbrHuj/tryp0vu99mfTSs
Q/YFn3awsj4a53qmd8iSaSkAv9sFlmXUf6gHYGed42rDlDusj8Hamzci/e/S
8hgsLt7m1DvuLuuBYj1aXKLlQKjXe26wTrjH+lixPkdP0vJYeE0923OT333W
x4n1mmeh5XHoMkqzOX5kJevjxfoNL6Tl8XjRqb6Ps/dfrE8Q69lQSMsTcL46
c+cVnwesB4n1zT9Gy0HILuz+xeqAKtaDxXqnnaHlYGSsys0NDa1mfaJY/0EV
tDwRWoP/kPbqh6yHiP1QU0UjBN6eFQdK0x7x8yFif2S/pS0hcG779YilSx4z
Hyr2S0rnGuZDUW9tPjJmKy2Hiv3T1/0J86G4sO+nUa2KaClM7Kf74U+ZD8Ou
TNfiY5W0HIagQ6+KQrxqmQ+D/+4/XZBBS+EYvuXY4j0XaIRj8PsLsqnbM+bD
0e+HZWPamWhLOFQLp2e9zaOlCHSeHfVxajONCLw/vckJE+qYj4CNRnXmzDra
EoHGyFZuhx/QUiTqxj1c1n/Yc+YjUe1zvsZ5MS1Hiv27pJy2ROLK5+v2zHF5
wXyU2M9VM2lE4XhH3fSrJ2g5SuzvAKd65qOQ99bNfdgkWooW+33bYRrRMD94
+WJ1u5fMR4v9b6OlLdFYftGy/3k+LcWIedDZNzAfg3mHl84JS6blGDEff+yj
LTGYvC1i+F77V8zHinlRaWjEImZFz6b2B2k5VszP/PavmY/F6G+rCqeaaClO
zNOtozTi4K7Lm3+2yxvm48R8eaXRljj0Cprn61ZKS/Fi3jb0bWQ+Hu08Alp+
v4CW48X8vblJW+LRpOpwovqrJuYTxDzGrqGRgIf21zMC62g5QczngYlvmU9A
xYstgTtyaClRzGtnh3fMJ6L4Vqq97RRaThTzO6uEtiRi3ymPUn3/ZuaTxDyX
L6WRhF/zWiwvekTLSSgNnzF7XV/JV8knIfBcpeuqaFpKRlFgTMWpTBrJGFlU
kokCWk5Ggbe356tHH/LJGHZwT/XTHi2YV2Pv4M/Wu4bQKjXcclaP++9CGmpk
ubZujMqn1Wq4vL8dBFfRshqbuj+Oz+j2kWKzGl3XJjm8C6YtavzsePH3/IW0
VQ3HZX6pOftpSYNltvt7VD6kVRrYLeh3LuFvLdlfg0WN69M/jaDVGjTPavul
wxJa1mBu7fzbfkdoswYvJ9etOPyctmgwvVKPGZ/bsL8Gj5Ov1k5OpiUtrh3d
fSd1Ga3SivtTfi4NLU7unLOlqZxWa8X9yv81LWuR/9Mow/c9WrG/Vty/ynxp
ixZb/uXg2tVIW7XiealZSks6rDJdqcraTat0yEw4PLa2jIYO6aHm3zwaaLVO
PF/nd7Nlfx3e366+LvahzTroBk0821ZPW3SI6zNkYGQmbdUpz+NdtKRHYPt3
T+9doFV6+NjcCxlQT0OPYa9P7k1zbs3+evH8LvCmZT1c7q2a0UJLm/XoWvFN
eWAGbdHD8Wz80JW/0Va9eN5fPUdLBjTv79PQ8zmtMuBltn2M0cmO/Q2oMT89
mDOCVhvE+0F9Mi0bcGPpoblei2izAWXzN934Lou2GHB65ndeZ0ppq0G8r33y
jJaM2J8U9C6208fsb0RO+ODkXz1pGLE10MlSlUirjdjg3aQatJCWjfhx8J0F
s7fTZiOWuBbfPVJCW4yQu2f72T6lrUbxfhjU0Z79TZhqO6vVGg9aZYKhMdZ4
M56GCfG1I0/2lmm1CWGVvftN2UrLJoy9brck7xRtNgHna6pfP6YtJngcLxvn
69iG/U3K++hQWkpBn90bHS7E0qoUdN+yINUpnUYKOv5iPJe0mVanwP6H8V9u
L6blFLRYOGhlzUPanIKGbzo9G9Legf1T8GTKm9B57rQ1BZXq23lF0bQ0CTcj
j3dsM8/B93+y061V
      "]]}, {}}, {{}, 
    {RGBColor[0.24720000000000014`, 0.24, 0.6], PointBox[CompressedData["
1:eJw11XlYzAkcx/GfY7fWscQuoYdZRzm2hNzSdwkRm0pW5/zmHsKu1u3Bz1rn
PitJjtUxhDZhY1VuI4lcWTlShilC0uGaNmNqre+nnqen59W7z3xn5o/mG+WP
AZqmgiBc/vj9/0/+qvGav2xpVLrXFsIv6JOLG23PToqF27J/2wp/hR4HO2K/
DXZie+2AZeyrO+Hu7JW74J7s0ATYmS1PgnvTSU/P4YUXDXBfdvVu+Fu2azLs
xl65F3Znl+2DB7DlKfBAduWf8CB29AHYgz3uIDyY3eYwPAT7v+ChZEpMdHpZ
nQ4PY78+Ag9nNxyFR7A7HINHsodmwKPYykzYk709Cx7Nvnsc9mJ3OwkTO+oU
W4LzT6N/xx5yFh1OOYc+hub6+X25MMaIzrYEnkcfy71LNjrb8gwWvLmfuIDO
tkTnoI/jHnkRnW3xzUUfz33AJXS2pctl9AncW+ahsy0NsODDve4KOttSexV9
InfrNXS2pdkN9EncHfLRJ1FWioud1eEmuu8nn0qAJba13z/ok7mfhiW21e8W
+hTuZbDEti4vQP+eu+NtdLY1Axb8uAfdQWdba2FhKvf4u2xiW8few9+jV8JG
9J2F2Ptz97mPPdtaB0voh4qwR1cVYx/A3ekB9mzrPVhCjzNhjx70EPtAvP5H
2LOtJlgKpAcFsatNAWbs2bYsWJj2yeFdS7Bn29bAEttUBRvRZ5RiH8T7HJjY
NvfH2LNNibARvfUT7KfzfgVMbFs1LLFNyjLs0e/Cwg+8n/wUe7YtG5bYphHP
sEf/GxZm8N7tOfZsWyossU3O5dij74WFYN73fIF9MDosBWNfgT16KiyE4P5L
7EPw/GApBM+/Env0bFgIxeuvwj4U7w8sheL9q8YevRoWwqiqPO+8m1cN9uz6
tbDEnpsPG8OopOHss0VfvMI+nArfJvt07A9TOOWXr0/NnAZL4ZT7cE6L6Uth
YzidKQiIfJfUuI+gY5eHXtt6sXEfQWlnnFw9Khr3EbT7aJNNBW1fYx9BO1Ke
VkUNgQU5Rcdf9WsXBpOc1sakpx9ZBUtyWr42zsE/pXEv//T5WXOtcS9S5E/y
gujXsEwkpcbbo7/jGzyeSMEhfeJueMKiSFM//vubo4IlkSZ4v5neegNsEOnj
p13WwcOwUSQPtzOOk2/DZpH69dizpKIOFhTU3XFd0cZub3FfQZ1azx7ZdxxM
Cmrb1D8+bxYsKsiudrBNvxmWFFRf0TnCPhM2KOidueFsSjFsVNDLO0+6TWjy
DvcVVHolT3rqDAtKun/ucMmaybBMSTePxY7pFQWTki6lLk7O2Q6LSjqbGN5c
fQaWlJQRO0bT7DFsUNLB9S65e+wtuK+k5OWtXMa4wWYl/RH1al1JICyoKEZ3
97m0BJapaH3YqYmyJJhUtMLfcOBcDiyqaMH4NS3lL2BJRbNHzppd36YW91Wk
cve7njAYNqoopJeHm2cobFaRf+dO0Q8kWFCTT5v66mX7YZmavJo/ntrlGkxq
GlJ36cjJV7CoJteqg+1COv6L+2rq+Tjm57pRsEFNnQsX3t6hhI1qcrgeOnjY
etisJvts2nbvECxoqCGzV+3CAlimIUtaixkd6mDSUKWh+nhG1zrc19CTuNud
grxhSUPFG08sfTsTNmjo1srE4tho2KihvPmrRw3KgM0aMs7UJ9wqggUtZUZM
qZ8nvMd9LR0KHCh3cIZJS3t9OhrTfWFRS7s8P8imzoMlLW0ZWLKqehts0NIG
l9zSTadho5Ykp7SxbqWwWUuLHDbvvW5nxX0dzf18wWdzXGGZjjTWYG2rQJh0
FFoz+lLaYljUUUBZj96+ibCko4lF9hteXIANOqL8yvIN5bBRR0Nzbk3q0+YD
7uvI7URW2mUPWNBTr8PxrfQhsExPTsmr5thJMOmp/Q7tjf37YFFPLX737T/+
Kizpqckv7pvLamCDnmoXfv3q1w423NdTVeR7/56jYLOeysRHRy8oYGEmPQjK
aa9aZ6P/AEgDEbA=
      "]]}, {}}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  Method->{},
  PlotRange->{{0, 200.}, {0, 233345.38465182378`}},
  PlotRangeClipping->True,
  PlotRangePadding->{{4., 4.}, {4666.907693036475, 
   4666.907693036475}}]], "Output",
 CellChangeTimes->{{3.876386544532989*^9, 3.876386615160928*^9}, 
   3.8763867513889647`*^9, 3.8763888749554234`*^9, {3.876392779359069*^9, 
   3.876392802703229*^9}, 3.876545897644453*^9, 3.8765468778294477`*^9, 
   3.876547012270105*^9, {3.876547218452455*^9, 3.87654728063571*^9}, 
   3.876556256678741*^9, 3.876556718152462*^9, 3.876560112492235*^9, 
   3.876560172024962*^9, {3.876560752113492*^9, 3.876560793047056*^9}, 
   3.8765609490264387`*^9, 3.876564978348179*^9, {3.876565089364338*^9, 
   3.876565142129472*^9}, 3.876572676364606*^9, 3.8769972626654167`*^9}],

Cell[BoxData[
 GraphicsBox[{{}, 
   {RGBColor[0.24720000000000014`, 0.24, 0.6], PointBox[CompressedData["
1:eJw9lAs41PkehyeXRJEOih6ksso1l6aVxMed0DDuDGbG3M4m1XY73Ri126oI
nVa6KBueIrLaZI92zew5anOwUdvGqXA6td03ts4+bunsefb3/XueeTyvz7zP
O37/3zMLpRuEch0ej3fj99f/f//xM+z/1sSpJ1931J/9AeYX/xJd5vqe8Qyk
VZZb9Pfp4A82RfYC22CnNn3G5vD33ep6IHAGY0sUyb/SbFMYMbaGeuVEv4Xb
LMZ2uF14unZWpTHjRXD+sC88v92EsT2Cv3xzJKR6NmMH2Ny/amhtb8p4KSbq
X2auExE7Qej48PK7XGIXrLgY5yIqJnZDbrh2rKaE2B1+L/yelX5K7AFNpGXV
lo3EngiJacuLiCP2wprh6pxvlhEvh4+bctTLgJiPyHV6Rc3N9HlX4OZi77aY
JOIP8aBY+O/et/T/eSNwq15UQTHxSrxd5GK1fymxD5aMZPrXael8VqEv4/mC
2mRiXxjol7WZjdB5rkZ2s05v+kFiP9zxP97z0J7YH2Oe5/nPtTMZA52a7Hd1
GYzVQKmxc/naKXpeAXjXxE9+c4axOgCPM99uuxhEeyCOCqq/K3pmyPZAaEUN
k+WljHlB6B59Y3h7Fe1BOGyfe9/nKd2PYEyI9QK7yhirg9E+568H1GG0h2Bn
mJOxYNyA7SE4WhVW5drImBeK9PINk6YK2kMhX281/NKW9jB4G/1w9XL/dLaH
YUJ5c0VaGWNeOFxFKydvxdMejrAy7yITC9ojcMaxMZt3l91vdQSak2X8kyfo
vq9BdX7Dhd5M2tegQ5jYWeZAeyR8s9ffevCLHtsjobHNizjxNWNeFEIWNBlq
99Eehf3vbAOCYmiPxlj091stbWmPRoVP1GfLX+myfS0i7TSNhd8yVq9F0tma
fXOLaRfgyOOMke8ltAuQ7WaIY3zaYzDcrGOyy4gxYrBFx1rw0ZAOe38MHEMP
rxG1MNbGwMnyWnl4MWNeLMK9ZixcomKMWMjrZs8bDiA/FoWLD42etCY/Fubz
qw4uHp3GfCE8+3bczfuRMYQo1XzpXNfEWC2EqOi/ZqeKGWuFcNx7IlWQQ34c
Nl1apd8WTX4cZvluFj50JT8O9y5UtFwxIT8OkfmNA/xhHvPjsW6DdXXqLcaI
h+k3PeZ2zYzV8bj31Dc7v5yxNh6ac5v/tHM3+QlwGE85oS8hPwHefs80C0LJ
T4B5315FuzP5CThrlnDs9RzyEyGt+OK9zwj7fkUioh5mDB3VMFYn4qZszuXH
hxlrE/Em8nWFQyZ9HydhjknFVIw7+Ulo3tS0KUOH/CQMP4rvDL8zxfwklFhN
WBjVMuYlw6+PX3JqD2Mk4/q0yrpxIWN1Mi5cHdRf4Eh+Ml50tm3U45GfAu9l
yT/X3H3H/BT8OJlSONXIWJ2CsUupTYYHGGtTkDM7oPPvUsa8VBw/8+3WD1aT
n4rFvnc2ulqSn4ro+l1nfnozyfxUmDavnm7Tw5iXhqihsY7JesZIwwvBObeP
DzJWpyFn1LFsh4r8NFTXeOXNDiNfhEovj0/cHcgXQZAdlNqvT74IqmlBPdN/
nmC+CH7eJXv5g4x56ajJ9VHYjDNGOkpTD1VPtyA/HW0r2pePLKN+OpYs3WTy
YA31M1BzV+74Tzn1M1B4pHXP12ryM/A8r0P3/EnyM5A86ny5/Ar5mdg/1fLZ
wV7yM3H6pz3b97wkPxOa1mu7NxnQ+WfC7IJxqWIRnb8YI0kGF0V0/nZinLMv
645LpuchxrGA7meRmxmLxci+/2oy+DA9HzGC5PJp/rWMK8X46knnyMp26omx
d9vDG/xBxkNiSGpL8z3HqS/Bduenc90t2H2yk0An1eoTN3e6jxIE/rqr3TWS
sViCoe6C264Kup8S2Db017vlM66UwMpkaZT7KbqvEhSYRjV4tjAekiBWvqqL
f4vurxS2SYKbskHqS6HbNfNWOjGk0ElZ/CqRWCyFwaF5djHEaim2xK//KIK4
Ugqb4frOQGKtFJuvSYJ8iYekuB49v5dPzMvC4XHllmVcPwvt1oMujlw/C/O6
544v4vpZcLpz/p4118/CkeMevXO5fhb+oRfzL1Oun4WWj2t+M+L6WZBU3PlA
n+vLMHRgv/L9APVliN0mvjpGDBlkN4wXviUWy/DppE/5L8RqGao+L7N7Rlwp
w1Hj3r/9h1grg/71U1kDxEMyvB5rsOkn5slR0dj65DbXl8Pmws7vfuD6cvhL
ztV1cH05Jp5PnG3n+nJEjy2v03B9ObZ3mGpbub4cl6wEj5q5vhw+tk1zm7i+
Aq3bxhLrub4CVvYvq89xfQUqc9fyqri+AmKXJ6rTXF8BI8vcgeNcX4HflPfF
n3N9BQq3dr0u4foKDLZbHSrk+krUTpZ4FnB9JXIG3j/ex/WVGPHj1+RxfSUG
pngbd3F9Ja7oLg/dzvWV+NW9cslmrq9EksTSfAPXV6JAvHHmOq6vQtej3cZK
rq/CeIHZ/Cyur8LZ60YemVxfhcThlcI0rq+Ch8uOPUlcX4WdOV9ciuP6Kmzf
Ujgi4PoqKDTmPlFc/884VmxRGP47/w8OkLuo
     "]]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesLabel->{None, None},
  AxesOrigin->{0, 1.500856221334229},
  Method->{},
  PlotRange->{{0, 200.}, {1.500856221334229, 4.141258146630319}},
  PlotRangeClipping->True,
  PlotRangePadding->{{4., 4.}, {0.052808038505921806`, 
   0.052808038505921806`}}]], "Output",
 CellChangeTimes->{{3.876386544532989*^9, 3.876386615160928*^9}, 
   3.8763867513889647`*^9, 3.8763888749554234`*^9, {3.876392779359069*^9, 
   3.876392802703229*^9}, 3.876545897644453*^9, 3.8765468778294477`*^9, 
   3.876547012270105*^9, {3.876547218452455*^9, 3.87654728063571*^9}, 
   3.876556256678741*^9, 3.876556718152462*^9, 3.876560112492235*^9, 
   3.876560172024962*^9, {3.876560752113492*^9, 3.876560793047056*^9}, 
   3.8765609490264387`*^9, 3.876564978348179*^9, {3.876565089364338*^9, 
   3.876565142129472*^9}, 3.876572676364606*^9, 3.876997262674312*^9}],

Cell[BoxData["1.500856221334229`"], "Output",
 CellChangeTimes->{{3.876386544532989*^9, 3.876386615160928*^9}, 
   3.8763867513889647`*^9, 3.8763888749554234`*^9, {3.876392779359069*^9, 
   3.876392802703229*^9}, 3.876545897644453*^9, 3.8765468778294477`*^9, 
   3.876547012270105*^9, {3.876547218452455*^9, 3.87654728063571*^9}, 
   3.876556256678741*^9, 3.876556718152462*^9, 3.876560112492235*^9, 
   3.876560172024962*^9, {3.876560752113492*^9, 3.876560793047056*^9}, 
   3.8765609490264387`*^9, 3.876564978348179*^9, {3.876565089364338*^9, 
   3.876565142129472*^9}, 3.876572676364606*^9, 3.876997262683505*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", "1", "}"}], "}"}]], "Output",
 CellChangeTimes->{{3.876386544532989*^9, 3.876386615160928*^9}, 
   3.8763867513889647`*^9, 3.8763888749554234`*^9, {3.876392779359069*^9, 
   3.876392802703229*^9}, 3.876545897644453*^9, 3.8765468778294477`*^9, 
   3.876547012270105*^9, {3.876547218452455*^9, 3.87654728063571*^9}, 
   3.876556256678741*^9, 3.876556718152462*^9, 3.876560112492235*^9, 
   3.876560172024962*^9, {3.876560752113492*^9, 3.876560793047056*^9}, 
   3.8765609490264387`*^9, 3.876564978348179*^9, {3.876565089364338*^9, 
   3.876565142129472*^9}, 3.876572676364606*^9, 3.876997262684637*^9}],

Cell[BoxData["4.141258146630319`"], "Output",
 CellChangeTimes->{{3.876386544532989*^9, 3.876386615160928*^9}, 
   3.8763867513889647`*^9, 3.8763888749554234`*^9, {3.876392779359069*^9, 
   3.876392802703229*^9}, 3.876545897644453*^9, 3.8765468778294477`*^9, 
   3.876547012270105*^9, {3.876547218452455*^9, 3.87654728063571*^9}, 
   3.876556256678741*^9, 3.876556718152462*^9, 3.876560112492235*^9, 
   3.876560172024962*^9, {3.876560752113492*^9, 3.876560793047056*^9}, 
   3.8765609490264387`*^9, 3.876564978348179*^9, {3.876565089364338*^9, 
   3.876565142129472*^9}, 3.876572676364606*^9, 3.876997262685747*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", "15", "}"}], "}"}]], "Output",
 CellChangeTimes->{{3.876386544532989*^9, 3.876386615160928*^9}, 
   3.8763867513889647`*^9, 3.8763888749554234`*^9, {3.876392779359069*^9, 
   3.876392802703229*^9}, 3.876545897644453*^9, 3.8765468778294477`*^9, 
   3.876547012270105*^9, {3.876547218452455*^9, 3.87654728063571*^9}, 
   3.876556256678741*^9, 3.876556718152462*^9, 3.876560112492235*^9, 
   3.876560172024962*^9, {3.876560752113492*^9, 3.876560793047056*^9}, 
   3.8765609490264387`*^9, 3.876564978348179*^9, {3.876565089364338*^9, 
   3.876565142129472*^9}, 3.876572676364606*^9, 3.876997262686842*^9}],

Cell[BoxData["1.9292684051551754`"], "Output",
 CellChangeTimes->{{3.876386544532989*^9, 3.876386615160928*^9}, 
   3.8763867513889647`*^9, 3.8763888749554234`*^9, {3.876392779359069*^9, 
   3.876392802703229*^9}, 3.876545897644453*^9, 3.8765468778294477`*^9, 
   3.876547012270105*^9, {3.876547218452455*^9, 3.87654728063571*^9}, 
   3.876556256678741*^9, 3.876556718152462*^9, 3.876560112492235*^9, 
   3.876560172024962*^9, {3.876560752113492*^9, 3.876560793047056*^9}, 
   3.8765609490264387`*^9, 3.876564978348179*^9, {3.876565089364338*^9, 
   3.876565142129472*^9}, 3.876572676364606*^9, 3.876997262687944*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "3", ",", "4", ",", "5", ",", "6", ",", "7", ",", "8", ",", "9", ",", "10", 
   ",", "11", ",", "12", ",", "13", ",", "14", ",", "15", ",", "16", ",", 
   "17", ",", "18", ",", "19", ",", "20", ",", "21", ",", "22", ",", "23", 
   ",", "24", ",", "25", ",", "26", ",", "27", ",", "28", ",", "29", ",", 
   "30", ",", "31", ",", "32", ",", "33", ",", "34", ",", "35", ",", "36", 
   ",", "37", ",", "38", ",", "39", ",", "40", ",", "41", ",", "42", ",", 
   "43", ",", "44", ",", "45", ",", "46", ",", "47", ",", "48", ",", "49", 
   ",", "50", ",", "51", ",", "52", ",", "53", ",", "54", ",", "55", ",", 
   "56", ",", "57", ",", "58", ",", "59", ",", "60", ",", "61", ",", "62", 
   ",", "63", ",", "64", ",", "65", ",", "66", ",", "67", ",", "68", ",", 
   "69", ",", "70", ",", "71", ",", "72", ",", "73", ",", "74", ",", "75", 
   ",", "76", ",", "77", ",", "78", ",", "79", ",", "80", ",", "81", ",", 
   "82", ",", "83", ",", "84", ",", "85", ",", "86", ",", "87"}], 
  "}"}]], "Output",
 CellChangeTimes->{{3.876386544532989*^9, 3.876386615160928*^9}, 
   3.8763867513889647`*^9, 3.8763888749554234`*^9, {3.876392779359069*^9, 
   3.876392802703229*^9}, 3.876545897644453*^9, 3.8765468778294477`*^9, 
   3.876547012270105*^9, {3.876547218452455*^9, 3.87654728063571*^9}, 
   3.876556256678741*^9, 3.876556718152462*^9, 3.876560112492235*^9, 
   3.876560172024962*^9, {3.876560752113492*^9, 3.876560793047056*^9}, 
   3.8765609490264387`*^9, 3.876564978348179*^9, {3.876565089364338*^9, 
   3.876565142129472*^9}, 3.876572676364606*^9, 3.876997262689302*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "5", ",", "6", ",", "7", ",", "8", ",", "9", ",", "10", ",", "11", ",", 
   "12", ",", "13", ",", "14", ",", "15", ",", "16", ",", "17", ",", "18", 
   ",", "19", ",", "20", ",", "21", ",", "22", ",", "23", ",", "24", ",", 
   "25", ",", "26", ",", "27", ",", "28", ",", "29", ",", "30", ",", "31", 
   ",", "32", ",", "33", ",", "34", ",", "35", ",", "36", ",", "37", ",", 
   "38", ",", "39", ",", "40", ",", "41"}], "}"}]], "Output",
 CellChangeTimes->{{3.876386544532989*^9, 3.876386615160928*^9}, 
   3.8763867513889647`*^9, 3.8763888749554234`*^9, {3.876392779359069*^9, 
   3.876392802703229*^9}, 3.876545897644453*^9, 3.8765468778294477`*^9, 
   3.876547012270105*^9, {3.876547218452455*^9, 3.87654728063571*^9}, 
   3.876556256678741*^9, 3.876556718152462*^9, 3.876560112492235*^9, 
   3.876560172024962*^9, {3.876560752113492*^9, 3.876560793047056*^9}, 
   3.8765609490264387`*^9, 3.876564978348179*^9, {3.876565089364338*^9, 
   3.876565142129472*^9}, 3.876572676364606*^9, 3.876997262690847*^9}]
}, Open  ]]
}, Open  ]]
},
WindowSize->{960, 1028},
WindowMargins->{{-10, Automatic}, {Automatic, -10}},
FrontEndVersion->"9.0 for Linux x86 (64-bit) (February 7, 2013)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[557, 20, 1845, 48, 129, "Code"],
Cell[CellGroupData[{
Cell[2427, 72, 3006, 79, 357, "Code"],
Cell[5436, 153, 499, 10, 32, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5972, 168, 904, 19, 90, "Code"],
Cell[6879, 189, 459, 6, 32, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[7375, 200, 122, 1, 80, "Section"],
Cell[CellGroupData[{
Cell[7522, 205, 110, 1, 44, "Subsection"],
Cell[CellGroupData[{
Cell[7657, 210, 2344, 47, 90, "Code"],
Cell[10004, 259, 787, 19, 70, "Output"],
Cell[10794, 280, 761, 18, 70, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[11604, 304, 111, 1, 44, "Subsection"],
Cell[CellGroupData[{
Cell[11740, 309, 2554, 50, 90, "Code"],
Cell[14297, 361, 563, 11, 70, "Output"],
Cell[14863, 374, 565, 11, 70, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[15477, 391, 104, 1, 44, "Subsection"],
Cell[CellGroupData[{
Cell[15606, 396, 2631, 72, 394, "Code"],
Cell[18240, 470, 950, 19, 70, "Output"],
Cell[19193, 491, 998, 24, 70, "Output"],
Cell[20194, 517, 701, 14, 70, "Output"],
Cell[20898, 533, 510, 7, 70, "Output"],
Cell[21411, 542, 528, 7, 70, "Output"],
Cell[21942, 551, 530, 8, 70, "Output"],
Cell[22475, 561, 511, 7, 70, "Output"],
Cell[22989, 570, 1350, 20, 70, "Output"],
Cell[24342, 592, 2385, 43, 70, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[26776, 641, 102, 1, 44, "Subsection"],
Cell[CellGroupData[{
Cell[26903, 646, 1027, 27, 128, "Code"],
Cell[27933, 675, 2835, 52, 70, "Output"]
}, Open  ]],
Cell[30783, 730, 93, 1, 32, "Input"]
}, Open  ]]
}, Closed]],
Cell[CellGroupData[{
Cell[30925, 737, 170, 2, 50, "Section"],
Cell[CellGroupData[{
Cell[31120, 743, 112, 1, 44, "Subsection"],
Cell[CellGroupData[{
Cell[31257, 748, 2707, 59, 90, "Code"],
Cell[33967, 809, 861, 20, 37, "Output"],
Cell[34831, 831, 833, 19, 37, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[35713, 856, 111, 1, 44, "Subsection"],
Cell[CellGroupData[{
Cell[35849, 861, 2688, 54, 90, "Code"],
Cell[38540, 917, 730, 17, 32, "Output"],
Cell[39273, 936, 733, 17, 32, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[40055, 959, 104, 1, 44, "Subsection"],
Cell[CellGroupData[{
Cell[40184, 964, 2752, 75, 394, "Code"],
Cell[42939, 1041, 1126, 28, 37, "Output"],
Cell[44068, 1071, 814, 18, 32, "Output"],
Cell[44885, 1091, 676, 14, 32, "Output"],
Cell[45564, 1107, 482, 7, 32, "Output"],
Cell[46049, 1116, 498, 7, 32, "Output"],
Cell[46550, 1125, 503, 8, 32, "Output"],
Cell[47056, 1135, 498, 7, 32, "Output"],
Cell[47557, 1144, 1218, 19, 77, "Output"],
Cell[48778, 1165, 2656, 47, 242, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[51483, 1218, 103, 1, 44, "Subsection"],
Cell[CellGroupData[{
Cell[51611, 1223, 1035, 27, 128, "Code"],
Cell[52649, 1252, 3341, 60, 227, "Output"]
}, Open  ]]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[56051, 1319, 94, 1, 80, "Section"],
Cell[CellGroupData[{
Cell[56170, 1324, 1840, 46, 356, "Code"],
Cell[58013, 1372, 5652, 98, 227, "Output"],
Cell[63668, 1472, 3636, 63, 241, "Output"],
Cell[67307, 1537, 616, 8, 32, "Output"],
Cell[67926, 1547, 644, 10, 32, "Output"],
Cell[68573, 1559, 616, 8, 32, "Output"],
Cell[69192, 1569, 645, 10, 32, "Output"],
Cell[69840, 1581, 617, 8, 32, "Output"],
Cell[70460, 1591, 1606, 24, 99, "Output"],
Cell[72069, 1617, 1049, 16, 55, "Output"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)