// Scell.fe
// Triply periodic minimal surface S from Fischer and Koch

// Programmer: Ken Brakke, brakke@susqu.edu, http://www.susqu.edu

/* Commands:
   gogo - typical evolution
   setcolor - to color one side yellow, as in my web page.

   To turn off showing all the edges in the graphics display,
      hit the "e" key in the graphics window.
*/

// 180 degree rotations about respective edges
view_transform_generators 12
swap_colors   0.  0.  1.  0.25    0. -1.  0.  0.25    1.  0.  0. -0.25    0  0  0  1   // a:
swap_colors  -1.  0.  0. -0.25    0.  0. -1. -0.25    0. -1.  0. -0.25    0  0  0  1   // b: 
swap_colors  -1.  0.  0. -0.25    0.  0.  1.  0.25    0.  1.  0. -0.25    0  0  0  1   // c: 
swap_colors   0. -1.  0.  0.25   -1.  0.  0.  0.25    0.  0. -1.  0.25    0  0  0  1   
swap_colors   0.  1.  0. -0.25    1.  0.  0.  0.25    0.  0. -1.  0.25    0  0  0  1   
swap_colors   0.  0. -1. -0.25    0. -1.  0. -0.25   -1.  0.  0. -0.25    0  0  0  1   
swap_colors   0.  0.  1. -0.25    0. -1.  0. -0.25    1.  0.  0.  0.25    0  0  0  1   
swap_colors  -1.  0.  0.  0.25    0.  0. -1.  0.25    0. -1.  0.  0.25    0  0  0  1   
swap_colors  -1.  0.  0.  0.25    0.  0.  1. -0.25    0.  1.  0.  0.25    0  0  0  1   
swap_colors   0. -1.  0. -0.25   -1.  0.  0. -0.25    0.  0. -1. -0.25    0  0  0  1   
swap_colors   0.  1.  0.  0.25    1.  0.  0. -0.25    0.  0. -1. -0.25    0  0  0  1   
swap_colors   0.  0. -1.  0.25    0. -1.  0.  0.25   -1.  0.  0.  0.25    0  0  0  1   

vertices
1  1/4 1/8 0  fixed
2  -1/8 1/8 -3/8  fixed
3  -1/8 0 -1/4  fixed
4  -1/8 3/8 1/8  fixed
5  0 1/4 1/8    fixed
6  -3/8 -1/8 1/8  fixed
7  -1/4 -1/8 0  fixed
8  1/8 -1/8 3/8  fixed
9  1/8 0 1/4  fixed
10  1/8 -3/8 -1/8  fixed 
11  0 -1/4 -1/8  fixed
12  3/8 1/8 -1/8  fixed

edges
1   1 2 fixed
2   2 3 fixed
3   3 4 fixed
4   4 5 fixed
5   5 6 fixed
6   6 7 fixed
7   7 8 fixed
8   8 9 fixed
9   9 10 fixed
10  10 11 fixed
11  11 12 fixed
12  12 1  fixed

faces
1   1 2 3 4 5 6 7 8 9 10 11 12

read

hessian_normal

gogo := { refine edge where valence==1; g 5; u; r; g 10; r; g 10; r; g 10; u; V;
          hessian; hessian;
        }

mgen := { printf "{\n";
          foreach edge ee where fixed do
           printf "rot[%g,%g,%g,%g,%g,%g],\n",ee.vertex[1].x,ee.vertex[1].y,
           ee.vertex[1].z,ee.vertex[2].x,ee.vertex[2].y,ee.vertex[2].z;
          printf "}\n";
        }

setcolor := { set facet backcolor yellow }


/* Commands:
   gogo - typical evolution
   setcolor - to color one side yellow, as in my web page.

   To turn off showing all the edges in the graphics display,
      hit the "e" key in the graphics window.
*/