// fe\test.fe: Dump of structure.
// double bubble
// Total energy: 3.84832028287475
SOAPFILM

LINEAR

// SCALE: 0.402401883419685

PARAMETER shift = 1.44

PARAMETER radius = 2.37

TOTAL_TIME 614.665

CONSTRAINT_TOLERANCE: 1e-012

VIEW_MATRIX 
 0.666666666666667   0.000000000000000   0.000000000000000  -0.333333333333333  
 0.000000000000000   0.666666666666667   0.000000000000000  -0.333333333333333  
 0.000000000000000   0.000000000000000   0.666666666666667  -0.333333333333333  
 0.000000000000000   0.000000000000000   0.000000000000000   1.000000000000000  

//potential engery
quantity poten energy global_method facet_vector_integral
vector_integrand
q1: 0
q2: 0
q3: 0.5*z^2-2*z

boundary 1 parameters 1 // right circle
x1: (shift*abs(cos(p1)) + sqrt(-(shift*sin(p1))^2+radius^2))*cos(p1)
x2: (shift*abs(cos(p1)) + sqrt(-(shift*sin(p1))^2+radius^2))*sin(p1)
x3: 0

vertices
1 0 0 0 
2 4*pi/8    boundary 1 fixed
3 5*pi/8    boundary 1 fixed
4 6*pi/8    boundary 1 fixed
5 7*pi/8    boundary 1 fixed
6 8*pi/8    boundary 1 fixed
7 9*pi/8    boundary 1 fixed
8 10*pi/8   boundary 1 fixed
9 11*pi/8   boundary 1 fixed
10 12*pi/8  boundary 1 fixed
11 13*pi/8  boundary 1 fixed
12 14*pi/8  boundary 1 fixed
13 15*pi/8  boundary 1 fixed
14 16*pi/8  boundary 1 fixed
15 17*pi/8  boundary 1 fixed
16 18*pi/8  boundary 1 fixed
17 19*pi/8  boundary 1 fixed
18 -shift 0 0 
19 shift 0 0 

edges /* given by endpoints*/
1 18 1
2 18 2
3 18 3
4 18 4
5 18 5
6 18 6
7 18 7
8 18 8
9 18 9
10 18 10
11 19 1
12 19 10
13 19 11
14 19 12
15 19 13
16 19 14 
17 19 15
18 19 16
19 19 17
20 19 2
21 1 2
22 2 3 boundary 1 fixed
23 3 4 boundary 1 fixed
24 4 5 boundary 1 fixed
25 5 6 boundary 1 fixed
26 6 7 boundary 1 fixed
27 7 8 boundary 1 fixed
28 8 9 boundary 1 fixed
29 9 10 boundary 1 fixed
30 10 1 
31 10 11 boundary 1 fixed
32 11 12 boundary 1 fixed
33 12 13 boundary 1 fixed
34 13 14 boundary 1 fixed
35 14 15 boundary 1 fixed
36 15 16 boundary 1 fixed
37 16 17 boundary 1 fixed
38 17 2 boundary 1 fixed


faces /* oriented edge loop*/
1 1 21 -2 
2 2 22 -3
3 3 23 -4 
4 4 24 -5 
5 5 25 -6 
6 6 26 -7
7 7 27 -8
8 8 28 -9
9 9 29 -10
10 10 30 -1
11 11 -30 -12
12 12 31 -13
13 13 32 -14
14 14 33 -15
15 15 34 -16
16 16 35 -17
17 17 36 -18
18 18 37 -19
19 19 38 -20
20 20 -21 -11


bodies
1 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20

read
// Typical evolution
gogo := {r; r; g 500; u; V; r; u; V; g 500; u; V; r; u; V; g 500; hessian; hessian; }