(* 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[      7962,        282]
NotebookOptionsPosition[      7121,        251]
NotebookOutlinePosition[      7464,        266]
CellTagsIndexPosition[      7421,        263]
WindowFrame->Normal
ContainsDynamic->False*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["The Singly Periodic Scherk Surfaces", "Title"],

Cell[BoxData[
 RowBox[{"<<", "Own`Mesh`"}]], "Input"],

Cell["\<\
We use a level curve parametrization for the singly periodic Scherk surface\
\>", "Text",
 CellChangeTimes->{{3.3881396932307367`*^9, 3.3881397063831778`*^9}}],

Cell[CellGroupData[{

Cell["Weierstrass Data", "Section",
 CellChangeTimes->{{3.388139683857562*^9, 3.38813968717456*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"dh", "[", "x_", "]"}], ":=", "1"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"G", "[", "x_", "]"}], ":=", 
   SqrtBox[
    RowBox[{"Cot", "[", "x", "]"}]]}], ";"}]}], "Input"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"\[Omega]1", "[", "z_", "]"}], ":=", " ", 
   RowBox[{
    RowBox[{"1", "/", "2"}], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"1", "/", 
       RowBox[{"G", "[", "z", "]"}]}], "-", 
      RowBox[{"G", "[", "z", "]"}]}], ")"}], 
    RowBox[{"dh", "[", "z", "]"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Omega]2", "[", "z_", "]"}], ":=", 
   RowBox[{
    RowBox[{"1", "/", "2"}], " ", "I", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"1", "/", 
       RowBox[{"G", "[", "z", "]"}]}], "+", 
      RowBox[{"G", "[", "z", "]"}]}], ")"}], 
    RowBox[{"dh", "[", "z", "]"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Omega]3", "[", "z_", "]"}], ":=", " ", 
   RowBox[{"dh", "[", "z", "]"}]}], ";"}]}], "Input"],

Cell["Explicit Integral", "Text",
 CellChangeTimes->{{3.388139722086543*^9, 3.38813972563868*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"f", "[", "z_", "]"}], ":=", 
   RowBox[{"Re", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{
       FractionBox["1", 
        RowBox[{"2", " ", 
         SqrtBox["2"]}]], 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Pi]"}], "+", 
         RowBox[{"Log", "[", 
          RowBox[{
           RowBox[{"-", "1"}], "+", 
           RowBox[{
            SqrtBox["2"], " ", 
            SqrtBox[
             RowBox[{"Cot", "[", "z", "]"}]]}], "-", 
           RowBox[{"Cot", "[", "z", "]"}]}], "]"}], "-", 
         RowBox[{"Log", "[", 
          RowBox[{"1", "+", 
           RowBox[{
            SqrtBox["2"], " ", 
            SqrtBox[
             RowBox[{"Cot", "[", "z", "]"}]]}], "+", 
           RowBox[{"Cot", "[", "z", "]"}]}], "]"}]}], ")"}]}], ",", 
      RowBox[{
       FractionBox["1", 
        SqrtBox["2"]], 
       RowBox[{"(", 
        RowBox[{"\[ImaginaryI]", " ", 
         RowBox[{"(", 
          RowBox[{"\[Pi]", "+", 
           RowBox[{"ArcTan", "[", 
            RowBox[{"1", "-", 
             RowBox[{
              SqrtBox["2"], " ", 
              SqrtBox[
               RowBox[{"Cot", "[", "z", "]"}]]}]}], "]"}], "-", 
           RowBox[{"ArcTan", "[", 
            RowBox[{"1", "+", 
             RowBox[{
              SqrtBox["2"], " ", 
              SqrtBox[
               RowBox[{"Cot", "[", "z", "]"}]]}]}], "]"}]}], ")"}]}], ")"}]}],
       ",", "z"}], "}"}], "]"}]}], ";"}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"n", "[", "z_", "]"}], ":=", 
  RowBox[{"StereographicProjection", "[", 
   RowBox[{"G", "[", "z", "]"}], "]"}]}]], "Input"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"y1", "=", "3"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"nr1", "=", "5"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"nt", "=", "5"}], ";"}]}], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"\[Epsilon]", "=", "0.000001"}], ";"}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"domain1", "=", 
   RowBox[{"RectangularDomain", "[", 
    RowBox[{
     RowBox[{"NRange", "[", 
      RowBox[{"\[Epsilon]", ",", 
       RowBox[{"Pi", "/", "2"}], ",", "nr1"}], "]"}], ",", 
     RowBox[{"NRange", "[", 
      RowBox[{"0", ",", "y1", ",", "nt"}], "]"}]}], "]"}]}], ";"}]], "Input"],

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"}], "}"}], ",", "domain1"}], "]"}]}], 
  ";"}]], "Input"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"fr2", "=", 
   RowBox[{"MeshReflect", "[", 
    RowBox[{"fr1", ",", 
     RowBox[{"Plane", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "1", ",", "0"}], "}"}], ",", "0"}], "]"}]}], 
    "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fr3", "=", 
   RowBox[{"MeshJoin", "[", 
    RowBox[{"fr1", ",", "fr2"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fr4", "=", 
   RowBox[{"MeshReflect", "[", 
    RowBox[{"fr3", ",", 
     RowBox[{"Plane", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", "0"}], "]"}]}], 
    "]"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"fr5", "=", 
   RowBox[{"MeshJoin", "[", 
    RowBox[{"fr3", ",", "fr4"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fr6", "=", 
   RowBox[{"MeshRotate", "[", 
    RowBox[{"fr5", ",", 
     RowBox[{"StraightLine", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"1", ",", "1", ",", "0"}], "}"}]}], "]"}]}], "]"}]}], 
  ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"fr7", "=", 
    RowBox[{"MeshJoin", "[", 
     RowBox[{"fr5", ",", "fr6"}], "]"}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fr8", "=", 
   RowBox[{"MeshTranslate", "[", 
    RowBox[{"fr7", ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "Pi"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fr9", "=", 
   RowBox[{"MeshTranslate", "[", 
    RowBox[{"fr7", ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", 
       RowBox[{"-", "Pi"}]}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"fr", "=", 
    RowBox[{"MeshJoin", "[", 
     RowBox[{"fr7", ",", "fr8", ",", "fr9"}], "]"}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]"}], "Input"],

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"gr", "=", 
   RowBox[{"Mesh3DToGraphics3D", "[", "fr", "]"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.388139751899102*^9, 3.388139753417307*^9}}]
}, Open  ]]
}, Open  ]]
},
WindowSize->{935, 802},
WindowMargins->{{Automatic, 416}, {Automatic, 120}},
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, 52, 0, 76, "Title"],
Cell[645, 25, 53, 1, 27, "Input"],
Cell[701, 28, 169, 3, 26, "Text"],
Cell[CellGroupData[{
Cell[895, 35, 100, 1, 67, "Section"],
Cell[998, 38, 248, 9, 48, "Input"],
Cell[1249, 49, 828, 26, 58, "Input"],
Cell[2080, 77, 98, 1, 26, "Text"],
Cell[2181, 80, 1536, 49, 92, "Input"],
Cell[3720, 131, 158, 4, 27, "Input"],
Cell[3881, 137, 209, 6, 58, "Input"],
Cell[4093, 145, 83, 2, 27, "Input"],
Cell[4179, 149, 331, 9, 27, "Input"],
Cell[4513, 160, 441, 15, 27, "Input"],
Cell[4957, 177, 1939, 64, 193, "Input"],
Cell[6899, 243, 194, 4, 27, "Input"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)