(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 6.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 14042, 475] NotebookOptionsPosition[ 12442, 420] NotebookOutlinePosition[ 12785, 435] CellTagsIndexPosition[ 12742, 432] WindowFrame->Normal ContainsDynamic->False*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["H -Surface of H.A. Schwarz", "Title", CellChangeTimes->{{3.3916904944746227`*^9, 3.39169049910471*^9}, { 3.391690786906435*^9, 3.391690793018215*^9}}], Cell[BoxData[ RowBox[{"<<", "Own`Mesh`"}]], "Input"], Cell[CellGroupData[{ Cell["Weierstrass Data", "Section", CellChangeTimes->{{3.391690511500511*^9, 3.3916905157767572`*^9}}], Cell["Choose lattice parameter", "Text"], Cell[BoxData[ RowBox[{ RowBox[{"\[Tau]", "=", RowBox[{"2.", "I"}]}], ";"}]], "Input", CellChangeTimes->{{3.391690523525725*^9, 3.391690524825197*^9}, 3.3918803918852177`*^9, {3.391880428221692*^9, 3.391880428467313*^9}, 3.391880552160636*^9}], Cell[BoxData[{ RowBox[{ RowBox[{"dh", "[", "z_", "]"}], ":=", "1"}], "\n", RowBox[{ RowBox[{"\[Rho]", "=", "1."}], ";"}], "\n", RowBox[{ RowBox[{"a0", "=", RowBox[{"1", "/", "4."}]}], ";"}], "\n", RowBox[{ RowBox[{"G2", "[", "z_", "]"}], ":=", RowBox[{"\[Rho]", " ", RowBox[{ RowBox[{"\[Theta]11", "[", RowBox[{ RowBox[{"z", "-", "a0"}], ",", "\[Tau]"}], "]"}], "/", RowBox[{"\[Theta]11", "[", RowBox[{ RowBox[{"z", "+", "a0"}], ",", "\[Tau]"}], "]"}]}]}]}]}], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"G", "[", "z_", "]"}], ":=", RowBox[{ RowBox[{"G2", "[", "z", "]"}], "^", RowBox[{"(", RowBox[{"2", "/", "3"}], ")"}]}]}]], "Input"], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"phi1", "[", "z_", "]"}], ":=", RowBox[{ RowBox[{"G", "[", "z", " ", "]"}], " ", RowBox[{"dh", "[", "z", "]"}]}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"phi2", "[", "z_", "]"}], ":=", RowBox[{ RowBox[{"1", "/", RowBox[{"G", "[", "z", "]"}]}], RowBox[{"dh", "[", "z", "]"}]}]}], ";"}]}], "Input"], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"\[Omega]1", "[", "z_", "]"}], " ", ":=", RowBox[{ RowBox[{"(", " ", RowBox[{ RowBox[{"phi2", "[", "z", "]"}], "-", RowBox[{"phi1", "[", "z", "]"}]}], " ", ")"}], "/", "2"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"\[Omega]2", "[", "z_", "]"}], " ", ":=", RowBox[{"I", RowBox[{ RowBox[{"(", " ", RowBox[{ RowBox[{"phi2", "[", "z", "]"}], "+", RowBox[{"phi1", "[", "z", "]"}]}], ")"}], "/", "2"}]}]}], ";"}], "\n", RowBox[{ RowBox[{"\[Omega]3", "[", "z_", "]"}], ":=", RowBox[{"dh", "[", "z", "]"}]}]}], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"cp", "=", RowBox[{".25", "+", RowBox[{"\[Tau]", "/", "4."}]}]}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"f0", "[", "w_", "]"}], ":=", RowBox[{"Re", "[", RowBox[{"{", "\n", "\t\t", RowBox[{ RowBox[{"NIntegrate", "[", RowBox[{ RowBox[{"\[Omega]1", "[", "z", "]"}], ",", RowBox[{"{", RowBox[{"z", ",", "cp", ",", "w"}], "}"}]}], "]"}], ",", RowBox[{"NIntegrate", "[", RowBox[{ RowBox[{"\[Omega]2", "[", "z", "]"}], ",", RowBox[{"{", RowBox[{"z", ",", "cp", ",", "w"}], "}"}]}], "]"}], ",", "\n", "\t", "w"}], "}"}], "]"}]}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"\[Delta]0", "=", RowBox[{"f0", "[", ".5", "]"}]}]], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", "0.2969972241227937`"}], ",", "0.38107394043563547`", ",", "0.5`"}], "}"}]], "Output", CellChangeTimes->{ 3.3895268773870564`*^9, {3.3916905544767303`*^9, 3.391690559065477*^9}, 3.391880395786188*^9, 3.3918804300517893`*^9, {3.391880498550592*^9, 3.391880503619698*^9}, {3.391880534024354*^9, 3.391880554662176*^9}}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"f", "[", "w_", "]"}], ":=", RowBox[{ RowBox[{"Re", "[", RowBox[{"{", "\n", "\t\t", RowBox[{ RowBox[{"NIntegrate", "[", RowBox[{ RowBox[{"\[Omega]1", "[", "z", "]"}], ",", RowBox[{"{", RowBox[{"z", ",", "cp", ",", "w"}], "}"}]}], "]"}], ",", RowBox[{"NIntegrate", "[", RowBox[{ RowBox[{"\[Omega]2", "[", "z", "]"}], ",", RowBox[{"{", RowBox[{"z", ",", "cp", ",", "w"}], "}"}]}], "]"}], ",", "\n", "\t", "w"}], "}"}], "]"}], "-", "\[Delta]0"}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"n", "[", "w_", "]"}], ":=", RowBox[{"StereographicProjection", "[", RowBox[{"G", "[", "w", "]"}], "]"}]}]], "Input"] }, Open ]], Cell[CellGroupData[{ Cell["Parametrization", "Section", CellChangeTimes->{{3.39169075054737*^9, 3.3916907545386143`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"\[Epsilon]", "=", RowBox[{"10", "^", RowBox[{"-", "7"}]}]}], ";"}]], "Input", CellChangeTimes->{{3.38952686934953*^9, 3.3895268711376057`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"y1", "=", RowBox[{ RowBox[{"Im", "[", "\[Tau]", "]"}], "/", "2."}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"nx", "=", "4"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"ny", "=", "5"}], ";"}]}], "Input", CellChangeTimes->{ 3.389527132523203*^9, {3.389528569474218*^9, 3.389528571253669*^9}, 3.391690570830476*^9, {3.391690766083724*^9, 3.391690772555951*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"domain", "=", RowBox[{"RectangularDomain", "[", RowBox[{ RowBox[{"Union", "[", RowBox[{ RowBox[{"NRange", "[", RowBox[{"0", ",", "a0", ",", "nx"}], "]"}], ",", RowBox[{"NRange", "[", RowBox[{"a0", ",", ".5", ",", "nx"}], "]"}], ",", RowBox[{"SameTest", "\[Rule]", "AlmostEqual"}]}], "]"}], ",", RowBox[{"NRange", "[", RowBox[{"0", ",", "y1", ",", "ny"}], "]"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.389526888524867*^9, 3.389526928192543*^9}, { 3.391690573867118*^9, 3.391690575779105*^9}}] }, Open ]], Cell[CellGroupData[{ Cell["Surface Patch Creation", "Section"], Cell[BoxData[ RowBox[{ RowBox[{"fr1", "=", RowBox[{"MeshPlot3D", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"f", "[", RowBox[{"x", "+", RowBox[{"I", " ", "y"}]}], "]"}], ",", RowBox[{"n", "[", RowBox[{"x", "+", RowBox[{"I", " ", "y"}]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"x", ",", "y"}], "}"}], ",", "domain"}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.389526950924272*^9, 3.389526954584331*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"c0", "=", RowBox[{ RowBox[{"f", "[", RowBox[{"\[Tau]", "/", "2"}], "]"}], ".", RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", RowBox[{"30.", " ", "Degree"}], "]"}], ",", RowBox[{"Sin", "[", RowBox[{ RowBox[{"-", "30"}], " ", "Degree"}], "]"}], ",", "0"}], "}"}]}]}], ";"}]], "Input", CellChangeTimes->{3.3916906406626177`*^9}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"len", "=", RowBox[{"2", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"f", "[", RowBox[{"\[Tau]", "/", "2"}], "]"}], "-", RowBox[{"f", "[", "a0", "]"}]}], ")"}], ".", RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", RowBox[{"30.", " ", "Degree"}], "]"}], ",", RowBox[{"Sin", "[", RowBox[{ RowBox[{"-", "30"}], " ", "Degree"}], "]"}], ",", "0"}], "}"}]}]}]}]], "Input"], Cell[BoxData["2.524566353920675`"], "Output", CellChangeTimes->{ 3.3895272307259808`*^9, 3.389528580572431*^9, 3.391690653489088*^9, { 3.391880402970707*^9, 3.391880431553011*^9}, {3.391880535513935*^9, 3.3918805563305283`*^9}}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"fr2", "=", RowBox[{"MeshJoin", "[", RowBox[{"fr1", ",", RowBox[{"MeshReflect", "[", RowBox[{"fr1", ",", RowBox[{"Plane", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "1", ",", "0"}], "}"}], ",", "0"}], "]"}]}], "]"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.389526970329463*^9, 3.389526982552957*^9}, 3.389527157452401*^9, 3.391690644702442*^9}], Cell[BoxData[ RowBox[{ RowBox[{"fr3", "=", RowBox[{"MeshJoin", "[", RowBox[{"fr2", ",", RowBox[{"MeshReflect", "[", RowBox[{"fr2", ",", RowBox[{"Plane", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", "1"}], "}"}], ",", RowBox[{"2", RowBox[{ RowBox[{"f", "[", "a0", "]"}], "[", RowBox[{"[", "3", "]"}], "]"}]}]}], "]"}]}], "]"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.3895269925910473`*^9, 3.389527010530532*^9}, 3.389527165332665*^9, 3.391690645652646*^9}], Cell[BoxData[ RowBox[{"\n", RowBox[{ RowBox[{ RowBox[{"fr4", "=", RowBox[{"MeshJoin", "[", RowBox[{"fr3", ",", RowBox[{"MeshReflect", "[", RowBox[{"fr3", ",", RowBox[{"Plane", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", RowBox[{"30.", " ", "Degree"}], "]"}], ",", RowBox[{"Sin", "[", RowBox[{"30", " ", "Degree"}], "]"}], ",", "0"}], "}"}], ",", "c0"}], "]"}]}], "]"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", "\n", RowBox[{ RowBox[{"fr5", "=", RowBox[{"MeshJoin", "[", RowBox[{"fr4", ",", RowBox[{"MeshReflect", "[", RowBox[{"fr4", ",", RowBox[{"Plane", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", RowBox[{"30", " ", "Degree"}], "]"}], ",", RowBox[{"Sin", "[", RowBox[{ RowBox[{"-", "30"}], " ", "Degree"}], "]"}], ",", "0"}], "}"}], ",", "c0"}], "]"}]}], "]"}]}], "]"}]}], ";"}]}]}]], "Input", CellChangeTimes->{{3.3895270303044*^9, 3.389527085580772*^9}, { 3.39169063642782*^9, 3.391690651776967*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"fr6", "=", RowBox[{"MeshJoin", "[", RowBox[{"fr5", ",", RowBox[{"MeshTranslate", "[", RowBox[{"fr5", ",", RowBox[{"len", RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", RowBox[{"30.", " ", "Degree"}], "]"}], ",", RowBox[{"Sin", "[", RowBox[{"30", " ", "Degree"}], "]"}], ",", "0"}], "}"}]}]}], "]"}], ",", RowBox[{"MeshTranslate", "[", RowBox[{"fr5", ",", RowBox[{"len", RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", RowBox[{"30.", " ", "Degree"}], "]"}], ",", RowBox[{"Sin", "[", RowBox[{ RowBox[{"-", "30"}], " ", "Degree"}], "]"}], ",", "0"}], "}"}]}]}], "]"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.3895272083120527`*^9, 3.389527220239347*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"fr7", "=", RowBox[{"MeshJoin", "[", RowBox[{"fr6", ",", RowBox[{"MeshTranslate", "[", RowBox[{"fr6", ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", "1"}], "}"}]}], "]"}], ",", RowBox[{"MeshTranslate", "[", RowBox[{"fr6", ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", RowBox[{"-", "1"}]}], "}"}]}], "]"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.389527271468606*^9, 3.389527282326833*^9}}], Cell[BoxData[ RowBox[{"Show", "[", RowBox[{ RowBox[{"gr", "=", RowBox[{"Mesh3DToGraphics3D", "[", "fr7", "]"}]}], ",", RowBox[{"Axes", "\[Rule]", "True"}]}], "]"}]], "Input", CellChangeTimes->{{3.3895271008517227`*^9, 3.389527103464707*^9}, 3.3895272381092033`*^9, {3.389527298934924*^9, 3.389527325891736*^9}, { 3.391690680273333*^9, 3.391690696818735*^9}, {3.391880448614911*^9, 3.391880475769706*^9}, 3.391880540952873*^9}], Cell[BoxData[""], "Input", CellChangeTimes->{{3.391880408037641*^9, 3.391880476127285*^9}, { 3.391880535972831*^9, 3.3918805721045427`*^9}}] }, Open ]] }, Open ]] }, WindowSize->{1005, 946}, WindowMargins->{{106, Automatic}, {Automatic, 12}}, FrontEndVersion->"6.0 for Mac OS X PowerPC (32-bit) (April 20, 2007)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[590, 23, 159, 2, 76, "Title"], Cell[752, 27, 53, 1, 27, "Input"], Cell[CellGroupData[{ Cell[830, 32, 103, 1, 67, "Section"], Cell[936, 35, 40, 0, 26, "Text"], Cell[979, 37, 258, 6, 27, "Input"], Cell[1240, 45, 528, 17, 73, "Input"], Cell[1771, 64, 181, 6, 27, "Input"], Cell[1955, 72, 384, 13, 43, "Input"], Cell[2342, 87, 643, 22, 58, "Input"], Cell[2988, 111, 125, 4, 27, "Input"], Cell[3116, 117, 552, 16, 58, "Input"], Cell[CellGroupData[{ Cell[3693, 137, 87, 2, 27, "Input"], Cell[3783, 141, 387, 8, 27, "Output"] }, Open ]], Cell[4185, 152, 596, 17, 58, "Input"], Cell[4784, 171, 158, 4, 27, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[4979, 180, 101, 1, 67, "Section"], Cell[5083, 183, 188, 5, 27, "Input"], Cell[5274, 190, 435, 12, 58, "Input"], Cell[5712, 204, 614, 15, 27, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[6363, 224, 41, 0, 67, "Section"], Cell[6407, 226, 506, 16, 27, "Input"], Cell[6916, 244, 422, 14, 27, "Input"], Cell[CellGroupData[{ Cell[7363, 262, 470, 16, 27, "Input"], Cell[7836, 280, 238, 4, 27, "Output"] }, Open ]], Cell[8089, 287, 460, 13, 27, "Input"], Cell[8552, 302, 589, 17, 27, "Input"], Cell[9144, 321, 1258, 36, 73, "Input"], Cell[10405, 359, 892, 26, 43, "Input"], Cell[11300, 387, 511, 14, 27, "Input"], Cell[11814, 403, 454, 9, 27, "Input"], Cell[12271, 414, 143, 2, 27, "Input"] }, Open ]] }, Open ]] } ] *) (* End of internal cache information *)