//genus-3 Lawson surface
// 45-degree quarter-slice
// simplified vertices 11/26/03

quantity smc energy method sq_mean_curvature global

CONSTRAINT 1
function: x2 = 0.0

CONSTRAINT 2
function: x1 = x2

CONSTRAINT 3
function: x3 = 0.0

view_transform_generators 4  // number of matrices below!

// mirror x-axis
-1,0,0,0
0,1,0,0
0,0,1,0
0,0,0,1

// mirror y-axis
1,0,0,0
0,-1,0,0
0,0,1,0
0,0,0,1

// mirror z-axis
1,0,0,0
0,1,0,0
0,0,-1,0
0,0,0,1

// mirror on x=y plane
0,1,0,0
1,0,0,0
0,0,1,0
0,0,0,1


vertices
1	0 0 8      CONSTRAINT 1 2
2	3.5 0 7      CONSTRAINT 1
3	2.3 2.3 7      CONSTRAINT 2
4	4 2 4  
5	5 0 2      CONSTRAINT 1
6	5 1.4 1.4  
7	5 2 0      CONSTRAINT 3       
8	5.5 3 0    CONSTRAINT 3
9       4.5 4.5 0  CONSTRAINT 2 3

20	4.3 0 1    CONSTRAINT 1
21	4.3 0.7 0.7 
22	4.3 1 0    CONSTRAINT 3

11	0 0 2.0      CONSTRAINT 1 2
12	2 0 1.7      CONSTRAINT 1
13	1.4 1.4 1.7  CONSTRAINT 2
14	3 1.5 1.5  
15	3 0 1.5      CONSTRAINT 1
16	3 1.4 1.2 
17	3 1.5 0      CONSTRAINT 3      
18	3 2.0 0  CONSTRAINT 3
19      2.5 2.5 0  CONSTRAINT 2 3

30	3.7 0 1    CONSTRAINT 1
31	3.7 0.7 0.7 
32	3.7 1 0    CONSTRAINT 3


edges
1	1 2   CONSTRAINT 1
2	1 3   CONSTRAINT 2
3	2 3
4	2 4
5	3 4
6	2 5   CONSTRAINT 1   
7	4 5
8	4 9   
9	3 9   CONSTRAINT 2
10	4 6
11	4 8   
12	6 8
13	5 6
14      5 20  CONSTRAINT 1
15      20 6
16      6 21
17      21 7
18      6 7
19      22 7  CONSTRAINT 3
20      20 21
21      21 22
22      7 8   CONSTRAINT 3
23      8 9   CONSTRAINT 3

31	11 12 CONSTRAINT 1
32	11 13 CONSTRAINT 2
33	12 13
34	12 14
35	13 14
36	12 15 CONSTRAINT 1   
37	14 15
38	14 19   
39	13 19 CONSTRAINT 2
40	14 16
41	14 18   
42	16 18
43	15 16
44      15 30 CONSTRAINT 1
45      30 16
46      16 31
47      31 17
48      16 17
49      32 17 CONSTRAINT 3
50      30 31
51      31 32
52      17 18 CONSTRAINT 3
53      18 19 CONSTRAINT 3

60	20 30  CONSTRAINT 1
61 	21 31
62 	22 32 CONSTRAINT 3


faces
1	1 3 -2
2	4 -5 -3
3	6 -7 -4
4	7 13 -10
5	5 8 -9
6	11 23 -8
7       10 12 -11
8	18 22 -12
9       14 15 -13
10	20 -16 -15
11      16 17 -18
12	21 19 -17

20	-20 60 50 -61
21	-21 61 51 -62

31	32 -33 -31
32	33 35 -34
33	34 37 -36
34	40 -43 -37
35	39 -38 -35
36	38 -53 -41
37       41 -42 -40
38	42 -52 -48
39       43 -45 -44
40	45 46 -50
41      48 -47 -46
42	47 -49 -51

bodies
1	 1 2 3 4 5 6 7 8 9 10 11 12  20 21  31 32 33 34 35 36 37 38 39 40 41 42

read
set facet tension 0

sx := { transform_expr "bcd"; }

gogo := { g5; r; g55; V; g55; u; U; g55; {w 0.02; l 1.0; u; g99}10 }

gg := { w 0.03; l 0.8; u; g99; }