(* 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[ 10064, 341] NotebookOptionsPosition[ 9050, 303] NotebookOutlinePosition[ 9413, 319] CellTagsIndexPosition[ 9370, 316] WindowFrame->Normal ContainsDynamic->False*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["Schwarz CLP Surfaces", "Title", CellChangeTimes->{3.391456687701777*^9}], Cell[BoxData[ RowBox[{"<<", "Own`Mesh`"}]], "Input"], Cell[CellGroupData[{ Cell["Weierstrass Representation", "Section", CellChangeTimes->{{3.391456261539761*^9, 3.391456268535554*^9}}], Cell["Choose lattice parameter:", "Text", CellChangeTimes->{{3.3914561868120937`*^9, 3.391456205753146*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"\[Tau]", "=", RowBox[{"1.", "I"}]}], ";"}]], "Input", CellChangeTimes->{{3.391456209363646*^9, 3.391456210568697*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"a", "=", ".25"}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"G", "[", "z_", "]"}], ":=", RowBox[{ RowBox[{"Sqrt", "[", "I", "]"}], RowBox[{ RowBox[{"\[Theta]11", "[", RowBox[{ RowBox[{"z", "-", "a"}], ",", "\[Tau]"}], "]"}], "^", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "/", "2."}], ")"}]}], RowBox[{ RowBox[{"\[Theta]11", "[", RowBox[{ RowBox[{"z", "+", "a"}], ",", "\[Tau]"}], "]"}], "^", RowBox[{"(", RowBox[{"1", "/", "2."}], ")"}]}], RowBox[{ RowBox[{"\[Theta]11", "[", RowBox[{ RowBox[{"z", "-", "a", "+", RowBox[{"\[Tau]", "/", "2"}]}], ",", "\[Tau]"}], "]"}], "^", RowBox[{"(", RowBox[{ RowBox[{"+", "1"}], "/", "2."}], ")"}]}], RowBox[{ RowBox[{"\[Theta]11", "[", RowBox[{ RowBox[{"z", "+", "a", "+", RowBox[{"\[Tau]", "/", "2"}]}], ",", "\[Tau]"}], "]"}], "^", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "/", "2."}], ")"}]}]}]}]], "Input", CellChangeTimes->{3.391456165541319*^9}], Cell[BoxData[ RowBox[{ RowBox[{"dh", "[", "z_", "]"}], ":=", "1"}]], "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", CellChangeTimes->{{3.3914561672344847`*^9, 3.391456174938344*^9}}], 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[{ RowBox[{"\[Omega]3", "[", "z_", "]"}], ":=", RowBox[{"dh", "[", "z", "]"}]}], ";"}]}], "Input", CellChangeTimes->{{3.391456161295723*^9, 3.391456184928474*^9}, { 3.391456227385852*^9, 3.391456278944255*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"\[Epsilon]", "=", RowBox[{"10", "^", RowBox[{"-", "12."}]}]}], ";"}], "\n", RowBox[{ RowBox[{"cp", "=", RowBox[{".25", "+", RowBox[{"\[Tau]", "/", "4."}]}]}], ";"}]}], "Input", CellChangeTimes->{{3.391456161295723*^9, 3.391456184928474*^9}, { 3.391456227385852*^9, 3.391456281745907*^9}}], 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"}], "}"}], "]"}]}], "\n", RowBox[{ RowBox[{ RowBox[{"\[Delta]0", "=", RowBox[{"f0", "[", "0", "]"}]}], ";"}], "\n"}], "\n", 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"}]}], "\n", RowBox[{ RowBox[{ RowBox[{"n", "[", "w_", "]"}], ":=", RowBox[{"StereographicProjection", "[", RowBox[{"G", "[", "w", "]"}], "]"}]}], ";"}]}], "Input", CellChangeTimes->{{3.391456161295723*^9, 3.391456184928474*^9}, { 3.391456227385852*^9, 3.391456256205233*^9}, 3.391456294144058*^9}] }, Open ]], Cell[CellGroupData[{ Cell["Parametrization", "Section", CellChangeTimes->{{3.391456297404413*^9, 3.391456301329183*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"y1", "=", RowBox[{ RowBox[{"Im", "[", "\[Tau]", "]"}], "/", "2."}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"nx", "=", "4"}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"ny", "=", "8"}], ";"}], "\[IndentingNewLine]"}], "\n", RowBox[{ RowBox[{"domain", "=", RowBox[{"RectangularDomain", "[", RowBox[{ RowBox[{"Union", "[", RowBox[{ RowBox[{"NRange", "[", RowBox[{"0", ",", "a", ",", "nx"}], "]"}], ",", RowBox[{"NRange", "[", RowBox[{"a", ",", ".5", ",", "nx"}], "]"}], ",", RowBox[{"SameTest", "\[Rule]", "AlmostEqual"}]}], "]"}], ",", RowBox[{"NRange", "[", RowBox[{"\[Epsilon]", ",", "y1", ",", "ny"}], "]"}]}], "]"}]}], ";"}]}], "Input", CellChangeTimes->{{3.391456161295723*^9, 3.391456184928474*^9}, { 3.391456227385852*^9, 3.391456256205233*^9}, {3.391456294144058*^9, 3.391456317499915*^9}, 3.3914563578898497`*^9, 3.391456573987247*^9, { 3.391456651075223*^9, 3.391456652887061*^9}}] }, Open ]], Cell[CellGroupData[{ Cell["Surface Generation", "Section", CellChangeTimes->{{3.3914563940393047`*^9, 3.3914563986326647`*^9}}], 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"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"fr2", "=", RowBox[{"MeshJoin", "[", RowBox[{"fr1", ",", RowBox[{"MeshReflect", "[", RowBox[{"fr1", ",", RowBox[{"Plane", "[", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", RowBox[{ RowBox[{"f", "[", RowBox[{".5", "+", RowBox[{"\[Tau]", "/", "2"}]}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], "]"}]}], "]"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"fr3", "=", RowBox[{"MeshJoin", "[", RowBox[{"fr2", ",", RowBox[{"MeshReflect", "[", RowBox[{"fr2", ",", RowBox[{"Plane", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "1", ",", "0"}], "}"}], ",", "0"}], "]"}]}], "]"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"fr4", "=", RowBox[{"MeshJoin", "[", RowBox[{"fr3", ",", RowBox[{"MeshReflect", "[", RowBox[{"fr3", ",", RowBox[{"Plane", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", "1"}], "}"}], ",", "0"}], "]"}]}], "]"}]}], "]"}]}], ";"}]}], "Input", CellChangeTimes->{{3.391456414791316*^9, 3.391456421455256*^9}, { 3.391456597781166*^9, 3.391456612682633*^9}}], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{ RowBox[{"gr", "=", RowBox[{"Mesh3DToGraphics3D", "[", "fr4", "]"}]}], ",", RowBox[{"Axes", "\[Rule]", "True"}], ",", RowBox[{"AxesLabel", "\[Rule]", RowBox[{"{", RowBox[{"x", ",", "y", ",", "z"}], "}"}]}]}], "]"}]}]], "Input", CellChangeTimes->{{3.391456414791316*^9, 3.391456421455256*^9}, { 3.391456597781166*^9, 3.391456612682633*^9}}] }, Open ]] }, Open ]] }, WindowSize->{1005, 946}, WindowMargins->{{89, Automatic}, {Automatic, 78}}, ShowSelection->True, 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, 79, 1, 76, "Title"], Cell[672, 26, 53, 1, 27, "Input"], Cell[CellGroupData[{ Cell[750, 31, 111, 1, 67, "Section"], Cell[864, 34, 109, 1, 26, "Text"], Cell[976, 37, 158, 4, 27, "Input"], Cell[1137, 43, 69, 2, 27, "Input"], Cell[1209, 47, 1030, 34, 27, "Input"], Cell[2242, 83, 80, 2, 27, "Input"], Cell[2325, 87, 452, 14, 43, "Input"], Cell[2780, 103, 778, 25, 58, "Input"], Cell[3561, 130, 349, 10, 43, "Input"], Cell[3913, 142, 1545, 44, 148, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[5495, 191, 100, 1, 67, "Section"], Cell[5598, 194, 1104, 30, 103, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[6739, 229, 107, 1, 67, "Section"], Cell[6849, 232, 1713, 54, 73, "Input"], Cell[8565, 288, 457, 11, 43, "Input"] }, Open ]] }, Open ]] } ] *) (* End of internal cache information *)