(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 11.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 158, 7] NotebookDataLength[ 237098, 5549] NotebookOptionsPosition[ 234768, 5513] NotebookOutlinePosition[ 235142, 5529] CellTagsIndexPosition[ 235099, 5526] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{ RowBox[{"F", "[", "x_", "]"}], ":=", RowBox[{ FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", RowBox[{"(", RowBox[{ RowBox[{"PolyLog", "[", RowBox[{"2", ",", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", SuperscriptBox["\[ExponentialE]", "x"]}]}], "]"}], "-", RowBox[{"PolyLog", "[", RowBox[{"2", ",", RowBox[{"\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", "x"]}]}], "]"}]}], ")"}]}]}]], "Input", CellChangeTimes->{{3.737196866673444*^9, 3.737196875196076*^9}}, CellLabel->"In[23]:=",ExpressionUUID->"748eabe5-c3d2-4801-b8b7-aa7a789097c7"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"G3", "[", RowBox[{"a_", ",", "b_"}], "]"}], ":=", RowBox[{ RowBox[{"F", "[", RowBox[{"a", "-", "b"}], "]"}], "+", RowBox[{"F", "[", RowBox[{"b", "-", "a"}], "]"}], "-", RowBox[{ FractionBox["\[Pi]", "2"], RowBox[{"(", RowBox[{"a", "+", "b"}], ")"}]}]}]}], ";"}]], "Input", CellChangeTimes->{ 3.737201965981352*^9, {3.737202002231657*^9, 3.737202033142387*^9}, { 3.737202126146511*^9, 3.7372021534299583`*^9}, {3.7378972597648673`*^9, 3.7378972820313244`*^9}, {3.737897408257219*^9, 3.7378974198940997`*^9}}, CellLabel->"In[22]:=",ExpressionUUID->"0ab29ce4-3405-44f4-bf4a-16ef730ff686"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"EnergySetup", "[", RowBox[{"n_", ",", "FacesCW_"}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{ "i", ",", "rhoV", ",", "rhoF", ",", "include", ",", "LineFaces", ",", "f", ",", "Faces", ",", "NM", ",", "nonlinefaces", ",", "radii", ",", "RADII"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"Faces", "=", RowBox[{"Sort", "/@", "FacesCW"}]}], ";", "\[IndentingNewLine]", RowBox[{"f", "=", RowBox[{"Length", "[", "Faces", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"LineFaces", "=", RowBox[{"{", "}"}]}], ";", RowBox[{"rhoV", "=", RowBox[{"Table", "[", RowBox[{ SubscriptBox["rV", "i"], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"n", "-", "2"}]}], "}"}]}], "]"}]}], ";", RowBox[{"rhoF", "=", RowBox[{"Table", "[", RowBox[{ SubscriptBox["rF", "i"], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"f", "-", "2"}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"include", " ", "=", " ", RowBox[{"{", "}"}]}], ";", "\[IndentingNewLine]", RowBox[{"nonlinefaces", "=", RowBox[{"{", "}"}]}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", "f"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"!", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"Faces", "[", RowBox[{"[", RowBox[{"i", ",", RowBox[{"-", "2"}]}], "]"}], "]"}], "\[Equal]", RowBox[{"n", "-", "1"}]}], "&&", RowBox[{ RowBox[{"Faces", "[", RowBox[{"[", RowBox[{"i", ",", RowBox[{"-", "1"}]}], "]"}], "]"}], "\[Equal]", "n"}]}], ")"}]}], ",", "\[IndentingNewLine]", RowBox[{"(*", " ", "then", " ", "*)"}], " ", RowBox[{ RowBox[{"AppendTo", "[", RowBox[{"include", ",", "i"}], "]"}], ";", RowBox[{"AppendTo", "[", RowBox[{"nonlinefaces", ",", RowBox[{"FacesCW", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], ",", "\[IndentingNewLine]", RowBox[{"(*", " ", "else", " ", "*)"}], " ", RowBox[{"AppendTo", "[", RowBox[{"LineFaces", ",", "i"}], "]"}]}], "]"}]}], "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"NM", "=", RowBox[{"NMinimize", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"k", "=", "1"}], RowBox[{"f", "-", "2"}]], RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], RowBox[{"Length", "[", RowBox[{"Faces", "[", RowBox[{"[", RowBox[{"include", "[", RowBox[{"[", "k", "]"}], "]"}], "]"}], "]"}], "]"}]], RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Faces", "[", RowBox[{"[", RowBox[{ RowBox[{"include", "[", RowBox[{"[", "k", "]"}], "]"}], ",", "i"}], "]"}], "]"}], "<", RowBox[{"n", "-", "1"}]}], ",", RowBox[{"G3", "[", RowBox[{ SubscriptBox["rV", RowBox[{"Faces", "[", RowBox[{"[", RowBox[{ RowBox[{"include", "[", RowBox[{"[", "k", "]"}], "]"}], ",", "i"}], "]"}], "]"}]], ",", SubscriptBox["rF", "k"]}], "]"}], ",", "0"}], "]"}]}]}], "+", RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], RowBox[{"n", "-", "2"}]], RowBox[{ SubscriptBox["rV", "i"], "\[Pi]", " ", RowBox[{"If", "[", RowBox[{ RowBox[{"MemberQ", "[", RowBox[{ RowBox[{"Flatten", "@", RowBox[{"Table", "[", RowBox[{ RowBox[{"Faces", "[", RowBox[{"[", RowBox[{"LineFaces", "[", RowBox[{"[", "m", "]"}], "]"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"m", ",", "1", ",", "2"}], "}"}]}], "]"}]}], ",", "i"}], "]"}], ",", RowBox[{"(*", " ", "then", " ", "*)"}], " ", "1", ",", RowBox[{"(*", " ", "else", " ", "*)"}], " ", "2"}], "]"}]}]}], "+", RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"j", "=", "1"}], RowBox[{"f", "-", "2"}]], RowBox[{ SubscriptBox["rF", "j"], "\[Pi]", " ", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"MemberQ", "[", RowBox[{ RowBox[{"Faces", "[", RowBox[{"[", RowBox[{"include", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}], "]"}], ",", RowBox[{"n", "-", "1"}]}], "]"}], "||", RowBox[{"MemberQ", "[", RowBox[{ RowBox[{"Faces", "[", RowBox[{"[", RowBox[{"include", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}], "]"}], ",", "n"}], "]"}]}], ",", RowBox[{"(*", " ", "then", " ", "*)"}], " ", "1", ",", RowBox[{"(*", " ", "else", " ", "*)"}], " ", "2"}], "]"}]}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], RowBox[{"n", "-", "2"}]], SubscriptBox["rV", "i"]}], "+", RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"j", "=", "1"}], RowBox[{"f", "-", "2"}]], SubscriptBox["rF", "j"]}]}], "\[Equal]", "0"}]}], "}"}], ",", RowBox[{"Join", "[", RowBox[{"rhoV", ",", "rhoF"}], "]"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"RADII", "=", RowBox[{"Exp", "[", RowBox[{"Re", "[", RowBox[{ RowBox[{"Join", "[", RowBox[{"rhoV", ",", "rhoF"}], "]"}], "/.", RowBox[{"(", RowBox[{"NM", "[", RowBox[{"[", "2", "]"}], "]"}], ")"}]}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"RADII", "/=", RowBox[{"RADII", "[", RowBox[{"[", "1", "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"radii", "=", RowBox[{"<|", "|>"}]}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", RowBox[{"n", "-", "2"}]}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"radii", "[", "i", "]"}], "=", RowBox[{"RADII", "[", RowBox[{"[", "i", "]"}], "]"}]}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", RowBox[{"Length", "[", "nonlinefaces", "]"}]}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"radii", "[", RowBox[{"nonlinefaces", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], "=", RowBox[{"RADII", "[", RowBox[{"[", RowBox[{"n", "-", "2", "+", "i"}], "]"}], "]"}]}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"radii", ",", "LineFaces"}], "}"}]}]}], "\[IndentingNewLine]", "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.7372045119396753`*^9, 3.7372046545187597`*^9}, { 3.737204685631832*^9, 3.7372047308167467`*^9}, {3.737204773625032*^9, 3.7372048110040503`*^9}, {3.737204842780097*^9, 3.737204924806987*^9}, { 3.7372050053122187`*^9, 3.737205008911607*^9}, {3.73720504639081*^9, 3.7372051214979887`*^9}, {3.7372051736233463`*^9, 3.737205179673382*^9}, { 3.737205261713398*^9, 3.7372052618754*^9}, {3.737205306100531*^9, 3.737205373527009*^9}, {3.737205491627762*^9, 3.737205535178033*^9}, { 3.737206922016953*^9, 3.737206989324239*^9}, {3.737207032865409*^9, 3.737207073209803*^9}, {3.737284915566703*^9, 3.737284955928123*^9}, { 3.737371291071519*^9, 3.7373713419977837`*^9}, {3.737468221809287*^9, 3.737468222498644*^9}, 3.737468268451545*^9, {3.7374683559250183`*^9, 3.737468513179357*^9}, {3.737468554130938*^9, 3.7374685561663094`*^9}, { 3.7374686045167637`*^9, 3.737468621354309*^9}, {3.737468670066434*^9, 3.737468700374711*^9}, {3.73746873598547*^9, 3.737468789149178*^9}, { 3.737469118685132*^9, 3.737469121422579*^9}, {3.737891244695794*^9, 3.737891247299219*^9}, {3.737995135349449*^9, 3.73799513678258*^9}}, CellLabel->"In[21]:=",ExpressionUUID->"72e544ec-21b3-4511-8714-78e40bf96f8f"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"Boundary", "[", RowBox[{"n_", ",", "FacesCW_", ",", "LineFaces_", ",", "radii_", ",", RowBox[{"\[Epsilon]_:", "0"}]}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{ "i", ",", "j", ",", "faces", ",", "vertices", ",", "lastvisited", ",", "d", ",", "h", ",", "H", ",", "breaker", ",", "unvisited", ",", "f", ",", "endcircle", ",", "ComputeEndcircle", ",", "ComputeEndcircle2", ",", "delete"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"unvisited", " ", "=", "FacesCW"}], ";", RowBox[{"(*", RowBox[{ "Keeps", " ", "track", " ", "of", " ", "the", " ", "faces", " ", "we", " ", RowBox[{"haven", "'"}], "t", " ", "visited"}], "*)"}], "\[IndentingNewLine]", RowBox[{"f", "=", RowBox[{"Length", "[", "FacesCW", "]"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"faces", "=", RowBox[{"<|", "|>"}]}], ";", " ", RowBox[{"(*", " ", RowBox[{"centers", " ", "of", " ", "face", " ", "circles", " ", RowBox[{"(", RowBox[{ RowBox[{"or", " ", "normals"}], ",", " ", RowBox[{"for", " ", "lines"}]}], ")"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", "f"}], ",", RowBox[{"i", "++"}], ",", RowBox[{ RowBox[{"faces", "[", RowBox[{"FacesCW", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], "=", RowBox[{"{", RowBox[{ RowBox[{"-", "\[Infinity]"}], ",", RowBox[{"-", "\[Infinity]"}]}], "}"}]}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"vertices", "=", RowBox[{"<|", "|>"}]}], ";", " ", RowBox[{"(*", " ", RowBox[{"centers", " ", "of", " ", "vertex", " ", "circles", " ", RowBox[{"(", RowBox[{ RowBox[{"or", " ", "normals"}], ",", " ", RowBox[{"for", " ", "lines"}]}], ")"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", RowBox[{"n", "-", "2"}]}], ",", RowBox[{"i", "++"}], ",", RowBox[{ RowBox[{"vertices", "[", "i", "]"}], "=", RowBox[{"{", RowBox[{ RowBox[{"-", "\[Infinity]"}], ",", RowBox[{"-", "\[Infinity]"}]}], "}"}]}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{ "Radii", " ", "is", " ", "an", " ", "assoc", " ", "from", " ", "v"}], "\[Rule]", RowBox[{ RowBox[{"rad", " ", "and", " ", "facelist"}], "\[Rule]", "rad"}]}], " ", "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ "Vertex", " ", "circles", " ", "going", " ", "vertically", " ", "upwards", " ", "along", " ", "first", " ", "FaceLine"}], "*)"}], "\[IndentingNewLine]", RowBox[{"h", "=", "0"}], ";", " ", RowBox[{"(*", " ", RowBox[{"height", " ", "along", " ", "the", " ", "left", " ", "line"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", RowBox[{ RowBox[{"Length", "[", RowBox[{"FacesCW", "[", RowBox[{"[", RowBox[{"LineFaces", "[", RowBox[{"[", "1", "]"}], "]"}], "]"}], "]"}], "]"}], "-", "2"}]}], RowBox[{"(*", RowBox[{ "Number", " ", "of", " ", "vertex", " ", "circles", " ", "that", " ", "lie", " ", "along", " ", "first", " ", "LineFace"}], "*)"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"vertices", "[", RowBox[{"FacesCW", "[", RowBox[{"[", RowBox[{ RowBox[{"LineFaces", "[", RowBox[{"[", "1", "]"}], "]"}], ",", "i"}], "]"}], "]"}], "]"}], "=", RowBox[{"{", RowBox[{"0", ",", RowBox[{"h", "+", RowBox[{"radii", "[", RowBox[{"FacesCW", "[", RowBox[{"[", RowBox[{ RowBox[{"LineFaces", "[", RowBox[{"[", "1", "]"}], "]"}], ",", "i"}], "]"}], "]"}], "]"}]}]}], "}"}]}], ";", RowBox[{"(*", RowBox[{ "Set", " ", "the", " ", "center", " ", "of", " ", "this", " ", "vertex", " ", "circle"}], "*)"}], "\[IndentingNewLine]", RowBox[{"h", "+=", RowBox[{"2", RowBox[{"radii", "[", RowBox[{"[", RowBox[{"FacesCW", "[", RowBox[{"[", RowBox[{ RowBox[{"LineFaces", "[", RowBox[{"[", "1", "]"}], "]"}], ",", "i"}], "]"}], "]"}], "]"}], "]"}]}]}]}]}], "]"}], ";", RowBox[{"(*", RowBox[{ "Keeps", " ", "track", " ", "of", " ", "our", " ", "cumulative", " ", "height", " ", "as", " ", "we", " ", "climb", " ", "up", " ", "the", " ", "first", " ", "FaceLine"}], "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"unvisited", "=", RowBox[{"Delete", "[", RowBox[{"unvisited", ",", RowBox[{"Position", "[", RowBox[{"unvisited", ",", RowBox[{"FacesCW", "[", RowBox[{"[", RowBox[{"LineFaces", "[", RowBox[{"[", "1", "]"}], "]"}], "]"}], "]"}]}], "]"}]}], "]"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"d", "=", "0"}], ";", RowBox[{"(*", RowBox[{"Horizontal", " ", "distance"}], "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"ComputeEndcircle", "[", RowBox[{"prev_", ",", "face_"}], "]"}], ":=", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ "Returns", " ", "which", " ", "vertex", " ", "circle", " ", "we", " ", "calculated", " ", "last", " ", "based", " ", "on", " ", "previous", " ", "vertex", " ", "circle", " ", "and", " ", "face"}], "*)"}], "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"face", "[", RowBox[{"[", RowBox[{ RowBox[{"Mod", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"Position", "[", RowBox[{"face", ",", RowBox[{"n", "-", "1"}]}], "]"}], "[", RowBox[{"[", RowBox[{"1", ",", "1"}], "]"}], "]"}], "-", "2"}], ",", RowBox[{"Length", "[", "face", "]"}]}], "]"}], "+", "1"}], "]"}], "]"}], "\[Equal]", "prev"}], ",", "\[IndentingNewLine]", RowBox[{"(*", "then", "*)"}], RowBox[{"face", "[", RowBox[{"[", RowBox[{ RowBox[{"Mod", "[", RowBox[{ RowBox[{ RowBox[{"Position", "[", RowBox[{"face", ",", RowBox[{"n", "-", "1"}]}], "]"}], "[", RowBox[{"[", RowBox[{"1", ",", "1"}], "]"}], "]"}], ",", RowBox[{"Length", "[", "face", "]"}]}], "]"}], "+", "1"}], "]"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"(*", "else", "*)"}], " ", RowBox[{"face", "[", RowBox[{"[", RowBox[{ RowBox[{"Mod", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"Position", "[", RowBox[{"face", ",", RowBox[{"n", "-", "1"}]}], "]"}], "[", RowBox[{"[", RowBox[{"1", ",", "1"}], "]"}], "]"}], "-", "2"}], ",", RowBox[{"Length", "[", "face", "]"}]}], "]"}], "+", "1"}], "]"}], "]"}]}], "]"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"endcircle", "=", RowBox[{"ComputeEndcircle", "[", RowBox[{"n", ",", RowBox[{"FacesCW", "[", RowBox[{"[", RowBox[{"LineFaces", "[", RowBox[{"[", "1", "]"}], "]"}], "]"}], "]"}]}], "]"}]}], ";", RowBox[{"(*", RowBox[{ RowBox[{ "Previous", " ", "vertex", " ", "circle", " ", "is", " ", "bottom", " ", "line"}], ",", " ", RowBox[{"face", " ", "is", " ", "first", " ", "FaceLine"}]}], "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ "Face", " ", "circles", " ", "going", " ", "horizontally", " ", "across", " ", "upper", " ", "vertex", " ", "line"}], "*)"}], "\[IndentingNewLine]", RowBox[{"breaker", "=", "True"}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "1"}], ",", RowBox[{"j", "\[LessEqual]", "f"}], " ", ",", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", RowBox[{"Length", "[", "unvisited", "]"}]}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"MemberQ", "[", RowBox[{ RowBox[{"unvisited", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"n", "-", "1"}]}], "]"}], ",", " ", RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{"determine", " ", "if", " ", "this", " ", "face"}], "-", RowBox[{"circle", " ", "contains", " ", "n"}], "-", "1"}], ",", " ", RowBox[{ "so", " ", "it", " ", "is", " ", "on", " ", "the", " ", "upper", " ", "vertex", " ", "line"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"MemberQ", "[", RowBox[{ RowBox[{"unvisited", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"ComputeEndcircle", "[", RowBox[{"endcircle", ",", RowBox[{"unvisited", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"determine", " ", "if", " ", "this", " ", "face"}], "-", RowBox[{ "circle", " ", "also", " ", "contains", " ", "the", " ", "vertex", " ", "right", " ", "before", " ", "n"}], "-", RowBox[{ "1", " ", "to", " ", "ensure", " ", "we", " ", "are", " ", "going", " ", "horizontally", " ", "across", " ", "in", " ", "the", " ", "right", " ", "order"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"unvisited", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", RowBox[{"FacesCW", "[", RowBox[{"[", RowBox[{"LineFaces", "[", RowBox[{"[", "2", "]"}], "]"}], "]"}], "]"}]}], ",", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"If", " ", RowBox[{"we", "'"}], "ve", " ", "reached", " ", "the", " ", "second", " ", "faceline"}], ",", " ", RowBox[{"then", " ", "do", " ", RowBox[{"nothing", "--"}], "find", " ", "the", " ", "other", " ", "face", " ", "circle", " ", "containing", " ", "endcircle"}]}], "*)"}], "\[IndentingNewLine]", ",", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ "else", " ", "if", " ", "the", " ", "face", " ", "circle", " ", "is", " ", "not", " ", "a", " ", "line"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"faces", "[", RowBox[{"unvisited", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], "=", RowBox[{"{", RowBox[{ RowBox[{"d", "+", RowBox[{"radii", "[", RowBox[{"unvisited", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}], ",", "h"}], "}"}]}], ";", " ", RowBox[{"(*", RowBox[{ "Set", " ", "the", " ", "center", " ", "of", " ", "this", " ", "face", " ", "circle"}], "*)"}], "\[IndentingNewLine]", RowBox[{"d", "+=", RowBox[{"2", RowBox[{"radii", "[", RowBox[{"unvisited", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}]}], ";", " ", RowBox[{"(*", " ", RowBox[{ "update", " ", "how", " ", "far", " ", "along", " ", "the", " ", "top", " ", "line", " ", "we", " ", "are"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"endcircle", "=", RowBox[{"ComputeEndcircle", "[", RowBox[{"endcircle", ",", RowBox[{"unvisited", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], ";", RowBox[{"(*", RowBox[{"Update", " ", "endcircle"}], "*)"}], "\[IndentingNewLine]", RowBox[{"unvisited", "=", RowBox[{"Delete", "[", RowBox[{"unvisited", ",", "i"}], "]"}]}], ";"}]}], "]"}]}], " ", RowBox[{"(*", " ", RowBox[{ "remove", " ", "that", " ", "face", " ", "from", " ", "the", " ", "unvisited", " ", "list"}], " ", "*)"}], "\[IndentingNewLine]", "]"}]}], "\[IndentingNewLine]", "]"}]}], "\[IndentingNewLine]", "]"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{ "Vertex", " ", "circles", " ", "along", " ", "the", " ", "second", " ", "FaceLine"}], ",", " ", RowBox[{"climbing", " ", "vertically", " ", "upwards"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"h", "=", "0"}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", RowBox[{ RowBox[{"Length", "[", RowBox[{"FacesCW", "[", RowBox[{"[", RowBox[{"LineFaces", "[", RowBox[{"[", "2", "]"}], "]"}], "]"}], "]"}], "]"}], "-", "2"}]}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"vertices", "[", RowBox[{"FacesCW", "[", RowBox[{"[", RowBox[{ RowBox[{"LineFaces", "[", RowBox[{"[", "2", "]"}], "]"}], ",", "i"}], "]"}], "]"}], "]"}], "=", RowBox[{"{", RowBox[{"d", ",", RowBox[{"h", "+", RowBox[{"radii", "[", RowBox[{"FacesCW", "[", RowBox[{"[", RowBox[{ RowBox[{"LineFaces", "[", RowBox[{"[", "2", "]"}], "]"}], ",", "i"}], "]"}], "]"}], "]"}]}]}], "}"}]}], ";", RowBox[{"h", "+=", RowBox[{"2", RowBox[{"radii", "[", RowBox[{"[", RowBox[{"FacesCW", "[", RowBox[{"[", RowBox[{ RowBox[{"LineFaces", "[", RowBox[{"[", "2", "]"}], "]"}], ",", "i"}], "]"}], "]"}], "]"}], "]"}]}]}]}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"H", "=", "h"}], ";", "\[IndentingNewLine]", RowBox[{"h", "=", "0"}], ";", RowBox[{"(*", RowBox[{"Reset", " ", "height"}], "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"ComputeEndcircle2", "[", RowBox[{"prev_", ",", "face_"}], "]"}], ":=", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"face", "[", RowBox[{"[", RowBox[{ RowBox[{"Mod", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"Position", "[", RowBox[{"face", ",", "n"}], "]"}], "[", RowBox[{"[", RowBox[{"1", ",", "1"}], "]"}], "]"}], "-", "2"}], ",", RowBox[{"Length", "[", "face", "]"}]}], "]"}], "+", "1"}], "]"}], "]"}], "\[Equal]", "prev"}], ",", "\[IndentingNewLine]", RowBox[{"(*", "then", "*)"}], RowBox[{"face", "[", RowBox[{"[", RowBox[{ RowBox[{"Mod", "[", RowBox[{ RowBox[{ RowBox[{"Position", "[", RowBox[{"face", ",", "n"}], "]"}], "[", RowBox[{"[", RowBox[{"1", ",", "1"}], "]"}], "]"}], ",", RowBox[{"Length", "[", "face", "]"}]}], "]"}], "+", "1"}], "]"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"(*", "else", "*)"}], " ", RowBox[{"face", "[", RowBox[{"[", RowBox[{ RowBox[{"Mod", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"Position", "[", RowBox[{"face", ",", "n"}], "]"}], "[", RowBox[{"[", RowBox[{"1", ",", "1"}], "]"}], "]"}], "-", "2"}], ",", RowBox[{"Length", "[", "face", "]"}]}], "]"}], "+", "1"}], "]"}], "]"}]}], "]"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"endcircle", "=", RowBox[{"ComputeEndcircle2", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", RowBox[{"FacesCW", "[", RowBox[{"[", RowBox[{"LineFaces", "[", RowBox[{"[", "1", "]"}], "]"}], "]"}], "]"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{ "Face", " ", "circles", " ", "along", " ", "lower", " ", "vertex", " ", "line"}], ",", " ", RowBox[{ "moving", " ", "horizontally", " ", "from", " ", "origin", " ", "to", " ", "second", " ", "FaceLine"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"d", "=", "0"}], ";", RowBox[{"(*", RowBox[{"Reset", " ", "distance"}], "*)"}], "\[IndentingNewLine]", RowBox[{"breaker", "=", "True"}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "1"}], ",", RowBox[{"j", "\[LessEqual]", "f"}], ",", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", RowBox[{"Length", "[", "unvisited", "]"}]}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"MemberQ", "[", RowBox[{ RowBox[{"unvisited", "[", RowBox[{"[", "i", "]"}], "]"}], ",", "n"}], "]"}], ",", " ", RowBox[{"(*", " ", RowBox[{ RowBox[{"determine", " ", "if", " ", "this", " ", "face"}], "-", RowBox[{"circle", " ", "contains", " ", "n"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"MemberQ", "[", RowBox[{ RowBox[{"unvisited", "[", RowBox[{"[", "i", "]"}], "]"}], ",", "endcircle"}], RowBox[{"(*", " ", RowBox[{"ComputeEndcircle2", "[", RowBox[{"endcircle", ",", RowBox[{"unvisited", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}], " ", "*)"}], "]"}], ",", RowBox[{"(*", " ", RowBox[{ RowBox[{"determine", " ", "if", " ", "this", " ", "face"}], "-", RowBox[{ "cirlce", " ", "also", " ", "contains", " ", "the", " ", "vertex", " ", "right", " ", "before", " ", "n"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"unvisited", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", RowBox[{"FacesCW", "[", RowBox[{"[", RowBox[{"LineFaces", "[", RowBox[{"[", "2", "]"}], "]"}], "]"}], "]"}]}], ",", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"If", " ", RowBox[{"we", "'"}], "ve", " ", "reached", " ", "the", " ", "second", " ", "faceline"}], ",", " ", RowBox[{"then", " ", "do", " ", "nothing"}]}], "*)"}], "\[IndentingNewLine]", ",", "\[IndentingNewLine]", RowBox[{"(*", " ", "else", " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"faces", "[", RowBox[{"unvisited", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], "=", RowBox[{"{", RowBox[{ RowBox[{"d", "+", RowBox[{"radii", "[", RowBox[{"unvisited", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}], ",", "h"}], "}"}]}], ";", " ", RowBox[{"(*", " ", RowBox[{"then", " ", RowBox[{"we", "'"}], "ll", " ", "make", " ", "this", " ", RowBox[{"face", "'"}], "s", " ", "center", " ", "be", " ", "where", " ", "we", " ", "know", " ", "it", " ", "is"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"d", "+=", RowBox[{"2", RowBox[{"radii", "[", RowBox[{"unvisited", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}]}], ";", " ", RowBox[{"(*", " ", RowBox[{ "update", " ", "how", " ", "far", " ", "along", " ", "the", " ", "bottom", " ", "line", " ", "we", " ", "are"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"endcircle", "=", RowBox[{"ComputeEndcircle2", "[", RowBox[{"endcircle", ",", RowBox[{"unvisited", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"unvisited", "=", RowBox[{"Delete", "[", RowBox[{"unvisited", ",", "i"}], "]"}]}], ";"}]}], "]"}]}], " ", RowBox[{"(*", " ", RowBox[{ "remove", " ", "that", " ", "face", " ", "from", " ", "the", " ", "unvisited", " ", "list"}], " ", "*)"}], "\[IndentingNewLine]", "]"}]}], "\[IndentingNewLine]", "]"}]}], "\[IndentingNewLine]", "]"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ "Set", " ", "the", " ", "normals", " ", "of", " ", "the", " ", "vertex", " ", "lines", " ", "and", " ", "face", " ", "lines"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"faces", "[", RowBox[{"FacesCW", "[", RowBox[{"[", RowBox[{"LineFaces", "[", RowBox[{"[", "2", "]"}], "]"}], "]"}], "]"}], "]"}], "=", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"faces", "[", RowBox[{"FacesCW", "[", RowBox[{"[", RowBox[{"LineFaces", "[", RowBox[{"[", "1", "]"}], "]"}], "]"}], "]"}], "]"}], "=", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "0"}], "}"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"vertices", "[", "n", "]"}], "=", RowBox[{"{", RowBox[{"0", ",", RowBox[{"-", "1"}]}], "}"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"vertices", "[", RowBox[{"n", "-", "1"}], "]"}], "=", RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"ThrowOrderError", "[", RowBox[{ "faces", ",", "vertices", ",", "radii", ",", "d", ",", "H", ",", "n", ",", RowBox[{"FacesCW", "[", RowBox[{"[", "LineFaces", "]"}], "]"}], ",", "0"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"faces", "[", RowBox[{"FacesCW", "[", RowBox[{"[", RowBox[{"LineFaces", "[", RowBox[{"[", "1", "]"}], "]"}], "]"}], "]"}], "]"}], "*=", RowBox[{"-", "1"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"faces", "[", RowBox[{"FacesCW", "[", RowBox[{"[", RowBox[{"LineFaces", "[", RowBox[{"[", "2", "]"}], "]"}], "]"}], "]"}], "]"}], "*=", RowBox[{"-", "1"}]}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", "n"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"i", "\[Equal]", RowBox[{"n", "-", "1"}]}], "||", RowBox[{"i", "\[Equal]", "n"}]}], ",", "\[IndentingNewLine]", ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{ RowBox[{"vertices", "[", "i", "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], "-", "d"}], "]"}], "\[LessEqual]", "\[Epsilon]"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"vertices", "[", "i", "]"}], "=", RowBox[{"{", RowBox[{"0", ",", RowBox[{ RowBox[{"vertices", "[", "i", "]"}], "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}]}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{ RowBox[{"vertices", "[", "i", "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], "-", "0"}], "]"}], "\[LessEqual]", "\[Epsilon]"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"vertices", "[", "i", "]"}], "=", RowBox[{"{", RowBox[{"d", ",", RowBox[{ RowBox[{"vertices", "[", "i", "]"}], "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}]}]}], "]"}]}], "]"}], ";"}]}], "\[IndentingNewLine]", "]"}]}], "]"}]}]}], "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"Return", " ", "the", " ", RowBox[{"assoc", "'"}], "s", " ", "of", " ", "the", " ", "faces", " ", "to", " ", "their", " ", "centers", " ", "and", " ", "vertices", " ", "to", " ", "their", " ", "centers"}], "*)"}], "\[IndentingNewLine]", RowBox[{"{", RowBox[{"vertices", ",", "faces", ",", "d", ",", "H"}], "}"}]}]}], "\[IndentingNewLine]", "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.737219030683413*^9, 3.737219068499385*^9}, { 3.737219129031519*^9, 3.737219157331717*^9}, {3.737219214776095*^9, 3.737219232499251*^9}, {3.737219530428797*^9, 3.737219737723432*^9}, { 3.737219973681823*^9, 3.737219981058498*^9}, {3.7372852522441998`*^9, 3.737285273052194*^9}, {3.737285379557729*^9, 3.7372853808850813`*^9}, { 3.737285506667365*^9, 3.737285515803762*^9}, {3.7372856044249487`*^9, 3.737285632285976*^9}, {3.7372857006889677`*^9, 3.7372857800137577`*^9}, { 3.73728584887117*^9, 3.737285857200169*^9}, {3.737286074306966*^9, 3.7372860829425373`*^9}, {3.7372861272847137`*^9, 3.737286147655734*^9}, { 3.737286180526772*^9, 3.737286183012919*^9}, {3.737286465705851*^9, 3.7372864888059*^9}, {3.7372865215606728`*^9, 3.737286606697309*^9}, { 3.737286645970365*^9, 3.737286648042437*^9}, {3.737286686140306*^9, 3.7372867642027607`*^9}, {3.737287269374832*^9, 3.737287342750917*^9}, { 3.737287384822444*^9, 3.7372877639349957`*^9}, {3.737287798168841*^9, 3.737287798698641*^9}, {3.737287874201215*^9, 3.737287886754616*^9}, { 3.73728794855105*^9, 3.73728794906709*^9}, {3.737287984087804*^9, 3.737287984258391*^9}, {3.737288167935055*^9, 3.7372881688958406`*^9}, { 3.737288229464828*^9, 3.7372884829941387`*^9}, {3.7372889651024017`*^9, 3.737289028710766*^9}, {3.73728913478118*^9, 3.7372891349574013`*^9}, { 3.737289176488153*^9, 3.7372892189711437`*^9}, {3.737289296434997*^9, 3.737289307920455*^9}, {3.737289377154367*^9, 3.7372893780167713`*^9}, { 3.73728946731621*^9, 3.7372894685890083`*^9}, {3.737289519964464*^9, 3.7372895646104527`*^9}, {3.737289627498472*^9, 3.7372896504948683`*^9}, { 3.737289698936853*^9, 3.7372897006060762`*^9}, {3.737289761143084*^9, 3.737289790713337*^9}, {3.737289852057372*^9, 3.737289856144655*^9}, { 3.7372899307846327`*^9, 3.737289954524984*^9}, {3.737290111472128*^9, 3.737290137729048*^9}, {3.7372903101204*^9, 3.7372903238682327`*^9}, { 3.737290363021036*^9, 3.737290395639821*^9}, {3.737290434642714*^9, 3.737290436215727*^9}, {3.737290481178443*^9, 3.737290524071663*^9}, { 3.737290651768982*^9, 3.7372906712458363`*^9}, {3.737290704533799*^9, 3.7372907061133537`*^9}, {3.737290744976763*^9, 3.737290766875083*^9}, { 3.737290999240056*^9, 3.737291019271243*^9}, {3.737291057840662*^9, 3.737291080456534*^9}, {3.737291128953987*^9, 3.737291140451013*^9}, { 3.737291260286092*^9, 3.7372914176080303`*^9}, {3.7372914768350573`*^9, 3.737291536874617*^9}, {3.7372915673943977`*^9, 3.73729157857458*^9}, { 3.73729459361259*^9, 3.737294594755714*^9}, {3.737294841256043*^9, 3.737294872151122*^9}, {3.73729499894667*^9, 3.737295013697496*^9}, { 3.737295117321186*^9, 3.737295146177869*^9}, {3.737295198600103*^9, 3.737295201561144*^9}, {3.737295332216126*^9, 3.7372953659228153`*^9}, { 3.7372958344031277`*^9, 3.737296025548025*^9}, {3.737302057273785*^9, 3.737302228939927*^9}, 3.7373022837982073`*^9, {3.737302343644175*^9, 3.737302375974662*^9}, 3.737302421339328*^9, {3.737302477569655*^9, 3.737302565791008*^9}, {3.737302777019927*^9, 3.737302823909573*^9}, { 3.737303123387793*^9, 3.737303214311923*^9}, {3.7373032530403357`*^9, 3.7373032965502157`*^9}, {3.737303417213832*^9, 3.737303518053563*^9}, { 3.737303590633803*^9, 3.737303600566893*^9}, {3.737303631996339*^9, 3.737303736869626*^9}, {3.737303779146243*^9, 3.737303825046517*^9}, { 3.737304130063211*^9, 3.737304181026523*^9}, {3.7373042470907288`*^9, 3.737304349277087*^9}, {3.7373043812177143`*^9, 3.737304407351365*^9}, { 3.737304581411367*^9, 3.737304616831106*^9}, {3.7373046620344*^9, 3.7373047068712873`*^9}, {3.737304771496949*^9, 3.7373048444644938`*^9}, { 3.737305306432076*^9, 3.7373053158581667`*^9}, {3.7373053546618013`*^9, 3.737305355635768*^9}, 3.7373054429132357`*^9, {3.737306098271888*^9, 3.737306106600956*^9}, {3.737306417696192*^9, 3.737306447095158*^9}, { 3.7373064781213903`*^9, 3.737306483985478*^9}, {3.737306573567268*^9, 3.737306587589746*^9}, {3.737306720838073*^9, 3.73730674819695*^9}, { 3.737306852180376*^9, 3.737306861658111*^9}, {3.737306929352565*^9, 3.7373069458296127`*^9}, {3.737306985973919*^9, 3.73730703321769*^9}, { 3.737307104020858*^9, 3.737307104380198*^9}, {3.7373072016479588`*^9, 3.737307219827284*^9}, {3.737307283259153*^9, 3.737307449938692*^9}, { 3.737307507268676*^9, 3.737307531627266*^9}, {3.737307624863483*^9, 3.737307630803578*^9}, 3.7373692366355057`*^9, {3.737378768428591*^9, 3.737378772696554*^9}, {3.737456548739272*^9, 3.737456676468967*^9}, { 3.737456732055797*^9, 3.737456785246457*^9}, {3.737456880379624*^9, 3.73745693261576*^9}, {3.7374569872793217`*^9, 3.737457002918715*^9}, { 3.737457034354722*^9, 3.7374570868457327`*^9}, {3.7374573446089993`*^9, 3.737457359698818*^9}, {3.737457621684537*^9, 3.737457654757286*^9}, { 3.7374577893447638`*^9, 3.7374578203273153`*^9}, 3.73745786787053*^9, { 3.737458414579453*^9, 3.7374584346434393`*^9}, {3.73745848723744*^9, 3.737458515616708*^9}, {3.737458663428769*^9, 3.73745867379788*^9}, { 3.737458742757172*^9, 3.7374588188135767`*^9}, {3.7374589524311037`*^9, 3.737458970370059*^9}, {3.737459004685546*^9, 3.737459017734614*^9}, { 3.737469143102393*^9, 3.737469152438889*^9}, {3.7374723221463118`*^9, 3.7374723517302732`*^9}, {3.737472517403573*^9, 3.73747252044235*^9}, { 3.737472566196344*^9, 3.737472599816885*^9}, {3.737472901109223*^9, 3.737472904411367*^9}, {3.737472959283648*^9, 3.737472959747048*^9}, { 3.737472997133112*^9, 3.737473040080291*^9}, {3.7374732109292183`*^9, 3.73747321567245*^9}, {3.737474220810277*^9, 3.737474243746545*^9}, { 3.737720030799748*^9, 3.7377200552173643`*^9}, {3.737720087337412*^9, 3.7377201185616283`*^9}, {3.737720170734948*^9, 3.737720176149166*^9}, { 3.7377202115121927`*^9, 3.7377202116963243`*^9}, {3.7377202481305037`*^9, 3.737720253754922*^9}, {3.73772028641499*^9, 3.7377203344757833`*^9}, { 3.737891223277732*^9, 3.7378912238571*^9}, {3.737906075544756*^9, 3.7379060758948555`*^9}, {3.737991347709439*^9, 3.7379913481113663`*^9}, { 3.7379951711905804`*^9, 3.7379951715466285`*^9}}, CellLabel->"In[20]:=",ExpressionUUID->"34558bf5-30de-4efd-a85d-1e42927d2f41"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{ RowBox[{"FindRest", "[", RowBox[{ RowBox[{"{", RowBox[{"x_", ",", "y_"}], "}"}], ",", "r_", ",", "face_", ",", "vertices_", ",", "radii_", ",", "lines_"}], "]"}], ":=", RowBox[{"(*", " ", RowBox[{ "face", " ", "is", " ", "all", " ", "the", " ", "vertices", " ", "that", " ", "the", " ", "face", " ", "intersects"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"{", RowBox[{"x_", ",", "y_"}], "}"}], ",", RowBox[{"r", ":", " ", RowBox[{ "known", " ", "center", " ", "of", " ", "opposite", " ", "kind", " ", "and", " ", "its", " ", "radius", "\[IndentingNewLine]", "face"}], ":", " ", RowBox[{ "all", " ", "the", " ", "circles", " ", "of", " ", "its", " ", "kind", " ", "that", " ", "pass", " ", "thru", " ", "that", " ", "opposite", " ", "known", " ", "circle", "\[IndentingNewLine]", RowBox[{"vertices", ":", " ", RowBox[{ "all", " ", "the", " ", "circles", " ", "of", " ", "its", " ", "kind", " ", "that", " ", "are", " ", "tangent", " ", "to", " ", "the", " ", "one", " ", RowBox[{"we", "'"}], "e", " ", "looking", " ", "for"}]}]}]}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"verts", "=", "vertices"}], ",", "i", ",", "j", ",", "\[Phi]", ",", "\[Theta]", ",", "\[Omega]", ",", RowBox[{"face2", "=", "face"}], ",", "b", ",", RowBox[{"b2", "=", "True"}]}], "}"}], ",", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"verts", "[", RowBox[{"face2", "[", RowBox[{"[", "1", "]"}], "]"}], "]"}], "\[Equal]", RowBox[{"{", RowBox[{ RowBox[{"-", "\[Infinity]"}], ",", RowBox[{"-", "\[Infinity]"}]}], "}"}]}], "||", RowBox[{"MemberQ", "[", RowBox[{"lines", ",", RowBox[{"face2", "[", RowBox[{"[", "1", "]"}], "]"}]}], "]"}]}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "2"}], ",", RowBox[{ RowBox[{"i", "\[LessEqual]", RowBox[{"Length", "[", "face2", "]"}]}], "&&", "b2"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"verts", "[", RowBox[{"face2", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], "\[NotEqual]", RowBox[{"{", RowBox[{ RowBox[{"-", "\[Infinity]"}], ",", RowBox[{"-", "\[Infinity]"}]}], "}"}]}], "&&", RowBox[{"!", RowBox[{"MemberQ", "[", RowBox[{"lines", ",", RowBox[{"face2", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"face2", "=", RowBox[{"Join", "[", RowBox[{ RowBox[{"face2", "[", RowBox[{"[", RowBox[{"i", ";;"}], "]"}], "]"}], ",", RowBox[{"face2", "[", RowBox[{"[", RowBox[{";;", RowBox[{"i", "-", "1"}]}], "]"}], "]"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"b2", "=", "False"}], ";"}]}], "\[IndentingNewLine]", "]"}]}], "]"}]}], "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"b", "=", "True"}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "2"}], ",", RowBox[{ RowBox[{"i", "\[LessEqual]", RowBox[{"Length", "[", "face2", "]"}]}], "&&", "b"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"MemberQ", "[", RowBox[{"lines", ",", RowBox[{"face2", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"AppendTo", "[", RowBox[{"face2", ",", RowBox[{"face2", "[", RowBox[{"[", "1", "]"}], "]"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", RowBox[{ RowBox[{"Length", "[", "face2", "]"}], "-", "1"}]}], ",", RowBox[{"j", ">", "i"}], ",", RowBox[{"j", "--"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"verts", "[", RowBox[{"face2", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}], "\[Equal]", RowBox[{"{", RowBox[{ RowBox[{"-", "\[Infinity]"}], ",", RowBox[{"-", "\[Infinity]"}]}], "}"}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Phi]", "=", RowBox[{"ArcTan", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"verts", "[", RowBox[{"face2", "[", RowBox[{"[", RowBox[{"j", "+", "1"}], "]"}], "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], "-", "x"}], ",", RowBox[{ RowBox[{ RowBox[{"verts", "[", RowBox[{"face2", "[", RowBox[{"[", RowBox[{"j", "+", "1"}], "]"}], "]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}], "-", "y"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"\[Theta]", "=", RowBox[{"ArcTan", "[", RowBox[{ RowBox[{"radii", "[", RowBox[{"face2", "[", RowBox[{"[", RowBox[{"j", "+", "1"}], "]"}], "]"}], "]"}], "/", "r"}], "]"}]}], ";", RowBox[{"(*", " ", RowBox[{"j", "+", RowBox[{ "1", " ", "refers", " ", "to", " ", "previously", " ", RowBox[{"computed", " ", "/", " ", "known"}]}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"\[Omega]", "=", RowBox[{"ArcTan", "[", RowBox[{ RowBox[{"radii", "[", RowBox[{"face2", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}], "/", "r"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"verts", "[", RowBox[{"face2", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}], "=", RowBox[{"{", RowBox[{ RowBox[{"x", "+", RowBox[{ SqrtBox[ RowBox[{ SuperscriptBox["r", "2"], "+", SuperscriptBox[ RowBox[{"radii", "[", RowBox[{"face2", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}], "2"]}]], RowBox[{"Cos", "[", RowBox[{"\[Phi]", "-", "\[Theta]", "-", "\[Omega]"}], "]"}]}]}], ",", RowBox[{"y", "+", RowBox[{ SqrtBox[ RowBox[{ SuperscriptBox["r", "2"], "+", SuperscriptBox[ RowBox[{"radii", "[", RowBox[{"face2", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}], "2"]}]], RowBox[{"Sin", "[", RowBox[{"\[Phi]", "-", "\[Theta]", "-", "\[Omega]"}], "]"}]}]}]}], "}"}]}]}]}], "\[IndentingNewLine]", "]"}]}], "\[IndentingNewLine]", "]"}], ";", RowBox[{"b", "=", "False"}], ";", "Break"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"verts", "[", RowBox[{"face2", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], "\[Equal]", RowBox[{"{", RowBox[{ RowBox[{"-", "\[Infinity]"}], ",", RowBox[{"-", "\[Infinity]"}]}], "}"}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Phi]", "=", RowBox[{"ArcTan", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"verts", "[", RowBox[{"face2", "[", RowBox[{"[", RowBox[{"i", "-", "1"}], "]"}], "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], "-", "x"}], ",", RowBox[{ RowBox[{ RowBox[{"verts", "[", RowBox[{"face2", "[", RowBox[{"[", RowBox[{"i", "-", "1"}], "]"}], "]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}], "-", "y"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"\[Theta]", "=", RowBox[{"ArcTan", "[", RowBox[{ RowBox[{"radii", "[", RowBox[{"face2", "[", RowBox[{"[", RowBox[{"i", "-", "1"}], "]"}], "]"}], "]"}], "/", "r"}], "]"}]}], ";", RowBox[{"(*", " ", RowBox[{"i", "-", RowBox[{ "1", " ", "refers", " ", "to", " ", "previously", " ", RowBox[{"computed", " ", "/", " ", "known"}]}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"\[Omega]", "=", RowBox[{"ArcTan", "[", RowBox[{ RowBox[{"radii", "[", RowBox[{"face2", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], "/", "r"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"verts", "[", RowBox[{"face2", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], "=", RowBox[{"{", RowBox[{ RowBox[{"x", "+", RowBox[{ SqrtBox[ RowBox[{ SuperscriptBox["r", "2"], "+", SuperscriptBox[ RowBox[{"radii", "[", RowBox[{"face2", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], "2"]}]], RowBox[{"Cos", "[", RowBox[{"\[Phi]", "+", "\[Theta]", "+", "\[Omega]"}], "]"}]}]}], ",", RowBox[{"y", "+", RowBox[{ SqrtBox[ RowBox[{ SuperscriptBox["r", "2"], "+", SuperscriptBox[ RowBox[{"radii", "[", RowBox[{"face2", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], "2"]}]], RowBox[{"Sin", "[", RowBox[{"\[Phi]", "+", "\[Theta]", "+", "\[Omega]"}], "]"}]}]}]}], "}"}]}]}]}], "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";"}]}], "]"}], ";", "\[IndentingNewLine]", "verts"}]}], "\[IndentingNewLine]", "]"}]}], ";"}], "\[IndentingNewLine]"}]], "Input", CellChangeTimes->CompressedData[" 1:eJwd030s1HEcB/BfylGhk+tmultSuqz0uKhUSsyVHk/IlWr0YDJEU5e7DCc0 xE8RihBdiLsJV1mblsgqD7XuGnk4VNfJUw/S1Pq9v3989vrn8/lsn733WRwY LjplRFHUbqagpUVmqkRr2NG6wv0a7HQYuQXPRnmXwOv+3EpokxSqhF+2f38I lZsLGqBOPdYCk/t9u2FwUvkgPLRzgQFSUv4kVJm++QOL5xwwucQoqxYT51v7 mUGnzw+JDds9raDUoOPDKKO/ofC6yjQCpi1nxcGy28OJcBEtS4NU/lAWDGfb 58Nkn7W3oLjOuQwq3aOI5W6FatgYW/ACjnq2vyLzWroT0q+tNfC4W2sXbNHv 0sGevSnERzZheii6T8+IYbR0KjSFuoxG4nBt8xGomNtOFA28CyJ9m35JoTwj RQ7NFV40XHfKIhveiPxBvLxzpADy7vNL4Pkuk2pYvLRQBV/lmb2FNQFUNzT1 2DYIs/VaImvAVw+nDyYRU/7uG4WRvMeT8MJklbGM0VbQSdwaqOmNZbyaOasP zlJXfILdnDaiE+fpnnjk4K8lardZ1EDbc5wGmDplPgrr+DQrgXHEksODD35X Ewd/S4TQ7n3FMbhJ4BcEKcd1wVB+xicUFj35FA63rF9zAXI/512CPar5V6BP 2+l00n+6loZ3J1xyoNgvpgDuize+A1ewKAVkf/xeAQUlC5TwnOZ9ExwLbP4I n2kbB+BGrnoY9ndMjcEjs4t2p+GuLuF+yB2mg+CcBMcQyBp/GQYpm3lSyLZ1 yYIeFvJiyAt9Vg4V3hNKuKiqvgmK3AJa4HOxuo3st/PvgBnCsSEYwa00wD5b q5mFyCvagQ1POKa332HcvLKmAzYFeWlgvGU6UajzHoScEGocrr4pnYbub4X/ YH1tolER/qU3h8h2dHWGkqyFG+FRL9VBqBt6RPxl3HMYTn/7EwATv+VGwH/H FkZC694l0fCrRC6BveuzY6H96ityOFsWdg+WunYr4A8zWStcddmlDZaXxq2J /sD888pcYp/m5DIpY7PxBgEUXHT7qWBMtooi/gegIJ/I "], CellLabel->"In[19]:=",ExpressionUUID->"1669407f-19e4-41f4-83d2-404c097a4b8c"], Cell[BoxData[ RowBox[{ RowBox[{"Cull", "[", "assoc_", "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"vis", "=", RowBox[{"{", "}"}]}], ",", RowBox[{"unv", "=", RowBox[{"{", "}"}]}], ",", "i"}], "}"}], ",", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", RowBox[{"Length", "[", RowBox[{"Keys", "[", "assoc", "]"}], "]"}]}], ",", RowBox[{"i", "++"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"assoc", "[", RowBox[{ RowBox[{"Keys", "[", "assoc", "]"}], "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], "\[Equal]", RowBox[{"{", RowBox[{ RowBox[{"-", "\[Infinity]"}], ",", RowBox[{"-", "\[Infinity]"}]}], "}"}]}], ",", RowBox[{"AppendTo", "[", RowBox[{"unv", ",", RowBox[{ RowBox[{"Keys", "[", "assoc", "]"}], "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}], ",", RowBox[{"AppendTo", "[", RowBox[{"vis", ",", RowBox[{ RowBox[{"Keys", "[", "assoc", "]"}], "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "]"}]}], "]"}], ";", RowBox[{"{", RowBox[{"vis", ",", "unv"}], "}"}]}]}], "]"}]}]], "Input", CellChangeTimes->{{3.737378883847493*^9, 3.7373790220179577`*^9}}, CellLabel->"In[18]:=",ExpressionUUID->"0ca661cb-0814-4baa-8045-ccb6be11fd34"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{ RowBox[{"Interior", "[", RowBox[{ "n_", ",", "FacesCW_", ",", "LineFaces_", ",", "radii_", ",", "vertices_", ",", "faces_", ",", "d_", ",", "h_"}], "]"}], ":=", "\[IndentingNewLine]", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Vunv", "=", RowBox[{"{", "}"}]}], ",", RowBox[{"Funv", "=", RowBox[{"{", "}"}]}], ",", RowBox[{"Vvis", "=", RowBox[{"{", "}"}]}], ",", RowBox[{"Fvis", "=", RowBox[{"{", "}"}]}], ",", "i", ",", "j", ",", "f", ",", RowBox[{"verts", "=", "vertices"}], ",", "verts2", ",", RowBox[{"fcs", "=", "faces"}], ",", "fcs2", ",", "fi", ",", "vi", ",", "x", ",", "y", ",", "k", ",", "b"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"f", "=", RowBox[{"Length", "[", "FacesCW", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", "n"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"vertices", "[", "i", "]"}], "\[Equal]", RowBox[{"{", RowBox[{ RowBox[{"-", "\[Infinity]"}], ",", RowBox[{"-", "\[Infinity]"}]}], "}"}]}], ",", RowBox[{"AppendTo", "[", RowBox[{"Vunv", ",", "i"}], "]"}], ",", RowBox[{"AppendTo", "[", RowBox[{"Vvis", ",", "i"}], "]"}]}], "]"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", "f"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"faces", "[", RowBox[{"FacesCW", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], "\[Equal]", RowBox[{"{", RowBox[{ RowBox[{"-", "\[Infinity]"}], ",", RowBox[{"-", "\[Infinity]"}]}], "}"}]}], ",", RowBox[{"AppendTo", "[", RowBox[{"Funv", ",", RowBox[{"FacesCW", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}], ",", RowBox[{"AppendTo", "[", RowBox[{"Fvis", ",", RowBox[{"FacesCW", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "]"}]}], "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "1"}], ",", RowBox[{ RowBox[{"j", "\[LessEqual]", RowBox[{"n", "+", RowBox[{"Length", "[", "FacesCW", "]"}]}]}], "&&", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"Length", "[", "Vunv", "]"}], "+", RowBox[{"Length", "[", "Funv", "]"}]}], ">", "0"}], ")"}]}], ",", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", RowBox[{"Length", "[", "Fvis", "]"}]}], ",", RowBox[{"i", "++"}], ",", RowBox[{"(*", " ", RowBox[{ "for", " ", "a", " ", "known", " ", "face", " ", "circle"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"Length", "[", RowBox[{"Intersection", "[", RowBox[{"Vvis", ",", RowBox[{"Fvis", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}], "]"}], "*", RowBox[{"Length", "[", RowBox[{"Intersection", "[", RowBox[{"Vunv", ",", RowBox[{"Fvis", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}], "]"}]}], ">", "0"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"fi", "=", RowBox[{"Fvis", "[", RowBox[{"[", "i", "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"{", RowBox[{"x", ",", "y"}], "}"}], "=", RowBox[{"fcs", "[", "fi", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"verts2", "=", RowBox[{"FindRest", "[", RowBox[{ RowBox[{"{", RowBox[{"x", ",", "y"}], "}"}], ",", RowBox[{"radii", "[", "fi", "]"}], ",", "fi", ",", "verts", ",", "radii", ",", RowBox[{"{", RowBox[{ RowBox[{"n", "-", "1"}], ",", "n"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"ThrowOrderError", "[", RowBox[{ "fcs", ",", "verts2", ",", "radii", ",", "d", ",", "h", ",", "n", ",", RowBox[{"{", RowBox[{ RowBox[{"FacesCW", "[", RowBox[{"[", RowBox[{"LineFaces", "[", RowBox[{"[", "1", "]"}], "]"}], "]"}], "]"}], ",", RowBox[{"FacesCW", "[", RowBox[{"[", RowBox[{"LineFaces", "[", RowBox[{"[", "2", "]"}], "]"}], "]"}], "]"}]}], "}"}], ",", "0"}], "]"}], ",", RowBox[{"verts", "=", RowBox[{"FindRest", "[", RowBox[{ RowBox[{"{", RowBox[{"x", ",", "y"}], "}"}], ",", RowBox[{"radii", "[", "fi", "]"}], ",", RowBox[{"Reverse", "[", "fi", "]"}], ",", "verts", ",", "radii", ",", RowBox[{"{", RowBox[{ RowBox[{"n", "-", "1"}], ",", "n"}], "}"}]}], "]"}]}], ",", "\[IndentingNewLine]", RowBox[{"verts", "=", "verts2"}]}], "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"{", RowBox[{"Vvis", ",", "Vunv"}], "}"}], "=", RowBox[{"Cull", "[", "verts", "]"}]}], ";"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", RowBox[{"Length", "[", "Vvis", "]"}]}], ",", RowBox[{"i", "++"}], ",", " ", RowBox[{"(*", " ", RowBox[{ "for", " ", "a", " ", "known", " ", "vertex", " ", "circle"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"b", "=", "False"}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"k", "=", "1"}], ",", RowBox[{"k", "\[LessEqual]", RowBox[{"Length", "[", "Fvis", "]"}]}], ",", RowBox[{"k", "++"}], ",", RowBox[{"b", "=", RowBox[{"b", "||", RowBox[{"If", "[", RowBox[{ RowBox[{"MemberQ", "[", RowBox[{ RowBox[{"Fvis", "[", RowBox[{"[", "k", "]"}], "]"}], ",", RowBox[{"Vvis", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}], ",", "True", ",", "False"}], "]"}]}]}]}], "]"}], ";", "b"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"b", "=", "False"}], ";", RowBox[{"For", "[", RowBox[{ RowBox[{"k", "=", "1"}], ",", RowBox[{"k", "\[LessEqual]", RowBox[{"Length", "[", "Funv", "]"}]}], ",", RowBox[{"k", "++"}], ",", RowBox[{"b", "=", RowBox[{"b", "||", RowBox[{"If", "[", RowBox[{ RowBox[{"MemberQ", "[", RowBox[{ RowBox[{"Funv", "[", RowBox[{"[", "k", "]"}], "]"}], ",", RowBox[{"Vvis", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}], ",", "True", ",", "False"}], "]"}]}]}]}], "]"}], ";", "b"}], ",", "\[IndentingNewLine]", RowBox[{"(*", " ", "then", " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"vi", "=", RowBox[{"Vvis", "[", RowBox[{"[", "i", "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"{", RowBox[{"x", ",", "y"}], "}"}], "=", RowBox[{"verts", "[", "vi", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"fcs2", "=", RowBox[{"FindRest", "[", RowBox[{ RowBox[{"{", RowBox[{"x", ",", "y"}], "}"}], ",", RowBox[{"radii", "[", "vi", "]"}], ",", RowBox[{"FacesList", "[", RowBox[{"vi", ",", "FacesCW"}], "]"}], ",", "fcs", ",", "radii", ",", RowBox[{"FacesCW", "[", RowBox[{"[", "LineFaces", "]"}], "]"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"ThrowOrderError", "[", RowBox[{ "fcs2", ",", "verts", ",", "radii", ",", "d", ",", "h", ",", "n", ",", RowBox[{"FacesCW", "[", RowBox[{"[", "LineFaces", "]"}], "]"}], ",", "0"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"fcs", "=", RowBox[{"FindRest", "[", RowBox[{ RowBox[{"{", RowBox[{"x", ",", "y"}], "}"}], ",", RowBox[{"radii", "[", "vi", "]"}], ",", "\[IndentingNewLine]", RowBox[{"Reverse", "[", RowBox[{"FacesList", "[", RowBox[{"vi", ",", "FacesCW"}], "]"}], "]"}], ",", "fcs", ",", "radii", ",", RowBox[{"FacesCW", "[", RowBox[{"[", "LineFaces", "]"}], "]"}]}], "]"}]}], ",", "\[IndentingNewLine]", RowBox[{"fcs", "=", "fcs2"}]}], "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"{", RowBox[{"Fvis", ",", "Funv"}], "}"}], "=", RowBox[{"Cull", "[", "fcs", "]"}]}], ";"}]}], "\[IndentingNewLine]", "]"}]}], "\[IndentingNewLine]", "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"verts", ",", "fcs"}], "}"}]}]}], "\[IndentingNewLine]", "]"}]}], ";"}], "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{{3.737371589519445*^9, 3.7373715939907503`*^9}, { 3.7373717937672253`*^9, 3.737371867174656*^9}, {3.737371907831069*^9, 3.73737217884262*^9}, {3.7373722186269712`*^9, 3.737372303076326*^9}, { 3.737372765799387*^9, 3.737372795651403*^9}, {3.737372858097958*^9, 3.73737290560703*^9}, {3.7373729483086967`*^9, 3.737372962791078*^9}, { 3.737376949244026*^9, 3.737377199857114*^9}, {3.737378821576833*^9, 3.737378856709831*^9}, {3.737379029954533*^9, 3.737379249070065*^9}, { 3.737379302218759*^9, 3.737379442833221*^9}, {3.737379504802657*^9, 3.737379571504554*^9}, {3.73737972110927*^9, 3.737379741421777*^9}, { 3.737380015082492*^9, 3.737380027586205*^9}, {3.737380096483328*^9, 3.737380123433268*^9}, {3.737380170681796*^9, 3.7373801708619337`*^9}, { 3.737380320948838*^9, 3.7373804262726393`*^9}, {3.73738047325154*^9, 3.737380485443796*^9}, {3.737380518592564*^9, 3.737380537424242*^9}, { 3.737380576718589*^9, 3.737380596012712*^9}, {3.7373859949503098`*^9, 3.7373860041894217`*^9}, {3.7373893401960897`*^9, 3.737389340516891*^9}, { 3.737390296785729*^9, 3.7373904116517477`*^9}, {3.7373904724338512`*^9, 3.737390578698831*^9}, {3.737390622302457*^9, 3.737390636926055*^9}, 3.737390871387472*^9, {3.737391458456609*^9, 3.7373915337313957`*^9}, { 3.7373916762920856`*^9, 3.7373916866307707`*^9}, {3.737391755289661*^9, 3.737391760387761*^9}, {3.737392424160388*^9, 3.737392446839201*^9}, { 3.7373926496569033`*^9, 3.7373926770265007`*^9}, {3.737415538959766*^9, 3.737415550423221*^9}, {3.7374157811082497`*^9, 3.737415783008522*^9}, { 3.737416229137042*^9, 3.73741626934859*^9}, {3.737416301320304*^9, 3.7374164268381977`*^9}, {3.737416488079246*^9, 3.737416506542727*^9}, { 3.737418543755864*^9, 3.737418614467947*^9}, {3.73745208819372*^9, 3.7374521079510527`*^9}, {3.737452248507903*^9, 3.737452260807358*^9}, 3.737452337464182*^9, {3.737452547756901*^9, 3.737452552175274*^9}, { 3.737452597327281*^9, 3.737452623027234*^9}, {3.737452951467626*^9, 3.7374529538565903`*^9}, {3.7374529931356688`*^9, 3.73745300450865*^9}, { 3.737453253031118*^9, 3.7374532838306513`*^9}, {3.737453448755361*^9, 3.7374535016621037`*^9}, {3.737453540364457*^9, 3.7374535737880087`*^9}, { 3.7374536292349997`*^9, 3.737453661316655*^9}, {3.737453692604967*^9, 3.737453705446292*^9}, {3.7374538230331917`*^9, 3.737453850633659*^9}, { 3.7374540473125362`*^9, 3.73745406845465*^9}, {3.737454116980228*^9, 3.737454141500236*^9}, {3.737454665545776*^9, 3.737454752199181*^9}, { 3.737454823087352*^9, 3.737454936342786*^9}, {3.73745497887639*^9, 3.7374549822716417`*^9}, {3.737455062923647*^9, 3.737455067024728*^9}, { 3.7374592709534693`*^9, 3.737459282809087*^9}, {3.737459718460556*^9, 3.7374598283067007`*^9}, {3.737459871128312*^9, 3.737459943923081*^9}, { 3.737460424909642*^9, 3.737460468107106*^9}, {3.737460531451869*^9, 3.7374605342064867`*^9}, {3.737460570370933*^9, 3.73746065914788*^9}, { 3.737469228334271*^9, 3.7374692297752113`*^9}, {3.7374759421424723`*^9, 3.737475959220735*^9}, {3.737475996982561*^9, 3.737476043275124*^9}, { 3.7374760741261263`*^9, 3.7374762044096413`*^9}, {3.737476237384112*^9, 3.737476257932431*^9}, {3.737476290722767*^9, 3.7374763178872833`*^9}, { 3.737476633241235*^9, 3.73747668584684*^9}, {3.737476775880768*^9, 3.737476807271324*^9}, {3.737476934280177*^9, 3.7374769365813828`*^9}, { 3.737476994438554*^9, 3.7374770257193117`*^9}, {3.737478615746154*^9, 3.737478664034655*^9}, {3.737478700737122*^9, 3.737478703803566*^9}, { 3.737891193943988*^9, 3.7378911945027657`*^9}, {3.7379034272313843`*^9, 3.737903448515432*^9}, 3.737903546475888*^9}, CellLabel->"In[17]:=",ExpressionUUID->"8eb99bd7-ecdf-447b-9254-91418df0046a"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"ThrowOrderError", "[", RowBox[{ "fcs2_", ",", "verts_", ",", "radii_", ",", "d_", ",", "h_", ",", "n_", ",", "fls_", ",", "\[Epsilon]_"}], "]"}], RowBox[{"(*", " ", RowBox[{"fcs2", " ", "is", " ", "proposed", " ", "fcs", " ", RowBox[{"association", ".", " ", "fls"}], " ", "is", " ", "facelines"}], " ", "*)"}], ":=", "\[IndentingNewLine]", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"i", ",", "j", ",", "b", ",", "b2", ",", "K"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"K", "=", RowBox[{"{", "}"}]}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", RowBox[{"Length", "[", RowBox[{"Keys", "[", "fcs2", "]"}], "]"}]}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"!", RowBox[{"MemberQ", "[", RowBox[{"fls", ",", RowBox[{ RowBox[{"Keys", "[", "fcs2", "]"}], "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], ",", RowBox[{"AppendTo", "[", RowBox[{"K", ",", RowBox[{ RowBox[{"Keys", "[", "fcs2", "]"}], "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "]"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"b", "=", "False"}], ";", " ", RowBox[{"(*", " ", "\"\\"", " ", "*)"}], "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{ RowBox[{"i", "\[LessEqual]", RowBox[{"Length", "[", "K", "]"}]}], "&&", RowBox[{"!", "b"}]}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"b2", "=", "False"}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"-", "\[Infinity]"}], "<", RowBox[{ RowBox[{"fcs2", "[", RowBox[{"K", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], "<", "0"}], "\[IndentingNewLine]", "||", RowBox[{ RowBox[{"-", "\[Infinity]"}], "<", RowBox[{ RowBox[{"fcs2", "[", RowBox[{"K", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}], "<", "0"}], "\[IndentingNewLine]", "||", RowBox[{ RowBox[{ RowBox[{"fcs2", "[", RowBox[{"K", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], ">", "d"}], "\[IndentingNewLine]", "||", RowBox[{ RowBox[{ RowBox[{"fcs2", "[", RowBox[{"K", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}], ">", "h"}], "\[IndentingNewLine]", "||", RowBox[{"(", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "1"}], ",", RowBox[{"j", "\[LessEqual]", RowBox[{"Length", "[", RowBox[{"K", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}], ",", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"K", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "\[NotEqual]", RowBox[{"n", "-", "1"}]}], "&&", RowBox[{ RowBox[{"K", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "\[NotEqual]", "n"}], "&&", RowBox[{ RowBox[{"verts", "[", RowBox[{"K", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "]"}], "\[NotEqual]", RowBox[{"{", RowBox[{ RowBox[{"-", "\[Infinity]"}], ",", RowBox[{"-", "\[Infinity]"}]}], "}"}]}], "&&", RowBox[{ RowBox[{"fcs2", "[", RowBox[{"K", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], "\[NotEqual]", RowBox[{"{", RowBox[{ RowBox[{"-", "\[Infinity]"}], ",", RowBox[{"-", "\[Infinity]"}]}], "}"}]}]}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"EuclideanDistance", "[", RowBox[{ RowBox[{"verts", "[", RowBox[{"K", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "]"}], ",", RowBox[{"fcs2", "[", RowBox[{"K", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}], "]"}], "\[GreaterEqual]", RowBox[{ RowBox[{"(", RowBox[{"1", "-", "\[Epsilon]"}], ")"}], RowBox[{"(", RowBox[{ RowBox[{"radii", "[", RowBox[{"K", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "]"}], "+", RowBox[{"radii", "[", RowBox[{"K", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}], ")"}]}]}], ",", RowBox[{ RowBox[{"b2", "=", "True"}], ";", "Break"}]}], "]"}]}], "]"}]}], "]"}], ";", "b2"}], ")"}]}]}], ",", "\[IndentingNewLine]", RowBox[{"b", "=", "True"}], ",", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "1"}], ",", RowBox[{ RowBox[{"j", "\[LessEqual]", RowBox[{"Length", "[", "K", "]"}]}], "&&", RowBox[{"!", "b"}]}], ",", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"i", "\[NotEqual]", "j"}], "&&", RowBox[{ RowBox[{"fcs2", "[", RowBox[{"K", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], "\[NotEqual]", RowBox[{"{", RowBox[{ RowBox[{"-", "\[Infinity]"}], ",", RowBox[{"-", "\[Infinity]"}]}], "}"}]}], "&&", RowBox[{ RowBox[{"fcs2", "[", RowBox[{"K", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}], "\[NotEqual]", RowBox[{"{", RowBox[{ RowBox[{"-", "\[Infinity]"}], ",", RowBox[{"-", "\[Infinity]"}]}], "}"}]}], "&&", "\[IndentingNewLine]", RowBox[{ RowBox[{"EuclideanDistance", "[", RowBox[{ RowBox[{"fcs2", "[", RowBox[{"K", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], ",", RowBox[{"fcs2", "[", RowBox[{"K", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}]}], "]"}], "<", RowBox[{ RowBox[{"(", RowBox[{"1", "+", "\[Epsilon]"}], ")"}], RowBox[{"(", RowBox[{ RowBox[{"radii", "[", RowBox[{"fcs2", "[", RowBox[{"K", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], "]"}], "+", RowBox[{"radii", "[", RowBox[{"fcs2", "[", RowBox[{"K", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}], "]"}]}], ")"}]}]}]}], "\[IndentingNewLine]", ",", RowBox[{"(*", " ", RowBox[{"then", ",", " ", "bad"}], " ", "*)"}], " ", RowBox[{"b", "=", "True"}]}], "]"}]}], "]"}]}], "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{ RowBox[{"i", "\[LessEqual]", RowBox[{"Length", "[", RowBox[{"fls", "[", RowBox[{"[", "1", "]"}], "]"}], "]"}]}], "&&", RowBox[{"!", "b"}]}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"verts", "[", RowBox[{"fls", "[", RowBox[{"[", RowBox[{"1", ",", "i"}], "]"}], "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], "\[NotEqual]", RowBox[{ FractionBox["d", "2"], "+", RowBox[{ RowBox[{ RowBox[{"fcs2", "[", RowBox[{"fls", "[", RowBox[{"[", "1", "]"}], "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], FractionBox["d", "2"]}]}]}], "&&", RowBox[{ RowBox[{"fls", "[", RowBox[{"[", RowBox[{"1", ",", "i"}], "]"}], "]"}], "<", RowBox[{"n", "-", "1"}]}]}], ",", RowBox[{"b", "=", "True"}]}], "]"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{ RowBox[{"i", "\[LessEqual]", RowBox[{"Length", "[", RowBox[{"fls", "[", RowBox[{"[", "2", "]"}], "]"}], "]"}]}], "&&", RowBox[{"!", "b"}]}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"verts", "[", RowBox[{"fls", "[", RowBox[{"[", RowBox[{"2", ",", "i"}], "]"}], "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], "\[NotEqual]", RowBox[{ FractionBox["d", "2"], "+", RowBox[{ RowBox[{ RowBox[{"fcs2", "[", RowBox[{"fls", "[", RowBox[{"[", "2", "]"}], "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], FractionBox["d", "2"]}]}]}], "&&", RowBox[{ RowBox[{"fls", "[", RowBox[{"[", RowBox[{"2", ",", "i"}], "]"}], "]"}], "<", RowBox[{"n", "-", "1"}]}]}], ",", RowBox[{"b", "=", "True"}]}], "]"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", "b"}]}], "\[IndentingNewLine]", "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.7373860513513117`*^9, 3.737386324009047*^9}, { 3.7373863565686207`*^9, 3.737386763058362*^9}, {3.737387011692689*^9, 3.7373874131979523`*^9}, {3.737387444285528*^9, 3.737387446417802*^9}, { 3.737387501724284*^9, 3.7373875162781897`*^9}, {3.7373875865809803`*^9, 3.7373876073787117`*^9}, {3.737387674890192*^9, 3.737387682382772*^9}, { 3.737387740454358*^9, 3.7373877582071867`*^9}, {3.737387870977488*^9, 3.737387912943714*^9}, {3.737388020284473*^9, 3.7373880216503363`*^9}, { 3.7373881133311872`*^9, 3.7373881231093473`*^9}, {3.737388218003031*^9, 3.737388251234787*^9}, {3.737388306336886*^9, 3.7373884325997267`*^9}, { 3.73738859353605*^9, 3.737388596101646*^9}, {3.7373886515625*^9, 3.737388731115857*^9}, {3.737388764603668*^9, 3.7373887825179358`*^9}, { 3.737388827395173*^9, 3.73738885534306*^9}, {3.7373889569580507`*^9, 3.7373891513974247`*^9}, {3.737389185971114*^9, 3.737389229300188*^9}, { 3.7373892871568117`*^9, 3.73738930336942*^9}, {3.737451960169162*^9, 3.737451962438613*^9}, {3.737452188533687*^9, 3.7374522140811863`*^9}, { 3.7374523846778708`*^9, 3.737452411688586*^9}, {3.7374524423850327`*^9, 3.7374525059763308`*^9}, {3.7374527627193413`*^9, 3.7374528272523108`*^9}, {3.7374529127507477`*^9, 3.737452943447597*^9}, { 3.737453021653625*^9, 3.737453033081242*^9}, 3.737453194007431*^9, { 3.737457949590412*^9, 3.737458106215234*^9}, {3.737458155680525*^9, 3.737458160666041*^9}, {3.7374582249762278`*^9, 3.737458233452886*^9}, 3.7379035853422227`*^9, {3.7379952037395506`*^9, 3.7379952046351557`*^9}}, CellLabel->"In[16]:=",ExpressionUUID->"df3b60f9-0301-4a00-901e-3f42ad1f6dcc"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"FacesList", "[", RowBox[{"v_", ",", "Faces_"}], "]"}], ":=", "\[IndentingNewLine]", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"ContainV", "=", RowBox[{"{", "}"}]}], ",", "i", ",", RowBox[{"L", "=", "2"}], ",", RowBox[{"p", "=", "3"}], ",", "temp"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", RowBox[{"Length", "[", "Faces", "]"}]}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"MemberQ", "[", RowBox[{ RowBox[{"Faces", "[", RowBox[{"[", "i", "]"}], "]"}], ",", "v"}], "]"}], ",", RowBox[{"AppendTo", "[", RowBox[{"ContainV", ",", RowBox[{"Faces", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "]"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"While", "[", RowBox[{ RowBox[{"L", "<", RowBox[{"Length", "[", "ContainV", "]"}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"While", "[", RowBox[{ RowBox[{"p", "\[LessEqual]", RowBox[{"Length", "[", "ContainV", "]"}]}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Length", "[", RowBox[{"Intersection", "[", RowBox[{ RowBox[{"ContainV", "[", RowBox[{"[", RowBox[{"L", "-", "1"}], "]"}], "]"}], ",", RowBox[{"ContainV", "[", RowBox[{"[", "p", "]"}], "]"}]}], "]"}], "]"}], "\[Equal]", "2"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"temp", "=", RowBox[{"ContainV", "[", RowBox[{"[", "p", "]"}], "]"}]}], ";", RowBox[{ RowBox[{"ContainV", "[", RowBox[{"[", "p", "]"}], "]"}], "=", RowBox[{"ContainV", "[", RowBox[{"[", "L", "]"}], "]"}]}], ";", RowBox[{ RowBox[{"ContainV", "[", RowBox[{"[", "L", "]"}], "]"}], "=", "temp"}], ";", RowBox[{"p", "=", RowBox[{"1", "+", RowBox[{"Length", "[", "ContainV", "]"}]}]}]}], ",", "\[IndentingNewLine]", RowBox[{"p", "++"}]}], "]"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", RowBox[{"L", "++"}], ";", "\[IndentingNewLine]", RowBox[{"p", "=", RowBox[{"L", "+", "1"}]}], ";"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", "ContainV"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.7373893963188553`*^9, 3.737389499461423*^9}, { 3.737389625681715*^9, 3.7373896502694798`*^9}, {3.7373897037433033`*^9, 3.737389735684783*^9}, {3.737389775354559*^9, 3.737389877602042*^9}, { 3.737389938671054*^9, 3.737389955021085*^9}, {3.737390098064752*^9, 3.7373901898598967`*^9}, {3.7373902333090754`*^9, 3.737390260404229*^9}, { 3.737477725750886*^9, 3.7374777259423532`*^9}, {3.737477772140749*^9, 3.737477809678876*^9}, {3.737477865110086*^9, 3.7374778881276903`*^9}, { 3.7374781429692793`*^9, 3.737478158525076*^9}, {3.737478266335372*^9, 3.737478275141066*^9}, {3.7374783969519587`*^9, 3.7374784008097754`*^9}, { 3.737903596846489*^9, 3.737903597250413*^9}}, CellLabel->"In[15]:=",ExpressionUUID->"94e58c8f-e9bc-4bc5-9f72-1e7d65b723b4"], Cell[BoxData[ RowBox[{ RowBox[{"Q", "=", RowBox[{"(", GridBox[{ {"0", RowBox[{"1", "/", "2"}], "0", "0"}, { RowBox[{"1", "/", "2"}], "0", "0", "0"}, {"0", "0", RowBox[{"-", "1"}], "0"}, {"0", "0", "0", RowBox[{"-", "1"}]} }], ")"}]}], ";"}]], "Input", CellLabel->"In[14]:=",ExpressionUUID->"e0269727-13f8-41e1-9bed-e3d214f9ade5"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"InvCoords", "[", RowBox[{"z_", ",", "r_", ",", RowBox[{"d_:", "0"}]}], "]"}], ":=", RowBox[{"If", "[", RowBox[{ RowBox[{"r", "\[Equal]", "\[Infinity]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"2", "d"}], ",", "0", ",", "z"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ FractionBox[ SuperscriptBox[ RowBox[{"Norm", "[", "z", "]"}], "2"], "r"], "-", "r"}], ",", RowBox[{"1", "/", "r"}], ",", RowBox[{"z", "/", "r"}]}], "}"}]}], "]"}]}], ";"}]], "Input", CellLabel->"In[13]:=",ExpressionUUID->"335fc47c-61a3-4d42-b8dc-1a6e429b3dbb"], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"Gram", "[", "SC_", "]"}], ":=", RowBox[{"SC", ".", "Q", ".", RowBox[{"Transpose", "[", "SC", "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"SuperCluster", "[", RowBox[{ "vertices_", ",", "faces_", ",", "radii_", ",", "fls_", ",", "d_", ",", "h_", ",", "n_"}], "]"}], ":=", "\[IndentingNewLine]", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"sc", "=", RowBox[{"{", "}"}]}], ",", "i"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", RowBox[{"n", "-", "2"}]}], ",", RowBox[{"i", "++"}], ",", RowBox[{"AppendTo", "[", RowBox[{"sc", ",", RowBox[{"Flatten", "[", RowBox[{"InvCoords", "[", RowBox[{ RowBox[{"vertices", "[", "i", "]"}], ",", RowBox[{"radii", "[", "i", "]"}]}], "]"}], "]"}]}], "]"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"AppendTo", "[", RowBox[{"sc", ",", RowBox[{"Flatten", "[", RowBox[{"InvCoords", "[", RowBox[{ RowBox[{"vertices", "[", RowBox[{"n", "-", "1"}], "]"}], ",", "\[Infinity]", ",", RowBox[{ FractionBox["h", "2"], "+", RowBox[{ RowBox[{ RowBox[{"vertices", "[", RowBox[{"n", "-", "1"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}], FractionBox["h", "2"]}]}]}], "]"}], "]"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"AppendTo", "[", RowBox[{"sc", ",", RowBox[{"Flatten", "[", RowBox[{"InvCoords", "[", RowBox[{ RowBox[{"vertices", "[", "n", "]"}], ",", "\[Infinity]", ",", RowBox[{ FractionBox["h", "2"], "+", RowBox[{ RowBox[{ RowBox[{"vertices", "[", "n", "]"}], "[", RowBox[{"[", "2", "]"}], "]"}], FractionBox["h", "2"]}]}]}], "]"}], "]"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", RowBox[{"Length", "[", RowBox[{"Keys", "[", "faces", "]"}], "]"}]}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"!", RowBox[{"MemberQ", "[", RowBox[{"fls", ",", RowBox[{ RowBox[{"Keys", "[", "faces", "]"}], "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], ",", "\[IndentingNewLine]", RowBox[{"AppendTo", "[", RowBox[{"sc", ",", RowBox[{"Flatten", "[", RowBox[{"InvCoords", "[", RowBox[{ RowBox[{"faces", "[", RowBox[{ RowBox[{"Keys", "[", "faces", "]"}], "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], ",", RowBox[{"radii", "[", RowBox[{ RowBox[{"Keys", "[", "faces", "]"}], "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}], "]"}], "]"}]}], "]"}]}], "\[IndentingNewLine]", "]"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"AppendTo", "[", RowBox[{"sc", ",", RowBox[{"Flatten", "[", RowBox[{"InvCoords", "[", RowBox[{ RowBox[{"faces", "[", RowBox[{"fls", "[", RowBox[{"[", "1", "]"}], "]"}], "]"}], ",", "\[Infinity]", ",", RowBox[{ FractionBox["d", "2"], "+", RowBox[{ RowBox[{ RowBox[{"faces", "[", RowBox[{"fls", "[", RowBox[{"[", "1", "]"}], "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], FractionBox["d", "2"]}]}]}], "]"}], "]"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"AppendTo", "[", RowBox[{"sc", ",", RowBox[{"Flatten", "[", RowBox[{"InvCoords", "[", RowBox[{ RowBox[{"faces", "[", RowBox[{"fls", "[", RowBox[{"[", "2", "]"}], "]"}], "]"}], ",", "\[Infinity]", ",", RowBox[{ FractionBox["d", "2"], "+", RowBox[{ RowBox[{ RowBox[{"faces", "[", RowBox[{"fls", "[", RowBox[{"[", "2", "]"}], "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], FractionBox["d", "2"]}]}]}], "]"}], "]"}]}], "]"}], ";", "\[IndentingNewLine]", "sc"}]}], "\[IndentingNewLine]", "]"}]}], ";"}]}], "Input", CellLabel->"In[11]:=",ExpressionUUID->"8f83b82e-c1db-40d9-ac4c-a723737f7704"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"DrawCircle", "[", RowBox[{"{", RowBox[{"bhat_", ",", "b_", ",", "bx_", ",", "by_"}], "}"}], "]"}], ":=", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Abs", "[", "b", "]"}], "\[GreaterEqual]", " ", ".01"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"bx", ",", "by"}], "}"}], "/", "b"}], ",", RowBox[{"1", "/", "b"}]}], "]"}], ",", RowBox[{"InfiniteLine", "[", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"bx", ",", "by"}], "}"}], RowBox[{"bhat", "/", "2"}]}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "by"}], ",", "bx"}], "}"}]}], "]"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.7383666121648498`*^9, 3.7383666304241414`*^9}, { 3.738366663056327*^9, 3.738366710683964*^9}, {3.7383669102892475`*^9, 3.73836691462473*^9}}, CellLabel->"In[10]:=",ExpressionUUID->"d37cc6ff-e0b5-488e-9c86-c606b2039282"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"Packing", "[", RowBox[{"n_", ",", "FacesCW_", ",", RowBox[{"\[Epsilon]_:", "0"}]}], "]"}], ":=", "\[IndentingNewLine]", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{ "radii", ",", "LineFaces", ",", "vertices", ",", "faces", ",", "d", ",", "h"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"radii", ",", "LineFaces"}], "}"}], "=", RowBox[{"EnergySetup", "[", RowBox[{"n", ",", "FacesCW"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"{", RowBox[{"vertices", ",", "faces", ",", "d", ",", "h"}], "}"}], "=", RowBox[{"Boundary", "[", RowBox[{ "n", ",", "FacesCW", ",", "LineFaces", ",", "radii", ",", "\[Epsilon]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"{", RowBox[{"vertices", ",", "faces"}], "}"}], "=", RowBox[{"Interior", "[", RowBox[{ "n", ",", "FacesCW", ",", "LineFaces", ",", "radii", ",", "vertices", ",", "faces", ",", "d", ",", "h"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"{", RowBox[{ "vertices", ",", "faces", ",", "radii", ",", "d", ",", "h", ",", "LineFaces"}], "}"}]}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.737469007369687*^9, 3.737469107407995*^9}, { 3.73746916343927*^9, 3.737469194858486*^9}, {3.737469239984439*^9, 3.737469246965393*^9}, {3.737469315964088*^9, 3.737469336634994*^9}, { 3.737469428057992*^9, 3.737469432696295*^9}, {3.737469465014307*^9, 3.737469481114959*^9}, {3.7374695687460957`*^9, 3.7374696091341476`*^9}, { 3.737470455200612*^9, 3.7374704661283417`*^9}, {3.737478361306628*^9, 3.7374783619053802`*^9}, {3.737720455890616*^9, 3.737720462414033*^9}}, CellLabel->"In[9]:=",ExpressionUUID->"9a3eb7b4-a69e-4770-9009-32f6122d8145"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"Render", "[", RowBox[{ "radii_", ",", "fcs_", ",", "verts_", ",", "fls_", ",", "n_", ",", "d_", ",", "h_"}], "]"}], ":=", RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{"Blue", ",", "Thick", ",", "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"i", "\[Equal]", RowBox[{"n", "-", "1"}]}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", RowBox[{"Max", "[", RowBox[{"Values", "[", "radii", "]"}], "]"}]}], ",", "h"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"d", "+", RowBox[{"Max", "[", RowBox[{"Values", "[", "radii", "]"}], "]"}]}], ",", "h"}], "}"}]}], "}"}], "]"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{"i", "\[Equal]", "n"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", RowBox[{"Max", "[", RowBox[{"Values", "[", "radii", "]"}], "]"}]}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"d", "+", RowBox[{"Max", "[", RowBox[{"Values", "[", "radii", "]"}], "]"}]}], ",", "0"}], "}"}]}], "}"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Circle", "[", RowBox[{ RowBox[{"verts", "[", RowBox[{ RowBox[{"Keys", "[", "verts", "]"}], "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], ",", RowBox[{"radii", "[", RowBox[{ RowBox[{"Keys", "[", "verts", "]"}], "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}], "]"}]}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}], ",", "\[IndentingNewLine]", "Red", ",", "Dashed", ",", "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"Keys", "[", "fcs", "]"}], "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", RowBox[{"fls", "[", RowBox[{"[", "1", "]"}], "]"}]}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", RowBox[{"-", RowBox[{"Max", "[", RowBox[{"Values", "[", "radii", "]"}], "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", RowBox[{"h", "+", RowBox[{"Max", "[", RowBox[{"Values", "[", "radii", "]"}], "]"}]}]}], "}"}]}], "}"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"Keys", "[", "fcs", "]"}], "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", RowBox[{"fls", "[", RowBox[{"[", "2", "]"}], "]"}]}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"d", ",", RowBox[{"-", RowBox[{"Max", "[", RowBox[{"Values", "[", "radii", "]"}], "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"d", ",", RowBox[{"h", "+", RowBox[{"Max", "[", RowBox[{"Values", "[", "radii", "]"}], "]"}]}]}], "}"}]}], "}"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Circle", "[", RowBox[{ RowBox[{"fcs", "[", RowBox[{ RowBox[{"Keys", "[", "fcs", "]"}], "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], ",", RowBox[{"radii", "[", RowBox[{ RowBox[{"Keys", "[", "fcs", "]"}], "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}], "]"}]}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"Length", "[", RowBox[{"Keys", "[", "fcs", "]"}], "]"}]}], "}"}]}], "]"}]}], "\[IndentingNewLine]", "}"}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.7373911026555433`*^9, 3.737391276476068*^9}, { 3.737391324439508*^9, 3.737391325806115*^9}, {3.7373913655803022`*^9, 3.73739137356105*^9}, {3.737460834501174*^9, 3.7374608791258917`*^9}, { 3.737460909495472*^9, 3.737460979640925*^9}, {3.737469858560081*^9, 3.737469918056663*^9}, {3.73746995156046*^9, 3.7374701969102716`*^9}, { 3.737470230362206*^9, 3.737470299948655*^9}, 3.7378928141832247`*^9}, CellLabel->"In[8]:=",ExpressionUUID->"affbea28-f93e-449d-873a-65b52e624ef9"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{ RowBox[{"SendToInt", "[", RowBox[{"\[Epsilon]_:", "0"}], "]"}], "[", "y_", "]"}], ":=", RowBox[{"With", "[", RowBox[{ RowBox[{"{", RowBox[{"x", "=", RowBox[{"Simplify", "[", "y", "]"}]}], "}"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Abs", "[", RowBox[{"x", "-", RowBox[{"Floor", "[", "x", "]"}]}], "]"}], "\[LessEqual]", "\[Epsilon]"}], " ", ",", RowBox[{"Floor", "[", "x", "]"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Abs", "[", RowBox[{"x", "-", RowBox[{"Ceiling", "[", "x", "]"}]}], "]"}], "\[LessEqual]", " ", "\[Epsilon]"}], " ", ",", RowBox[{"Ceiling", "[", "x", "]"}], ",", "x"}], "]"}], ",", "x"}], "]"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.7384110559425135`*^9, 3.7384110603257904`*^9}, { 3.7384111195493803`*^9, 3.7384111627798147`*^9}, {3.73887183730313*^9, 3.738871862607346*^9}}, CellLabel->"In[69]:=",ExpressionUUID->"15495297-a588-42f2-9830-ff3a5a9d76a4"], Cell[BoxData[{ RowBox[{ RowBox[{"R", "[", "vhat_", "]"}], ":=", RowBox[{ RowBox[{"IdentityMatrix", "[", "4", "]"}], "+", RowBox[{"2", RowBox[{"Q", ".", RowBox[{"Transpose", "[", RowBox[{"{", "vhat", "}"}], "]"}], ".", RowBox[{"{", "vhat", "}"}]}]}]}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"BMat", "[", "n_", "]"}], ":=", RowBox[{"Table", "[", RowBox[{ RowBox[{"b", "[", RowBox[{"i", ",", "j"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"Bend", "[", RowBox[{"V_", ",", "m_"}], "]"}], "[", "i_", "]"}], ":=", "\[IndentingNewLine]", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{ "sol", ",", "assignments", ",", "rest", ",", "l", ",", "others1", ",", "others2"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"sol", "=", RowBox[{ RowBox[{"Solve", "[", RowBox[{ RowBox[{ RowBox[{"BMat", "[", "m", "]"}], ".", RowBox[{"(", RowBox[{"V", "[", RowBox[{"[", RowBox[{";;", "m"}], "]"}], "]"}], ")"}]}], "\[Equal]", RowBox[{ RowBox[{"(", RowBox[{"V", "[", RowBox[{"[", RowBox[{";;", "m"}], "]"}], "]"}], ")"}], ".", RowBox[{"R", "[", RowBox[{"V", "[", RowBox[{"[", RowBox[{"m", "+", "i"}], "]"}], "]"}], "]"}]}]}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"assignments", "=", RowBox[{"Association", "[", "sol", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"rest", "=", RowBox[{"Complement", "[", RowBox[{ RowBox[{"Flatten", "[", RowBox[{"BMat", "[", "m", "]"}], "]"}], ",", RowBox[{"Keys", "[", "assignments", "]"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"l", "=", RowBox[{"Length", "[", "rest", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"others1", "=", RowBox[{ RowBox[{"Solve", "[", RowBox[{"rest", "\[Equal]", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "l"}], "}"}]}], "]"}]}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"others2", "=", RowBox[{"<|", "|>"}]}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "1"}], ",", RowBox[{"j", "\[LessEqual]", "l"}], ",", RowBox[{"j", "++"}], ",", RowBox[{ RowBox[{"others2", "[", RowBox[{"rest", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}], "=", SubscriptBox["a", "j"]}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"others2", "=", RowBox[{"Normal", "[", "others2", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Column", "[", RowBox[{"{", "\[IndentingNewLine]", RowBox[{ RowBox[{"MatrixForm", "[", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"BMat", "[", "m", "]"}], "/.", "assignments"}], ")"}], "/.", "others2"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"MatrixForm", "[", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"BMat", "[", "m", "]"}], "/.", "assignments"}], ")"}], "/.", "others1"}], "]"}]}], "\[IndentingNewLine]", "}"}], "]"}]}]}], "\[IndentingNewLine]", "]"}]}], ";"}]}], "Input", CellChangeTimes->{{3.7378245147847567`*^9, 3.737824531315222*^9}, { 3.737995674554207*^9, 3.7379956992022705`*^9}, 3.7380623568729105`*^9}, CellLabel->"In[4]:=",ExpressionUUID->"cd35cdcc-79cb-420e-8b24-2c510cf42e01"], Cell[BoxData[ RowBox[{ RowBox[{"Everything", "[", RowBox[{"n_", ",", "FacesCW_", ",", RowBox[{"\[Epsilon]_:", "0"}], ",", "steps_", ",", RowBox[{"sizemax_:", "\[Infinity]"}]}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"vertices", ",", "faces", ",", RowBox[{"f", "=", RowBox[{"Length", "[", "FacesCW", "]"}]}], ",", "radii", ",", "d", ",", "h", ",", "LineFaces", ",", "Vnye", ",", "Vout", ",", RowBox[{"bends", "=", RowBox[{"{", "}"}]}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{ "vertices", ",", "faces", ",", "radii", ",", "d", ",", "h", ",", "LineFaces"}], "}"}], "=", RowBox[{"Packing", "[", RowBox[{"n", ",", "FacesCW", ",", "\[Epsilon]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Vnye", "=", RowBox[{ RowBox[{"SendToInt", "[", "\[Epsilon]", "]"}], "//@", RowBox[{"SuperCluster", "[", RowBox[{"vertices", ",", "faces", ",", "radii", ",", RowBox[{"FacesCW", "[", RowBox[{"[", "LineFaces", "]"}], "]"}], ",", "d", ",", "h", ",", "n"}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Vout", "=", RowBox[{"MatrixForm", "[", RowBox[{"RootApproximant", "[", "Vnye", "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"bends", "=", RowBox[{ RowBox[{"Bend", "[", RowBox[{ RowBox[{"Vout", "[", RowBox[{"[", "1", "]"}], "]"}], ",", "n"}], "]"}], "/@", RowBox[{"Range", "@", "f"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Column", "[", RowBox[{"{", "\[IndentingNewLine]", RowBox[{ RowBox[{"Render", "[", RowBox[{"radii", ",", "faces", ",", "vertices", ",", RowBox[{"FacesCW", "[", RowBox[{"[", "LineFaces", "]"}], "]"}], ",", "n", ",", "d", ",", "h"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Graphics", "[", RowBox[{ RowBox[{"{", RowBox[{"Thick", ",", " ", RowBox[{"DrawCircle", "/@", RowBox[{"Orbit", "[", RowBox[{ RowBox[{"Vnye", "[", RowBox[{"[", RowBox[{ RowBox[{"-", "f"}], ";;"}], "]"}], "]"}], ",", RowBox[{"Vnye", "[", RowBox[{"[", RowBox[{";;", "n"}], "]"}], "]"}], ",", "steps", ",", "sizemax"}], "]"}]}], ",", "Red", ",", RowBox[{"DrawCircle", "/@", RowBox[{"Vnye", "[", RowBox[{"[", RowBox[{ RowBox[{"-", "f"}], ";;"}], "]"}], "]"}]}], ",", " ", "Blue", ",", " ", RowBox[{"DrawCircle", "/@", RowBox[{"Vnye", "[", RowBox[{"[", RowBox[{";;", "n"}], "]"}], "]"}]}]}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "d"}], ",", RowBox[{"2", " ", "d"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "1.2"}], " ", RowBox[{"Max", "[", RowBox[{"Values", "[", "radii", "]"}], "]"}]}], ",", RowBox[{"h", "+", RowBox[{"1.2", " ", RowBox[{"Max", "[", RowBox[{"Values", "[", "radii", "]"}], "]"}]}]}]}], "}"}]}], "}"}]}], ",", RowBox[{"PlotRangeClipping", "->", "True"}]}], "]"}], ",", "\[IndentingNewLine]", "Vout", ",", "\[IndentingNewLine]", "bends", ",", "\[IndentingNewLine]", RowBox[{"MatrixForm", "[", RowBox[{"RootApproximant", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"SendToInt", "[", "\[Epsilon]", "]"}], "//@", RowBox[{"Gram", "[", "Vnye", "]"}]}], "\[IndentingNewLine]", "]"}], "]"}]}], "\[IndentingNewLine]", "}"}], "]"}]}]}], "]"}]}]], "Input",\ CellChangeTimes->CompressedData[" 1:eJwdy38s1HEcx/ET/bo1h0qGWoc4M86P/GixuvKrzrQ5P/Iz0TFjyK8hP2pZ LmuhlI7ralwhkuqY8qPD+XVUWml+niJtnKn86seovq/PH+899t77/WRGJPjw N9FoNDdqYOcTi8bqERWnZr6jFZqFeHVCbQfJB8hSk1rXUBbmrtjB4dC9QbCd XhEC2beKEuGNxy3pcCAg3pA5quLoxYr2wYkRyVsjyqXsj0QTe36pH+WjyMU7 cEg2djiN8v3OS8ehdOSWN5QrhnnwjUzrLNm7mESJNCMwm7KOn0d0fpDaWkXp vucJMU3njBIalEo+QWV6yAK524gXoZ5vg6+CUkvL2A/Ss5Kj4G2FTzTs7HfM gCzNmWzYXyjNg4KtM0RaZd9VGK8UFUCp2XwZdLa4J4bnhNznMMYxoRmWCFPb SM/r6YC6c7m9sLhKQBT+WVfC+QL6FMy3r/4K3ZXcBWgZJv4BJTlLxMFuxhqc vlJO9JPZqA1Qhp/OIwarfhlAxewhQyiObdkP5eopTNikEcWG6qFxtjAs4KET 6V+yOdC2eJILMzWiT8E6YbA/bObFB0K6rnFk2JiKY9UoioYlK5oJMOfygURI kxemwFouPw1ayz9nwiH9ZaKoQdFznpJ5obUfWm3whyDLLvIdPD05tJpEeTDm 6FoS6c3XYXUtm5jiEaqXTCmYriemZ8g94OAxU0840Gu9ABlzgu+wVOGikULp L6wkilyV+rBMPk3MG84yhpF1LUSjegYLlivUzcmff58l9BUZs6HMxMoFhrsu HYEOzyzcyd+bCG8YNTjhA1lBDB48UdsVAm8HPwiDXzbbxEPnkw5E0U+auJ/y 4bftxG12kw1w1+5pYjtjJNN3XMXRYXtmQ+uNnHzY3qYtgKtT6qXw5ooTUfdp bTmscRu9C/9NeL0gfdE7YnnF605Iv/+7G84MKMahzavWCXjRM24RXr+WSYzI bVqGszK1NfjXcouOH+Ua15Co2tHAg6aclgD4HwSTeW4= "], CellLabel->"In[75]:=",ExpressionUUID->"b49300f9-4eb0-4e74-b00e-238ad2cc61d5"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"Orbit", "[", RowBox[{"co_", ",", "cl_", ",", "steps_", ",", RowBox[{"sizemax_:", "\[Infinity]"}]}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"i", ",", "j", ",", "k", ",", RowBox[{"orbit", "=", "cl"}], ",", "L"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", "steps"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"L", "=", RowBox[{"Length", "[", "orbit", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "1"}], ",", RowBox[{"j", "\[LessEqual]", RowBox[{"Length", "[", "co", "]"}]}], ",", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"k", "=", "1"}], ",", RowBox[{"k", "\[LessEqual]", "L"}], ",", RowBox[{"k", "++"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ "orbit", "\[LeftDoubleBracket]", "k", "\[RightDoubleBracket]"}], ".", RowBox[{"R", "[", RowBox[{ "co", "\[LeftDoubleBracket]", "j", "\[RightDoubleBracket]"}], "]"}]}], ")"}], "\[LeftDoubleBracket]", "2", "\[RightDoubleBracket]"}], "\[LessEqual]", "sizemax"}], "&&", RowBox[{"!", RowBox[{"MemberQ", "[", RowBox[{ RowBox[{"Round", "[", RowBox[{ RowBox[{"N", "[", "orbit", "]"}], ",", ".01"}], "]"}], ",", RowBox[{"Round", "[", RowBox[{ RowBox[{"N", "[", RowBox[{ RowBox[{ "orbit", "\[LeftDoubleBracket]", "k", "\[RightDoubleBracket]"}], ".", RowBox[{"R", "[", RowBox[{ "co", "\[LeftDoubleBracket]", "j", "\[RightDoubleBracket]"}], "]"}]}], "]"}], ",", ".01"}], "]"}]}], RowBox[{"(*", RowBox[{"new", " ", "thing"}], "*)"}], "]"}]}]}], ",", RowBox[{"AppendTo", "[", RowBox[{"orbit", ",", RowBox[{ RowBox[{ "orbit", "\[LeftDoubleBracket]", "k", "\[RightDoubleBracket]"}], ".", RowBox[{"R", "[", RowBox[{ "co", "\[LeftDoubleBracket]", "j", "\[RightDoubleBracket]"}], "]"}]}]}], RowBox[{"(*", "newthing", "*)"}], "]"}]}], "]"}]}], "\[IndentingNewLine]", "]"}]}], "\[IndentingNewLine]", "]"}]}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", "orbit"}]}], "\[IndentingNewLine]", "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.738062389806858*^9, 3.738062427087122*^9}, { 3.7380625273360224`*^9, 3.7380625840383434`*^9}, {3.738064608037428*^9, 3.738064608439394*^9}, {3.738408767299783*^9, 3.738408799676605*^9}, { 3.738586345119252*^9, 3.738586346798085*^9}, {3.7388709170245705`*^9, 3.7388709187843447`*^9}}, CellLabel->"In[54]:=",ExpressionUUID->"b43a5c15-ff8f-4841-a157-1d63bba9b7c4"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Everything", "[", RowBox[{"7", ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1", ",", "7"}], "}"}], ",", RowBox[{"{", RowBox[{"3", ",", "2", ",", "7"}], "}"}], ",", RowBox[{"{", RowBox[{"4", ",", "3", ",", "7"}], "}"}], ",", RowBox[{"{", RowBox[{"5", ",", "4", ",", "7"}], "}"}], ",", RowBox[{"{", RowBox[{"5", ",", "7", ",", "6"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "6", ",", "7"}], "}"}], ",", RowBox[{"{", RowBox[{"5", ",", "4", ",", "3", ",", "2", ",", "1", ",", "6"}], "}"}]}], "}"}], ",", RowBox[{"1", "/", "100"}], ",", "5", ",", "35"}], "]"}]], "Input", CellChangeTimes->{{3.737720388591385*^9, 3.7377204147139273`*^9}, { 3.73799104566405*^9, 3.737991046219566*^9}, {3.7379911247406125`*^9, 3.737991126272558*^9}, 3.737992714004908*^9, {3.7383399270780163`*^9, 3.738339928102619*^9}, {3.7383399833318253`*^9, 3.738340023435418*^9}, { 3.7383669397920876`*^9, 3.73836699772031*^9}, {3.7383676488210897`*^9, 3.738367650275752*^9}, {3.738411438453149*^9, 3.738411446817583*^9}, { 3.7384118416165533`*^9, 3.7384118445254*^9}, {3.738586039916703*^9, 3.738586105678995*^9}, {3.73887251798787*^9, 3.738872520627311*^9}}, CellLabel->"In[76]:=",ExpressionUUID->"8c0d00eb-3a38-4d87-9e21-66e5a04e235e"], Cell[BoxData[ TagBox[GridBox[{ { GraphicsBox[ {RGBColor[0, 0, 1], Thickness[ Large], {CircleBox[{3.464101615137755, 1.}], CircleBox[{2.309401076758503, 0.3333333333333335}, 0.3333333333333333], CircleBox[{1.732050807568877, 0.24999999999999994}, 0.25], CircleBox[{1.1547005383792512, 0.3333333333333332}, 0.33333333333333337], CircleBox[{0, 1.0000000000000002}, 1.0000000000000002], LineBox[{{-1.7320508075688774`, 2.}, {5.196152422706632, 2.}}], LineBox[{{-1.7320508075688774`, 0}, {5.196152422706632, 0}}]}, {RGBColor[1, 0, 0], Dashing[{Small, Small}], CircleBox[{2.8867513459481287, 0}, 0.5773502691896256], CircleBox[{2.02072594216369, 0}, 0.28867513459481287], CircleBox[{1.4433756729740643, 0}, 0.28867513459481287], CircleBox[{0.5773502691896257, 0}, 0.5773502691896257], LineBox[{{0, -1.7320508075688774`}, {0, 3.7320508075688776`}}], LineBox[{{3.4641016151377544`, -1.7320508075688774`}, { 3.4641016151377544`, 3.7320508075688776`}}], CircleBox[{1.7320508075688774, 2.0000000000000004}, 1.7320508075688774]}}]}, { GraphicsBox[ {Thickness[ Large], {CircleBox[{3.464101615137755, 1}], CircleBox[ NCache[{2.309401076758503, Rational[1, 3]}, {2.309401076758503, 0.3333333333333333}], NCache[Rational[1, 3], 0.3333333333333333]], CircleBox[ NCache[{1.732050807568877, Rational[1, 4]}, {1.732050807568877, 0.25}], NCache[Rational[1, 4], 0.25]], CircleBox[ NCache[{1.154700538379251, Rational[1, 3]}, {1.154700538379251, 0.3333333333333333}], NCache[Rational[1, 3], 0.3333333333333333]], CircleBox[{0, 1}], TagBox[ConicHullRegionBox[{0, 2}, {{-1, 0}}], "InfiniteLine"], TagBox[ConicHullRegionBox[{0, 0}, {{1, 0}}], "InfiniteLine"], CircleBox[{2.5980762113533196, 0.06249999999999983}, 0.06249999999999983], CircleBox[{2.694301256218253, 0.037037037037036986}, 0.037037037037036986], CircleBox[{2.7712812921102032, 0.03999999999999999}, 0.03999999999999999], CircleBox[{2.8867513459481287, 0.0833333333333333}, 0.0833333333333333], CircleBox[{2.0784609690826508, 0.04}, 0.04], CircleBox[{1.9245008972987527, 0.037037037037037}, 0.037037037037037], CircleBox[{1.539600717839004, 0.037037037037037035}, 0.037037037037037035], CircleBox[{1.3856406460551016, 0.04}, 0.04], CircleBox[{0.6928203230275508, 0.03999999999999999}, 0.03999999999999999], CircleBox[{0.7698003589195008, 0.03703703703703703}, 0.03703703703703703], CircleBox[{0.8660254037844383, 0.06249999999999996}, 0.06249999999999996], CircleBox[{0.5773502691896257, 0.08333333333333331}, 0.08333333333333331], CircleBox[{1.7320508075688772, 1.2500000000000004}, 0.7499999999999999], CircleBox[{-3.464101615137755, 1.}], CircleBox[{-2.309401076758503, 0.3333333333333333}, 0.3333333333333333], CircleBox[{-1.732050807568877, 0.25}, 0.25], CircleBox[{-1.154700538379251, 0.3333333333333333}, 0.3333333333333333], CircleBox[{4.618802153517006, 0.3333333333333333}, 0.3333333333333333], CircleBox[{5.196152422706632, 0.25}, 0.25], CircleBox[{5.773502691896258, 0.3333333333333333}, 0.3333333333333333], CircleBox[{6.928203230275509, 1.}], CircleBox[{2.721794126179663, 0.17857142857142835}, 0.10714285714285696], CircleBox[{3.0792014356780073, 0.03703703703703708}, 0.03703703703703708], CircleBox[{1.9985201625794748, 0.09615384615384615}, 0.05769230769230766], CircleBox[{1.4655814525582807, 0.09615384615384606}, 0.05769230769230761], CircleBox[{0.7423074889580902, 0.1785714285714285}, 0.10714285714285705], CircleBox[{0.3849001794597505, 0.03703703703703704}, 0.03703703703703704], CircleBox[{2.3094010767585065, 0.7083333333333309}, 0.04166666666666662], CircleBox[{2.424871130596426, 0.8500000000000005}, 0.04999999999999994], CircleBox[{1.9944827481096135, 0.5151515151515167}, 0.030303030303030262], CircleBox[{1.880512305360495, 0.4857142857142879}, 0.0285714285714285], CircleBox[{1.5835893097772604, 0.4857142857142842}, 0.028571428571428595], CircleBox[{1.4696188670281383, 0.5151515151515159}, 0.03030303030303029], CircleBox[{1.1547005383792508, 0.7083333333333338}, 0.0416666666666666], CircleBox[{1.0392304845413265, 0.8500000000000003}, 0.049999999999999975], CircleBox[{1.1547005383792515, 1.8888888888888893}, 0.1111111111111111], CircleBox[{1.0939268258329748, 1.736842105263158}, 0.05263157894736841], CircleBox[{1.0392304845413263, 1.65}, 0.04999999999999999], CircleBox[{0.9447549859466605, 1.5454545454545454}, 0.09090909090909091], CircleBox[{2.519346629191095, 1.5454545454545459}, 0.09090909090909095], CircleBox[{2.4248711305964274, 1.6500000000000004}, 0.04999999999999999], CircleBox[{2.3701747893047798, 1.7368421052631582}, 0.05263157894736843], CircleBox[{2.309401076758504, 1.8888888888888904}, 0.11111111111111117], CircleBox[{-2.5980762113533196, 0.06249999999999983}, 0.06249999999999983], CircleBox[{-2.694301256218253, 0.037037037037036986}, 0.037037037037036986], CircleBox[{-2.7712812921102032, 0.03999999999999999}, 0.03999999999999999], CircleBox[{-2.8867513459481287, 0.0833333333333333}, 0.0833333333333333], CircleBox[{-2.0784609690826508, 0.04}, 0.04], CircleBox[{-1.9245008972987527, 0.037037037037037}, 0.037037037037037], CircleBox[{-1.539600717839004, 0.037037037037037035}, 0.037037037037037035], CircleBox[{-1.3856406460551016, 0.04}, 0.04], CircleBox[{-0.6928203230275508, 0.03999999999999999}, 0.03999999999999999], CircleBox[{-0.7698003589195008, 0.03703703703703703}, 0.03703703703703703], CircleBox[{-0.8660254037844383, 0.06249999999999996}, 0.06249999999999996], CircleBox[{-0.5773502691896257, 0.08333333333333331}, 0.08333333333333331], CircleBox[{-1.7320508075688772, 1.2500000000000004}, 0.7499999999999999], CircleBox[{-4.618802153517006, 0.3333333333333333}, 0.3333333333333333], CircleBox[{-5.196152422706632, 0.25}, 0.25], CircleBox[{-5.773502691896258, 0.3333333333333333}, 0.3333333333333333], CircleBox[{-6.928203230275509, 1.}], CircleBox[{4.330127018922189, 0.06249999999999983}, 0.06249999999999983], CircleBox[{4.233901974057255, 0.037037037037036986}, 0.037037037037036986], CircleBox[{4.156921938165306, 0.03999999999999999}, 0.03999999999999999], CircleBox[{4.041451884327381, 0.0833333333333333}, 0.0833333333333333], CircleBox[{4.8497422611928585, 0.04}, 0.04], CircleBox[{5.003702332976756, 0.037037037037037}, 0.037037037037037], CircleBox[{5.388602512436504, 0.037037037037037035}, 0.037037037037037035], CircleBox[{5.5425625842204065, 0.04}, 0.04], CircleBox[{6.235382907247959, 0.03999999999999999}, 0.03999999999999999], CircleBox[{6.1584028713560075, 0.03703703703703703}, 0.03703703703703703], CircleBox[{6.06217782649107, 0.06249999999999996}, 0.06249999999999996], CircleBox[{6.350852961085882, 0.08333333333333331}, 0.08333333333333331], CircleBox[{5.196152422706632, 1.2500000000000004}, 0.7499999999999999], CircleBox[{10.392304845413264, 1.}], CircleBox[{9.237604307034012, 0.3333333333333333}, 0.3333333333333333], CircleBox[{8.660254037844386, 0.25}, 0.25], CircleBox[{8.08290376865476, 0.3333333333333333}, 0.3333333333333333], CircleBox[{3.008298771040681, 0.06578947368421045}, 0.03947368421052626], CircleBox[{0.45580284409707295, 0.06578947368421052}, 0.0394736842105263], CircleBox[{2.4248711305964257, 0.7250000000000002}, 0.07499999999999978], CircleBox[{1.9485571585149875, 0.45312499999999845}, 0.046875], CircleBox[{1.5155444566227678, 0.4531250000000012}, 0.0468749999999999], CircleBox[{1.0392304845413265, 0.7250000000000003}, 0.0749999999999999], CircleBox[{0.8660254037844385, 1.8125000000000002}, 0.18749999999999997], CircleBox[{1.3856406460551018, 1.9600000000000006}, 0.04000000000000001], CircleBox[{2.5980762113533156, 1.8125000000000002}, 0.18749999999999997], CircleBox[{2.078460969082652, 1.9600000000000006}, 0.039999999999999994], CircleBox[{-2.721794126179663, 0.17857142857142835}, 0.10714285714285696], CircleBox[{-3.0792014356780073, 0.03703703703703708}, 0.03703703703703708], CircleBox[{-1.9985201625794748, 0.09615384615384615}, 0.05769230769230766], CircleBox[{-1.4655814525582807, 0.09615384615384606}, 0.05769230769230761], CircleBox[{-0.7423074889580902, 0.1785714285714285}, 0.10714285714285705], CircleBox[{-0.3849001794597505, 0.03703703703703704}, 0.03703703703703704], CircleBox[{-2.3094010767585065, 0.7083333333333309}, 0.04166666666666662], CircleBox[{-2.424871130596426, 0.8500000000000005}, 0.04999999999999994], CircleBox[{-1.9944827481096135, 0.5151515151515167}, 0.030303030303030262], CircleBox[{-1.880512305360495, 0.4857142857142879}, 0.0285714285714285], CircleBox[{-1.5835893097772604, 0.4857142857142842}, 0.028571428571428595], CircleBox[{-1.4696188670281383, 0.5151515151515159}, 0.03030303030303029], CircleBox[{-1.1547005383792508, 0.7083333333333338}, 0.0416666666666666], CircleBox[{-1.0392304845413265, 0.8500000000000003}, 0.049999999999999975], CircleBox[{-1.1547005383792515, 1.8888888888888893}, 0.1111111111111111], CircleBox[{-1.0939268258329748, 1.736842105263158}, 0.05263157894736841], CircleBox[{-1.0392304845413263, 1.65}, 0.04999999999999999], CircleBox[{-0.9447549859466605, 1.5454545454545454}, 0.09090909090909091], CircleBox[{-2.519346629191095, 1.5454545454545459}, 0.09090909090909095], CircleBox[{-2.4248711305964274, 1.6500000000000004}, 0.04999999999999999], CircleBox[{-2.3701747893047798, 1.7368421052631582}, 0.05263157894736843], CircleBox[{-2.309401076758504, 1.8888888888888904}, 0.11111111111111117], CircleBox[{-4.330127018922189, 0.06249999999999983}, 0.06249999999999983], CircleBox[{-4.233901974057255, 0.037037037037036986}, 0.037037037037036986], CircleBox[{-4.156921938165306, 0.03999999999999999}, 0.03999999999999999], CircleBox[{-4.041451884327381, 0.0833333333333333}, 0.0833333333333333], CircleBox[{-4.8497422611928585, 0.04}, 0.04], CircleBox[{-5.003702332976756, 0.037037037037037}, 0.037037037037037], CircleBox[{-5.388602512436504, 0.037037037037037035}, 0.037037037037037035], CircleBox[{-5.5425625842204065, 0.04}, 0.04], CircleBox[{-6.235382907247959, 0.03999999999999999}, 0.03999999999999999], CircleBox[{-6.1584028713560075, 0.03703703703703703}, 0.03703703703703703], CircleBox[{-6.06217782649107, 0.06249999999999996}, 0.06249999999999996], CircleBox[{-6.350852961085882, 0.08333333333333331}, 0.08333333333333331], CircleBox[{-5.196152422706632, 1.2500000000000004}, 0.7499999999999999], CircleBox[{-10.392304845413264, 1.}], CircleBox[{-9.237604307034012, 0.3333333333333333}, 0.3333333333333333], CircleBox[{-8.660254037844386, 0.25}, 0.25], CircleBox[{-8.08290376865476, 0.3333333333333333}, 0.3333333333333333], CircleBox[{4.2064091040958465, 0.17857142857142835}, 0.10714285714285696], CircleBox[{3.849001794597502, 0.03703703703703708}, 0.03703703703703708], CircleBox[{4.9296830676960335, 0.09615384615384615}, 0.05769230769230766], CircleBox[{5.462621777717228, 0.09615384615384606}, 0.05769230769230761], CircleBox[{6.1858957413174185, 0.1785714285714285}, 0.10714285714285705], CircleBox[{6.543303050815759, 0.03703703703703704}, 0.03703703703703704], CircleBox[{4.618802153517003, 0.7083333333333309}, 0.04166666666666662], CircleBox[{4.503332099679082, 0.8500000000000005}, 0.04999999999999994], CircleBox[{4.933720482165896, 0.5151515151515167}, 0.030303030303030262], CircleBox[{5.047690924915014, 0.4857142857142879}, 0.0285714285714285], CircleBox[{5.344613920498248, 0.4857142857142842}, 0.028571428571428595], CircleBox[{5.458584363247371, 0.5151515151515159}, 0.03030303030303029], CircleBox[{5.773502691896257, 0.7083333333333338}, 0.0416666666666666], CircleBox[{5.888972745734183, 0.8500000000000003}, 0.049999999999999975], CircleBox[{5.773502691896256, 1.8888888888888893}, 0.1111111111111111], CircleBox[{5.834276404442534, 1.736842105263158}, 0.05263157894736841], CircleBox[{5.888972745734182, 1.65}, 0.04999999999999999], CircleBox[{5.983448244328849, 1.5454545454545454}, 0.09090909090909091], CircleBox[{4.4088566010844135, 1.5454545454545459}, 0.09090909090909095], CircleBox[{4.503332099679081, 1.6500000000000004}, 0.04999999999999999], CircleBox[{4.5580284409707295, 1.7368421052631582}, 0.05263157894736843], CircleBox[{4.618802153517005, 1.8888888888888904}, 0.11111111111111117], CircleBox[{9.52627944162883, 0.06249999999999983}, 0.06249999999999983], CircleBox[{9.622504486493762, 0.037037037037036986}, 0.037037037037036986], CircleBox[{9.699484522385712, 0.03999999999999999}, 0.03999999999999999], CircleBox[{9.814954576223638, 0.0833333333333333}, 0.0833333333333333], CircleBox[{9.006664199358159, 0.04}, 0.04], CircleBox[{8.852704127574262, 0.037037037037037}, 0.037037037037037], CircleBox[{8.467803948114511, 0.037037037037037035}, 0.037037037037037035], CircleBox[{8.313843876330612, 0.04}, 0.04], CircleBox[{7.6210235533030595, 0.03999999999999999}, 0.03999999999999999], CircleBox[{7.698003589195009, 0.03703703703703703}, 0.03703703703703703], CircleBox[{7.794228634059946, 0.06249999999999996}, 0.06249999999999996], CircleBox[{7.505553499465134, 0.08333333333333331}, 0.08333333333333331], CircleBox[{8.660254037844386, 1.2500000000000004}, 0.7499999999999999], CircleBox[{11.547005383792516, 0.3333333333333333}, 0.3333333333333333], CircleBox[{12.12435565298214, 0.25}, 0.25], CircleBox[{12.701705922171769, 0.3333333333333333}, 0.3333333333333333], CircleBox[{13.856406460551018, 1.}], CircleBox[{2.6768057935155403, 0.3295454545454501}, 0.03409090909090852], CircleBox[{0.7872958216222171, 0.3295454545454546}, 0.03409090909090903], CircleBox[{0.9447549859466604, 1.4204545454545459}, 0.03409090909090906], CircleBox[{0.692820323027551, 1.9599999999999997}, 0.039999999999999994], CircleBox[{0.6415002990995844, 1.8148148148148158}, 0.037037037037037035], CircleBox[{1.299038105676658, 1.9531250000000002}, 0.04687499999999999], CircleBox[{2.5193466291910944, 1.4204545454545459}, 0.03409090909090901], CircleBox[{2.822601316038171, 1.814814814814816}, 0.03703703703703707], CircleBox[{2.7712812921102055, 1.9600000000000013}, 0.0400000000000001], CircleBox[{2.165063509461097, 1.9531250000000004}, 0.046875], CircleBox[{-3.008298771040681, 0.06578947368421045}, 0.03947368421052626], CircleBox[{-0.45580284409707295, 0.06578947368421052}, 0.0394736842105263], CircleBox[{-2.4248711305964257, 0.7250000000000002}, 0.07499999999999978], CircleBox[{-1.9485571585149875, 0.45312499999999845}, 0.046875], CircleBox[{-1.5155444566227678, 0.4531250000000012}, 0.0468749999999999], CircleBox[{-1.0392304845413265, 0.7250000000000003}, 0.0749999999999999], CircleBox[{-0.8660254037844385, 1.8125000000000002}, 0.18749999999999997], CircleBox[{-1.3856406460551018, 1.9600000000000006}, 0.04000000000000001], CircleBox[{-2.5980762113533156, 1.8125000000000002}, 0.18749999999999997], CircleBox[{-2.078460969082652, 1.9600000000000006}, 0.039999999999999994], CircleBox[{-4.2064091040958465, 0.17857142857142835}, 0.10714285714285696], CircleBox[{-3.849001794597502, 0.03703703703703708}, 0.03703703703703708], CircleBox[{-4.9296830676960335, 0.09615384615384615}, 0.05769230769230766], CircleBox[{-5.462621777717228, 0.09615384615384606}, 0.05769230769230761], CircleBox[{-6.1858957413174185, 0.1785714285714285}, 0.10714285714285705], CircleBox[{-6.543303050815759, 0.03703703703703704}, 0.03703703703703704], CircleBox[{-4.618802153517003, 0.7083333333333309}, 0.04166666666666662], CircleBox[{-4.503332099679082, 0.8500000000000005}, 0.04999999999999994], CircleBox[{-4.933720482165896, 0.5151515151515167}, 0.030303030303030262], CircleBox[{-5.047690924915014, 0.4857142857142879}, 0.0285714285714285], CircleBox[{-5.344613920498248, 0.4857142857142842}, 0.028571428571428595], CircleBox[{-5.458584363247371, 0.5151515151515159}, 0.03030303030303029], CircleBox[{-5.773502691896257, 0.7083333333333338}, 0.0416666666666666], CircleBox[{-5.888972745734183, 0.8500000000000003}, 0.049999999999999975], CircleBox[{-5.773502691896256, 1.8888888888888893}, 0.1111111111111111], CircleBox[{-5.834276404442534, 1.736842105263158}, 0.05263157894736841], CircleBox[{-5.888972745734182, 1.65}, 0.04999999999999999], CircleBox[{-5.983448244328849, 1.5454545454545454}, 0.09090909090909091], CircleBox[{-4.4088566010844135, 1.5454545454545459}, 0.09090909090909095], CircleBox[{-4.503332099679081, 1.6500000000000004}, 0.04999999999999999], CircleBox[{-4.5580284409707295, 1.7368421052631582}, 0.05263157894736843], CircleBox[{-4.618802153517005, 1.8888888888888904}, 0.11111111111111117], CircleBox[{-9.52627944162883, 0.06249999999999983}, 0.06249999999999983], CircleBox[{-9.622504486493762, 0.037037037037036986}, 0.037037037037036986], CircleBox[{-9.699484522385712, 0.03999999999999999}, 0.03999999999999999], CircleBox[{-9.814954576223638, 0.0833333333333333}, 0.0833333333333333], CircleBox[{-9.006664199358159, 0.04}, 0.04], CircleBox[{-8.852704127574262, 0.037037037037037}, 0.037037037037037], CircleBox[{-8.467803948114511, 0.037037037037037035}, 0.037037037037037035], CircleBox[{-8.313843876330612, 0.04}, 0.04], CircleBox[{-7.6210235533030595, 0.03999999999999999}, 0.03999999999999999], CircleBox[{-7.698003589195009, 0.03703703703703703}, 0.03703703703703703], CircleBox[{-7.794228634059946, 0.06249999999999996}, 0.06249999999999996], CircleBox[{-7.505553499465134, 0.08333333333333331}, 0.08333333333333331], CircleBox[{-8.660254037844386, 1.2500000000000004}, 0.7499999999999999], CircleBox[{-11.547005383792516, 0.3333333333333333}, 0.3333333333333333], CircleBox[{-12.12435565298214, 0.25}, 0.25], CircleBox[{-12.701705922171769, 0.3333333333333333}, 0.3333333333333333], CircleBox[{-13.856406460551018, 1.}], CircleBox[{3.9199044592348273, 0.06578947368421045}, 0.03947368421052626], CircleBox[{6.472400386178436, 0.06578947368421052}, 0.0394736842105263], CircleBox[{4.503332099679082, 0.7250000000000002}, 0.07499999999999978], CircleBox[{4.979646071760522, 0.45312499999999845}, 0.046875], CircleBox[{5.4126587736527405, 0.4531250000000012}, 0.0468749999999999], CircleBox[{5.888972745734181, 0.7250000000000003}, 0.0749999999999999], CircleBox[{6.062177826491071, 1.8125000000000002}, 0.18749999999999997], CircleBox[{5.542562584220407, 1.9600000000000006}, 0.04000000000000001], CircleBox[{4.330127018922193, 1.8125000000000002}, 0.18749999999999997], CircleBox[{4.849742261192858, 1.9600000000000006}, 0.039999999999999994], CircleBox[{9.649997356455172, 0.17857142857142835}, 0.10714285714285696], CircleBox[{10.007404665953516, 0.03703703703703708}, 0.03703703703703708], CircleBox[{8.926723392854983, 0.09615384615384615}, 0.05769230769230766], CircleBox[{8.393784682833788, 0.09615384615384606}, 0.05769230769230761], CircleBox[{7.670510719233598, 0.1785714285714285}, 0.10714285714285705], CircleBox[{7.31310340973526, 0.03703703703703704}, 0.03703703703703704], CircleBox[{9.237604307034017, 0.7083333333333309}, 0.04166666666666662], CircleBox[{9.353074360871934, 0.8500000000000005}, 0.04999999999999994], CircleBox[{8.922685978385122, 0.5151515151515167}, 0.030303030303030262], CircleBox[{8.808715535636004, 0.4857142857142879}, 0.0285714285714285], CircleBox[{8.51179254005277, 0.4857142857142842}, 0.028571428571428595], CircleBox[{8.397822097303646, 0.5151515151515159}, 0.03030303030303029], CircleBox[{8.08290376865476, 0.7083333333333338}, 0.0416666666666666], CircleBox[{7.967433714816835, 0.8500000000000003}, 0.049999999999999975], CircleBox[{8.082903768654761, 1.8888888888888893}, 0.1111111111111111], CircleBox[{8.022130056108484, 1.736842105263158}, 0.05263157894736841], CircleBox[{7.967433714816835, 1.65}, 0.04999999999999999], CircleBox[{7.87295821622217, 1.5454545454545454}, 0.09090909090909091], CircleBox[{9.447549859466603, 1.5454545454545459}, 0.09090909090909095], CircleBox[{9.353074360871936, 1.6500000000000004}, 0.04999999999999999], CircleBox[{9.298378019580287, 1.7368421052631582}, 0.05263157894736843], CircleBox[{9.237604307034012, 1.8888888888888904}, 0.11111111111111117], CircleBox[{11.258330249197698, 0.06249999999999983}, 0.06249999999999983], CircleBox[{11.162105204332763, 0.037037037037036986}, 0.037037037037036986], CircleBox[{11.085125168440815, 0.03999999999999999}, 0.03999999999999999], CircleBox[{10.96965511460289, 0.0833333333333333}, 0.0833333333333333], CircleBox[{11.777945491468367, 0.04}, 0.04], CircleBox[{11.931905563252265, 0.037037037037037}, 0.037037037037037], CircleBox[{12.316805742712011, 0.037037037037037035}, 0.037037037037037035], CircleBox[{12.470765814495918, 0.04}, 0.04], CircleBox[{13.163586137523467, 0.03999999999999999}, 0.03999999999999999], CircleBox[{13.086606101631517, 0.03703703703703703}, 0.03703703703703703], CircleBox[{12.990381056766578, 0.06249999999999996}, 0.06249999999999996], CircleBox[{13.279056191361388, 0.08333333333333331}, 0.08333333333333331], CircleBox[{12.124355652982143, 1.2500000000000004}, 0.7499999999999999], CircleBox[{17.32050807568877, 1.}], CircleBox[{16.16580753730952, 0.3333333333333333}, 0.3333333333333333], CircleBox[{15.588457268119894, 0.25}, 0.25], CircleBox[{15.011106998930268, 0.3333333333333333}, 0.3333333333333333], CircleBox[{0.5773502691896261, 1.9166666666666665}, 0.08333333333333331], CircleBox[{1.0392304845413265, 1.97}, 0.029999999999999992], CircleBox[{2.4248711305964292, 1.9700000000000004}, 0.03000000000000001], CircleBox[{2.8867513459481278, 1.9166666666666667}, 0.08333333333333327], CircleBox[{-2.6768057935155403, 0.3295454545454501}, 0.03409090909090852], CircleBox[{-0.7872958216222171, 0.3295454545454546}, 0.03409090909090903], CircleBox[{-0.9447549859466604, 1.4204545454545459}, 0.03409090909090906], CircleBox[{-0.692820323027551, 1.9599999999999997}, 0.039999999999999994], CircleBox[{-0.6415002990995844, 1.8148148148148158}, 0.037037037037037035], CircleBox[{-1.299038105676658, 1.9531250000000002}, 0.04687499999999999], CircleBox[{-2.5193466291910944, 1.4204545454545459}, 0.03409090909090901], CircleBox[{-2.822601316038171, 1.814814814814816}, 0.03703703703703707], CircleBox[{-2.7712812921102055, 1.9600000000000013}, 0.0400000000000001], CircleBox[{-2.165063509461097, 1.9531250000000004}, 0.046875], CircleBox[{-3.9199044592348273, 0.06578947368421045}, 0.03947368421052626], CircleBox[{-6.472400386178436, 0.06578947368421052}, 0.0394736842105263], CircleBox[{-4.503332099679082, 0.7250000000000002}, 0.07499999999999978], CircleBox[{-4.979646071760522, 0.45312499999999845}, 0.046875], CircleBox[{-5.4126587736527405, 0.4531250000000012}, 0.0468749999999999], CircleBox[{-5.888972745734181, 0.7250000000000003}, 0.0749999999999999], CircleBox[{-6.062177826491071, 1.8125000000000002}, 0.18749999999999997], CircleBox[{-5.542562584220407, 1.9600000000000006}, 0.04000000000000001], CircleBox[{-4.330127018922193, 1.8125000000000002}, 0.18749999999999997], CircleBox[{-4.849742261192858, 1.9600000000000006}, 0.039999999999999994], CircleBox[{-9.649997356455172, 0.17857142857142835}, 0.10714285714285696], CircleBox[{-10.007404665953516, 0.03703703703703708}, 0.03703703703703708], CircleBox[{-8.926723392854983, 0.09615384615384615}, 0.05769230769230766], CircleBox[{-8.393784682833788, 0.09615384615384606}, 0.05769230769230761], CircleBox[{-7.670510719233598, 0.1785714285714285}, 0.10714285714285705], CircleBox[{-7.31310340973526, 0.03703703703703704}, 0.03703703703703704], CircleBox[{-9.237604307034017, 0.7083333333333309}, 0.04166666666666662], CircleBox[{-9.353074360871934, 0.8500000000000005}, 0.04999999999999994], CircleBox[{-8.922685978385122, 0.5151515151515167}, 0.030303030303030262], CircleBox[{-8.808715535636004, 0.4857142857142879}, 0.0285714285714285], CircleBox[{-8.51179254005277, 0.4857142857142842}, 0.028571428571428595], CircleBox[{-8.397822097303646, 0.5151515151515159}, 0.03030303030303029], CircleBox[{-8.08290376865476, 0.7083333333333338}, 0.0416666666666666], CircleBox[{-7.967433714816835, 0.8500000000000003}, 0.049999999999999975], CircleBox[{-8.082903768654761, 1.8888888888888893}, 0.1111111111111111], CircleBox[{-8.022130056108484, 1.736842105263158}, 0.05263157894736841], CircleBox[{-7.967433714816835, 1.65}, 0.04999999999999999], CircleBox[{-7.87295821622217, 1.5454545454545454}, 0.09090909090909091], CircleBox[{-9.447549859466603, 1.5454545454545459}, 0.09090909090909095], CircleBox[{-9.353074360871936, 1.6500000000000004}, 0.04999999999999999], CircleBox[{-9.298378019580287, 1.7368421052631582}, 0.05263157894736843], CircleBox[{-9.237604307034012, 1.8888888888888904}, 0.11111111111111117], CircleBox[{-11.258330249197698, 0.06249999999999983}, 0.06249999999999983], CircleBox[{-11.162105204332763, 0.037037037037036986}, 0.037037037037036986], CircleBox[{-11.085125168440815, 0.03999999999999999}, 0.03999999999999999], CircleBox[{-10.96965511460289, 0.0833333333333333}, 0.0833333333333333], CircleBox[{-11.777945491468367, 0.04}, 0.04], CircleBox[{-11.931905563252265, 0.037037037037037}, 0.037037037037037], CircleBox[{-12.316805742712011, 0.037037037037037035}, 0.037037037037037035], CircleBox[{-12.470765814495918, 0.04}, 0.04], CircleBox[{-13.163586137523467, 0.03999999999999999}, 0.03999999999999999], CircleBox[{-13.086606101631517, 0.03703703703703703}, 0.03703703703703703], CircleBox[{-12.990381056766578, 0.06249999999999996}, 0.06249999999999996], CircleBox[{-13.279056191361388, 0.08333333333333331}, 0.08333333333333331], CircleBox[{-12.124355652982143, 1.2500000000000004}, 0.7499999999999999], CircleBox[{-17.32050807568877, 1.}], CircleBox[{-16.16580753730952, 0.3333333333333333}, 0.3333333333333333], CircleBox[{-15.588457268119894, 0.25}, 0.25], CircleBox[{-15.011106998930268, 0.3333333333333333}, 0.3333333333333333], CircleBox[{4.251397436759969, 0.3295454545454501}, 0.03409090909090852], CircleBox[{6.140907408653292, 0.3295454545454546}, 0.03409090909090903], CircleBox[{5.983448244328849, 1.4204545454545459}, 0.03409090909090906], CircleBox[{6.235382907247957, 1.9599999999999997}, 0.039999999999999994], CircleBox[{6.286702931175924, 1.8148148148148158}, 0.037037037037037035], CircleBox[{5.629165124598851, 1.9531250000000002}, 0.04687499999999999], CircleBox[{4.408856601084415, 1.4204545454545459}, 0.03409090909090901], CircleBox[{4.105601914237338, 1.814814814814816}, 0.03703703703703707], CircleBox[{4.156921938165303, 1.9600000000000013}, 0.0400000000000001], CircleBox[{4.7631397208144115, 1.9531250000000004}, 0.046875], CircleBox[{9.936502001316189, 0.06578947368421045}, 0.03947368421052626], CircleBox[{7.384006074372581, 0.06578947368421052}, 0.0394736842105263], CircleBox[{9.353074360871934, 0.7250000000000002}, 0.07499999999999978], CircleBox[{8.876760388790496, 0.45312499999999845}, 0.046875], CircleBox[{8.443747686898277, 0.4531250000000012}, 0.0468749999999999], CircleBox[{7.967433714816835, 0.7250000000000003}, 0.0749999999999999], CircleBox[{7.794228634059947, 1.8125000000000002}, 0.18749999999999997], CircleBox[{8.31384387633061, 1.9600000000000006}, 0.04000000000000001], CircleBox[{9.526279441628825, 1.8125000000000002}, 0.18749999999999997], CircleBox[{9.00666419935816, 1.9600000000000006}, 0.039999999999999994], CircleBox[{11.134612334371356, 0.17857142857142835}, 0.10714285714285696], CircleBox[{10.77720502487301, 0.03703703703703708}, 0.03703703703703708], CircleBox[{11.857886297971543, 0.09615384615384615}, 0.05769230769230766], CircleBox[{12.390825007992737, 0.09615384615384606}, 0.05769230769230761], CircleBox[{13.114098971592927, 0.1785714285714285}, 0.10714285714285705], CircleBox[{13.471506281091267, 0.03703703703703704}, 0.03703703703703704], CircleBox[{11.547005383792511, 0.7083333333333309}, 0.04166666666666662], CircleBox[{11.431535329954592, 0.8500000000000005}, 0.04999999999999994], CircleBox[{11.861923712441405, 0.5151515151515167}, 0.030303030303030262], CircleBox[{11.975894155190522, 0.4857142857142879}, 0.0285714285714285], CircleBox[{12.272817150773758, 0.4857142857142842}, 0.028571428571428595], CircleBox[{12.386787593522879, 0.5151515151515159}, 0.03030303030303029], CircleBox[{12.701705922171765, 0.7083333333333338}, 0.0416666666666666], CircleBox[{12.817175976009693, 0.8500000000000003}, 0.049999999999999975], CircleBox[{12.701705922171765, 1.8888888888888893}, 0.1111111111111111], CircleBox[{12.762479634718042, 1.736842105263158}, 0.05263157894736841], CircleBox[{12.81717597600969, 1.65}, 0.04999999999999999], CircleBox[{12.911651474604358, 1.5454545454545454}, 0.09090909090909091], CircleBox[{11.337059831359921, 1.5454545454545459}, 0.09090909090909095], CircleBox[{11.43153532995459, 1.6500000000000004}, 0.04999999999999999], CircleBox[{11.48623167124624, 1.7368421052631582}, 0.05263157894736843], CircleBox[{11.547005383792513, 1.8888888888888904}, 0.11111111111111117], CircleBox[{16.45448267190434, 0.06249999999999983}, 0.06249999999999983], CircleBox[{16.55070771676927, 0.037037037037036986}, 0.037037037037036986], CircleBox[{16.62768775266122, 0.03999999999999999}, 0.03999999999999999], CircleBox[{16.74315780649915, 0.0833333333333333}, 0.0833333333333333], CircleBox[{15.93486742963367, 0.04}, 0.04], CircleBox[{15.78090735784977, 0.037037037037037}, 0.037037037037037], CircleBox[{15.39600717839002, 0.037037037037037035}, 0.037037037037037035], CircleBox[{15.242047106606119, 0.04}, 0.04], CircleBox[{14.54922678357857, 0.03999999999999999}, 0.03999999999999999], CircleBox[{14.626206819470518, 0.03703703703703703}, 0.03703703703703703], CircleBox[{14.722431864335455, 0.06249999999999996}, 0.06249999999999996], CircleBox[{14.433756729740644, 0.08333333333333331}, 0.08333333333333331], CircleBox[{15.588457268119894, 1.2500000000000004}, 0.7499999999999999], CircleBox[{18.475208614068023, 0.3333333333333333}, 0.3333333333333333], CircleBox[{19.05255888325765, 0.25}, 0.25], CircleBox[{19.629909152447276, 0.3333333333333333}, 0.3333333333333333], CircleBox[{20.784609690826528, 1.}]}, {RGBColor[1, 0, 0], CircleBox[{2.886751345948128, 0.}, 0.5773502691896256], CircleBox[{2.02072594216369, 0.}, 0.28867513459481287], CircleBox[{1.4433756729740643, 0.}, 0.28867513459481287], CircleBox[{0.5773502691896257, 0.}, 0.5773502691896257], CircleBox[{1.7320508075688774, 2.0000000000000004}, 1.7320508075688774], TagBox[ConicHullRegionBox[{0, 0}, {{0, -1}}], "InfiniteLine"], TagBox[ConicHullRegionBox[{3.4641016151377544`, 0.}, {{0, 1}}], "InfiniteLine"]}, {RGBColor[0, 0, 1], CircleBox[{3.464101615137755, 1}], CircleBox[ NCache[{2.309401076758503, Rational[1, 3]}, {2.309401076758503, 0.3333333333333333}], NCache[Rational[1, 3], 0.3333333333333333]], CircleBox[ NCache[{1.732050807568877, Rational[1, 4]}, {1.732050807568877, 0.25}], NCache[Rational[1, 4], 0.25]], CircleBox[ NCache[{1.154700538379251, Rational[1, 3]}, {1.154700538379251, 0.3333333333333333}], NCache[Rational[1, 3], 0.3333333333333333]], CircleBox[{0, 1}], TagBox[ConicHullRegionBox[{0, 2}, {{-1, 0}}], "InfiniteLine"], TagBox[ConicHullRegionBox[{0, 0}, {{1, 0}}], "InfiniteLine"]}}, PlotRange->{{-3.4641016151377544`, 6.928203230275509}, {-2.078460969082653, 4.078460969082653}}, PlotRangeClipping->True]}, { TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"12", "1", RowBox[{"2", " ", SqrtBox["3"]}], "1"}, {"16", "3", RowBox[{"4", " ", SqrtBox["3"]}], "1"}, {"12", "4", RowBox[{"4", " ", SqrtBox["3"]}], "1"}, {"4", "3", RowBox[{"2", " ", SqrtBox["3"]}], "1"}, {"0", "1", "0", "1"}, {"4", "0", "0", "1"}, {"0", "0", "0", RowBox[{"-", "1"}]}, { RowBox[{"8", " ", SqrtBox["3"]}], SqrtBox["3"], "5", "0"}, { RowBox[{"8", " ", SqrtBox["3"]}], RowBox[{"2", " ", SqrtBox["3"]}], "7", "0"}, { RowBox[{"4", " ", SqrtBox["3"]}], RowBox[{"2", " ", SqrtBox["3"]}], "5", "0"}, {"0", SqrtBox["3"], "1", "0"}, { FractionBox["4", SqrtBox["3"]], FractionBox["1", SqrtBox["3"]], "1", FractionBox["2", SqrtBox["3"]]}, {"0", "0", RowBox[{"-", "1"}], "0"}, { RowBox[{"4", " ", SqrtBox["3"]}], "0", "1", "0"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]}, { RowBox[{"{", RowBox[{ TagBox[GridBox[{ { TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { SubscriptBox["a", "1"], SubscriptBox["a", "2"], SubscriptBox["a", "3"], RowBox[{"1", "-", SubscriptBox["a", "1"], "-", RowBox[{"2", " ", SubscriptBox["a", "2"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "3"]}]}], RowBox[{ RowBox[{"-", "2"}], "+", RowBox[{"2", " ", SubscriptBox["a", "1"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "2"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "3"]}]}], RowBox[{"2", "-", RowBox[{"2", " ", SubscriptBox["a", "1"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "2"]}], "-", SubscriptBox["a", "3"]}], "0"}, { SubscriptBox["a", "4"], SubscriptBox["a", "5"], SubscriptBox["a", "6"], RowBox[{"2", "-", SubscriptBox["a", "4"], "-", RowBox[{"2", " ", SubscriptBox["a", "5"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "6"]}]}], RowBox[{ RowBox[{"-", "3"}], "+", RowBox[{"2", " ", SubscriptBox["a", "4"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "5"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "6"]}]}], RowBox[{"2", "-", RowBox[{"2", " ", SubscriptBox["a", "4"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "5"]}], "-", SubscriptBox["a", "6"]}], "0"}, { SubscriptBox["a", "7"], SubscriptBox["a", "8"], SubscriptBox["a", "9"], RowBox[{"12", "-", SubscriptBox["a", "7"], "-", RowBox[{"2", " ", SubscriptBox["a", "8"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "9"]}]}], RowBox[{ RowBox[{"-", "20"}], "+", RowBox[{"2", " ", SubscriptBox["a", "7"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "8"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "9"]}]}], RowBox[{"15", "-", RowBox[{"2", " ", SubscriptBox["a", "7"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "8"]}], "-", SubscriptBox["a", "9"]}], "6"}, { SubscriptBox["a", "10"], SubscriptBox["a", "11"], SubscriptBox["a", "12"], RowBox[{"21", "-", SubscriptBox["a", "10"], "-", RowBox[{"2", " ", SubscriptBox["a", "11"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "12"]}]}], RowBox[{ RowBox[{"-", "36"}], "+", RowBox[{"2", " ", SubscriptBox["a", "10"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "11"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "12"]}]}], RowBox[{"28", "-", RowBox[{"2", " ", SubscriptBox["a", "10"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "11"]}], "-", SubscriptBox["a", "12"]}], "12"}, { SubscriptBox["a", "13"], SubscriptBox["a", "14"], SubscriptBox["a", "15"], RowBox[{"20", "-", SubscriptBox["a", "13"], "-", RowBox[{"2", " ", SubscriptBox["a", "14"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "15"]}]}], RowBox[{ RowBox[{"-", "35"}], "+", RowBox[{"2", " ", SubscriptBox["a", "13"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "14"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "15"]}]}], RowBox[{"28", "-", RowBox[{"2", " ", SubscriptBox["a", "13"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "14"]}], "-", SubscriptBox["a", "15"]}], "12"}, { SubscriptBox["a", "16"], SubscriptBox["a", "17"], SubscriptBox["a", "18"], RowBox[{"10", "-", SubscriptBox["a", "16"], "-", RowBox[{"2", " ", SubscriptBox["a", "17"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "18"]}]}], RowBox[{ RowBox[{"-", "18"}], "+", RowBox[{"2", " ", SubscriptBox["a", "16"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "17"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "18"]}]}], RowBox[{"15", "-", RowBox[{"2", " ", SubscriptBox["a", "16"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "17"]}], "-", SubscriptBox["a", "18"]}], "6"}, { SubscriptBox["a", "19"], SubscriptBox["a", "20"], SubscriptBox["a", "21"], RowBox[{ RowBox[{"-", SubscriptBox["a", "19"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "20"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "21"]}]}], RowBox[{ RowBox[{"2", " ", SubscriptBox["a", "19"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "20"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "21"]}]}], RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", SubscriptBox["a", "19"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "20"]}], "-", SubscriptBox["a", "21"]}], "1"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]}, { TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0", "0", "0", "1", RowBox[{"-", "2"}], "2", "0"}, {"0", "0", "0", "2", RowBox[{"-", "3"}], "2", "0"}, {"0", "0", "0", "12", RowBox[{"-", "20"}], "15", "6"}, {"0", "0", "0", "21", RowBox[{"-", "36"}], "28", "12"}, {"0", "0", "0", "20", RowBox[{"-", "35"}], "28", "12"}, {"0", "0", "0", "10", RowBox[{"-", "18"}], "15", "6"}, {"0", "0", "0", "0", "0", "0", "1"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]} }, DefaultBaseStyle->"Column", GridBoxAlignment->{"Columns" -> {{Left}}}, GridBoxItemSize->{ "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}], "Column"], ",", TagBox[GridBox[{ { TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { SubscriptBox["a", "1"], SubscriptBox["a", "2"], SubscriptBox["a", "3"], RowBox[{"15", "-", SubscriptBox["a", "1"], "-", RowBox[{"2", " ", SubscriptBox["a", "2"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "3"]}]}], RowBox[{ RowBox[{"-", "20"}], "+", RowBox[{"2", " ", SubscriptBox["a", "1"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "2"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "3"]}]}], RowBox[{"12", "-", RowBox[{"2", " ", SubscriptBox["a", "1"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "2"]}], "-", SubscriptBox["a", "3"]}], "6"}, { SubscriptBox["a", "4"], SubscriptBox["a", "5"], SubscriptBox["a", "6"], RowBox[{"2", "-", SubscriptBox["a", "4"], "-", RowBox[{"2", " ", SubscriptBox["a", "5"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "6"]}]}], RowBox[{ RowBox[{"-", "3"}], "+", RowBox[{"2", " ", SubscriptBox["a", "4"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "5"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "6"]}]}], RowBox[{"2", "-", RowBox[{"2", " ", SubscriptBox["a", "4"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "5"]}], "-", SubscriptBox["a", "6"]}], "0"}, { SubscriptBox["a", "7"], SubscriptBox["a", "8"], SubscriptBox["a", "9"], RowBox[{"2", "-", SubscriptBox["a", "7"], "-", RowBox[{"2", " ", SubscriptBox["a", "8"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "9"]}]}], RowBox[{ RowBox[{"-", "2"}], "+", RowBox[{"2", " ", SubscriptBox["a", "7"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "8"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "9"]}]}], RowBox[{"1", "-", RowBox[{"2", " ", SubscriptBox["a", "7"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "8"]}], "-", SubscriptBox["a", "9"]}], "0"}, { SubscriptBox["a", "10"], SubscriptBox["a", "11"], SubscriptBox["a", "12"], RowBox[{"15", "-", SubscriptBox["a", "10"], "-", RowBox[{"2", " ", SubscriptBox["a", "11"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "12"]}]}], RowBox[{ RowBox[{"-", "18"}], "+", RowBox[{"2", " ", SubscriptBox["a", "10"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "11"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "12"]}]}], RowBox[{"10", "-", RowBox[{"2", " ", SubscriptBox["a", "10"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "11"]}], "-", SubscriptBox["a", "12"]}], "6"}, { SubscriptBox["a", "13"], SubscriptBox["a", "14"], SubscriptBox["a", "15"], RowBox[{"28", "-", SubscriptBox["a", "13"], "-", RowBox[{"2", " ", SubscriptBox["a", "14"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "15"]}]}], RowBox[{ RowBox[{"-", "35"}], "+", RowBox[{"2", " ", SubscriptBox["a", "13"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "14"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "15"]}]}], RowBox[{"20", "-", RowBox[{"2", " ", SubscriptBox["a", "13"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "14"]}], "-", SubscriptBox["a", "15"]}], "12"}, { SubscriptBox["a", "16"], SubscriptBox["a", "17"], SubscriptBox["a", "18"], RowBox[{"28", "-", SubscriptBox["a", "16"], "-", RowBox[{"2", " ", SubscriptBox["a", "17"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "18"]}]}], RowBox[{ RowBox[{"-", "36"}], "+", RowBox[{"2", " ", SubscriptBox["a", "16"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "17"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "18"]}]}], RowBox[{"21", "-", RowBox[{"2", " ", SubscriptBox["a", "16"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "17"]}], "-", SubscriptBox["a", "18"]}], "12"}, { SubscriptBox["a", "19"], SubscriptBox["a", "20"], SubscriptBox["a", "21"], RowBox[{ RowBox[{"-", SubscriptBox["a", "19"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "20"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "21"]}]}], RowBox[{ RowBox[{"2", " ", SubscriptBox["a", "19"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "20"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "21"]}]}], RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", SubscriptBox["a", "19"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "20"]}], "-", SubscriptBox["a", "21"]}], "1"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]}, { TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0", "0", "0", "15", RowBox[{"-", "20"}], "12", "6"}, {"0", "0", "0", "2", RowBox[{"-", "3"}], "2", "0"}, {"0", "0", "0", "2", RowBox[{"-", "2"}], "1", "0"}, {"0", "0", "0", "15", RowBox[{"-", "18"}], "10", "6"}, {"0", "0", "0", "28", RowBox[{"-", "35"}], "20", "12"}, {"0", "0", "0", "28", RowBox[{"-", "36"}], "21", "12"}, {"0", "0", "0", "0", "0", "0", "1"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]} }, DefaultBaseStyle->"Column", GridBoxAlignment->{"Columns" -> {{Left}}}, GridBoxItemSize->{ "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}], "Column"], ",", TagBox[GridBox[{ { TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { SubscriptBox["a", "1"], SubscriptBox["a", "2"], SubscriptBox["a", "3"], RowBox[{"21", "-", SubscriptBox["a", "1"], "-", RowBox[{"2", " ", SubscriptBox["a", "2"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "3"]}]}], RowBox[{ RowBox[{"-", "14"}], "+", RowBox[{"2", " ", SubscriptBox["a", "1"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "2"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "3"]}]}], RowBox[{"6", "-", RowBox[{"2", " ", SubscriptBox["a", "1"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "2"]}], "-", SubscriptBox["a", "3"]}], "12"}, { SubscriptBox["a", "4"], SubscriptBox["a", "5"], SubscriptBox["a", "6"], RowBox[{"12", "-", SubscriptBox["a", "4"], "-", RowBox[{"2", " ", SubscriptBox["a", "5"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "6"]}]}], RowBox[{ RowBox[{"-", "9"}], "+", RowBox[{"2", " ", SubscriptBox["a", "4"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "5"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "6"]}]}], RowBox[{"4", "-", RowBox[{"2", " ", SubscriptBox["a", "4"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "5"]}], "-", SubscriptBox["a", "6"]}], "6"}, { SubscriptBox["a", "7"], SubscriptBox["a", "8"], SubscriptBox["a", "9"], RowBox[{"2", "-", SubscriptBox["a", "7"], "-", RowBox[{"2", " ", SubscriptBox["a", "8"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "9"]}]}], RowBox[{ RowBox[{"-", "2"}], "+", RowBox[{"2", " ", SubscriptBox["a", "7"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "8"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "9"]}]}], RowBox[{"1", "-", RowBox[{"2", " ", SubscriptBox["a", "7"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "8"]}], "-", SubscriptBox["a", "9"]}], "0"}, { SubscriptBox["a", "10"], SubscriptBox["a", "11"], SubscriptBox["a", "12"], RowBox[{"1", "-", SubscriptBox["a", "10"], "-", RowBox[{"2", " ", SubscriptBox["a", "11"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "12"]}]}], RowBox[{ RowBox[{"2", " ", SubscriptBox["a", "10"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "11"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "12"]}]}], RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", SubscriptBox["a", "10"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "11"]}], "-", SubscriptBox["a", "12"]}], "0"}, { SubscriptBox["a", "13"], SubscriptBox["a", "14"], SubscriptBox["a", "15"], RowBox[{"10", "-", SubscriptBox["a", "13"], "-", RowBox[{"2", " ", SubscriptBox["a", "14"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "15"]}]}], RowBox[{ RowBox[{"-", "5"}], "+", RowBox[{"2", " ", SubscriptBox["a", "13"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "14"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "15"]}]}], RowBox[{"2", "-", RowBox[{"2", " ", SubscriptBox["a", "13"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "14"]}], "-", SubscriptBox["a", "15"]}], "6"}, { SubscriptBox["a", "16"], SubscriptBox["a", "17"], SubscriptBox["a", "18"], RowBox[{"20", "-", SubscriptBox["a", "16"], "-", RowBox[{"2", " ", SubscriptBox["a", "17"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "18"]}]}], RowBox[{ RowBox[{"-", "12"}], "+", RowBox[{"2", " ", SubscriptBox["a", "16"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "17"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "18"]}]}], RowBox[{"5", "-", RowBox[{"2", " ", SubscriptBox["a", "16"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "17"]}], "-", SubscriptBox["a", "18"]}], "12"}, { SubscriptBox["a", "19"], SubscriptBox["a", "20"], SubscriptBox["a", "21"], RowBox[{ RowBox[{"-", SubscriptBox["a", "19"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "20"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "21"]}]}], RowBox[{ RowBox[{"2", " ", SubscriptBox["a", "19"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "20"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "21"]}]}], RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", SubscriptBox["a", "19"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "20"]}], "-", SubscriptBox["a", "21"]}], "1"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]}, { TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0", "0", "0", "21", RowBox[{"-", "14"}], "6", "12"}, {"0", "0", "0", "12", RowBox[{"-", "9"}], "4", "6"}, {"0", "0", "0", "2", RowBox[{"-", "2"}], "1", "0"}, {"0", "0", "0", "1", "0", "0", "0"}, {"0", "0", "0", "10", RowBox[{"-", "5"}], "2", "6"}, {"0", "0", "0", "20", RowBox[{"-", "12"}], "5", "12"}, {"0", "0", "0", "0", "0", "0", "1"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]} }, DefaultBaseStyle->"Column", GridBoxAlignment->{"Columns" -> {{Left}}}, GridBoxItemSize->{ "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}], "Column"], ",", TagBox[GridBox[{ { TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { SubscriptBox["a", "1"], SubscriptBox["a", "2"], SubscriptBox["a", "3"], RowBox[{"5", "-", SubscriptBox["a", "1"], "-", RowBox[{"2", " ", SubscriptBox["a", "2"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "3"]}]}], RowBox[{"10", "+", RowBox[{"2", " ", SubscriptBox["a", "1"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "2"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "3"]}]}], RowBox[{ RowBox[{"-", "2"}], "-", RowBox[{"2", " ", SubscriptBox["a", "1"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "2"]}], "-", SubscriptBox["a", "3"]}], "12"}, { SubscriptBox["a", "4"], SubscriptBox["a", "5"], SubscriptBox["a", "6"], RowBox[{"6", "-", SubscriptBox["a", "4"], "-", RowBox[{"2", " ", SubscriptBox["a", "5"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "6"]}]}], RowBox[{"9", "+", RowBox[{"2", " ", SubscriptBox["a", "4"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "5"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "6"]}]}], RowBox[{ RowBox[{"-", "2"}], "-", RowBox[{"2", " ", SubscriptBox["a", "4"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "5"]}], "-", SubscriptBox["a", "6"]}], "12"}, { SubscriptBox["a", "7"], SubscriptBox["a", "8"], SubscriptBox["a", "9"], RowBox[{"4", "-", SubscriptBox["a", "7"], "-", RowBox[{"2", " ", SubscriptBox["a", "8"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "9"]}]}], RowBox[{"4", "+", RowBox[{"2", " ", SubscriptBox["a", "7"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "8"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "9"]}]}], RowBox[{ RowBox[{"-", "1"}], "-", RowBox[{"2", " ", SubscriptBox["a", "7"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "8"]}], "-", SubscriptBox["a", "9"]}], "6"}, { SubscriptBox["a", "10"], SubscriptBox["a", "11"], SubscriptBox["a", "12"], RowBox[{"1", "-", SubscriptBox["a", "10"], "-", RowBox[{"2", " ", SubscriptBox["a", "11"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "12"]}]}], RowBox[{ RowBox[{"2", " ", SubscriptBox["a", "10"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "11"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "12"]}]}], RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", SubscriptBox["a", "10"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "11"]}], "-", SubscriptBox["a", "12"]}], "0"}, { SubscriptBox["a", "13"], SubscriptBox["a", "14"], SubscriptBox["a", "15"], RowBox[{ RowBox[{"-", SubscriptBox["a", "13"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "14"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "15"]}]}], RowBox[{"1", "+", RowBox[{"2", " ", SubscriptBox["a", "13"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "14"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "15"]}]}], RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", SubscriptBox["a", "13"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "14"]}], "-", SubscriptBox["a", "15"]}], "0"}, { SubscriptBox["a", "16"], SubscriptBox["a", "17"], SubscriptBox["a", "18"], RowBox[{"2", "-", SubscriptBox["a", "16"], "-", RowBox[{"2", " ", SubscriptBox["a", "17"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "18"]}]}], RowBox[{"6", "+", RowBox[{"2", " ", SubscriptBox["a", "16"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "17"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "18"]}]}], RowBox[{ RowBox[{"-", "1"}], "-", RowBox[{"2", " ", SubscriptBox["a", "16"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "17"]}], "-", SubscriptBox["a", "18"]}], "6"}, { SubscriptBox["a", "19"], SubscriptBox["a", "20"], SubscriptBox["a", "21"], RowBox[{ RowBox[{"-", SubscriptBox["a", "19"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "20"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "21"]}]}], RowBox[{ RowBox[{"2", " ", SubscriptBox["a", "19"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "20"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "21"]}]}], RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", SubscriptBox["a", "19"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "20"]}], "-", SubscriptBox["a", "21"]}], "1"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]}, { TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0", "0", "0", "5", "10", RowBox[{"-", "2"}], "12"}, {"0", "0", "0", "6", "9", RowBox[{"-", "2"}], "12"}, {"0", "0", "0", "4", "4", RowBox[{"-", "1"}], "6"}, {"0", "0", "0", "1", "0", "0", "0"}, {"0", "0", "0", "0", "1", "0", "0"}, {"0", "0", "0", "2", "6", RowBox[{"-", "1"}], "6"}, {"0", "0", "0", "0", "0", "0", "1"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]} }, DefaultBaseStyle->"Column", GridBoxAlignment->{"Columns" -> {{Left}}}, GridBoxItemSize->{ "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}], "Column"], ",", TagBox[GridBox[{ { TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { SubscriptBox["a", "1"], SubscriptBox["a", "2"], SubscriptBox["a", "3"], RowBox[{"1", "-", SubscriptBox["a", "1"], "-", RowBox[{"2", " ", SubscriptBox["a", "2"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "3"]}]}], RowBox[{ RowBox[{"-", "2"}], "+", RowBox[{"2", " ", SubscriptBox["a", "1"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "2"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "3"]}]}], RowBox[{"2", "-", RowBox[{"2", " ", SubscriptBox["a", "1"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "2"]}], "-", SubscriptBox["a", "3"]}], "0"}, { SubscriptBox["a", "4"], SubscriptBox["a", "5"], SubscriptBox["a", "6"], RowBox[{"2", "-", SubscriptBox["a", "4"], "-", RowBox[{"2", " ", SubscriptBox["a", "5"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "6"]}]}], RowBox[{ RowBox[{"-", "3"}], "+", RowBox[{"2", " ", SubscriptBox["a", "4"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "5"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "6"]}]}], RowBox[{"2", "-", RowBox[{"2", " ", SubscriptBox["a", "4"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "5"]}], "-", SubscriptBox["a", "6"]}], "0"}, { SubscriptBox["a", "7"], SubscriptBox["a", "8"], SubscriptBox["a", "9"], RowBox[{"2", "-", SubscriptBox["a", "7"], "-", RowBox[{"2", " ", SubscriptBox["a", "8"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "9"]}]}], RowBox[{ RowBox[{"-", "2"}], "+", RowBox[{"2", " ", SubscriptBox["a", "7"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "8"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "9"]}]}], RowBox[{"1", "-", RowBox[{"2", " ", SubscriptBox["a", "7"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "8"]}], "-", SubscriptBox["a", "9"]}], "0"}, { SubscriptBox["a", "10"], SubscriptBox["a", "11"], SubscriptBox["a", "12"], RowBox[{"1", "-", SubscriptBox["a", "10"], "-", RowBox[{"2", " ", SubscriptBox["a", "11"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "12"]}]}], RowBox[{ RowBox[{"2", " ", SubscriptBox["a", "10"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "11"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "12"]}]}], RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", SubscriptBox["a", "10"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "11"]}], "-", SubscriptBox["a", "12"]}], "0"}, { SubscriptBox["a", "13"], SubscriptBox["a", "14"], SubscriptBox["a", "15"], RowBox[{ RowBox[{"-", SubscriptBox["a", "13"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "14"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "15"]}]}], RowBox[{"1", "+", RowBox[{"2", " ", SubscriptBox["a", "13"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "14"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "15"]}]}], RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", SubscriptBox["a", "13"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "14"]}], "-", SubscriptBox["a", "15"]}], "0"}, { SubscriptBox["a", "16"], SubscriptBox["a", "17"], SubscriptBox["a", "18"], RowBox[{ RowBox[{"-", SubscriptBox["a", "16"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "17"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "18"]}]}], RowBox[{ RowBox[{"2", " ", SubscriptBox["a", "16"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "17"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "18"]}]}], RowBox[{"1", "-", RowBox[{"2", " ", SubscriptBox["a", "16"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "17"]}], "-", SubscriptBox["a", "18"]}], "0"}, { SubscriptBox["a", "19"], SubscriptBox["a", "20"], SubscriptBox["a", "21"], RowBox[{ FractionBox["2", "3"], "-", SubscriptBox["a", "19"], "-", RowBox[{"2", " ", SubscriptBox["a", "20"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "21"]}]}], RowBox[{ RowBox[{"-", FractionBox["2", "3"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "19"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "20"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "21"]}]}], RowBox[{ FractionBox["2", "3"], "-", RowBox[{"2", " ", SubscriptBox["a", "19"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "20"]}], "-", SubscriptBox["a", "21"]}], RowBox[{"-", "1"}]} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]}, { TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0", "0", "0", "1", RowBox[{"-", "2"}], "2", "0"}, {"0", "0", "0", "2", RowBox[{"-", "3"}], "2", "0"}, {"0", "0", "0", "2", RowBox[{"-", "2"}], "1", "0"}, {"0", "0", "0", "1", "0", "0", "0"}, {"0", "0", "0", "0", "1", "0", "0"}, {"0", "0", "0", "0", "0", "1", "0"}, {"0", "0", "0", FractionBox["2", "3"], RowBox[{"-", FractionBox["2", "3"]}], FractionBox["2", "3"], RowBox[{"-", "1"}]} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]} }, DefaultBaseStyle->"Column", GridBoxAlignment->{"Columns" -> {{Left}}}, GridBoxItemSize->{ "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}], "Column"], ",", TagBox[GridBox[{ { TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { SubscriptBox["a", "1"], SubscriptBox["a", "2"], SubscriptBox["a", "3"], RowBox[{ RowBox[{"-", "1"}], "-", SubscriptBox["a", "1"], "-", RowBox[{"2", " ", SubscriptBox["a", "2"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "3"]}]}], RowBox[{"4", "+", RowBox[{"2", " ", SubscriptBox["a", "1"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "2"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "3"]}]}], RowBox[{"4", "-", RowBox[{"2", " ", SubscriptBox["a", "1"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "2"]}], "-", SubscriptBox["a", "3"]}], "6"}, { SubscriptBox["a", "4"], SubscriptBox["a", "5"], SubscriptBox["a", "6"], RowBox[{ RowBox[{"-", "2"}], "-", SubscriptBox["a", "4"], "-", RowBox[{"2", " ", SubscriptBox["a", "5"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "6"]}]}], RowBox[{"9", "+", RowBox[{"2", " ", SubscriptBox["a", "4"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "5"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "6"]}]}], RowBox[{"6", "-", RowBox[{"2", " ", SubscriptBox["a", "4"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "5"]}], "-", SubscriptBox["a", "6"]}], "12"}, { SubscriptBox["a", "7"], SubscriptBox["a", "8"], SubscriptBox["a", "9"], RowBox[{ RowBox[{"-", "2"}], "-", SubscriptBox["a", "7"], "-", RowBox[{"2", " ", SubscriptBox["a", "8"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "9"]}]}], RowBox[{"10", "+", RowBox[{"2", " ", SubscriptBox["a", "7"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "8"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "9"]}]}], RowBox[{"5", "-", RowBox[{"2", " ", SubscriptBox["a", "7"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "8"]}], "-", SubscriptBox["a", "9"]}], "12"}, { SubscriptBox["a", "10"], SubscriptBox["a", "11"], SubscriptBox["a", "12"], RowBox[{ RowBox[{"-", "1"}], "-", SubscriptBox["a", "10"], "-", RowBox[{"2", " ", SubscriptBox["a", "11"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "12"]}]}], RowBox[{"6", "+", RowBox[{"2", " ", SubscriptBox["a", "10"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "11"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "12"]}]}], RowBox[{"2", "-", RowBox[{"2", " ", SubscriptBox["a", "10"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "11"]}], "-", SubscriptBox["a", "12"]}], "6"}, { SubscriptBox["a", "13"], SubscriptBox["a", "14"], SubscriptBox["a", "15"], RowBox[{ RowBox[{"-", SubscriptBox["a", "13"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "14"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "15"]}]}], RowBox[{"1", "+", RowBox[{"2", " ", SubscriptBox["a", "13"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "14"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "15"]}]}], RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", SubscriptBox["a", "13"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "14"]}], "-", SubscriptBox["a", "15"]}], "0"}, { SubscriptBox["a", "16"], SubscriptBox["a", "17"], SubscriptBox["a", "18"], RowBox[{ RowBox[{"-", SubscriptBox["a", "16"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "17"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "18"]}]}], RowBox[{ RowBox[{"2", " ", SubscriptBox["a", "16"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "17"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "18"]}]}], RowBox[{"1", "-", RowBox[{"2", " ", SubscriptBox["a", "16"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "17"]}], "-", SubscriptBox["a", "18"]}], "0"}, { SubscriptBox["a", "19"], SubscriptBox["a", "20"], SubscriptBox["a", "21"], RowBox[{ RowBox[{"-", SubscriptBox["a", "19"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "20"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "21"]}]}], RowBox[{ RowBox[{"2", " ", SubscriptBox["a", "19"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "20"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "21"]}]}], RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", SubscriptBox["a", "19"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "20"]}], "-", SubscriptBox["a", "21"]}], "1"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]}, { TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0", "0", "0", RowBox[{"-", "1"}], "4", "4", "6"}, {"0", "0", "0", RowBox[{"-", "2"}], "9", "6", "12"}, {"0", "0", "0", RowBox[{"-", "2"}], "10", "5", "12"}, {"0", "0", "0", RowBox[{"-", "1"}], "6", "2", "6"}, {"0", "0", "0", "0", "1", "0", "0"}, {"0", "0", "0", "0", "0", "1", "0"}, {"0", "0", "0", "0", "0", "0", "1"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]} }, DefaultBaseStyle->"Column", GridBoxAlignment->{"Columns" -> {{Left}}}, GridBoxItemSize->{ "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}], "Column"], ",", TagBox[GridBox[{ { TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { SubscriptBox["a", "1"], SubscriptBox["a", "2"], SubscriptBox["a", "3"], RowBox[{"1", "-", SubscriptBox["a", "1"], "-", RowBox[{"2", " ", SubscriptBox["a", "2"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "3"]}]}], RowBox[{ RowBox[{"-", "2"}], "+", RowBox[{"2", " ", SubscriptBox["a", "1"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "2"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "3"]}]}], RowBox[{"2", "-", RowBox[{"2", " ", SubscriptBox["a", "1"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "2"]}], "-", SubscriptBox["a", "3"]}], "0"}, { SubscriptBox["a", "4"], SubscriptBox["a", "5"], SubscriptBox["a", "6"], RowBox[{"4", "-", SubscriptBox["a", "4"], "-", RowBox[{"2", " ", SubscriptBox["a", "5"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "6"]}]}], RowBox[{ RowBox[{"-", "9"}], "+", RowBox[{"2", " ", SubscriptBox["a", "4"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "5"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "6"]}]}], RowBox[{"12", "-", RowBox[{"2", " ", SubscriptBox["a", "4"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "5"]}], "-", SubscriptBox["a", "6"]}], "6"}, { SubscriptBox["a", "7"], SubscriptBox["a", "8"], SubscriptBox["a", "9"], RowBox[{"6", "-", SubscriptBox["a", "7"], "-", RowBox[{"2", " ", SubscriptBox["a", "8"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "9"]}]}], RowBox[{ RowBox[{"-", "14"}], "+", RowBox[{"2", " ", SubscriptBox["a", "7"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "8"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "9"]}]}], RowBox[{"21", "-", RowBox[{"2", " ", SubscriptBox["a", "7"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "8"]}], "-", SubscriptBox["a", "9"]}], "12"}, { SubscriptBox["a", "10"], SubscriptBox["a", "11"], SubscriptBox["a", "12"], RowBox[{"5", "-", SubscriptBox["a", "10"], "-", RowBox[{"2", " ", SubscriptBox["a", "11"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "12"]}]}], RowBox[{ RowBox[{"-", "12"}], "+", RowBox[{"2", " ", SubscriptBox["a", "10"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "11"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "12"]}]}], RowBox[{"20", "-", RowBox[{"2", " ", SubscriptBox["a", "10"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "11"]}], "-", SubscriptBox["a", "12"]}], "12"}, { SubscriptBox["a", "13"], SubscriptBox["a", "14"], SubscriptBox["a", "15"], RowBox[{"2", "-", SubscriptBox["a", "13"], "-", RowBox[{"2", " ", SubscriptBox["a", "14"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "15"]}]}], RowBox[{ RowBox[{"-", "5"}], "+", RowBox[{"2", " ", SubscriptBox["a", "13"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "14"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "15"]}]}], RowBox[{"10", "-", RowBox[{"2", " ", SubscriptBox["a", "13"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "14"]}], "-", SubscriptBox["a", "15"]}], "6"}, { SubscriptBox["a", "16"], SubscriptBox["a", "17"], SubscriptBox["a", "18"], RowBox[{ RowBox[{"-", SubscriptBox["a", "16"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "17"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "18"]}]}], RowBox[{ RowBox[{"2", " ", SubscriptBox["a", "16"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "17"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "18"]}]}], RowBox[{"1", "-", RowBox[{"2", " ", SubscriptBox["a", "16"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "17"]}], "-", SubscriptBox["a", "18"]}], "0"}, { SubscriptBox["a", "19"], SubscriptBox["a", "20"], SubscriptBox["a", "21"], RowBox[{ RowBox[{"-", SubscriptBox["a", "19"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "20"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "21"]}]}], RowBox[{ RowBox[{"2", " ", SubscriptBox["a", "19"]}], "+", RowBox[{"3", " ", SubscriptBox["a", "20"]}], "+", RowBox[{"2", " ", SubscriptBox["a", "21"]}]}], RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", SubscriptBox["a", "19"]}], "-", RowBox[{"2", " ", SubscriptBox["a", "20"]}], "-", SubscriptBox["a", "21"]}], "1"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]}, { TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0", "0", "0", "1", RowBox[{"-", "2"}], "2", "0"}, {"0", "0", "0", "4", RowBox[{"-", "9"}], "12", "6"}, {"0", "0", "0", "6", RowBox[{"-", "14"}], "21", "12"}, {"0", "0", "0", "5", RowBox[{"-", "12"}], "20", "12"}, {"0", "0", "0", "2", RowBox[{"-", "5"}], "10", "6"}, {"0", "0", "0", "0", "0", "1", "0"}, {"0", "0", "0", "0", "0", "0", "1"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]} }, DefaultBaseStyle->"Column", GridBoxAlignment->{"Columns" -> {{Left}}}, GridBoxItemSize->{ "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}], "Column"]}], "}"}]}, { TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{"-", "1"}], "1", "5", "7", "5", "1", "1", "0", RowBox[{"2", " ", SqrtBox["3"]}], RowBox[{"4", " ", SqrtBox["3"]}], RowBox[{"4", " ", SqrtBox["3"]}], "0", RowBox[{"2", " ", SqrtBox["3"]}], "0"}, {"1", RowBox[{"-", "1"}], "1", "5", "7", "5", "1", "0", "0", RowBox[{"2", " ", SqrtBox["3"]}], RowBox[{"4", " ", SqrtBox["3"]}], "0", RowBox[{"4", " ", SqrtBox["3"]}], RowBox[{"2", " ", SqrtBox["3"]}]}, {"5", "1", RowBox[{"-", "1"}], "1", "5", "7", "1", RowBox[{"2", " ", SqrtBox["3"]}], "0", "0", RowBox[{"2", " ", SqrtBox["3"]}], "0", RowBox[{"4", " ", SqrtBox["3"]}], RowBox[{"4", " ", SqrtBox["3"]}]}, {"7", "5", "1", RowBox[{"-", "1"}], "1", "5", "1", RowBox[{"4", " ", SqrtBox["3"]}], RowBox[{"2", " ", SqrtBox["3"]}], "0", "0", "0", RowBox[{"2", " ", SqrtBox["3"]}], RowBox[{"4", " ", SqrtBox["3"]}]}, {"5", "7", "5", "1", RowBox[{"-", "1"}], "1", "1", RowBox[{"4", " ", SqrtBox["3"]}], RowBox[{"4", " ", SqrtBox["3"]}], RowBox[{"2", " ", SqrtBox["3"]}], "0", "0", "0", RowBox[{"2", " ", SqrtBox["3"]}]}, {"1", "5", "7", "5", "1", RowBox[{"-", "1"}], "1", RowBox[{"2", " ", SqrtBox["3"]}], RowBox[{"4", " ", SqrtBox["3"]}], RowBox[{"4", " ", SqrtBox["3"]}], RowBox[{"2", " ", SqrtBox["3"]}], "0", "0", "0"}, {"1", "1", "1", "1", "1", "1", RowBox[{"-", "1"}], "0", "0", "0", "0", FractionBox["2", SqrtBox["3"]], "0", "0"}, {"0", "0", RowBox[{"2", " ", SqrtBox["3"]}], RowBox[{"4", " ", SqrtBox["3"]}], RowBox[{"4", " ", SqrtBox["3"]}], RowBox[{"2", " ", SqrtBox["3"]}], "0", RowBox[{"-", "1"}], "1", "5", "7", "1", "5", "1"}, { RowBox[{"2", " ", SqrtBox["3"]}], "0", "0", RowBox[{"2", " ", SqrtBox["3"]}], RowBox[{"4", " ", SqrtBox["3"]}], RowBox[{"4", " ", SqrtBox["3"]}], "0", "1", RowBox[{"-", "1"}], "1", "5", "1", "7", "5"}, { RowBox[{"4", " ", SqrtBox["3"]}], RowBox[{"2", " ", SqrtBox["3"]}], "0", "0", RowBox[{"2", " ", SqrtBox["3"]}], RowBox[{"4", " ", SqrtBox["3"]}], "0", "5", "1", RowBox[{"-", "1"}], "1", "1", "5", "7"}, { RowBox[{"4", " ", SqrtBox["3"]}], RowBox[{"4", " ", SqrtBox["3"]}], RowBox[{"2", " ", SqrtBox["3"]}], "0", "0", RowBox[{"2", " ", SqrtBox["3"]}], "0", "7", "5", "1", RowBox[{"-", "1"}], "1", "1", "5"}, {"0", "0", "0", "0", "0", "0", FractionBox["2", SqrtBox["3"]], "1", "1", "1", "1", RowBox[{"-", "1"}], "1", "1"}, { RowBox[{"2", " ", SqrtBox["3"]}], RowBox[{"4", " ", SqrtBox["3"]}], RowBox[{"4", " ", SqrtBox["3"]}], RowBox[{"2", " ", SqrtBox["3"]}], "0", "0", "0", "5", "7", "5", "1", "1", RowBox[{"-", "1"}], "1"}, {"0", RowBox[{"2", " ", SqrtBox["3"]}], RowBox[{"4", " ", SqrtBox["3"]}], RowBox[{"4", " ", SqrtBox["3"]}], RowBox[{"2", " ", SqrtBox["3"]}], "0", "0", "1", "5", "7", "5", "1", "1", RowBox[{"-", "1"}]} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]} }, DefaultBaseStyle->"Column", GridBoxAlignment->{"Columns" -> {{Left}}}, GridBoxItemSize->{"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}], "Column"]], "Output", CellChangeTimes->{3.7388725601667094`*^9}, CellLabel->"Out[76]=",ExpressionUUID->"d341dffc-22e0-491e-9c00-a61d51a92e7c"] }, Open ]] }, WindowSize->{1060, 661}, WindowMargins->{{Automatic, -1}, {Automatic, 37}}, Magnification:>0.95 Inherited, FrontEndVersion->"11.3 for Microsoft Windows (64-bit) (March 6, 2018)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[558, 20, 665, 18, 50, "Input",ExpressionUUID->"748eabe5-c3d2-4801-b8b7-aa7a789097c7"], Cell[1226, 40, 696, 18, 47, "Input",ExpressionUUID->"0ab29ce4-3405-44f4-bf4a-16ef730ff686"], Cell[1925, 60, 9876, 235, 649, "Input",ExpressionUUID->"72e544ec-21b3-4511-8714-78e40bf96f8f"], Cell[11804, 297, 35697, 782, 1936, "Input",ExpressionUUID->"34558bf5-30de-4efd-a85d-1e42927d2f41"], Cell[47504, 1081, 13767, 317, 771, "Input",ExpressionUUID->"1669407f-19e4-41f4-83d2-404c097a4b8c"], Cell[61274, 1400, 1573, 43, 82, "Input",ExpressionUUID->"0ca661cb-0814-4baa-8045-ccb6be11fd34"], Cell[62850, 1445, 16049, 333, 766, "Input",ExpressionUUID->"8eb99bd7-ecdf-447b-9254-91418df0046a"], Cell[78902, 1780, 13583, 313, 776, "Input",ExpressionUUID->"df3b60f9-0301-4a00-901e-3f42ad1f6dcc"], Cell[92488, 2095, 3769, 87, 262, "Input",ExpressionUUID->"94e58c8f-e9bc-4bc5-9f72-1e7d65b723b4"], Cell[96260, 2184, 397, 13, 72, "Input",ExpressionUUID->"e0269727-13f8-41e1-9bed-e3d214f9ade5"], Cell[96660, 2199, 682, 20, 53, "Input",ExpressionUUID->"335fc47c-61a3-4d42-b8dc-1a6e429b3dbb"], Cell[97345, 2221, 4976, 131, 389, "Input",ExpressionUUID->"8f83b82e-c1db-40d9-ac4c-a723737f7704"], Cell[102324, 2354, 1032, 29, 27, "Input",ExpressionUUID->"d37cc6ff-e0b5-488e-9c86-c606b2039282"], Cell[103359, 2385, 1960, 44, 118, "Input",ExpressionUUID->"9a3eb7b4-a69e-4770-9009-32f6122d8145"], Cell[105322, 2431, 5279, 132, 172, "Input",ExpressionUUID->"affbea28-f93e-449d-873a-65b52e624ef9"], Cell[110604, 2565, 1181, 31, 46, "Input",ExpressionUUID->"15495297-a588-42f2-9830-ff3a5a9d76a4"], Cell[111788, 2598, 4112, 110, 336, "Input",ExpressionUUID->"cd35cdcc-79cb-420e-8b24-2c510cf42e01"], Cell[115903, 2710, 5326, 123, 316, "Input",ExpressionUUID->"b49300f9-4eb0-4e74-b00e-238ad2cc61d5"], Cell[121232, 2835, 3798, 88, 226, "Input",ExpressionUUID->"b43a5c15-ff8f-4841-a157-1d63bba9b7c4"], Cell[CellGroupData[{ Cell[125055, 2927, 1384, 29, 27, "Input",ExpressionUUID->"8c0d00eb-3a38-4d87-9e21-66e5a04e235e"], Cell[126442, 2958, 108310, 2552, 2624, "Output",ExpressionUUID->"d341dffc-22e0-491e-9c00-a61d51a92e7c"] }, Open ]] } ] *)