1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55
|
// quadm.fe
// quadrilateral boundary wire
// test of metric
INTEGRAL_ORDER 7
METRIC
/*
1+4*x1^2 4*x1*x2 2*x1
4*x1*x2 1+4*x2^2 2*x2
2*x1 2*x2 1
*/
1 0 0
0 1 0
0 0 1
constraint 1
formula: x1 = 0
constraint 2
formula: x1 = 1
constraint 3
formula: x2 = 0
constraint 4
formula: x2 = 1
constraint 5
formula: x1^2 + x2^2 + x3 = 0
vertices
1 0 0 0 fixed
2 1 0 -1 fixed
3 1 1 -2 fixed
4 0 1 -1 fixed
edges
1 1 2 fixed constraints 3,5
2 2 3 fixed constraints 2,5
3 3 4 fixed constraints 4,5
4 4 1 fixed constraints 1,5
5 2 4
faces
1 1 5 4
2 2 3 -5
read
// Typical evolution
gogo := { r; g 12; U; r; g 12; r; g 12; r; g 20; }
|