(* 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[ 11953, 420] NotebookOptionsPosition[ 10504, 370] NotebookOutlinePosition[ 10846, 385] CellTagsIndexPosition[ 10803, 382] WindowFrame->Normal ContainsDynamic->False*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["D-Surface of H.A. Schwarz", "Section", CellChangeTimes->{{3.391694560065063*^9, 3.391694571742551*^9}}], Cell[BoxData[ RowBox[{"<<", "Own`Mesh`"}]], "Input"] }, Open ]], Cell[CellGroupData[{ Cell["Weierstrass Data", "Section", CellChangeTimes->{{3.3916937632609367`*^9, 3.391693767898058*^9}}], Cell["\<\ Choose lattice parameter. The constant below gives the most symmetric surface\ \ \>", "Text", CellChangeTimes->{{3.3916945758267736`*^9, 3.3916946032869062`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"\[Tau]", "=", RowBox[{"0.3908504810515956`", "I"}]}], ";"}]], "Input"], 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[{"1", "/", "2"}], ")"}]}]}]], "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"] }, Open ]], Cell[CellGroupData[{ Cell["Parametrization", "Section", CellChangeTimes->{{3.39169389742931*^9, 3.391693899971046*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"\[Epsilon]", "=", RowBox[{"10", "^", RowBox[{"-", "6"}]}]}], ";"}]], "Input", CellChangeTimes->{{3.391693885571022*^9, 3.391693890979126*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"cp", "=", RowBox[{".25", "+", RowBox[{"\[Tau]", "/", "4."}]}]}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"f0", "[", "w_", "]"}], ":=", RowBox[{"Im", "[", 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.14177760196511296`"}], ",", "0.2490728786008459`", ",", "0"}], "}"}]], "Output", CellChangeTimes->{3.391693908638855*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"f", "[", "w_", "]"}], ":=", RowBox[{ RowBox[{"Im", "[", 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"], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"y1", "=", RowBox[{ RowBox[{"Im", "[", "\[Tau]", "]"}], "/", "2."}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"nx", "=", "4"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"ny", "=", "5"}], ";"}]}]}]], "Input", CellChangeTimes->{ 3.391693915669951*^9, {3.391693988102202*^9, 3.3916939903930407`*^9}, { 3.391694478008788*^9, 3.391694480114407*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"domain", "=", RowBox[{"RectangularDomain", "[", RowBox[{ RowBox[{"Union", "[", RowBox[{ RowBox[{"NRange", "[", RowBox[{"0", ",", "a0", ",", "nx"}], "]"}], ",", RowBox[{"NRange", "[", RowBox[{"a0", ",", ".5", ",", RowBox[{"2", "nx"}]}], "]"}], ",", RowBox[{"SameTest", "\[Rule]", "AlmostEqual"}]}], "]"}], ",", RowBox[{"NRange", "[", RowBox[{"0", ",", "y1", ",", "ny"}], "]"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.3916939187893553`*^9, 3.3916939747827873`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell["Surface Patch Creation", "Section"], Cell[BoxData[{ RowBox[{ RowBox[{"yv", "=", RowBox[{ RowBox[{"f", "[", "1", "]"}], "[", RowBox[{"[", "2", "]"}], "]"}]}], ";"}], "\n", 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.3916940999372177`*^9, 3.391694107864546*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"fr2", "=", RowBox[{"MeshJoin", "[", RowBox[{"fr1", ",", RowBox[{"MeshRotate", "[", RowBox[{"fr1", ",", RowBox[{"StraightLine", "[", RowBox[{ RowBox[{"{", RowBox[{"1", ",", RowBox[{"-", "1"}], ",", "y1"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "1", ",", "y1"}], "}"}]}], "]"}]}], "]"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.391694113198019*^9, 3.3916941594460373`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"fr3", "=", RowBox[{"MeshJoin", "[", RowBox[{"fr2", ",", RowBox[{"MeshRotate", "[", RowBox[{"fr2", ",", RowBox[{"StraightLine", "[", RowBox[{ RowBox[{"{", RowBox[{"yv", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"yv", ",", "0", ",", "1"}], "}"}]}], "]"}]}], "]"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.391694130614297*^9, 3.391694184754208*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"fr4", "=", RowBox[{"MeshJoin", "[", RowBox[{"fr3", ",", RowBox[{"MeshRotate", "[", RowBox[{"fr3", ",", RowBox[{"StraightLine", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"yv", ",", "0", ",", "0"}], "}"}]}], "]"}]}], "]"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.3916942680422373`*^9, 3.391694292099065*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"tr1", "=", RowBox[{"{", RowBox[{ RowBox[{"2", RowBox[{ RowBox[{"f", "[", "1", "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], ",", "0", ",", "0"}], "}"}]}], ";"}], "\n", RowBox[{ RowBox[{"fr5", "=", RowBox[{"MeshJoin", "[", RowBox[{"fr4", ",", RowBox[{"MeshTranslate", "[", RowBox[{"fr4", ",", "tr1"}], "]"}], ",", RowBox[{"MeshTranslate", "[", RowBox[{"fr4", ",", RowBox[{"-", "tr1"}]}], "]"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"tr2", "=", RowBox[{"{", RowBox[{"0", ",", RowBox[{"2", RowBox[{ RowBox[{"f", "[", "1", "]"}], "[", RowBox[{"[", "2", "]"}], "]"}]}], ",", "0"}], "}"}]}], ";"}], "\n", RowBox[{ RowBox[{"fr6", "=", RowBox[{"MeshJoin", "[", RowBox[{"fr5", ",", RowBox[{"MeshTranslate", "[", RowBox[{"fr5", ",", "tr2"}], "]"}], ",", RowBox[{"MeshTranslate", "[", RowBox[{"fr5", ",", RowBox[{"-", "tr2"}]}], "]"}]}], "]"}]}], ";"}]}], "Input", CellChangeTimes->{{3.391694366589098*^9, 3.391694406298986*^9}, 3.391694469919409*^9}], Cell[BoxData[ RowBox[{ RowBox[{"fr7", "=", RowBox[{"MeshJoin", "[", RowBox[{"fr6", ",", RowBox[{"MeshTranslate", "[", RowBox[{"fr6", ",", RowBox[{"4", RowBox[{"{", RowBox[{"0", ",", "0", ",", "y1"}], "}"}]}]}], "]"}], ",", RowBox[{"MeshTranslate", "[", RowBox[{"fr6", ",", RowBox[{ RowBox[{"-", "4"}], RowBox[{"{", RowBox[{"0", ",", "0", ",", "y1"}], "}"}]}]}], "]"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.3916944104918747`*^9, 3.391694418362138*^9}}], Cell[BoxData[ RowBox[{"Show", "[", RowBox[{ RowBox[{"gr", "=", RowBox[{"Mesh3DToGraphics3D", "[", "fr7", "]"}]}], ",", RowBox[{"Axes", "\[Rule]", "False"}], ",", RowBox[{"Boxed", "\[Rule]", "False"}]}], "]"}]], "Input", CellChangeTimes->{{3.391694431620379*^9, 3.391694440003141*^9}}], Cell[BoxData[""], "Input", CellChangeTimes->{{3.391694434122511*^9, 3.391694466251926*^9}}] }, Open ]] }, WindowSize->{1005, 946}, WindowMargins->{{80, Automatic}, {-6, Automatic}}, 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, 110, 1, 67, "Section"], Cell[703, 26, 53, 1, 27, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[793, 32, 103, 1, 67, "Section"], Cell[899, 35, 173, 4, 26, "Text"], Cell[1075, 41, 109, 3, 27, "Input"], Cell[1187, 46, 528, 17, 73, "Input"], Cell[1718, 65, 181, 6, 27, "Input"], Cell[1902, 73, 384, 13, 43, "Input"], Cell[2289, 88, 643, 22, 58, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[2969, 115, 99, 1, 67, "Section"], Cell[3071, 118, 187, 5, 27, "Input"], Cell[3261, 125, 125, 4, 27, "Input"], Cell[3389, 131, 552, 16, 58, "Input"], Cell[CellGroupData[{ Cell[3966, 151, 87, 2, 27, "Input"], Cell[4056, 155, 185, 5, 27, "Output"] }, Open ]], Cell[4256, 163, 596, 17, 58, "Input"], Cell[4855, 182, 158, 4, 27, "Input"], Cell[5016, 188, 478, 14, 73, "Input"], Cell[5497, 204, 594, 15, 27, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[6128, 224, 41, 0, 67, "Section"], Cell[6172, 226, 646, 21, 43, "Input"], Cell[6821, 249, 547, 16, 27, "Input"], Cell[7371, 267, 489, 14, 27, "Input"], Cell[7863, 283, 490, 14, 27, "Input"], Cell[8356, 299, 1159, 37, 73, "Input"], Cell[9518, 338, 567, 17, 27, "Input"], Cell[10088, 357, 305, 7, 27, "Input"], Cell[10396, 366, 92, 1, 27, "Input"] }, Open ]] } ] *) (* End of internal cache information *)