(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 12.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 158, 7] NotebookDataLength[ 49982, 1163] NotebookOptionsPosition[ 49389, 1145] NotebookOutlinePosition[ 49796, 1161] CellTagsIndexPosition[ 49753, 1158] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[TextData[{ StyleBox["Weighted ", FontColor->RGBColor[1, 0, 0]], StyleBox["k", FontSlant->"Italic", FontColor->RGBColor[1, 0, 0]], StyleBox["-means algorithm", FontColor->RGBColor[1, 0, 0]] }], "Subsubsection", CellChangeTimes->{{3.73112424028771*^9, 3.7311242546173754`*^9}, { 3.8106047392458076`*^9, 3.8106047456549196`*^9}, {3.8131286886441045`*^9, 3.8131286920600576`*^9}, 3.8209872570601444`*^9, 3.821011326978313*^9},ExpressionUUID->"e66d6e93-d6b2-4c5b-9154-\ 8959b5b959a5"], Cell[BoxData[ RowBox[{" ", RowBox[{ RowBox[{ RowBox[{ StyleBox["Input", FontColor->RGBColor[1, 0, 0]], StyleBox[":", FontColor->RGBColor[1, 0, 0]], StyleBox[" ", FontColor->RGBColor[0, 0, 1]], StyleBox[ RowBox[{ RowBox[{"data", " ", "set", " ", "\[ScriptCapitalA]"}], " ", "\[Subset]", " ", SuperscriptBox["\[DoubleStruckCapitalR]", "n"]}], FontColor->RGBColor[0, 0, 1]]}], StyleBox[",", FontColor->RGBColor[0, 0, 1]], StyleBox[" ", FontColor->RGBColor[0, 0, 1]], StyleBox[ RowBox[{"data", " ", "weights", " ", "w"}], FontColor->RGBColor[0, 0, 1]], StyleBox[",", FontColor->RGBColor[0, 0, 1]], StyleBox[" ", FontColor->RGBColor[0, 0, 1]], StyleBox[ RowBox[{ RowBox[{"initial", " ", "circle", " ", "centers", " ", "Circ"}], ";"}], FontColor->RGBColor[0, 0, 1]]}], StyleBox["\[IndentingNewLine]", FontColor->RGBColor[0, 0, 1]], RowBox[{ RowBox[{ StyleBox["Output", FontColor->RGBColor[1, 0, 0]], StyleBox[":", FontColor->RGBColor[1, 0, 0]], StyleBox[" ", FontColor->RGBColor[0, 0, 1]], StyleBox[ RowBox[{"partition", " ", "\[CapitalPi]"}], FontColor->RGBColor[0, 0, 1]]}], StyleBox["=", FontColor->RGBColor[0, 0, 1]], StyleBox[ RowBox[{ RowBox[{ RowBox[{"{", RowBox[{ SubscriptBox["\[Pi]", "1"], ",", "\[TripleDot]", ",", SubscriptBox["\[Pi]", "k"]}], "}"}], " ", "of", " ", "the", " ", "set", " ", "\[ScriptCapitalA]"}], " ", "\[Subset]", " ", SuperscriptBox["\[DoubleStruckCapitalR]", "n"]}], FontColor->RGBColor[0, 0, 1]]}], StyleBox["\[IndentingNewLine]", FontColor->RGBColor[0, 0, 1]], StyleBox[" ", FontColor->RGBColor[0, 0, 1]], RowBox[{ StyleBox[ RowBox[{"(", RowBox[{ RowBox[{ "associated", " ", "to", " ", "each", " ", "element", " ", "of", " ", "the", " ", "cluster", " ", SubscriptBox["\[Pi]", "j"]}], "\[Element]", RowBox[{ "\[CapitalPi]", " ", "is", " ", "the", " ", "corresponding", " ", "weight"}]}], ")"}], FontColor->RGBColor[0, 0, 1]], StyleBox[",", FontColor->RGBColor[0, 0, 1]], StyleBox[" ", FontColor->RGBColor[0, 0, 1]], RowBox[{ StyleBox["generalized", FontColor->RGBColor[0, 0, 1]], StyleBox[ StyleBox[ RowBox[{" ", " "}]], FontColor->RGBColor[0, 0, 1]], StyleBox["circle", FontColor->RGBColor[0, 0, 1]], StyleBox[" ", FontColor->RGBColor[0, 0, 1]], "centers", " ", SubscriptBox["Cir", "1"]}], StyleBox[",", FontColor->RGBColor[0, 0, 1]], StyleBox["\[TripleDot]", FontColor->RGBColor[0, 0, 1]], ",", SubscriptBox["Cir", "k"], StyleBox[",", FontColor->RGBColor[0, 0, 1]], StyleBox[" ", FontColor->RGBColor[0, 0, 1]], StyleBox[ RowBox[{ RowBox[{"number", " ", "of", " ", "iterations", " ", StyleBox["It", FontSlant->"Italic"]}], ";"}], FontColor->RGBColor[0, 0, 1]]}], StyleBox["\[IndentingNewLine]", FontColor->RGBColor[0, 0, 1]], RowBox[{ RowBox[{ StyleBox["Submodules", FontColor->RGBColor[1, 0, 0]], StyleBox[":", FontColor->RGBColor[1, 0, 0]], StyleBox[" ", FontColor->RGBColor[0, 0, 1]], StyleBox[ RowBox[{"DIRECT", "[", "]"}], FontColor->RGBColor[0, 0, 1]]}], StyleBox[",", FontColor->RGBColor[0, 0, 1]], StyleBox[" ", FontColor->RGBColor[0, 0, 1]], StyleBox[ RowBox[{"dist", "[", "]"}], FontColor->RGBColor[0, 0, 1]], StyleBox[",", FontColor->RGBColor[0, 0, 1]], StyleBox[" ", FontColor->RGBColor[0, 0, 1]], StyleBox[ RowBox[{"Dmn", "[", "]"}], FontColor->RGBColor[0, 0, 1]], StyleBox[",", FontColor->RGBColor[0, 0, 1]], StyleBox[" ", FontColor->RGBColor[0, 0, 1]], StyleBox[ RowBox[{ RowBox[{"WGLS", "[", "]"}], ";"}], FontColor->RGBColor[0, 0, 1]]}], StyleBox[" ", FontColor->RGBColor[0, 0, 1]], "\[IndentingNewLine]"}]}]], "Text", CellChangeTimes->{{3.7311239388490767`*^9, 3.7311241887665453`*^9}, { 3.731124262852821*^9, 3.7311242818390455`*^9}, {3.731124334688675*^9, 3.7311243931480875`*^9}, {3.731124464437002*^9, 3.7311245820590296`*^9}, { 3.731124626641944*^9, 3.7311246685996666`*^9}, 3.73112478187739*^9, 3.7311249716632557`*^9, 3.731125008229858*^9, {3.731131363164089*^9, 3.7311315307912507`*^9}, {3.8106047626816826`*^9, 3.810604763185858*^9}, { 3.810606445967278*^9, 3.810606530225013*^9}, {3.810606641466092*^9, 3.810606684833495*^9}, {3.810606749868683*^9, 3.810606752984932*^9}, { 3.8106068274532113`*^9, 3.810606866094061*^9}, {3.810606901083395*^9, 3.810607019637113*^9}, {3.810607054277118*^9, 3.8106071829172134`*^9}, { 3.8106072199097958`*^9, 3.8106072546405845`*^9}, 3.8106074697262*^9, { 3.8131282699269085`*^9, 3.813128270182798*^9}, {3.8131283217992973`*^9, 3.8131283433443713`*^9}, {3.8131285602764816`*^9, 3.813128631734671*^9}, { 3.813128716159551*^9, 3.8131287310664387`*^9}, {3.813128829720832*^9, 3.8131288437367177`*^9}, {3.8131288968861256`*^9, 3.8131288976456766`*^9}, {3.8208206004327154`*^9, 3.820820602339464*^9}, { 3.8208232348027654`*^9, 3.820823235970901*^9}, {3.8209872382550364`*^9, 3.820987241863156*^9}, {3.8209873121364293`*^9, 3.820987370577606*^9}, { 3.820987470767851*^9, 3.8209875036173415`*^9}, {3.8209875433950415`*^9, 3.8209876048601007`*^9}, {3.82098765240862*^9, 3.820987699634896*^9}, 3.8209877620138216`*^9, {3.820987851324072*^9, 3.820987889081056*^9}, { 3.820987919084846*^9, 3.8209879407206955`*^9}, {3.8209879711907043`*^9, 3.8209882843502274`*^9}, {3.8223173587958217`*^9, 3.822317368644164*^9}, { 3.8223174160400667`*^9, 3.8223174317795243`*^9}, {3.8224885876165414`*^9, 3.8224887155109468`*^9}}, FontSize->12,ExpressionUUID->"a16ffa28-7742-4592-bcc2-c1a697822085"], Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "Circ", " ", "are", " ", "initial", " ", "approximations", " ", "for", " ", "generalized", " ", "circles", "\[IndentingNewLine]", "It", " ", "is", " ", "used", " ", "DIRECT", " ", "algorithm", " ", "without", " ", "Compile"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"WKGC", "[", RowBox[{"w_", ",", "A_", ",", "Circ_", ",", "Iter_"}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"It", "=", "0"}], ",", RowBox[{"m", "=", RowBox[{"Length", "[", "w", "]"}]}], ",", RowBox[{"k", "=", RowBox[{"Length", "[", "Circ", "]"}]}], ",", "Kr", ",", "PI", ",", "tab", ",", "fc", ",", "min", ",", "imin", ",", "Pod", ",", "F0", ",", "F1", ",", "p", ",", "q", ",", "r", ",", "dom", ",", "min1", ",", "fCompile", ",", "FF", ",", RowBox[{"EPS", "=", "0.5"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"Kr", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Pod", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"w", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"A", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "m"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", "\t\t\t\t\t\t", RowBox[{"(*", " ", RowBox[{"Odredjivanje", " ", "particije"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"PI", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "k"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Do", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"tab", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"DGcir", "[", RowBox[{ RowBox[{"Pod", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", RowBox[{"Circ", "[", RowBox[{"[", RowBox[{"j", ",", "1"}], "]"}], "]"}], ",", RowBox[{"Circ", "[", RowBox[{"[", RowBox[{"j", ",", "2"}], "]"}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}]}], "]"}], "//", "N"}]}], ";", "\[IndentingNewLine]", RowBox[{"imin", "=", RowBox[{ RowBox[{"Ordering", "[", RowBox[{"tab", ",", "1"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", "imin", "]"}], "]"}], "=", RowBox[{"Append", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", "imin", "]"}], "]"}], ",", RowBox[{"Pod", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}]}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"i", ",", "m"}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"PI", "=", RowBox[{"DeleteCases", "[", RowBox[{"PI", ",", RowBox[{"{", "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"k", "=", RowBox[{"Length", "[", "PI", "]"}]}], ";", " ", "\[IndentingNewLine]", RowBox[{"F1", "=", RowBox[{"WGLS", "[", RowBox[{"PI", ",", "Circ"}], "]"}]}], ";", " ", "\[IndentingNewLine]", " ", RowBox[{"(*", " ", RowBox[{ RowBox[{"Calculation", " ", "of", " ", "circle"}], "-", "centers"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"Do", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"c0", "=", RowBox[{"Mean", "[", RowBox[{"WeightedData", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"j", ",", "All", ",", "2"}], "]"}], "]"}], ",", RowBox[{"PI", "[", RowBox[{"[", RowBox[{"j", ",", "All", ",", "1"}], "]"}], "]"}]}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Cov", "=", RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"j", ",", "s", ",", "1"}], "]"}], "]"}], RowBox[{"KroneckerProduct", "[", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"j", ",", "s", ",", "2"}], "]"}], "]"}], "-", "c0"}], ")"}], ",", RowBox[{"(", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"j", ",", "s", ",", "2"}], "]"}], "]"}], "-", "c0"}], ")"}]}], "]"}]}], ",", RowBox[{"{", RowBox[{"s", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}]}], "}"}]}], "]"}], "/", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"v", "=", RowBox[{ RowBox[{"Eigensystem", "[", "Cov", "]"}], "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"d0", "=", RowBox[{"2", RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{"Abs", "[", RowBox[{"Dot", "[", RowBox[{ RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"j", ",", "s", ",", "2"}], "]"}], "]"}], "-", "c0"}], ",", "v"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"s", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}]}], "}"}]}], "]"}], "/", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{"seg0", "=", RowBox[{"{", RowBox[{ RowBox[{"c0", "-", RowBox[{".5", "d0", " ", "v"}]}], ",", RowBox[{"c0", "+", RowBox[{".5", "d0", " ", "v"}]}]}], "}"}]}], ";", "\[IndentingNewLine]", RowBox[{"r0", "=", RowBox[{"Mean", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"dist", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"j", ",", "s", ",", "2"}], "]"}], "]"}], ",", "seg0"}], "]"}], ",", RowBox[{"{", RowBox[{"s", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}]}], "}"}]}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"Print", "[", RowBox[{"j", ",", "\"\<. appr: \>\"", ",", RowBox[{"{", RowBox[{"seg0", ",", "r0"}], "}"}]}], "]"}], ";"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"lam", "=", RowBox[{"d0", "/", "4"}]}], ";", "\[IndentingNewLine]", RowBox[{"dom", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"seg0", "[", RowBox[{"[", RowBox[{"1", ",", "1"}], "]"}], "]"}], "-", "lam"}], ",", RowBox[{ RowBox[{"seg0", "[", RowBox[{"[", RowBox[{"1", ",", "1"}], "]"}], "]"}], "+", "lam"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"seg0", "[", RowBox[{"[", RowBox[{"1", ",", "2"}], "]"}], "]"}], "-", "lam"}], ",", RowBox[{ RowBox[{"seg0", "[", RowBox[{"[", RowBox[{"1", ",", "2"}], "]"}], "]"}], "+", "lam"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"seg0", "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}], "-", "lam"}], ",", RowBox[{ RowBox[{"seg0", "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}], "+", "lam"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"seg0", "[", RowBox[{"[", RowBox[{"2", ",", "2"}], "]"}], "]"}], "-", "lam"}], ",", RowBox[{ RowBox[{"seg0", "[", RowBox[{"[", RowBox[{"2", ",", "2"}], "]"}], "]"}], "+", "lam"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"r0", "-", "lam"}], ",", RowBox[{"r0", "+", "lam"}]}], "}"}]}], "}"}]}], ";", RowBox[{"(*", " ", RowBox[{ RowBox[{"Print", "[", RowBox[{"j", ",", "\"\<. dom: \>\"", ",", "dom"}], "]"}], ";"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"distGC", "[", RowBox[{"xm_", ",", "ym_", ",", "xn_", ",", "yn_", ",", "r_"}], "]"}], ":=", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"Dmn", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"j", ",", "s", ",", "2"}], "]"}], "]"}], ",", "xm", ",", "ym", ",", "xn", ",", "yn"}], "]"}], "^", "2"}], "-", RowBox[{"r", "^", "2"}]}], ")"}], "^", "2"}], ",", RowBox[{"{", RowBox[{"s", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"FCompile", "=", RowBox[{"Compile", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"z1", ",", "_Real"}], "}"}], ",", RowBox[{"{", RowBox[{"z2", ",", "_Real"}], "}"}], ",", RowBox[{"{", RowBox[{"z3", ",", "_Real"}], "}"}], ",", RowBox[{"{", RowBox[{"z4", ",", "_Real"}], "}"}], ",", RowBox[{"{", RowBox[{"z5", ",", "_Real"}], "}"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"Evaluate", "[", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"Dmn", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"j", ",", "s", ",", "2"}], "]"}], "]"}], ",", "z1", ",", "z2", ",", "z3", ",", "z4"}], "]"}], "^", "2"}], "-", RowBox[{"z5", "^", "2"}]}], ")"}], "^", "2"}], ",", RowBox[{"{", RowBox[{"s", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}]}], "}"}]}], "]"}], "]"}]}], "]"}]}], ";", "\t", "\[IndentingNewLine]", RowBox[{ RowBox[{"distGC", "[", RowBox[{"z1_", ",", "z2_", ",", "z3_", ",", "z4_", ",", "z5_"}], "]"}], ":=", RowBox[{"Apply", "[", RowBox[{"FCompile", ",", RowBox[{"{", RowBox[{"z1", ",", "z2", ",", "z3", ",", "z4", ",", "z5"}], "}"}]}], "]"}]}], ";"}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"Print", "[", RowBox[{"distGC", "[", RowBox[{"j", ",", RowBox[{"seg0", "[", RowBox[{"[", RowBox[{"1", ",", "1"}], "]"}], "]"}], ",", RowBox[{"seg0", "[", RowBox[{"[", RowBox[{"1", ",", "2"}], "]"}], "]"}], ",", RowBox[{"seg0", "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}], ",", RowBox[{"seg0", "[", RowBox[{"[", RowBox[{"2", ",", "2"}], "]"}], "]"}], ",", "r0"}], "]"}], "]"}], ";"}], "*)"}], "\[IndentingNewLine]", RowBox[{"tDIR", "=", RowBox[{"Timing", "[", RowBox[{"min", "=", RowBox[{"DIRECT", "[", RowBox[{"distGC", ",", "dom", ",", RowBox[{"DMin", "\[Rule]", RowBox[{"10.", "^", RowBox[{"(", RowBox[{"-", "4"}], ")"}]}]}], ",", RowBox[{"MaxIter", "\[Rule]", "30"}], ",", " ", RowBox[{"Eps", "\[Rule]", RowBox[{"10", "^", RowBox[{"(", RowBox[{"-", "4"}], ")"}]}]}]}], "]"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"Kr", "[", RowBox[{"[", "j", "]"}], "]"}], "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"min", "[", RowBox[{"[", RowBox[{"1", ",", "1"}], "]"}], "]"}], ",", RowBox[{"min", "[", RowBox[{"[", RowBox[{"1", ",", "2"}], "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"min", "[", RowBox[{"[", RowBox[{"1", ",", "3"}], "]"}], "]"}], ",", RowBox[{"min", "[", RowBox[{"[", RowBox[{"1", ",", "4"}], "]"}], "]"}]}], "}"}]}], "}"}], ",", RowBox[{"min", "[", RowBox[{"[", RowBox[{"1", ",", "5"}], "]"}], "]"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", ",", RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"F0", "=", RowBox[{"WGLS", "[", RowBox[{"PI", ",", "Kr"}], "]"}]}], ";", " ", "\[IndentingNewLine]", RowBox[{"While", "[", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{ RowBox[{"Chop", "[", RowBox[{"F0", "-", "F1"}], "]"}], "\[NotEqual]", "0"}], "&&", RowBox[{"It", "<", "Iter"}]}], ","}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"Abs", "[", RowBox[{"F0", "-", "F1"}], "]"}], "/", "F0"}], ">", "EPS"}], "&&", RowBox[{"It", "<", "Iter"}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"F0", "=", "F1"}], ";", RowBox[{"It", "=", RowBox[{"It", "+", "1"}]}], ";", "\[IndentingNewLine]", RowBox[{"PI", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "k"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Do", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"tab", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"DGcir", "[", RowBox[{ RowBox[{"Pod", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", RowBox[{"Kr", "[", RowBox[{"[", RowBox[{"j", ",", "1"}], "]"}], "]"}], ",", RowBox[{"Kr", "[", RowBox[{"[", RowBox[{"j", ",", "2"}], "]"}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}]}], "]"}], "//", "N"}]}], ";", "\[IndentingNewLine]", RowBox[{"imin", "=", RowBox[{ RowBox[{"Ordering", "[", RowBox[{"tab", ",", "1"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", "imin", "]"}], "]"}], "=", RowBox[{"Append", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", "imin", "]"}], "]"}], ",", RowBox[{"Pod", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}]}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"i", ",", "m"}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"PI", "=", RowBox[{"DeleteCases", "[", RowBox[{"PI", ",", RowBox[{"{", "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"k", "=", RowBox[{"Length", "[", "PI", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Kr", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Do", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"c0", "=", RowBox[{"Mean", "[", RowBox[{"WeightedData", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"j", ",", "All", ",", "2"}], "]"}], "]"}], ",", RowBox[{"PI", "[", RowBox[{"[", RowBox[{"j", ",", "All", ",", "1"}], "]"}], "]"}]}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Cov", "=", RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"j", ",", "s", ",", "1"}], "]"}], "]"}], RowBox[{"KroneckerProduct", "[", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"j", ",", "s", ",", "2"}], "]"}], "]"}], "-", "c0"}], ")"}], ",", RowBox[{"(", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"j", ",", "s", ",", "2"}], "]"}], "]"}], "-", "c0"}], ")"}]}], "]"}]}], ",", RowBox[{"{", RowBox[{"s", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}]}], "}"}]}], "]"}], "/", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"v", "=", RowBox[{ RowBox[{"Eigensystem", "[", "Cov", "]"}], "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"d0", "=", RowBox[{"2", RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{"Abs", "[", RowBox[{"Dot", "[", RowBox[{ RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"j", ",", "s", ",", "2"}], "]"}], "]"}], "-", "c0"}], ",", "v"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"s", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}]}], "}"}]}], "]"}], "/", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{"seg0", "=", RowBox[{"{", RowBox[{ RowBox[{"c0", "-", RowBox[{".5", "d0", " ", "v"}]}], ",", RowBox[{"c0", "+", RowBox[{".5", "d0", " ", "v"}]}]}], "}"}]}], ";", "\[IndentingNewLine]", RowBox[{"r0", "=", RowBox[{"Mean", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"dist", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"j", ",", "s", ",", "2"}], "]"}], "]"}], ",", "seg0"}], "]"}], ",", RowBox[{"{", RowBox[{"s", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}]}], "}"}]}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"Print", "[", RowBox[{"j", ",", "\"\<. appr: \>\"", ",", RowBox[{"{", RowBox[{"seg0", ",", "r0"}], "}"}]}], "]"}], ";"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"lam", "=", RowBox[{"d0", "/", "4"}]}], ";", "\[IndentingNewLine]", RowBox[{"dom", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"seg0", "[", RowBox[{"[", RowBox[{"1", ",", "1"}], "]"}], "]"}], "-", "lam"}], ",", RowBox[{ RowBox[{"seg0", "[", RowBox[{"[", RowBox[{"1", ",", "1"}], "]"}], "]"}], "+", "lam"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"seg0", "[", RowBox[{"[", RowBox[{"1", ",", "2"}], "]"}], "]"}], "-", "lam"}], ",", RowBox[{ RowBox[{"seg0", "[", RowBox[{"[", RowBox[{"1", ",", "2"}], "]"}], "]"}], "+", "lam"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"seg0", "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}], "-", "lam"}], ",", RowBox[{ RowBox[{"seg0", "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}], "+", "lam"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"seg0", "[", RowBox[{"[", RowBox[{"2", ",", "2"}], "]"}], "]"}], "-", "lam"}], ",", RowBox[{ RowBox[{"seg0", "[", RowBox[{"[", RowBox[{"2", ",", "2"}], "]"}], "]"}], "+", "lam"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"r0", "-", "lam"}], ",", RowBox[{"r0", "+", "lam"}]}], "}"}]}], "}"}]}], ";", RowBox[{"(*", " ", RowBox[{ RowBox[{"Print", "[", RowBox[{"j", ",", "\"\<. dom: \>\"", ",", "dom"}], "]"}], ";"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"distGC", "[", RowBox[{"xm_", ",", "ym_", ",", "xn_", ",", "yn_", ",", "r_"}], "]"}], ":=", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"Dmn", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"j", ",", "s", ",", "2"}], "]"}], "]"}], ",", "xm", ",", "ym", ",", "xn", ",", "yn"}], "]"}], "^", "2"}], "-", RowBox[{"r", "^", "2"}]}], ")"}], "^", "2"}], ",", RowBox[{"{", RowBox[{"s", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"FCompile", "=", RowBox[{"Compile", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"z1", ",", "_Real"}], "}"}], ",", RowBox[{"{", RowBox[{"z2", ",", "_Real"}], "}"}], ",", RowBox[{"{", RowBox[{"z3", ",", "_Real"}], "}"}], ",", RowBox[{"{", RowBox[{"z4", ",", "_Real"}], "}"}], ",", RowBox[{"{", RowBox[{"z5", ",", "_Real"}], "}"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"Evaluate", "[", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"Dmn", "[", RowBox[{ RowBox[{"PI", "[", RowBox[{"[", RowBox[{"j", ",", "s", ",", "2"}], "]"}], "]"}], ",", "z1", ",", "z2", ",", "z3", ",", "z4"}], "]"}], "^", "2"}], "-", RowBox[{"z5", "^", "2"}]}], ")"}], "^", "2"}], ",", RowBox[{"{", RowBox[{"s", ",", RowBox[{"Length", "[", RowBox[{"PI", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}]}], "}"}]}], "]"}], "]"}]}], "]"}]}], ";", "\t", "\[IndentingNewLine]", "distGCz1_"}], ",", "z2_", ",", "z3_", ",", "z4_", ",", "z5_"}], "]"}], ":=", RowBox[{"Apply", "[", RowBox[{"FCompile", ",", RowBox[{"{", RowBox[{"z1", ",", "z2", ",", "z3", ",", "z4", ",", "z5"}], "}"}]}], "]"}]}], ";"}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"Print", "[", RowBox[{"distGC", "[", RowBox[{"j", ",", RowBox[{"seg0", "[", RowBox[{"[", RowBox[{"1", ",", "1"}], "]"}], "]"}], ",", RowBox[{"seg0", "[", RowBox[{"[", RowBox[{"1", ",", "2"}], "]"}], "]"}], ",", RowBox[{"seg0", "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}], ",", RowBox[{"seg0", "[", RowBox[{"[", RowBox[{"2", ",", "2"}], "]"}], "]"}], ",", "r0"}], "]"}], "]"}], ";"}], "*)"}], "\[IndentingNewLine]", RowBox[{"tDIR", "=", RowBox[{"Timing", "[", RowBox[{"min", "=", RowBox[{"DIRECT", "[", RowBox[{"distGC", ",", "dom", ",", RowBox[{"DMin", "\[Rule]", RowBox[{"10.", "^", RowBox[{"(", RowBox[{"-", "4"}], ")"}]}]}], ",", RowBox[{"MaxIter", "\[Rule]", "15"}], ",", " ", RowBox[{"Eps", "\[Rule]", RowBox[{"10", "^", RowBox[{"(", RowBox[{"-", "4"}], ")"}]}]}]}], "]"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"Kr", "[", RowBox[{"[", "j", "]"}], "]"}], "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"min", "[", RowBox[{"[", RowBox[{"1", ",", "1"}], "]"}], "]"}], ",", RowBox[{"min", "[", RowBox[{"[", RowBox[{"1", ",", "2"}], "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"min", "[", RowBox[{"[", RowBox[{"1", ",", "3"}], "]"}], "]"}], ",", RowBox[{"min", "[", RowBox[{"[", RowBox[{"1", ",", "4"}], "]"}], "]"}]}], "}"}]}], "}"}], ",", RowBox[{"min", "[", RowBox[{"[", RowBox[{"1", ",", "5"}], "]"}], "]"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"Print", "[", RowBox[{"Kr", "[", RowBox[{"[", "j", "]"}], "]"}], "]"}], ";"}], "*)"}], "\[IndentingNewLine]", ",", RowBox[{"{", RowBox[{"j", ",", "k"}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"F1", "=", RowBox[{"WGLS", "[", RowBox[{"PI", ",", "Kr"}], "]"}]}], ";"}]}], " ", "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"PI", ",", "Kr", ",", "F1", ",", "It"}], "}"}]}]}], "\[IndentingNewLine]", "]"}]}]}]], "Input", CellChangeTimes->{{3.525616360604874*^9, 3.5256163644892807`*^9}, 3.5256170782217345`*^9, 3.5262058902050314`*^9, 3.5262060212488317`*^9, 3.526206146804221*^9, {3.5262146495319996`*^9, 3.526214660301216*^9}, { 3.526268235603675*^9, 3.5262682603453183`*^9}, {3.526271497881405*^9, 3.526271504183816*^9}, {3.526272345602494*^9, 3.526272371124139*^9}, { 3.526272761748825*^9, 3.5262727798448563`*^9}, {3.5262733652982845`*^9, 3.5262734414576187`*^9}, {3.526273604665105*^9, 3.5262736247735405`*^9}, 3.526289526950795*^9, {3.5712988897963448`*^9, 3.5712989095705805`*^9}, { 3.5713006259192133`*^9, 3.5713007138277683`*^9}, {3.5713007531946383`*^9, 3.571300771948871*^9}, {3.5713008158719487`*^9, 3.5713008213017583`*^9}, { 3.5713010538319693`*^9, 3.571301059058978*^9}, {3.5713026879974575`*^9, 3.5713027468205614`*^9}, {3.571302831529111*^9, 3.571302841748129*^9}, { 3.5713028969822264`*^9, 3.571302908216246*^9}, {3.571303208091176*^9, 3.571303209541979*^9}, {3.5713032974143333`*^9, 3.5713033007069397`*^9}, { 3.5713034979916883`*^9, 3.5713034993176904`*^9}, {3.571303581655036*^9, 3.5713035867260447`*^9}, {3.571304578139395*^9, 3.5713046065042453`*^9}, { 3.571304657245535*^9, 3.571304680884577*^9}, {3.571310017656273*^9, 3.571310032448099*^9}, {3.571310291879758*^9, 3.571310305546382*^9}, { 3.5713118990221987`*^9, 3.5713118994433994`*^9}, {3.5713119570019007`*^9, 3.5713119573451014`*^9}, 3.5713122470594134`*^9, 3.571365674705989*^9, { 3.57136576008634*^9, 3.57136582210345*^9}, 3.5713658564929104`*^9, { 3.571365903925994*^9, 3.5713659117592077`*^9}, {3.571365953060481*^9, 3.571365988119743*^9}, 3.571366114458166*^9, {3.5713663292159452`*^9, 3.5713664037264786`*^9}, {3.5713664475247555`*^9, 3.5713664652181873`*^9}, {3.5814910655797153`*^9, 3.581491069245722*^9}, { 3.5814911427368517`*^9, 3.581491161320484*^9}, {3.5814912211107903`*^9, 3.5814912292302055`*^9}, 3.581491340924204*^9, 3.581491489466466*^9, { 3.5814941706831613`*^9, 3.581494179358776*^9}, {3.5814942296632657`*^9, 3.581494273082543*^9}, {3.5814943407558627`*^9, 3.581494362584301*^9}, { 3.581494398689765*^9, 3.581494410485386*^9}, {3.5814951461976876`*^9, 3.5814952734253125`*^9}, {3.581495422466776*^9, 3.581495493897702*^9}, { 3.5814955450767927`*^9, 3.581495592931077*^9}, {3.5814959122566423`*^9, 3.581495914783847*^9}, 3.5814959583627234`*^9, {3.581495990769981*^9, 3.581496042323472*^9}, {3.581496166860293*^9, 3.5814961681404953`*^9}, { 3.581496211220171*^9, 3.581496268267872*^9}, {3.581496314563154*^9, 3.5814963205389643`*^9}, {3.5814963777708654`*^9, 3.581496378831668*^9}, { 3.5814966241231017`*^9, 3.581496841866887*^9}, {3.5814968760847483`*^9, 3.58149695638349*^9}, {3.5814969870435443`*^9, 3.581497021245805*^9}, { 3.581497070472892*^9, 3.5814971273021927`*^9}, {3.5814980924357*^9, 3.5814981186331463`*^9}, {3.581498195994083*^9, 3.5814981984598875`*^9}, { 3.5814988817003007`*^9, 3.5814989105195513`*^9}, {3.581499086589263*^9, 3.5814990933450747`*^9}, {3.58149912442633*^9, 3.5814991691208096`*^9}, { 3.5814992165694933`*^9, 3.5814992169760942`*^9}, {3.5814992662967815`*^9, 3.5814992682945843`*^9}, {3.581499308815856*^9, 3.581499316070869*^9}, { 3.5814993868243947`*^9, 3.5814993873079953`*^9}, {3.5814994617809267`*^9, 3.581499468567939*^9}, {3.5815106299000654`*^9, 3.5815106391372814`*^9}, { 3.5815110135057435`*^9, 3.5815110286759706`*^9}, {3.5815110831466665`*^9, 3.581511086064872*^9}, {3.581511151193387*^9, 3.5815111676700163`*^9}, { 3.5815112910278344`*^9, 3.5815114290888786`*^9}, 3.5815116061363916`*^9, 3.58151167332311*^9, {3.581511950196599*^9, 3.5815119640680237`*^9}, { 3.581512028259537*^9, 3.5815120550819845`*^9}, {3.58151212555471*^9, 3.5815121260539107`*^9}, {3.581512233530901*^9, 3.5815122389445114`*^9}, { 3.5815122978788157`*^9, 3.5815123561403184`*^9}, {3.5815537107726717`*^9, 3.5815537805705957`*^9}, {3.58155607793946*^9, 3.5815561234249415`*^9}, { 3.581556166614818*^9, 3.5815562609625854`*^9}, {3.5815566302322392`*^9, 3.5815566438384633`*^9}, {3.5815567567266626`*^9, 3.581556867089059*^9}, { 3.5815569383317842`*^9, 3.58155705465699*^9}, {3.5815572419219217`*^9, 3.5815572745671797`*^9}, {3.581557449847497*^9, 3.5815574603639154`*^9}, { 3.5815577815810843`*^9, 3.5815578343971777`*^9}, 3.581559273295124*^9, 3.5815594090433645`*^9, {3.581600453911972*^9, 3.5816004860006294`*^9}, { 3.5816574214652042`*^9, 3.581657437782833*^9}, {3.5816575291053934`*^9, 3.581657529744994*^9}, {3.581657718552126*^9, 3.5816577278497424`*^9}, { 3.583410946631775*^9, 3.583410952435985*^9}, {3.583470336209275*^9, 3.583470341796085*^9}, {3.583470594135929*^9, 3.5834705990811377`*^9}, { 3.5835652299728966`*^9, 3.583565241205916*^9}, {3.5842569059185286`*^9, 3.5842569160137463`*^9}, {3.5842575753765182`*^9, 3.58425757649972*^9}, { 3.7236229836409836`*^9, 3.7236229889170213`*^9}, 3.7236236940018296`*^9, { 3.7273412111531963`*^9, 3.727341235434759*^9}, 3.727341338991079*^9, { 3.727341580125145*^9, 3.7273416634025517`*^9}, {3.727359040773569*^9, 3.7273590510829735`*^9}, {3.727359131615052*^9, 3.72735913643587*^9}, 3.72748997770763*^9, {3.7275097914288425`*^9, 3.727509938239212*^9}, { 3.7275099750829835`*^9, 3.727509980551524*^9}, {3.7275332980145874`*^9, 3.7275333823096633`*^9}, {3.727533518902872*^9, 3.7275335418589*^9}, { 3.7275335943985844`*^9, 3.7275336037063446`*^9}, {3.7275336378400474`*^9, 3.727533722237434*^9}, {3.7275717530834756`*^9, 3.7275717639994974`*^9}, { 3.727571801095117*^9, 3.727571851914236*^9}, {3.7275719577626963`*^9, 3.727572009238515*^9}, {3.7275720456172357`*^9, 3.727572072194888*^9}, 3.727572164679779*^9, {3.7275722241960554`*^9, 3.7275722290238547`*^9}, { 3.7275723340601096`*^9, 3.727572345158622*^9}, {3.7275724023477106`*^9, 3.7275724552152233`*^9}, {3.727572628689478*^9, 3.7275726538403463`*^9}, { 3.7275727822102537`*^9, 3.727572870095913*^9}, {3.7275729047800913`*^9, 3.727572912322177*^9}, {3.7275729471678257`*^9, 3.7275729522041764`*^9}, { 3.7275730366557007`*^9, 3.7275730391092243`*^9}, 3.727573126793394*^9, { 3.7275734643488526`*^9, 3.7275734677588725`*^9}, {3.7275736102106075`*^9, 3.7275736160721927`*^9}, {3.727573711003586*^9, 3.727573746640335*^9}, { 3.727626606947974*^9, 3.727626622702861*^9}, {3.7276266601172953`*^9, 3.7276267031346703`*^9}, {3.727626975135805*^9, 3.7276270043915873`*^9}, { 3.727627151090649*^9, 3.7276271931735325`*^9}, {3.7276273591076612`*^9, 3.727627375422084*^9}, {3.7277458809910774`*^9, 3.7277458815692863`*^9}, { 3.727838582935392*^9, 3.727838594249132*^9}, {3.728573264024838*^9, 3.7285732984111276`*^9}, {3.7285738727197027`*^9, 3.7285738731417522`*^9}, {3.7285758808493805`*^9, 3.72857588165664*^9}, { 3.7286101443550544`*^9, 3.7286101692380047`*^9}, {3.72861276013577*^9, 3.7286127606499667`*^9}, {3.7286233690913043`*^9, 3.7286233699757605`*^9}, {3.7286657242510815`*^9, 3.728665764201351*^9}, 3.728672548312504*^9, {3.7289675285739627`*^9, 3.728967529597266*^9}, { 3.72896762505945*^9, 3.7289676508858595`*^9}, {3.7294358423890896`*^9, 3.7294358739079885`*^9}, {3.7294359206005864`*^9, 3.729435968386664*^9}, { 3.7294360320027275`*^9, 3.729436110167185*^9}, {3.729436183034219*^9, 3.7294361841593723`*^9}, {3.729436289669063*^9, 3.729436421370212*^9}, { 3.7294365105203924`*^9, 3.7294365493992896`*^9}, 3.729436608671262*^9, { 3.7294366904141045`*^9, 3.729436716260569*^9}, 3.7294374130048695`*^9, { 3.7295244532527256`*^9, 3.729524453510421*^9}, {3.7295244864316254`*^9, 3.7295245899464273`*^9}, {3.729524807114545*^9, 3.729524828444416*^9}, 3.72952486414518*^9, {3.729524899424323*^9, 3.7295249523886566`*^9}, { 3.729524987244117*^9, 3.729525094451185*^9}, 3.729525124873105*^9, { 3.7295308194951906`*^9, 3.729530851952281*^9}, {3.729557515175043*^9, 3.7295575157980776`*^9}, {3.72955977085687*^9, 3.7295597885300894`*^9}, { 3.729559828856864*^9, 3.729559854493911*^9}, {3.7298314522941947`*^9, 3.7298314947991133`*^9}, {3.7298320454051285`*^9, 3.729832086537828*^9}, { 3.729832558441718*^9, 3.7298325741164846`*^9}, {3.7394081256351185`*^9, 3.7394081514866257`*^9}, {3.739408337214862*^9, 3.739408353260339*^9}, 3.7404514115410013`*^9, 3.740451503140458*^9, {3.7404515384793625`*^9, 3.7404515390602827`*^9}, {3.7411358549186172`*^9, 3.741135876639139*^9}, 3.7411432916227493`*^9, {3.7411434180505857`*^9, 3.7411434534721136`*^9}, { 3.7411829487432814`*^9, 3.741182952354971*^9}, {3.7411830219634895`*^9, 3.741183051286788*^9}, {3.741183636969272*^9, 3.7411836481639905`*^9}, { 3.7411840747634277`*^9, 3.7411840858538322`*^9}, {3.74118415120074*^9, 3.741184156452607*^9}, {3.741184260727316*^9, 3.7411842652322807`*^9}, { 3.741184928164649*^9, 3.741184928586512*^9}, {3.7411850560332756`*^9, 3.7411850688625565`*^9}, 3.7411851501317053`*^9, {3.7411855593217106`*^9, 3.7411856847092953`*^9}, {3.741185722839398*^9, 3.7411857247471876`*^9}, { 3.7411857747193823`*^9, 3.741185857521716*^9}, {3.7411859324297123`*^9, 3.741185948579706*^9}, {3.7411860851863995`*^9, 3.7411861213201933`*^9}, { 3.7411863263649144`*^9, 3.7411863564406805`*^9}, {3.7412447802056575`*^9, 3.7412448072148743`*^9}, {3.741332297786372*^9, 3.741332303019165*^9}, { 3.741332566951375*^9, 3.741332567422721*^9}, {3.7413326009773645`*^9, 3.74133263160271*^9}, {3.7413328414357605`*^9, 3.741332865196314*^9}, { 3.741333952177419*^9, 3.7413339574539027`*^9}, {3.7413341290171137`*^9, 3.7413341296023893`*^9}, {3.7760494893864474`*^9, 3.776049492432605*^9}, { 3.822380484700011*^9, 3.8223805515550766`*^9}, {3.8223805861001577`*^9, 3.8223806085775275`*^9}, {3.822381684587931*^9, 3.8223816881667137`*^9}, { 3.8223821919468665`*^9, 3.82238219951668*^9}, {3.8223827353021526`*^9, 3.822382744425187*^9}, {3.8223831213845315`*^9, 3.8223831615909634`*^9}, { 3.8223975738730855`*^9, 3.8223975843639317`*^9}, {3.8223976259339476`*^9, 3.822397652322728*^9}, 3.82244894682337*^9, 3.822457736760254*^9, 3.822458856384435*^9, 3.822458891387639*^9, {3.8224887632176266`*^9, 3.8224888314405932`*^9}},ExpressionUUID->"96ebfb21-588b-4f59-bcd5-\ b127f17db049"] }, Open ]] }, WindowSize->{1002, 779}, WindowMargins->{{-1019, Automatic}, {Automatic, 196}}, TaggingRules->{"SlideshowSettings" -> {"Toolbar" -> True}}, FrontEndVersion->"12.0 for Microsoft Windows (64-bit) (April 8, 2019)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[580, 22, 507, 13, 44, "Subsubsection",ExpressionUUID->"e66d6e93-d6b2-4c5b-9154-8959b5b959a5"], Cell[1090, 37, 6121, 166, 122, "Text",ExpressionUUID->"a16ffa28-7742-4592-bcc2-c1a697822085"], Cell[7214, 205, 42159, 937, 1511, "Input",ExpressionUUID->"96ebfb21-588b-4f59-bcd5-b127f17db049"] }, Open ]] } ] *)