From b0bde507f71ed3c026c8a6557466e43a291ac322 Mon Sep 17 00:00:00 2001 From: Leila Ghaffari Date: Thu, 20 Jun 2024 10:47:57 -0700 Subject: [PATCH] add some cube meshes none of which is working --- examples/cube/cube_pec.json | 2 +- examples/cube/mesh/cube.geo | 42 ++ examples/cube/mesh/cube.msh | 489 ++++++++++++++++ examples/cube/mesh/cube1.geo | 69 +++ examples/cube/mesh/cube1.msh | 1048 ++++++++++++++++++++++++++++++++++ examples/cube/mesh/cube2.geo | 15 + examples/cube/mesh/cube2.msh | 724 +++++++++++++++++++++++ 7 files changed, 2388 insertions(+), 1 deletion(-) create mode 100644 examples/cube/mesh/cube.geo create mode 100644 examples/cube/mesh/cube.msh create mode 100644 examples/cube/mesh/cube1.geo create mode 100644 examples/cube/mesh/cube1.msh create mode 100644 examples/cube/mesh/cube2.geo create mode 100644 examples/cube/mesh/cube2.msh diff --git a/examples/cube/cube_pec.json b/examples/cube/cube_pec.json index 7c605a0eb..34cbefb24 100644 --- a/examples/cube/cube_pec.json +++ b/examples/cube/cube_pec.json @@ -7,7 +7,7 @@ }, "Model": { - "Mesh": "mesh/cube.mesh", + "Mesh": "mesh/cube.msh", "L0": 1.0e-2, // cm "Refinement": { diff --git a/examples/cube/mesh/cube.geo b/examples/cube/mesh/cube.geo new file mode 100644 index 000000000..2d6620c51 --- /dev/null +++ b/examples/cube/mesh/cube.geo @@ -0,0 +1,42 @@ +SetFactory("OpenCASCADE"); +gridsize = 100/3; + +// Face +// -- Points +Point(1) = {0,0,0,gridsize}; +Point(2) = {0,0,100,gridsize}; +Point(3) = {0,100,100,gridsize}; +Point(4) = {0,100,0,gridsize}; + +// -- Line +Line(5) = {1,2}; +Line(6) = {2,3}; +Line(7) = {3,4}; +Line(8) = {4,1}; +Line Loop(9) = {5,6,7,8}; + +// -- Surface +Plane Surface(10) = {9}; +Recombine Surface{10}; + +// Extrude to 3D +bodyExtrusion[] = +Extrude { 100,0,0 } +{ + Surface{10}; + Layers{100/gridsize}; + Recombine; +}; + +// Boundaries +Physical Surface("back") = {10}; +Physical Surface("right") = {bodyExtrusion[2]}; +Physical Surface("top") = {bodyExtrusion[3]}; +Physical Surface("left") = {bodyExtrusion[4]}; +Physical Surface("bottom") = {bodyExtrusion[5]}; +Physical Surface("front") = {bodyExtrusion[0]}; + +// Volume +Physical Volume("mesh") = {bodyExtrusion[1]}; +Mesh 3; +Coherence Mesh; \ No newline at end of file diff --git a/examples/cube/mesh/cube.msh b/examples/cube/mesh/cube.msh new file mode 100644 index 000000000..924724a46 --- /dev/null +++ b/examples/cube/mesh/cube.msh @@ -0,0 +1,489 @@ +$MeshFormat +4.1 0 8 +$EndMeshFormat +$PhysicalNames +7 +2 1 "back" +2 2 "right" +2 3 "top" +2 4 "left" +2 5 "bottom" +2 6 "front" +3 7 "mesh" +$EndPhysicalNames +$Entities +8 12 6 1 +1 0 0 0 0 +2 0 0 100 0 +3 0 100 100 0 +4 0 100 0 0 +5 100 0 0 0 +6 100 0 100 0 +7 100 100 100 0 +8 100 100 0 0 +5 -1e-07 -1e-07 -1.000000011686097e-07 1e-07 1e-07 100.0000001 0 2 1 -2 +6 -1e-07 -1.000000011686097e-07 99.99999990000001 1e-07 100.0000001 100.0000001 0 2 2 -3 +7 -1e-07 99.99999990000001 -1.000000011686097e-07 1e-07 100.0000001 100.0000001 0 2 3 -4 +8 -1e-07 -1.000000011686097e-07 -1e-07 1e-07 100.0000001 1e-07 0 2 4 -1 +9 -1.000000011686097e-07 -1e-07 -1e-07 100.0000001 1e-07 1e-07 0 2 1 -5 +10 -1.000000011686097e-07 -1e-07 99.99999990000001 100.0000001 1e-07 100.0000001 0 2 2 -6 +11 99.99999990000001 -1e-07 -1.000000011686097e-07 100.0000001 1e-07 100.0000001 0 2 5 -6 +12 -1.000000011686097e-07 99.99999990000001 99.99999990000001 100.0000001 100.0000001 100.0000001 0 2 3 -7 +13 99.99999990000001 -1.000000011686097e-07 99.99999990000001 100.0000001 100.0000001 100.0000001 0 2 6 -7 +14 -1.000000011686097e-07 99.99999990000001 -1e-07 100.0000001 100.0000001 1e-07 0 2 4 -8 +15 99.99999990000001 99.99999990000001 -1.000000011686097e-07 100.0000001 100.0000001 100.0000001 0 2 7 -8 +16 99.99999990000001 -1.000000011686097e-07 -1e-07 100.0000001 100.0000001 1e-07 0 2 8 -5 +10 -1e-07 -1.000000011686097e-07 -1.000000011686097e-07 1e-07 100.0000001 100.0000001 1 1 4 5 6 7 8 +11 -1.000000011686097e-07 -1e-07 -1.000000011686097e-07 100.0000001 1e-07 100.0000001 1 2 4 9 11 -10 -5 +12 -1.000000011686097e-07 -1.000000011686097e-07 99.99999990000001 100.0000001 100.0000001 100.0000001 1 3 4 10 13 -12 -6 +13 -1.000000011686097e-07 99.99999990000001 -1.000000011686097e-07 100.0000001 100.0000001 100.0000001 1 4 4 12 15 -14 -7 +14 -1.000000011686097e-07 -1.000000011686097e-07 -1e-07 100.0000001 100.0000001 1e-07 1 5 4 14 16 -9 -8 +15 99.99999990000001 -1.000000011686097e-07 -1.000000011686097e-07 100.0000001 100.0000001 100.0000001 1 6 4 11 13 15 16 +1 -1.000000011686097e-07 -1.000000011686097e-07 -1.000000011686097e-07 100.0000001 100.0000001 100.0000001 1 7 6 11 12 13 14 10 -15 +$EndEntities +$Nodes +27 124 1 124 +0 1 0 1 +1 +0 0 0 +0 2 0 1 +2 +0 0 100 +0 3 0 1 +3 +0 100 100 +0 4 0 1 +4 +0 100 0 +0 5 0 1 +5 +100 0 0 +0 6 0 1 +6 +100 0 100 +0 7 0 1 +7 +100 100 100 +0 8 0 1 +8 +100 100 0 +1 5 0 3 +9 +10 +11 +0 0 25 +0 0 50 +0 0 75 +1 6 0 3 +12 +13 +14 +0 25 100 +0 50 100 +0 75 100 +1 7 0 3 +15 +16 +17 +0 100 75 +0 100 50 +0 100 25 +1 8 0 3 +18 +19 +20 +0 75 0 +0 50 0 +0 25 0 +1 9 0 2 +21 +22 +33.33333333333333 0 0 +66.66666666666666 0 0 +1 10 0 2 +23 +24 +33.33333333333333 0 100 +66.66666666666666 0 100 +1 11 0 3 +25 +26 +27 +100 0 25 +100 0 50 +100 0 75 +1 12 0 2 +28 +29 +33.33333333333333 100 100 +66.66666666666666 100 100 +1 13 0 3 +30 +31 +32 +100 25 100 +100 50 100 +100 75 100 +1 14 0 2 +33 +34 +33.33333333333333 100 0 +66.66666666666666 100 0 +1 15 0 3 +35 +36 +37 +100 100 75 +100 100 50 +100 100 25 +1 16 0 3 +38 +39 +40 +100 75 0 +100 50 0 +100 25 0 +2 10 0 15 +41 +42 +43 +44 +45 +46 +47 +48 +49 +50 +51 +52 +53 +54 +55 +0 37.6059357266501 78.8307029219294 +0 21.43427655108268 40.72204128987809 +0 80.85941028877747 37.68069268138613 +0 58.64030531554702 72.54195134384774 +0 48.60942435472928 51.85657868617768 +0 27.99242978990447 60.98094056314455 +0 41.1359678245144 32.2450375555237 +0 76.92476999784536 62.68318576745698 +0 35.97820959953506 18.11422321226597 +0 59.57162545096222 24.85529455971059 +0 22.49664856597377 85.3474726464806 +0 84.64424817292736 77.62465917641187 +0 79.00151556333624 22.1850598725291 +0 67.43756762511735 44.08313791101965 +0 18.69294775670395 18.35388528705161 +2 11 0 6 +56 +57 +58 +59 +60 +61 +33.33333333333333 0 25 +66.66666666666666 0 25 +33.33333333333333 0 50 +66.66666666666666 0 50 +33.33333333333333 0 75 +66.66666666666666 0 75 +2 12 0 6 +62 +63 +64 +65 +66 +67 +33.33333333333333 25 100 +66.66666666666666 25 100 +33.33333333333333 50 100 +66.66666666666666 50 100 +33.33333333333333 75 100 +66.66666666666666 75 100 +2 13 0 6 +68 +69 +70 +71 +72 +73 +33.33333333333333 100 75 +66.66666666666666 100 75 +33.33333333333333 100 50 +66.66666666666666 100 50 +33.33333333333333 100 25 +66.66666666666666 100 25 +2 14 0 6 +74 +75 +76 +77 +78 +79 +33.33333333333333 75 0 +66.66666666666666 75 0 +33.33333333333333 50 0 +66.66666666666666 50 0 +33.33333333333333 25 0 +66.66666666666666 25 0 +2 15 0 15 +80 +81 +82 +83 +84 +85 +86 +87 +88 +89 +90 +91 +92 +93 +94 +100 37.6059357266501 78.8307029219294 +100 21.43427655108268 40.72204128987809 +100 80.85941028877747 37.68069268138613 +100 58.64030531554702 72.54195134384774 +100 48.60942435472928 51.85657868617768 +100 27.99242978990447 60.98094056314455 +100 41.1359678245144 32.2450375555237 +100 76.92476999784536 62.68318576745698 +100 35.97820959953506 18.11422321226597 +100 59.57162545096222 24.85529455971059 +100 22.49664856597377 85.3474726464806 +100 84.64424817292736 77.62465917641187 +100 79.00151556333624 22.1850598725291 +100 67.43756762511735 44.08313791101965 +100 18.69294775670395 18.35388528705161 +3 1 0 30 +95 +96 +97 +98 +99 +100 +101 +102 +103 +104 +105 +106 +107 +108 +109 +110 +111 +112 +113 +114 +115 +116 +117 +118 +119 +120 +121 +122 +123 +124 +33.33333333333333 37.6059357266501 78.8307029219294 +66.66666666666666 37.6059357266501 78.8307029219294 +33.33333333333333 21.43427655108268 40.72204128987809 +66.66666666666666 21.43427655108268 40.72204128987809 +33.33333333333333 80.85941028877747 37.68069268138613 +66.66666666666666 80.85941028877747 37.68069268138613 +33.33333333333333 58.64030531554702 72.54195134384774 +66.66666666666666 58.64030531554702 72.54195134384774 +33.33333333333333 48.60942435472928 51.85657868617768 +66.66666666666666 48.60942435472928 51.85657868617768 +33.33333333333333 27.99242978990447 60.98094056314455 +66.66666666666666 27.99242978990447 60.98094056314455 +33.33333333333333 41.1359678245144 32.2450375555237 +66.66666666666666 41.1359678245144 32.2450375555237 +33.33333333333333 76.92476999784536 62.68318576745698 +66.66666666666666 76.92476999784536 62.68318576745698 +33.33333333333333 35.97820959953506 18.11422321226597 +66.66666666666666 35.97820959953506 18.11422321226597 +33.33333333333333 59.57162545096222 24.85529455971059 +66.66666666666666 59.57162545096222 24.85529455971059 +33.33333333333333 22.49664856597377 85.3474726464806 +66.66666666666666 22.49664856597377 85.3474726464806 +33.33333333333333 84.64424817292736 77.62465917641187 +66.66666666666666 84.64424817292736 77.62465917641187 +33.33333333333333 79.00151556333624 22.1850598725291 +66.66666666666666 79.00151556333624 22.1850598725291 +33.33333333333333 67.43756762511735 44.08313791101965 +66.66666666666666 67.43756762511735 44.08313791101965 +33.33333333333333 18.69294775670395 18.35388528705161 +66.66666666666666 18.69294775670395 18.35388528705161 +$EndNodes +$Elements +7 158 1 158 +2 10 3 22 +1 47 42 46 45 +2 44 48 54 45 +3 49 47 50 19 +4 46 41 44 45 +5 14 44 41 13 +6 11 46 42 10 +7 14 3 15 52 +8 12 51 11 2 +9 18 53 17 4 +10 16 43 54 48 +11 54 43 53 50 +12 10 42 55 9 +13 45 54 50 47 +14 16 48 52 15 +15 16 17 53 43 +16 13 41 51 12 +17 47 49 55 42 +18 9 55 20 1 +19 14 52 48 44 +20 11 51 41 46 +21 19 50 53 18 +22 19 20 55 49 +2 11 3 12 +23 1 21 56 9 +24 21 22 57 56 +25 22 5 25 57 +26 9 56 58 10 +27 56 57 59 58 +28 57 25 26 59 +29 10 58 60 11 +30 58 59 61 60 +31 59 26 27 61 +32 11 60 23 2 +33 60 61 24 23 +34 61 27 6 24 +2 12 3 12 +35 2 23 62 12 +36 23 24 63 62 +37 24 6 30 63 +38 12 62 64 13 +39 62 63 65 64 +40 63 30 31 65 +41 13 64 66 14 +42 64 65 67 66 +43 65 31 32 67 +44 14 66 28 3 +45 66 67 29 28 +46 67 32 7 29 +2 13 3 12 +47 3 28 68 15 +48 28 29 69 68 +49 29 7 35 69 +50 15 68 70 16 +51 68 69 71 70 +52 69 35 36 71 +53 16 70 72 17 +54 70 71 73 72 +55 71 36 37 73 +56 17 72 33 4 +57 72 73 34 33 +58 73 37 8 34 +2 14 3 12 +59 4 33 74 18 +60 33 34 75 74 +61 34 8 38 75 +62 18 74 76 19 +63 74 75 77 76 +64 75 38 39 77 +65 19 76 78 20 +66 76 77 79 78 +67 77 39 40 79 +68 20 78 21 1 +69 78 79 22 21 +70 79 40 5 22 +2 15 3 22 +71 86 84 85 81 +72 83 84 93 87 +73 88 39 89 86 +74 85 84 83 80 +75 32 31 80 83 +76 27 26 81 85 +77 32 91 35 7 +78 30 6 27 90 +79 38 8 37 92 +80 36 87 93 82 +81 93 89 92 82 +82 26 25 94 81 +83 84 86 89 93 +84 36 35 91 87 +85 36 82 92 37 +86 31 30 90 80 +87 86 81 94 88 +88 25 5 40 94 +89 32 83 87 91 +90 27 85 80 90 +91 39 38 92 89 +92 39 88 94 40 +3 1 5 66 +93 47 45 46 42 107 103 105 97 +94 107 103 105 97 108 104 106 98 +95 108 104 106 98 86 84 85 81 +96 44 45 54 48 101 103 121 109 +97 101 103 121 109 102 104 122 110 +98 102 104 122 110 83 84 93 87 +99 49 19 50 47 111 76 113 107 +100 111 76 113 107 112 77 114 108 +101 112 77 114 108 88 39 89 86 +102 46 45 44 41 105 103 101 95 +103 105 103 101 95 106 104 102 96 +104 106 104 102 96 85 84 83 80 +105 14 13 41 44 66 64 95 101 +106 66 64 95 101 67 65 96 102 +107 67 65 96 102 32 31 80 83 +108 11 10 42 46 60 58 97 105 +109 60 58 97 105 61 59 98 106 +110 61 59 98 106 27 26 81 85 +111 14 52 15 3 66 117 68 28 +112 66 117 68 28 67 118 69 29 +113 67 118 69 29 32 91 35 7 +114 12 2 11 51 62 23 60 115 +115 62 23 60 115 63 24 61 116 +116 63 24 61 116 30 6 27 90 +117 18 4 17 53 74 33 72 119 +118 74 33 72 119 75 34 73 120 +119 75 34 73 120 38 8 37 92 +120 16 48 54 43 70 109 121 99 +121 70 109 121 99 71 110 122 100 +122 71 110 122 100 36 87 93 82 +123 54 50 53 43 121 113 119 99 +124 121 113 119 99 122 114 120 100 +125 122 114 120 100 93 89 92 82 +126 10 9 55 42 58 56 123 97 +127 58 56 123 97 59 57 124 98 +128 59 57 124 98 26 25 94 81 +129 45 47 50 54 103 107 113 121 +130 103 107 113 121 104 108 114 122 +131 104 108 114 122 84 86 89 93 +132 16 15 52 48 70 68 117 109 +133 70 68 117 109 71 69 118 110 +134 71 69 118 110 36 35 91 87 +135 16 43 53 17 70 99 119 72 +136 70 99 119 72 71 100 120 73 +137 71 100 120 73 36 82 92 37 +138 13 12 51 41 64 62 115 95 +139 64 62 115 95 65 63 116 96 +140 65 63 116 96 31 30 90 80 +141 47 42 55 49 107 97 123 111 +142 107 97 123 111 108 98 124 112 +143 108 98 124 112 86 81 94 88 +144 9 1 20 55 56 21 78 123 +145 56 21 78 123 57 22 79 124 +146 57 22 79 124 25 5 40 94 +147 14 44 48 52 66 101 109 117 +148 66 101 109 117 67 102 110 118 +149 67 102 110 118 32 83 87 91 +150 11 46 41 51 60 105 95 115 +151 60 105 95 115 61 106 96 116 +152 61 106 96 116 27 85 80 90 +153 19 18 53 50 76 74 119 113 +154 76 74 119 113 77 75 120 114 +155 77 75 120 114 39 38 92 89 +156 19 49 55 20 76 111 123 78 +157 76 111 123 78 77 112 124 79 +158 77 112 124 79 39 88 94 40 +$EndElements diff --git a/examples/cube/mesh/cube1.geo b/examples/cube/mesh/cube1.geo new file mode 100644 index 000000000..ec852011f --- /dev/null +++ b/examples/cube/mesh/cube1.geo @@ -0,0 +1,69 @@ +// Define a unit cube with 3 elements along each direction + +// Define the size of the cube and number of elements +lc = 1.0 / 3.0; // characteristic length (edge length of the smallest elements) + +// Define points +Point(1) = {0, 0, 0, lc}; +Point(2) = {1, 0, 0, lc}; +Point(3) = {1, 1, 0, lc}; +Point(4) = {0, 1, 0, lc}; +Point(5) = {0, 0, 1, lc}; +Point(6) = {1, 0, 1, lc}; +Point(7) = {1, 1, 1, lc}; +Point(8) = {0, 1, 1, lc}; + +// Define lines +Line(1) = {1, 2}; +Line(2) = {2, 3}; +Line(3) = {3, 4}; +Line(4) = {4, 1}; +Line(5) = {5, 6}; +Line(6) = {6, 7}; +Line(7) = {7, 8}; +Line(8) = {8, 5}; +Line(9) = {1, 5}; +Line(10) = {2, 6}; +Line(11) = {3, 7}; +Line(12) = {4, 8}; + +// Define surfaces +Line Loop(1) = {1, 2, 3, 4}; // bottom face +Plane Surface(1) = {1}; + +Line Loop(2) = {5, 6, 7, 8}; // top face +Plane Surface(2) = {2}; + +Line Loop(3) = {1, 10, -5, -9}; // front face +Plane Surface(3) = {3}; + +Line Loop(4) = {2, 11, -6, -10}; // back face +Plane Surface(4) = {4}; + +Line Loop(5) = {3, 12, -7, -11}; // right face +Plane Surface(5) = {5}; + +Line Loop(6) = {4, 9, -8, -12}; // left face +Plane Surface(6) = {6}; + +// Define volume +Surface Loop(1) = {1, 2, 3, 4, 5, 6}; +Volume(1) = {1}; + +// Define physical entities +Physical Volume(1) = {1}; // volume +Physical Surface(1) = {1}; // bottom face +Physical Surface(2) = {2}; // top face +Physical Surface(3) = {3}; // front face +Physical Surface(4) = {4}; // back face +Physical Surface(5) = {5}; // right face +Physical Surface(6) = {6}; // left face + +// Meshing options +Transfinite Line {1, 2, 3, 4, 5, 6, 7, 8} = 3; // divide each edge into 3 segments + +// Generate mesh +Mesh.RecombineAll = 1; // ensure elements are tetrahedra (if using 3D) + +// Save mesh +Mesh 3; diff --git a/examples/cube/mesh/cube1.msh b/examples/cube/mesh/cube1.msh new file mode 100644 index 000000000..a079ace36 --- /dev/null +++ b/examples/cube/mesh/cube1.msh @@ -0,0 +1,1048 @@ +$MeshFormat +4.1 0 8 +$EndMeshFormat +$Entities +8 12 6 1 +1 0 0 0 0 +2 1 0 0 0 +3 1 1 0 0 +4 0 1 0 0 +5 0 0 1 0 +6 1 0 1 0 +7 1 1 1 0 +8 0 1 1 0 +1 0 0 0 1 0 0 0 2 1 -2 +2 1 0 0 1 1 0 0 2 2 -3 +3 0 1 0 1 1 0 0 2 3 -4 +4 0 0 0 0 1 0 0 2 4 -1 +5 0 0 1 1 0 1 0 2 5 -6 +6 1 0 1 1 1 1 0 2 6 -7 +7 0 1 1 1 1 1 0 2 7 -8 +8 0 0 1 0 1 1 0 2 8 -5 +9 0 0 0 0 0 1 0 2 1 -5 +10 1 0 0 1 0 1 0 2 2 -6 +11 1 1 0 1 1 1 0 2 3 -7 +12 0 1 0 0 1 1 0 2 4 -8 +1 0 0 0 1 1 0 1 1 4 1 2 3 4 +2 0 0 1 1 1 1 1 2 4 5 6 7 8 +3 0 0 0 1 0 1 1 3 4 1 10 -5 -9 +4 1 0 0 1 1 1 1 4 4 2 11 -6 -10 +5 0 1 0 1 1 1 1 5 4 3 12 -7 -11 +6 0 0 0 0 1 1 1 6 4 4 9 -8 -12 +1 0 0 0 1 1 1 1 1 6 1 2 3 4 5 6 +$EndEntities +$Nodes +27 149 1 149 +0 1 0 1 +1 +0 0 0 +0 2 0 1 +2 +1 0 0 +0 3 0 1 +3 +1 1 0 +0 4 0 1 +4 +0 1 0 +0 5 0 1 +5 +0 0 1 +0 6 0 1 +6 +1 0 1 +0 7 0 1 +7 +1 1 1 +0 8 0 1 +8 +0 1 1 +1 1 0 1 +9 +0.4999999999986921 0 0 +1 2 0 1 +10 +1 0.4999999999986921 0 +1 3 0 1 +11 +0.5000000000020595 1 0 +1 4 0 1 +12 +0 0.5000000000020595 0 +1 5 0 1 +13 +0.4999999999986921 0 1 +1 6 0 1 +14 +1 0.4999999999986921 1 +1 7 0 1 +15 +0.5000000000020595 1 1 +1 8 0 1 +16 +0 0.5000000000020595 1 +1 9 0 3 +17 +18 +19 +0 0 0.2499999999994123 +0 0 0.499999999998694 +0 0 0.7499999999993415 +1 10 0 3 +20 +21 +22 +1 0 0.2499999999994123 +1 0 0.499999999998694 +1 0 0.7499999999993415 +1 11 0 3 +23 +24 +25 +1 1 0.2499999999994123 +1 1 0.499999999998694 +1 1 0.7499999999993415 +1 12 0 3 +26 +27 +28 +0 1 0.2499999999994123 +0 1 0.499999999998694 +0 1 0.7499999999993415 +2 1 0 3 +29 +30 +31 +0.500868160179623 0.5000407391735242 0 +0.3538373570722549 0.3535353995019713 0 +0.6494286120981281 0.6455331172534629 0 +2 2 0 3 +32 +33 +34 +0.500868160179623 0.5000407391735242 1 +0.3538373570722549 0.3535353995019713 1 +0.6494286120981281 0.6455331172534629 1 +2 3 0 5 +35 +36 +37 +38 +39 +0.7974353479321052 0 0.6877862819462278 +0.2129255667833261 0 0.6965629630906244 +0.8010066550213139 0 0.4317669806347234 +0.5024448049034224 0 0.425279600534871 +0.2049646025880945 0 0.4296239968590346 +2 4 0 5 +40 +41 +42 +43 +44 +1 0.7870744332158536 0.6965629630906426 +1 0.2025646520676803 0.6877862819463714 +1 0.19899334497806 0.431766980634766 +1 0.4975551950954128 0.4252796005348975 +1 0.7950353974112077 0.4296239968590537 +2 5 0 5 +45 +46 +47 +48 +49 +0.2025646520679871 1 0.6877862819462277 +0.7870744332169047 1 0.6965629630906213 +0.1989933449788684 1 0.4317669806347214 +0.497555195096912 1 0.4252796005348685 +0.7950353974120982 1 0.4296239968590294 +2 6 0 5 +50 +51 +52 +53 +54 +0 0.2025646520679871 0.6877862819462277 +0 0.7870744332169047 0.6965629630906213 +0 0.1989933449788684 0.4317669806347214 +0 0.497555195096912 0.4252796005348684 +0 0.7950353974120982 0.4296239968590294 +3 1 0 95 +55 +56 +57 +58 +59 +60 +61 +62 +63 +64 +65 +66 +67 +68 +69 +70 +71 +72 +73 +74 +75 +76 +77 +78 +79 +80 +81 +82 +83 +84 +85 +86 +87 +88 +89 +90 +91 +92 +93 +94 +95 +96 +97 +98 +99 +100 +101 +102 +103 +104 +105 +106 +107 +108 +109 +110 +111 +112 +113 +114 +115 +116 +117 +118 +119 +120 +121 +122 +123 +124 +125 +126 +127 +128 +129 +130 +131 +132 +133 +134 +135 +136 +137 +138 +139 +140 +141 +142 +143 +144 +145 +146 +147 +148 +149 +0.1674207824539155 0.5534795165914752 0.7429915539693416 +0.5281601325430061 0.5376671863817584 0.5142992676660743 +0.5190556788080727 0.1971770907608744 0.7589273910962548 +0.2882350156407151 0.2931432811575913 0.7104380230346685 +0.2904972297322184 0.7081449835140236 0.6702993452273934 +0.7871118179315179 0.3339268205178395 0.719492114327923 +0.9550152443978956 0.4394843237660069 0.2700490516878797 +0.6821309125896221 0.8503663956569449 0.7116967897694256 +0.3137521083884309 0.3123911352212271 0.3839768025215791 +0.7220693403136235 0.1558580984882748 0.8836336997540579 +0.7631603319939644 0.6957523037103721 0.3093540525728788 +0.6535867454148008 0.3330684769100799 0.3301515557342764 +0.748754507456579 0.5805647993193322 0.745508861527491 +0.28417168298705 0.7084565652651623 0.3238408343578692 +0.2088535057263311 0.7233035260120336 0.8307846316397498 +0.5600160247411565 0.4019860273593701 0.7458430984021528 +0.5345194475406101 0.6680805282488037 0.7484467121923934 +0.4570142440966071 0.7778265892709529 0.5329236523405105 +0.2516642202314041 0.5016862188557994 0.5488835054784834 +0.408006440335403 0.8496590015277079 0.7365904811259594 +0.7866926748504027 0.7196820971340713 0.5767528820233768 +0.801536499629657 0.4298882052265103 0.5498516366080943 +0.8663220186545487 0.4116923541339826 0.7622922960428037 +0.1779835591127313 0.8273804205401304 0.612023437700934 +0.8531414615501419 0.6216101467060552 0.8147599274467597 +0.6658672414115111 0.2311964620498285 0.5530459913298923 +0.1869163043682643 0.2023768598260871 0.5811440194110624 +0.4142736929522999 0.2014985029771879 0.5375931090015299 +0.345840577228217 0.5301806379249432 0.7773140589176472 +0.1786339262430983 0.3411751388026198 0.7972110206848291 +0.8416750874614651 0.1583249125386375 0.6327529833365979 +0.8546430426259828 0.251955365000736 0.2219351932182087 +0.3610208016285903 0.1745019172164392 0.8325016517600635 +0.1800974184628683 0.4783774081166433 0.2269164959921525 +0.8301229655144178 0.8295619408724958 0.7427191605051309 +0.5037990272849764 0.6322660577740594 0.2205947356610432 +0.4374881888016602 0.2164440013856706 0.2399325559912465 +0.6222814324814426 0.8193270839913476 0.8762190626027534 +0.6538429787521974 0.8965896282641519 0.79125864039266 +0.8754051881852166 0.5344514199288928 0.6691824641774401 +0.211666037014531 0.2197996052567972 0.259569218509684 +0.5365821485458884 0.7962282039816932 0.2501829072568148 +0.6191756130889289 0.287378068958063 0.8948969657154443 +0.646795436297198 0.8182938123347052 0.4775339906118063 +0.8212459260253788 0.8787064610771852 0.5916822669496049 +0.425816003044372 0.4268687235636786 0.2247652514703541 +0.5916348529112476 0.2511687836478189 0.1272032489533531 +0.7264371813150733 0.1852129686727517 0.776047065106618 +0.7799546458227461 0.5368824010674271 0.188912122718714 +0.1457886584882689 0.1455700139682676 0.7456771844586016 +0.1153617174611301 0.8833670400308884 0.7611469247891579 +0.2575112427516381 0.2302801256238736 0.1305811019527028 +0.7582761247403024 0.7576842824635533 0.1441293977478046 +0.7875741930694378 0.4113934641064198 0.1605472281843 +0.2136763793129695 0.5883940346693888 0.1810680198318228 +0.5955846110366332 0.213803704380484 0.03821458425470083 +0.4080774430106722 0.7884038577442863 0.1675041469809888 +0.2092822296178801 0.208901686346849 0.8671813527536747 +0.7799581698803135 0.7752140141976021 0.868608638809155 +0.771695226450963 0.3636782709296283 0.8586049863697229 +0.2228486366990566 0.5520996581520117 0.8890371618459596 +0.5909041384039783 0.1941205388598714 0.9075257166512823 +0.3921003127405192 0.7565003223613107 0.8587475042571067 +0.3550837435683837 0.1575557789960127 0.6378666401211325 +0.1036940953027152 0.08650980516090523 0.5942077659514047 +0.6502217019638834 0.1580263191036872 0.6362082157789555 +0.8930784375691317 0.06019444698603681 0.5978520749589592 +0.3144347516910219 0.1334751584738271 0.2608974981801674 +0.6670083160843355 0.1487298256826016 0.2377791768521491 +0.1782313916955046 0.1280469984108517 0.8616407407724914 +0.8086228164746284 0.09107613542257691 0.8599351555082412 +0.9044889565629333 0.6532612212424351 0.6754907566835836 +0.8958793023419975 0.8946436957629933 0.5947511151394069 +0.9050292974255363 0.3632933226513138 0.651180510006245 +0.9097420426522265 0.1003894992614351 0.5923883156447932 +0.9657950099387258 0.3205036597456594 0.2257714836199959 +0.818113776797283 0.6956118166859163 0.2136575008314129 +0.8573713999842805 0.1817442384590212 0.851657438099172 +0.8970462307473724 0.8109633935667312 0.8507006162770029 +0.6252197436343142 0.904981915694632 0.6406209408429453 +0.8818950578455461 0.9312626300648108 0.6200094969748375 +0.3412951833306019 0.8903769317552416 0.6293297155057672 +0.1116635838382449 0.9301034861562817 0.6159277738533527 +0.6981476481277674 0.8140961712989514 0.2137258993484745 +0.2690611438221781 0.8698132421906811 0.2408539561418523 +0.8087120750289876 0.9139362068653413 0.8524731490020848 +0.2052234056002745 0.8858471275144386 0.8472428810020498 +0.1175910527964432 0.6474302527619328 0.6338942575461223 +0.09316864834391006 0.8955274576572507 0.5940467399869215 +0.1580263191037795 0.3497782980364568 0.6362082157789544 +0.0670690556463219 0.1106750950075519 0.616027128358326 +0.1195184152919033 0.3001036158839015 0.2023632509322149 +0.06196808060673629 0.6981476481277674 0.2137258993484744 +0.1286074768946308 0.1756411630175116 0.8594465704863923 +0.09893387302519846 0.8017268185849988 0.842838680792444 +$EndNodes +$Elements +12 672 1 672 +2 1 3 6 +1 12 1 9 30 +2 10 3 11 31 +3 2 10 31 29 +4 4 12 30 29 +5 2 29 30 9 +6 4 29 31 11 +2 2 3 6 +7 16 5 13 33 +8 14 7 15 34 +9 6 14 34 32 +10 8 16 33 32 +11 6 32 33 13 +12 8 32 34 15 +2 3 2 4 +13 17 1 39 +14 37 2 20 +15 39 18 17 +16 20 21 37 +2 3 3 8 +17 13 36 39 38 +18 39 36 19 18 +19 13 38 37 35 +20 22 35 37 21 +21 38 39 1 9 +22 38 9 2 37 +23 13 5 19 36 +24 13 35 22 6 +2 4 2 4 +25 20 2 42 +26 44 3 23 +27 42 21 20 +28 23 24 44 +2 4 3 8 +29 14 43 44 40 +30 44 24 25 40 +31 14 41 42 43 +32 22 21 42 41 +33 43 42 2 10 +34 43 10 3 44 +35 14 6 22 41 +36 14 40 25 7 +2 5 2 4 +37 23 3 49 +38 47 4 26 +39 49 24 23 +40 26 27 47 +2 5 3 8 +41 15 46 49 48 +42 49 46 25 24 +43 15 48 47 45 +44 28 45 47 27 +45 48 49 3 11 +46 48 11 4 47 +47 15 7 25 46 +48 15 45 28 8 +2 6 2 4 +49 52 1 17 +50 26 4 54 +51 17 18 52 +52 54 27 26 +2 6 3 8 +53 16 51 54 53 +54 54 51 28 27 +55 16 53 52 50 +56 19 50 52 18 +57 53 12 1 52 +58 53 54 4 12 +59 16 50 19 5 +60 16 8 28 51 +3 1 4 568 +61 40 126 79 89 +62 132 125 6 64 +63 51 55 142 59 +64 40 75 126 89 +65 16 112 148 84 +66 58 73 70 82 +67 67 71 62 75 +68 56 70 73 82 +69 134 74 62 93 +70 56 63 68 100 +71 108 100 29 101 +72 68 72 56 90 +73 67 79 126 89 +74 71 72 62 75 +75 114 60 132 102 +76 108 66 100 101 +77 10 131 107 103 +78 60 76 128 85 +79 68 72 59 73 +80 14 34 113 79 +81 16 33 112 84 +82 56 68 63 73 +83 59 73 142 78 +84 51 69 59 105 +85 126 75 67 89 +86 56 73 63 82 +87 51 55 59 69 +88 29 108 90 100 +89 11 138 48 96 +90 63 66 56 82 +91 1 52 17 95 +92 48 139 11 96 +93 68 73 59 78 +94 72 90 68 96 +95 142 59 55 73 +96 59 72 68 78 +97 56 66 63 100 +98 55 58 144 73 +99 134 72 62 74 +100 13 112 33 87 +101 132 60 114 77 +102 31 107 11 96 +103 66 90 108 100 +104 108 90 66 103 +105 113 79 67 89 +106 56 68 90 100 +107 15 117 141 74 +108 56 71 67 75 +109 67 70 34 71 +110 66 80 56 82 +111 141 117 69 74 +112 144 58 55 84 +113 132 125 64 102 +114 117 69 59 83 +115 1 17 39 95 +116 56 80 70 82 +117 33 116 87 97 +118 114 77 67 79 +119 14 132 114 77 +120 114 34 67 70 +121 66 82 63 91 +122 33 16 115 84 +123 51 149 69 105 +124 56 90 72 98 +125 14 114 34 79 +126 58 63 73 82 +127 12 147 53 88 +128 62 71 67 92 +129 56 72 68 73 +130 56 72 71 75 +131 113 67 34 92 +132 9 122 38 91 +133 58 144 73 81 +134 61 66 65 103 +135 118 58 36 87 +136 53 146 12 88 +137 56 90 66 100 +138 70 73 58 83 +139 38 123 9 91 +140 56 71 59 83 +141 115 55 69 83 +142 59 69 55 83 +143 31 11 111 96 +144 109 29 90 100 +145 107 10 3 131 +146 53 55 144 73 +147 3 11 107 138 +148 32 34 70 71 +149 33 116 13 87 +150 60 67 114 77 +151 66 90 65 103 +152 1 9 106 122 +153 2 10 108 86 +154 53 142 55 73 +155 107 31 65 96 +156 12 30 109 88 +157 38 123 66 80 +158 106 12 1 146 +159 29 100 30 101 +160 51 149 55 69 +161 59 69 117 74 +162 114 67 34 79 +163 30 9 110 101 +164 114 32 34 70 +165 56 59 73 83 +166 31 29 108 90 +167 109 90 68 100 +168 117 59 71 83 +169 90 96 72 98 +170 65 66 61 76 +171 56 73 70 83 +172 114 67 60 70 +173 56 59 72 73 +174 33 112 58 87 +175 67 71 34 92 +176 62 74 71 92 +177 38 82 66 91 +178 33 115 83 84 +179 142 51 59 78 +180 11 107 138 96 +181 34 67 113 79 +182 57 87 116 97 +183 58 81 63 82 +184 62 72 71 74 +185 56 66 65 76 +186 74 92 62 93 +187 60 67 77 94 +188 139 48 68 96 +189 31 90 108 103 +190 12 4 147 109 +191 58 73 63 81 +192 67 75 62 89 +193 111 11 4 139 +194 65 66 56 90 +195 53 73 63 88 +196 35 22 85 102 +197 123 38 66 91 +198 60 114 97 102 +199 60 80 76 85 +200 107 65 138 96 +201 108 86 66 101 +202 55 73 59 83 +203 55 58 73 83 +204 33 70 58 83 +205 67 79 77 94 +206 91 100 66 101 +207 2 130 10 86 +208 117 71 32 83 +209 139 68 111 96 +210 56 59 71 72 +211 30 100 91 101 +212 67 75 126 94 +213 65 90 56 98 +214 53 68 73 88 +215 75 89 40 99 +216 35 125 22 102 +217 147 68 53 88 +218 109 68 88 100 +219 63 91 122 95 +220 122 9 106 91 +221 14 77 114 79 +222 112 58 87 104 +223 118 58 57 82 +224 48 72 68 96 +225 55 83 115 84 +226 140 133 113 89 +227 109 68 147 88 +228 31 90 65 96 +229 36 58 118 81 +230 36 87 58 104 +231 123 91 66 101 +232 67 126 79 94 +233 38 80 66 82 +234 60 97 70 102 +235 58 83 55 84 +236 132 60 41 102 +237 66 86 123 101 +238 118 57 58 87 +239 33 83 58 84 +240 133 113 7 140 +241 140 89 113 92 +242 60 114 70 97 +243 125 132 22 102 +244 60 70 67 76 +245 56 70 71 83 +246 62 89 75 99 +247 63 73 68 88 +248 56 70 67 71 +249 9 91 123 101 +250 112 84 58 104 +251 60 70 80 102 +252 140 89 92 93 +253 124 5 112 148 +254 57 58 70 82 +255 57 80 120 82 +256 117 71 59 74 +257 32 70 114 97 +258 33 58 112 84 +259 60 85 41 102 +260 62 75 72 98 +261 136 59 72 74 +262 38 120 80 82 +263 51 59 78 105 +264 60 76 67 94 +265 70 76 60 80 +266 56 76 70 80 +267 35 22 121 85 +268 129 22 41 85 +269 107 65 31 103 +270 62 92 89 93 +271 109 30 29 100 +272 65 96 90 98 +273 112 124 148 104 +274 56 66 76 80 +275 67 70 56 76 +276 129 121 22 85 +277 120 85 80 102 +278 135 127 89 99 +279 107 131 65 103 +280 57 70 80 82 +281 54 142 73 78 +282 25 135 127 89 +283 120 35 85 102 +284 52 73 144 81 +285 133 79 113 89 +286 136 72 59 78 +287 16 55 115 84 +288 63 88 68 100 +289 32 71 70 83 +290 117 34 32 71 +291 109 88 30 100 +292 118 57 38 82 +293 52 53 144 73 +294 57 70 97 102 +295 57 80 70 102 +296 38 57 120 82 +297 59 74 136 78 +298 58 81 36 104 +299 32 70 33 83 +300 54 142 53 73 +301 118 81 58 82 +302 56 75 67 76 +303 72 75 56 98 +304 65 75 56 76 +305 71 72 59 74 +306 31 111 29 90 +307 32 33 70 97 +308 46 135 25 89 +309 25 127 40 89 +310 63 95 88 100 +311 111 90 31 96 +312 128 41 60 85 +313 112 148 84 104 +314 116 97 64 102 +315 147 139 54 68 +316 42 76 66 80 +317 68 90 111 96 +318 52 63 73 81 +319 43 75 76 94 +320 54 73 68 78 +321 25 133 140 89 +322 127 40 89 99 +323 52 63 53 73 +324 65 56 75 98 +325 124 19 148 104 +326 144 81 58 84 +327 33 115 32 83 +328 46 135 89 99 +329 144 50 81 84 +330 7 133 140 25 +331 138 96 65 98 +332 43 75 65 76 +333 54 53 68 73 +334 12 109 147 88 +335 148 124 5 19 +336 43 65 61 76 +337 15 92 74 93 +338 116 57 97 102 +339 48 72 96 98 +340 57 120 80 102 +341 133 40 79 89 +342 68 72 47 78 +343 34 71 117 92 +344 80 85 60 102 +345 91 95 63 100 +346 31 65 90 103 +347 76 80 42 85 +348 71 74 117 92 +349 6 22 125 132 +350 106 122 91 95 +351 43 76 128 94 +352 42 128 76 85 +353 67 76 75 94 +354 48 47 68 72 +355 22 121 129 21 +356 116 57 13 87 +357 39 122 63 82 +358 25 40 133 89 +359 140 46 25 89 +360 37 123 66 86 +361 37 66 123 80 +362 27 143 137 28 +363 145 19 18 119 +364 138 48 96 98 +365 136 47 72 78 +366 11 139 111 96 +367 10 108 61 103 +368 25 127 135 24 +369 61 65 131 103 +370 43 131 61 65 +371 63 66 91 100 +372 38 37 123 80 +373 54 4 139 147 +374 48 136 47 72 +375 28 143 137 105 +376 47 139 48 68 +377 145 19 119 104 +378 39 38 122 82 +379 10 61 131 103 +380 54 53 147 68 +381 39 122 1 95 +382 28 137 45 105 +383 1 146 52 95 +384 122 106 1 95 +385 39 63 122 95 +386 46 62 134 99 +387 28 51 143 105 +388 145 50 19 104 +389 146 1 106 95 +390 52 63 81 95 +391 39 81 63 95 +392 50 81 84 104 +393 28 149 51 105 +394 36 119 19 104 +395 36 19 124 104 +396 50 84 148 104 +397 129 37 121 85 +398 45 141 28 105 +399 63 81 39 82 +400 10 61 108 86 +401 2 123 37 86 +402 15 117 74 92 +403 134 62 46 93 +404 148 19 50 104 +405 42 130 2 86 +406 143 78 137 105 +407 134 72 48 98 +408 43 42 128 76 +409 134 62 72 98 +410 46 89 62 99 +411 145 119 81 104 +412 58 84 81 104 +413 37 42 66 80 +414 47 54 139 68 +415 137 78 45 105 +416 37 80 120 85 +417 120 37 38 80 +418 10 130 61 86 +419 35 37 120 85 +420 128 42 41 85 +421 51 78 143 105 +422 137 143 47 78 +423 145 81 50 104 +424 126 40 44 75 +425 118 38 39 82 +426 41 42 129 85 +427 121 37 35 85 +428 37 66 42 86 +429 52 145 119 81 +430 36 81 119 104 +431 118 39 81 82 +432 143 47 27 137 +433 18 52 145 119 +434 44 65 43 75 +435 129 37 21 121 +436 4 47 54 139 +437 138 65 49 98 +438 75 98 62 99 +439 62 98 134 99 +440 54 143 51 78 +441 49 135 46 99 +442 142 54 51 78 +443 37 42 80 85 +444 46 134 49 99 +445 44 131 43 65 +446 36 39 119 81 +447 47 54 68 78 +448 136 45 47 78 +449 39 36 118 81 +450 65 75 49 98 +451 137 47 45 78 +452 49 98 75 99 +453 52 144 50 81 +454 52 50 145 81 +455 138 49 48 98 +456 52 81 39 95 +457 134 48 49 98 +458 49 75 44 99 +459 49 65 44 75 +460 134 98 49 99 +461 37 129 42 85 +462 39 17 52 95 +463 143 54 47 78 +464 52 119 39 81 +465 52 18 39 119 +466 47 27 54 143 +467 37 129 21 42 +468 49 24 44 23 +469 17 52 18 39 +470 27 54 26 47 +471 42 21 37 20 +472 113 7 14 133 +473 124 112 5 13 +474 47 4 54 26 +475 5 16 112 148 +476 113 15 7 140 +477 14 132 6 114 +478 117 141 8 15 +479 16 142 51 55 +480 51 149 16 55 +481 53 142 16 55 +482 16 144 53 55 +483 13 120 38 57 +484 13 38 118 57 +485 10 130 43 61 +486 43 131 10 61 +487 116 6 13 64 +488 6 125 13 64 +489 8 141 117 69 +490 8 117 32 69 +491 32 115 8 69 +492 14 40 126 79 +493 14 133 40 79 +494 14 113 133 79 +495 117 32 69 83 +496 32 115 69 83 +497 50 144 16 84 +498 144 55 16 84 +499 16 148 50 84 +500 36 13 118 87 +501 13 57 118 87 +502 15 140 113 92 +503 15 113 34 92 +504 34 117 15 92 +505 15 134 46 93 +506 46 140 15 93 +507 140 92 15 93 +508 134 15 74 93 +509 43 128 14 94 +510 128 77 14 94 +511 14 126 43 94 +512 14 77 79 94 +513 14 79 126 94 +514 32 114 6 97 +515 6 116 32 97 +516 6 64 116 97 +517 32 116 33 97 +518 110 29 30 101 +519 35 120 13 102 +520 120 57 13 102 +521 13 125 35 102 +522 13 64 125 102 +523 31 10 107 103 +524 31 108 10 103 +525 139 111 147 4 +526 139 147 111 68 +527 131 107 44 3 +528 131 44 107 65 +529 44 107 23 3 +530 23 107 49 3 +531 49 107 138 3 +532 49 138 107 65 +533 87 97 70 33 +534 70 97 87 57 +535 70 58 87 33 +536 87 58 70 57 +537 89 92 67 113 +538 67 92 89 62 +539 82 122 91 38 +540 82 91 122 63 +541 91 30 9 106 +542 91 9 30 101 +543 99 40 44 127 +544 99 44 40 75 +545 66 108 61 86 +546 61 108 66 103 +547 135 127 49 24 +548 135 49 127 99 +549 44 49 127 24 +550 44 127 49 99 +551 8 69 105 141 +552 105 69 8 149 +553 105 28 8 141 +554 8 28 105 149 +555 41 102 22 132 +556 22 102 41 85 +557 94 60 128 77 +558 94 128 60 76 +559 102 116 13 64 +560 102 13 116 57 +561 93 89 46 140 +562 93 46 89 62 +563 91 106 100 30 +564 91 100 106 95 +565 100 106 88 30 +566 100 88 106 95 +567 88 146 106 95 +568 53 146 95 52 +569 95 146 53 88 +570 95 63 53 52 +571 53 63 95 88 +572 126 43 75 44 +573 126 75 43 94 +574 130 43 76 42 +575 76 43 130 61 +576 76 66 86 42 +577 86 66 76 61 +578 130 76 86 42 +579 86 76 130 61 +580 132 102 97 114 +581 97 102 132 64 +582 97 6 132 114 +583 132 6 97 64 +584 8 149 115 16 +585 8 115 149 69 +586 55 115 149 16 +587 55 149 115 69 +588 15 136 141 45 +589 15 141 136 74 +590 141 136 105 45 +591 105 136 78 45 +592 72 136 134 48 +593 72 134 136 74 +594 15 134 136 48 +595 15 136 134 74 +596 60 128 132 41 +597 60 132 128 77 +598 14 132 128 41 +599 14 128 132 77 +600 2 42 37 20 +601 37 42 2 86 +602 112 104 36 124 +603 36 104 112 87 +604 36 13 112 124 +605 112 13 36 87 +606 147 111 109 4 +607 29 109 111 4 +608 29 111 109 90 +609 9 110 123 2 +610 9 123 110 101 +611 123 110 86 2 +612 123 86 110 101 +613 86 110 108 2 +614 86 108 110 101 +615 108 110 29 2 +616 108 29 110 101 +617 68 109 111 90 +618 68 111 109 147 +619 44 107 49 23 +620 49 107 44 65 +621 106 12 88 30 +622 106 88 12 146 +623 141 136 78 105 +624 78 136 141 74 +625 141 78 59 105 +626 59 78 141 74 +627 141 59 69 105 +628 69 59 141 74 +3 1 7 44 +629 12 1 9 30 106 +630 10 3 11 31 107 +631 2 10 31 29 108 +632 4 12 30 29 109 +633 2 29 30 9 110 +634 4 29 31 11 111 +635 13 5 16 33 112 +636 15 7 14 34 113 +637 34 14 6 32 114 +638 33 16 8 32 115 +639 33 32 6 13 116 +640 34 32 8 15 117 +641 39 36 13 38 118 +642 19 36 39 18 119 +643 37 38 13 35 120 +644 37 35 22 21 121 +645 1 39 38 9 122 +646 2 9 38 37 123 +647 19 5 13 36 124 +648 22 35 13 6 125 +649 44 43 14 40 126 +650 25 24 44 40 127 +651 42 41 14 43 128 +652 42 21 22 41 129 +653 2 42 43 10 130 +654 3 10 43 44 131 +655 22 6 14 41 132 +656 25 40 14 7 133 +657 49 46 15 48 134 +658 25 46 49 24 135 +659 47 48 15 45 136 +660 47 45 28 27 137 +661 3 49 48 11 138 +662 4 11 48 47 139 +663 25 7 15 46 140 +664 28 45 15 8 141 +665 54 51 16 53 142 +666 28 51 54 27 143 +667 52 53 16 50 144 +668 52 50 19 18 145 +669 1 12 53 52 146 +670 4 54 53 12 147 +671 19 50 16 5 148 +672 28 8 16 51 149 +$EndElements diff --git a/examples/cube/mesh/cube2.geo b/examples/cube/mesh/cube2.geo new file mode 100644 index 000000000..a67feb62a --- /dev/null +++ b/examples/cube/mesh/cube2.geo @@ -0,0 +1,15 @@ +lc = 0.3; +Point(1) = {0.0,0.0,0.0,lc}; +Point(2) = {1,0.0,0.0,lc}; +Point(3) = {1,1,0.0,lc}; +Point(4) = {0,1,0.0,lc}; +Line(1) = {4,3}; +Line(2) = {3,2}; +Line(3) = {2,1}; +Line(4) = {1,4}; +Line Loop(5) = {2,3,4,1}; +Plane Surface(6) = {5}; +tmp[] = Extrude {0,0.0,1} { + Surface{6}; +}; +Physical Volume(1) = tmp[1]; diff --git a/examples/cube/mesh/cube2.msh b/examples/cube/mesh/cube2.msh new file mode 100644 index 000000000..30c91f753 --- /dev/null +++ b/examples/cube/mesh/cube2.msh @@ -0,0 +1,724 @@ +$MeshFormat +4.1 0 8 +$EndMeshFormat +$Entities +8 12 6 1 +1 0 0 0 0 +2 1 0 0 0 +3 1 1 0 0 +4 0 1 0 0 +5 1 1 1 0 +6 1 0 1 0 +10 0 0 1 0 +14 0 1 1 0 +1 0 1 0 1 1 0 0 2 4 -3 +2 1 0 0 1 1 0 0 2 3 -2 +3 0 0 0 1 0 0 0 2 2 -1 +4 0 0 0 0 1 0 0 2 1 -4 +8 1 0 1 1 1 1 0 2 5 -6 +9 0 0 1 1 0 1 0 2 6 -10 +10 0 0 1 0 1 1 0 2 10 -14 +11 0 1 1 1 1 1 0 2 14 -5 +13 1 1 0 1 1 1 0 2 3 -5 +14 1 0 0 1 0 1 0 2 2 -6 +18 0 0 0 0 0 1 0 2 1 -10 +22 0 1 0 0 1 1 0 2 4 -14 +6 0 0 0 1 1 0 0 4 2 3 4 1 +15 1 0 0 1 1 1 0 4 2 14 -8 -13 +19 0 0 0 1 0 1 0 4 3 18 -9 -14 +23 0 0 0 0 1 1 0 4 4 22 -10 -18 +27 0 1 0 1 1 1 0 4 1 13 -11 -22 +28 0 0 1 1 1 1 0 4 8 9 10 11 +1 0 0 0 1 1 1 1 1 6 -6 28 15 19 23 27 +$EndEntities +$Nodes +27 141 1 141 +0 1 0 1 +1 +0 0 0 +0 2 0 1 +2 +1 0 0 +0 3 0 1 +3 +1 1 0 +0 4 0 1 +4 +0 1 0 +0 5 0 1 +5 +1 1 1 +0 6 0 1 +6 +1 0 1 +0 10 0 1 +7 +0 0 1 +0 14 0 1 +8 +0 1 1 +1 1 0 3 +9 +10 +11 +0.2499999999994109 1 0 +0.4999999999986918 1 0 +0.7499999999993401 1 0 +1 2 0 3 +12 +13 +14 +1 0.7500000000003477 0 +1 0.5000000000020616 0 +1 0.2500000000010419 0 +1 3 0 3 +15 +16 +17 +0.7500000000003477 0 0 +0.5000000000020616 0 0 +0.2500000000010419 0 0 +1 4 0 3 +18 +19 +20 +0 0.2499999999994109 0 +0 0.4999999999986918 0 +0 0.7499999999993401 0 +1 8 0 3 +21 +22 +23 +1 0.7500000000003477 1 +1 0.5000000000020616 1 +1 0.2500000000010419 1 +1 9 0 3 +24 +25 +26 +0.7500000000003477 0 1 +0.5000000000020616 0 1 +0.2500000000010419 0 1 +1 10 0 3 +27 +28 +29 +0 0.2499999999994109 1 +0 0.4999999999986918 1 +0 0.7499999999993401 1 +1 11 0 3 +30 +31 +32 +0.2499999999994109 1 1 +0.4999999999986918 1 1 +0.7499999999993401 1 1 +1 13 0 3 +33 +34 +35 +1 1 0.2499999999994109 +1 1 0.4999999999986918 +1 1 0.7499999999993401 +1 14 0 3 +36 +37 +38 +1 0 0.2499999999994109 +1 0 0.4999999999986918 +1 0 0.7499999999993401 +1 18 0 3 +39 +40 +41 +0 0 0.2499999999994109 +0 0 0.4999999999986918 +0 0 0.7499999999993401 +1 22 0 3 +42 +43 +44 +0 1 0.2499999999994109 +0 1 0.4999999999986918 +0 1 0.7499999999993401 +2 6 0 14 +45 +46 +47 +48 +49 +50 +51 +52 +53 +54 +55 +56 +57 +58 +0.7867687832230401 0.3640932128839357 0 +0.3753337416957953 0.2110423218192964 0 +0.2544869979106266 0.614600373335793 0 +0.6243689539557368 0.8080509614599706 0 +0.6215671767908427 0.1979210047255828 0 +0.50565027269992 0.4308090314147043 0 +0.2167990069479854 0.3712961061742477 0 +0.7478815024664029 0.6139924015723902 0 +0.3855439771726532 0.8112741566735273 0 +0.5035863408410679 0.6557453848912771 0 +0.8169872981072784 0.1830127018930939 0 +0.8169872981071886 0.8169872981074586 0 +0.1830127018931349 0.1830127018926979 0 +0.1780061950165381 0.8351749060017322 0 +2 15 0 15 +59 +60 +61 +62 +63 +64 +65 +66 +67 +68 +69 +70 +71 +72 +73 +1 0.3829082618847883 0.7742952263263159 +1 0.384668784278834 0.2109240761334862 +1 0.8023234404950805 0.6408098273465964 +1 0.6245279258008938 0.2113539510935045 +1 0.5030682640978668 0.4047588473860917 +1 0.23045971405121 0.3942368581047184 +1 0.7764715641531658 0.3930291397001214 +1 0.6409766963415464 0.8053554665151262 +1 0.1815045070665119 0.6312063377284696 +1 0.3773688473460519 0.5680660677565404 +1 0.816987297978409 0.8169872979781391 +1 0.1830127020221846 0.1830127026658178 +1 0.8169872981083395 0.1830127025356663 +1 0.1815079327373217 0.8184920672627829 +1 0.5936411810720301 0.6000404470685039 +2 19 0 15 +74 +75 +76 +77 +78 +79 +80 +81 +82 +83 +84 +85 +86 +87 +88 +0.3829082625566119 0 0.7742952259365105 +0.3846687842254827 0 0.2109240760307497 +0.8023234409646326 0 0.640809826816817 +0.6245279258568575 0 0.2113539509042247 +0.5030682642536575 0 0.404758847061386 +0.2304597140658443 0 0.3942368580436552 +0.7764715650546368 0 0.393029138952583 +0.6409766971335645 0 0.8053554660456254 +0.1815045071248842 0 0.6312063376901305 +0.3773688481401301 0 0.5680660675612745 +0.1830127020221846 0 0.1830127026658178 +0.8169872982141551 0 0.8169872982138852 +0.8169872981083395 0 0.1830127025356663 +0.1815079327373218 0 0.8184920672627828 +0.5936411833492412 0 0.6000404460828872 +2 23 0 15 +89 +90 +91 +92 +93 +94 +95 +96 +97 +98 +99 +100 +101 +102 +103 +0 0.7742952261041663 0.6170917381413945 +0 0.210924076326279 0.6153312163442729 +0 0.6408098271182379 0.1976765594203343 +0 0.2113539518277431 0.3754720741933714 +0 0.4047588472500019 0.4969317363650173 +0 0.3942368581045613 0.7695402865545821 +0 0.3930291394809293 0.2235284353964096 +0 0.8053554662084675 0.3590233032895219 +0 0.6312063377100441 0.8184954930814853 +0 0.5680660676375914 0.6226311530688312 +0 0.1830127026656377 0.816987297978098 +0 0.816987297978098 0.1830127026656377 +0 0.183012701891818 0.183012701891818 +0 0.8184920672683385 0.8184920672683385 +0 0.6000404461074688 0.4063588179957743 +2 27 0 15 +104 +105 +106 +107 +108 +109 +110 +111 +112 +113 +114 +115 +116 +117 +118 +0.7742952261041663 1 0.6170917381413945 +0.210924076326279 1 0.6153312163442729 +0.6408098271182379 1 0.1976765594203343 +0.2113539518277431 1 0.3754720741933714 +0.4047588472500019 1 0.4969317363650173 +0.3942368581045613 1 0.7695402865545821 +0.3930291394809293 1 0.2235284353964096 +0.8053554662084675 1 0.3590233032895219 +0.6312063377100441 1 0.8184954930814853 +0.5680660676375914 1 0.6226311530688312 +0.1830127026656377 1 0.816987297978098 +0.816987297978098 1 0.1830127026656377 +0.183012701891818 1 0.183012701891818 +0.8184920672683385 1 0.8184920672683385 +0.6000404461074688 1 0.4063588179957743 +2 28 0 14 +119 +120 +121 +122 +123 +124 +125 +126 +127 +128 +129 +130 +131 +132 +0.7867687832230401 0.3640932128839357 1 +0.3753337416957953 0.2110423218192964 1 +0.2544869979106266 0.614600373335793 1 +0.6243689539557368 0.8080509614599706 1 +0.6215671767908427 0.1979210047255828 1 +0.50565027269992 0.4308090314147043 1 +0.2167990069479854 0.3712961061742477 1 +0.7478815024664029 0.6139924015723902 1 +0.3855439771726532 0.8112741566735273 1 +0.5035863408410679 0.6557453848912771 1 +0.8169872981072784 0.1830127018930939 1 +0.8169872981071886 0.8169872981074586 1 +0.1830127018931349 0.1830127018926979 1 +0.1780061950165381 0.8351749060017322 1 +3 1 0 9 +133 +134 +135 +136 +137 +138 +139 +140 +141 +0.5216464791328072 0.5233310092201628 0.514005813484948 +0.7079341732272582 0.3031667650223043 0.6732481463868246 +0.3031667653511403 0.2920107740972996 0.6732481460614643 +0.3290644531198246 0.3183526393001203 0.3206671843227681 +0.6859802748358628 0.3118179121057249 0.3143348718675962 +0.3071929996135236 0.6879595812927457 0.3187642236686014 +0.3082899826046207 0.6981438521670125 0.6856486853762523 +0.6732481463617611 0.7198757996728931 0.6968332351749583 +0.6871018922896833 0.7081239524571338 0.3150440478978938 +$EndNodes +$Elements +1 375 1 375 +3 1 4 375 +1 93 98 135 139 +2 73 63 68 134 +3 83 88 78 135 +4 93 135 133 139 +5 133 135 93 136 +6 133 138 93 139 +7 133 136 93 138 +8 108 118 113 140 +9 134 135 133 136 +10 134 136 133 137 +11 63 68 134 137 +12 77 49 75 137 +13 54 138 50 141 +14 81 74 123 134 +15 83 135 78 136 +16 120 123 74 135 +17 48 53 54 141 +18 73 133 63 134 +19 133 139 128 140 +20 123 134 74 135 +21 52 60 62 137 +22 48 110 53 141 +23 50 138 133 141 +24 46 75 49 136 +25 133 137 50 141 +26 105 89 107 138 +27 128 133 124 139 +28 50 136 133 138 +29 63 62 60 137 +30 124 133 128 140 +31 108 113 139 140 +32 63 134 133 137 +33 73 63 133 141 +34 133 136 50 137 +35 66 73 59 140 +36 49 136 75 137 +37 53 138 54 141 +38 107 89 96 138 +39 133 138 108 141 +40 124 123 120 135 +41 110 48 106 141 +42 133 135 124 139 +43 92 90 79 135 +44 52 137 62 141 +45 110 138 53 141 +46 73 134 59 140 +47 133 134 73 140 +48 103 93 136 138 +49 77 75 78 137 +50 98 97 94 139 +51 73 133 140 141 +52 98 94 135 139 +53 63 62 137 141 +54 133 134 124 135 +55 63 137 133 141 +56 108 139 133 140 +57 124 134 133 140 +58 125 135 94 139 +59 124 134 123 135 +60 104 65 61 140 +61 105 89 138 139 +62 118 108 138 141 +63 128 127 122 140 +64 108 138 133 139 +65 81 88 74 134 +66 76 67 80 134 +67 105 107 108 138 +68 90 82 79 135 +69 52 45 60 137 +70 79 92 135 136 +71 74 134 88 135 +72 93 90 92 135 +73 119 126 59 134 +74 67 64 80 134 +75 73 61 65 140 +76 104 65 140 141 +77 104 111 65 141 +78 124 126 134 140 +79 47 91 95 138 +80 128 139 127 140 +81 109 127 139 140 +82 75 136 78 137 +83 47 95 136 138 +84 73 140 65 141 +85 80 134 64 137 +86 92 93 135 136 +87 49 50 46 136 +88 95 51 47 136 +89 105 138 108 139 +90 64 134 68 137 +91 50 137 52 141 +92 103 136 95 138 +93 59 73 68 134 +94 103 89 98 139 +95 74 88 83 135 +96 91 103 95 138 +97 83 79 135 136 +98 88 134 80 137 +99 67 68 64 134 +100 118 138 110 141 +101 49 50 136 137 +102 82 83 79 135 +103 96 89 103 138 +104 54 50 52 141 +105 103 138 89 139 +106 124 126 119 134 +107 47 136 50 138 +108 124 128 126 140 +109 103 93 95 136 +110 121 128 124 139 +111 109 139 113 140 +112 73 65 63 141 +113 54 47 50 138 +114 118 104 140 141 +115 118 108 110 138 +116 58 53 116 138 +117 113 118 104 140 +118 118 111 104 141 +119 94 93 98 135 +120 106 118 110 141 +121 78 88 80 137 +122 92 84 101 136 +123 63 64 68 137 +124 88 76 80 134 +125 112 109 113 140 +126 58 100 91 138 +127 84 75 46 136 +128 57 101 84 136 +129 122 127 31 112 +130 84 92 79 136 +131 132 97 102 139 +132 10 110 53 48 +133 46 57 84 136 +134 127 132 114 139 +135 49 75 16 77 +136 83 78 79 136 +137 92 40 79 90 +138 45 52 50 137 +139 49 77 86 137 +140 82 79 40 90 +141 100 58 116 138 +142 46 75 16 49 +143 108 113 109 139 +144 62 52 13 60 +145 10 110 48 106 +146 110 116 53 138 +147 56 62 71 141 +148 123 120 74 25 +149 80 70 86 137 +150 49 86 55 137 +151 51 50 47 136 +152 131 120 87 135 +153 80 67 76 37 +154 70 55 86 137 +155 25 74 123 81 +156 34 61 104 65 +157 59 22 66 126 +158 132 102 114 139 +159 114 102 105 139 +160 129 81 123 134 +161 64 70 80 137 +162 129 85 81 134 +163 109 31 127 112 +164 57 51 101 136 +165 65 104 34 111 +166 127 114 109 139 +167 96 116 107 138 +168 106 48 56 141 +169 87 120 74 135 +170 72 85 129 134 +171 37 80 67 64 +172 87 99 131 135 +173 58 91 47 138 +174 45 70 60 137 +175 105 102 89 139 +176 69 66 130 140 +177 91 95 19 47 +178 116 96 100 138 +179 82 90 99 135 +180 97 28 94 125 +181 95 101 51 136 +182 121 97 132 139 +183 82 99 87 135 +184 92 95 93 136 +185 70 45 55 137 +186 56 71 115 141 +187 52 62 56 141 +188 131 99 125 135 +189 117 61 69 140 +190 56 115 106 141 +191 45 13 52 60 +192 22 59 119 126 +193 120 125 124 135 +194 112 130 122 140 +195 72 129 119 134 +196 117 69 130 140 +197 94 125 99 135 +198 124 119 123 134 +199 65 62 63 141 +200 63 60 64 137 +201 95 92 101 136 +202 111 115 71 141 +203 117 104 61 140 +204 119 59 72 134 +205 78 80 77 137 +206 130 112 117 140 +207 79 78 75 136 +208 108 107 110 138 +209 66 126 130 140 +210 65 111 71 141 +211 84 79 75 136 +212 90 93 94 135 +213 28 97 121 125 +214 58 47 53 138 +215 95 51 19 47 +216 53 47 54 138 +217 61 73 66 140 +218 64 60 70 137 +219 69 61 66 140 +220 105 108 109 139 +221 89 97 98 139 +222 80 86 77 137 +223 48 54 52 141 +224 127 128 121 139 +225 127 121 132 139 +226 49 45 50 137 +227 110 107 116 138 +228 50 51 46 136 +229 129 123 119 134 +230 46 51 57 136 +231 100 96 91 138 +232 65 71 62 141 +233 102 97 89 139 +234 125 120 131 135 +235 85 76 81 134 +236 94 99 90 135 +237 86 80 36 70 +238 114 44 105 102 +239 84 75 17 46 +240 49 55 45 137 +241 67 59 68 134 +242 83 82 74 135 +243 80 64 36 70 +244 44 89 105 102 +245 122 126 128 140 +246 109 114 105 139 +247 48 52 56 141 +248 17 57 84 46 +249 72 59 67 134 +250 58 53 9 116 +251 87 74 82 135 +252 76 88 81 134 +253 91 96 103 138 +254 87 131 26 120 +255 30 132 114 127 +256 15 49 77 86 +257 92 39 84 79 +258 39 84 101 92 +259 26 74 87 120 +260 111 106 115 141 +261 15 49 86 55 +262 122 130 126 140 +263 12 71 56 62 +264 48 56 11 106 +265 109 30 114 127 +266 53 110 9 116 +267 58 20 91 100 +268 60 45 14 70 +269 112 113 104 140 +270 122 112 32 130 +271 29 102 132 97 +272 117 112 104 140 +273 52 12 56 62 +274 99 82 41 90 +275 116 42 96 107 +276 106 111 118 141 +277 47 20 91 58 +278 67 76 38 85 +279 70 45 14 55 +280 101 57 18 51 +281 56 115 11 106 +282 94 125 27 99 +283 117 35 61 104 +284 18 95 101 51 +285 131 27 125 99 +286 72 67 38 85 +287 112 32 130 117 +288 87 99 82 41 +289 72 129 23 119 +290 119 23 72 59 +291 69 117 35 61 +292 132 121 29 97 +293 123 24 81 129 +294 85 81 24 129 +295 100 116 42 96 +296 39 92 40 79 +297 17 75 16 46 +298 111 33 65 71 +299 64 37 80 36 +300 89 105 43 44 +301 120 26 74 25 +302 33 111 115 71 +303 52 13 12 62 +304 49 16 15 77 +305 35 61 104 34 +306 21 66 126 130 +307 71 115 3 56 +308 2 86 70 55 +309 1 57 101 84 +310 30 109 31 127 +311 6 72 85 129 +312 69 130 5 117 +313 87 7 131 99 +314 100 58 4 116 +315 21 66 130 69 +316 102 114 8 132 +317 43 96 42 107 +318 110 9 10 53 +319 48 11 10 106 +320 41 82 40 90 +321 20 91 19 47 +322 33 65 34 111 +323 45 14 13 60 +324 22 21 66 126 +325 28 94 125 27 +326 59 119 23 22 +327 51 18 95 19 +328 121 29 97 28 +329 38 76 67 37 +330 123 24 25 81 +331 31 122 112 32 +332 11 3 115 56 +333 86 2 15 55 +334 85 24 6 129 +335 116 4 100 42 +336 1 57 84 17 +337 70 14 2 55 +338 72 38 6 85 +339 130 32 5 117 +340 99 87 7 41 +341 8 44 114 102 +342 129 6 72 23 +343 5 69 117 35 +344 87 131 7 26 +345 101 57 1 18 +346 3 33 115 71 +347 7 27 131 99 +348 86 36 2 70 +349 84 101 1 39 +350 58 9 4 116 +351 71 3 12 56 +352 130 5 21 69 +353 132 114 8 30 +354 58 4 20 100 +355 8 102 132 29 +356 139 103 93 98 +357 139 93 103 138 +358 137 88 136 78 +359 137 136 88 134 +360 135 136 88 78 +361 135 88 136 134 +362 141 108 140 118 +363 141 140 108 133 +364 59 140 126 66 +365 126 140 59 134 +366 112 127 140 122 +367 112 140 127 109 +368 139 125 124 121 +369 139 124 125 135 +370 107 89 43 96 +371 43 89 107 105 +372 125 97 139 94 +373 125 139 97 121 +374 67 85 134 76 +375 134 85 67 72 +$EndElements