diff --git a/dune/mmdg/mmdg.hh b/dune/mmdg/mmdg.hh index e881d5e45fce9c20c5f45b69dccb2b95fabc3e67..4d58fd3c12ffb5349624ddf9df763cdf209b4104 100644 --- a/dune/mmdg/mmdg.hh +++ b/dune/mmdg/mmdg.hh @@ -133,7 +133,7 @@ private: //in total we use the basis given by the union of // { (phi_elem,i , 0) | elem in elements(gridView), i = 0,...,dim } //and - // { (0 , phi_iElem,i ) | iElem in elements(iGridView), i = 0,...,dim-1 } + // { (0 , phi_iElem,i) | iElem in elements(iGridView), i = 0,...,dim-1 } //thus the stiffness matrix A and the load vector b exhibit the following //block structure diff --git a/dune/mmdg/problems/mmdgproblem4.hh b/dune/mmdg/problems/mmdgproblem4.hh index 50777e829fe6f3ccc05b2fae67eefde119346c4e..61d19dc60b91ff31607991fcdbfd3110b199175c 100644 --- a/dune/mmdg/problems/mmdgproblem4.hh +++ b/dune/mmdg/problems/mmdgproblem4.hh @@ -32,7 +32,6 @@ public: //the exact solution on the interface at position pos Scalar exactInterfaceSolution (const Coordinate& pos) const { - //return 1.5; return pos[0] * pos[0] - pos[1] * pos[1] + 0.5 * d_; } @@ -51,7 +50,7 @@ public: //interface source term at position pos Scalar qInterface (const Coordinate& pos) const { - return 2.0 * d_; //d^2 + return 2.0 * d_; } //aperture d of the fracture at position pos @@ -76,7 +75,6 @@ public: permeability[i][i] = 1.0 / d_; } - //return identity matrix return permeability; } diff --git a/dune/mmdg/problems/mmdgproblem5.hh b/dune/mmdg/problems/mmdgproblem5.hh new file mode 100644 index 0000000000000000000000000000000000000000..0a3ad3d11a558330a941ee54917bc3f671e50b37 --- /dev/null +++ b/dune/mmdg/problems/mmdgproblem5.hh @@ -0,0 +1,128 @@ +#ifndef __DUNE_MMDG_MMDGPROBLEM5_HH__ +#define __DUNE_MMDG_MMDGPROBLEM5_HH__ + +#include <cmath> +#include <dune/mmdg/problems/coupleddgproblem.hh> + +//a coupled dg problem +template<class Coordinate, class Scalar = double> +class MMDGProblem5 : public CoupledDGProblem<Coordinate, Scalar> +{ +public: + static constexpr int dim = Coordinate::dimension; + using Base = CoupledDGProblem<Coordinate, Scalar>; + using Matrix = typename Base::Matrix; + + //constructor + MMDGProblem5 (const Scalar d) : d_(d) {} + + //the exact bulk solution at position pos + Scalar exactSolution (const Coordinate& pos) const + { + return (pos[1] > pos[0]) ? + 1.0 / (1.0 + pos[0] * pos[0]) : -1.0 / (1.0 + pos[1] * pos[1]); + } + + //the exact solution on the interface at position pos + Scalar exactInterfaceSolution (const Coordinate& pos) const + { + return 0.0; + } + + //indicates whether an exact solution is implemented for the problem + bool hasExactSolution () const + { + return true; + }; + + //source term at position pos + Scalar q (const Coordinate& pos) const + { + const Scalar x2Plus1 = pos[0] * pos[0] + 1.0; + const Scalar y2Plus1 = pos[1] * pos[1] + 1.0; + + return (pos[1] > pos[0]) ? + -4.0 * pos[0] / (x2Plus1 * x2Plus1) : 4.0 * pos[1] / (y2Plus1 * y2Plus1); + } + + //interface source term at position pos + Scalar qInterface (const Coordinate& pos) const + { + return 0.0; + } + + //aperture d of the fracture at position pos + Scalar aperture (const Coordinate& pos) const + { + return d_; + } + + //permeability tensor at position pos + Matrix K (const Coordinate& pos) const + { + Matrix permeability(0.0); + + for (int i = 0; i < dim; i++) + { + permeability[i][i] = pos[i] + 1.0 / pos[i]; + } + + return permeability; + } + + //permeability of the fracture in normal direction at position pos + Scalar Kperp (const Coordinate& pos) const + { + constexpr Scalar Kvalue = 1.0 / sqrt(2.0); + return Kvalue; + } + + //tangential permeability tensor of the interface at position pos + Matrix Kparallel (const Coordinate& pos) const + { + Matrix permeability(0.0); + Scalar y2Plus1_2 = 1.0 + pos[1] * pos[1]; + y2Plus1_2 *= y2Plus1_2; + + for (int i = 0; i < dim; i++) + { + permeability[i][i] = y2Plus1_2; + } + + return permeability; + } + + //returns the recommended quadrature order to compute an integral + //over x * boundary(x) + int quadratureOrderBoundary () const + { + return 10; + } + + //returns the recommended quadrature order to compute an integral + //over x * interfaceBoundary(x) + int quadratureOrderInterfaceBoundary () const + { //NOTE: only relevant for 3d, not implemented! + return 10; + } + + //returns the recommended quadrature order to compute an integral + //over x * qInterface(x) + int quadratureOrder_qInterface () const + { + return 10; + } + + //returns the recommended quadrature order to compute an integral + //over x * q(x) + int quadratureOrder_q () const + { + return 10; + } + +private: + Scalar d_; //aperture + +}; + +#endif diff --git a/src/grids/mmdg5.geo b/src/grids/mmdg5.geo new file mode 100644 index 0000000000000000000000000000000000000000..d13f9ab8cd101a5dd0289069232482975d1587ae --- /dev/null +++ b/src/grids/mmdg5.geo @@ -0,0 +1,27 @@ +//characteristic length +lc = 8e-2; + +//domain corners, domain = [0,1] x [0,1] +Point(1) = {1, 1, 0, lc}; +Point(2) = {2, 1, 0, lc}; +Point(3) = {2, 2, 0, lc}; +Point(4) = {1, 2, 0, lc}; + +//domain outline +Line(1) = {1, 2}; +Line(2) = {2, 3}; +Line(3) = {3, 4}; +Line(4) = {4, 1}; + +//interface outline +Line(5) = {1,3}; + +//curve loops +Curve Loop(1) = {1:4}; //domain boundary + +//surfaces +Plane Surface(1) = {1}; +Physical Surface(1) = {1}; + +Curve{5} In Surface{1}; +Physical Curve(2) = {5}; diff --git a/src/grids/mmdg5_1.msh b/src/grids/mmdg5_1.msh new file mode 100644 index 0000000000000000000000000000000000000000..f66babb60a94f99b0e3742f6338a8c4c320cc628 --- /dev/null +++ b/src/grids/mmdg5_1.msh @@ -0,0 +1,20 @@ +$MeshFormat +2.2 0 8 +$EndMeshFormat +$Nodes +5 +1 1 1 0 +2 2 1 0 +3 2 2 0 +4 1 2 0 +5 1.500000000000001 1.500000000000001 0 +$EndNodes +$Elements +6 +1 1 2 2 5 1 5 +2 1 2 2 5 5 3 +3 2 2 1 1 2 5 1 +4 2 2 1 1 1 5 4 +5 2 2 1 1 2 3 5 +6 2 2 1 1 5 3 4 +$EndElements diff --git a/src/grids/mmdg5_2e-1.msh b/src/grids/mmdg5_2e-1.msh new file mode 100644 index 0000000000000000000000000000000000000000..db679d1d2cc570bf4c336583d7a6dded54ed9359 --- /dev/null +++ b/src/grids/mmdg5_2e-1.msh @@ -0,0 +1,148 @@ +$MeshFormat +2.2 0 8 +$EndMeshFormat +$Nodes +51 +1 1 1 0 +2 2 1 0 +3 2 2 0 +4 1 2 0 +5 1.2 1 0 +6 1.4 1 0 +7 1.6 1 0 +8 1.8 1 0 +9 2 1.2 0 +10 2 1.4 0 +11 2 1.6 0 +12 2 1.8 0 +13 1.8 2 0 +14 1.6 2 0 +15 1.4 2 0 +16 1.2 2 0 +17 1 1.8 0 +18 1 1.6 0 +19 1 1.4 0 +20 1 1.2 0 +21 1.125 1.125 0 +22 1.25 1.25 0 +23 1.375 1.375 0 +24 1.5 1.5 0 +25 1.625 1.625 0 +26 1.75 1.75 0 +27 1.875 1.875 0 +28 1.703849158412359 1.293938654098996 0 +29 1.293938654098996 1.703849158412359 0 +30 1.174024904605511 1.509339147710085 0 +31 1.509339147710085 1.174024904605511 0 +32 1.487248935408783 1.818081513491189 0 +33 1.818081513491189 1.487248935408783 0 +34 1.479200251300387 1.654652792268168 0 +35 1.654652792268168 1.479200251300387 0 +36 1.832916666666667 1.167083333333333 0 +37 1.167083333333334 1.832916666666667 0 +38 1.878194428793538 1.693976527745575 0 +39 1.121134348255501 1.305733343267934 0 +40 1.693976527745575 1.878194428793538 0 +41 1.305733343267933 1.121134348255501 0 +42 1.347973506156167 1.521924870408076 0 +43 1.521924870408076 1.347973506156168 0 +44 1.870969467714043 1.309654184568222 0 +45 1.689220994557822 1.127009378407568 0 +46 1.127009378407568 1.689220994557822 0 +47 1.309654184568223 1.870969467714043 0 +48 1.241915030879185 1.383084969120815 0 +49 1.383084969120815 1.241915030879185 0 +50 1.616915030879185 1.758084969120815 0 +51 1.758084969120814 1.616915030879186 0 +$EndNodes +$Elements +88 +1 1 2 2 5 1 21 +2 1 2 2 5 21 22 +3 1 2 2 5 22 23 +4 1 2 2 5 23 24 +5 1 2 2 5 24 25 +6 1 2 2 5 25 26 +7 1 2 2 5 26 27 +8 1 2 2 5 27 3 +9 2 2 1 1 1 21 20 +10 2 2 1 1 12 3 27 +11 2 2 1 1 5 21 1 +12 2 2 1 1 27 3 13 +13 2 2 1 1 41 31 49 +14 2 2 1 1 40 32 50 +15 2 2 1 1 30 39 48 +16 2 2 1 1 33 38 51 +17 2 2 1 1 6 7 31 +18 2 2 1 1 18 19 30 +19 2 2 1 1 10 11 33 +20 2 2 1 1 14 15 32 +21 2 2 1 1 6 31 41 +22 2 2 1 1 14 32 40 +23 2 2 1 1 33 11 38 +24 2 2 1 1 30 19 39 +25 2 2 1 1 29 30 42 +26 2 2 1 1 31 28 43 +27 2 2 1 1 2 9 36 +28 2 2 1 1 4 17 37 +29 2 2 1 1 16 4 37 +30 2 2 1 1 8 2 36 +31 2 2 1 1 25 35 51 +32 2 2 1 1 34 25 50 +33 2 2 1 1 28 33 35 +34 2 2 1 1 32 29 34 +35 2 2 1 1 20 21 39 +36 2 2 1 1 27 13 40 +37 2 2 1 1 12 27 38 +38 2 2 1 1 21 5 41 +39 2 2 1 1 21 22 39 +40 2 2 1 1 22 21 41 +41 2 2 1 1 27 26 38 +42 2 2 1 1 26 27 40 +43 2 2 1 1 30 29 46 +44 2 2 1 1 33 28 44 +45 2 2 1 1 28 31 45 +46 2 2 1 1 29 32 47 +47 2 2 1 1 35 33 51 +48 2 2 1 1 32 34 50 +49 2 2 1 1 24 25 34 +50 2 2 1 1 25 24 35 +51 2 2 1 1 13 14 40 +52 2 2 1 1 5 6 41 +53 2 2 1 1 11 12 38 +54 2 2 1 1 19 20 39 +55 2 2 1 1 18 30 46 +56 2 2 1 1 31 7 45 +57 2 2 1 1 10 33 44 +58 2 2 1 1 32 15 47 +59 2 2 1 1 17 18 46 +60 2 2 1 1 9 10 44 +61 2 2 1 1 7 8 45 +62 2 2 1 1 15 16 47 +63 2 2 1 1 37 17 46 +64 2 2 1 1 8 36 45 +65 2 2 1 1 36 9 44 +66 2 2 1 1 16 37 47 +67 2 2 1 1 28 35 43 +68 2 2 1 1 34 29 42 +69 2 2 1 1 23 24 42 +70 2 2 1 1 24 23 43 +71 2 2 1 1 36 28 45 +72 2 2 1 1 29 37 46 +73 2 2 1 1 28 36 44 +74 2 2 1 1 37 29 47 +75 2 2 1 1 23 42 48 +76 2 2 1 1 43 23 49 +77 2 2 1 1 22 23 48 +78 2 2 1 1 23 22 49 +79 2 2 1 1 25 26 50 +80 2 2 1 1 26 25 51 +81 2 2 1 1 31 43 49 +82 2 2 1 1 42 30 48 +83 2 2 1 1 24 34 42 +84 2 2 1 1 35 24 43 +85 2 2 1 1 22 41 49 +86 2 2 1 1 39 22 48 +87 2 2 1 1 26 40 50 +88 2 2 1 1 38 26 51 +$EndElements diff --git a/src/grids/mmdg5_5e-1.msh b/src/grids/mmdg5_5e-1.msh new file mode 100644 index 0000000000000000000000000000000000000000..f5c5415e77a1946b7529c3e665534bc12d7ba246 --- /dev/null +++ b/src/grids/mmdg5_5e-1.msh @@ -0,0 +1,38 @@ +$MeshFormat +2.2 0 8 +$EndMeshFormat +$Nodes +12 +1 1 1 0 +2 2 1 0 +3 2 2 0 +4 1 2 0 +5 1.5 1 0 +6 2 1.5 0 +7 1.5 2 0 +8 1 1.5 0 +9 1.333333333333334 1.333333333333334 0 +10 1.666666666666669 1.666666666666669 0 +11 1.7 1.3 0 +12 1.3 1.7 0 +$EndNodes +$Elements +17 +1 1 2 2 5 1 9 +2 1 2 2 5 9 10 +3 1 2 2 5 10 3 +4 2 2 1 1 10 3 7 +5 2 2 1 1 6 3 10 +6 2 2 1 1 1 9 8 +7 2 2 1 1 5 9 1 +8 2 2 1 1 10 9 11 +9 2 2 1 1 9 10 12 +10 2 2 1 1 5 2 11 +11 2 2 1 1 4 8 12 +12 2 2 1 1 2 6 11 +13 2 2 1 1 7 4 12 +14 2 2 1 1 8 9 12 +15 2 2 1 1 9 5 11 +16 2 2 1 1 10 7 12 +17 2 2 1 1 6 10 11 +$EndElements diff --git a/src/grids/mmdg5_5e-2.msh b/src/grids/mmdg5_5e-2.msh new file mode 100644 index 0000000000000000000000000000000000000000..ef0f0c2641b6f8bd2c121e9a6b1a0f8eb48268c6 --- /dev/null +++ b/src/grids/mmdg5_5e-2.msh @@ -0,0 +1,1630 @@ +$MeshFormat +2.2 0 8 +$EndMeshFormat +$Nodes +558 +1 1 1 0 +2 2 1 0 +3 2 2 0 +4 1 2 0 +5 1.05 1 0 +6 1.1 1 0 +7 1.15 1 0 +8 1.2 1 0 +9 1.25 1 0 +10 1.3 1 0 +11 1.35 1 0 +12 1.4 1 0 +13 1.45 1 0 +14 1.5 1 0 +15 1.55 1 0 +16 1.6 1 0 +17 1.65 1 0 +18 1.7 1 0 +19 1.75 1 0 +20 1.8 1 0 +21 1.85 1 0 +22 1.9 1 0 +23 1.95 1 0 +24 2 1.05 0 +25 2 1.1 0 +26 2 1.15 0 +27 2 1.2 0 +28 2 1.25 0 +29 2 1.3 0 +30 2 1.35 0 +31 2 1.4 0 +32 2 1.45 0 +33 2 1.5 0 +34 2 1.55 0 +35 2 1.6 0 +36 2 1.65 0 +37 2 1.7 0 +38 2 1.75 0 +39 2 1.8 0 +40 2 1.85 0 +41 2 1.9 0 +42 2 1.95 0 +43 1.95 2 0 +44 1.9 2 0 +45 1.85 2 0 +46 1.8 2 0 +47 1.75 2 0 +48 1.7 2 0 +49 1.65 2 0 +50 1.6 2 0 +51 1.55 2 0 +52 1.5 2 0 +53 1.45 2 0 +54 1.4 2 0 +55 1.35 2 0 +56 1.3 2 0 +57 1.25 2 0 +58 1.2 2 0 +59 1.15 2 0 +60 1.1 2 0 +61 1.05 2 0 +62 1 1.95 0 +63 1 1.9 0 +64 1 1.85 0 +65 1 1.8 0 +66 1 1.75 0 +67 1 1.7 0 +68 1 1.65 0 +69 1 1.6 0 +70 1 1.55 0 +71 1 1.5 0 +72 1 1.45 0 +73 1 1.4 0 +74 1 1.35 0 +75 1 1.3 0 +76 1 1.25 0 +77 1 1.2 0 +78 1 1.15 0 +79 1 1.1 0 +80 1 1.05 0 +81 1.03448275862069 1.03448275862069 0 +82 1.068965517241379 1.068965517241379 0 +83 1.103448275862069 1.103448275862069 0 +84 1.137931034482759 1.137931034482759 0 +85 1.172413793103448 1.172413793103448 0 +86 1.206896551724138 1.206896551724138 0 +87 1.241379310344828 1.241379310344828 0 +88 1.275862068965517 1.275862068965517 0 +89 1.310344827586207 1.310344827586207 0 +90 1.344827586206897 1.344827586206897 0 +91 1.379310344827586 1.379310344827586 0 +92 1.413793103448276 1.413793103448276 0 +93 1.448275862068966 1.448275862068966 0 +94 1.482758620689655 1.482758620689655 0 +95 1.517241379310345 1.517241379310345 0 +96 1.551724137931034 1.551724137931034 0 +97 1.586206896551724 1.586206896551724 0 +98 1.620689655172414 1.620689655172414 0 +99 1.655172413793104 1.655172413793104 0 +100 1.689655172413793 1.689655172413793 0 +101 1.724137931034483 1.724137931034483 0 +102 1.758620689655173 1.758620689655173 0 +103 1.793103448275862 1.793103448275862 0 +104 1.827586206896552 1.827586206896552 0 +105 1.862068965517242 1.862068965517242 0 +106 1.896551724137931 1.896551724137931 0 +107 1.931034482758621 1.931034482758621 0 +108 1.96551724137931 1.96551724137931 0 +109 1.293908159595183 1.707173976843308 0 +110 1.707173976843305 1.293908159595184 0 +111 1.488455395827184 1.785564290506268 0 +112 1.785564290506269 1.488455395827183 0 +113 1.217626737629361 1.513527943507945 0 +114 1.513527943507946 1.217626737629362 0 +115 1.832216434209611 1.16742030246278 0 +116 1.167420302462781 1.832216434209612 0 +117 1.152978619675205 1.368922049152326 0 +118 1.368922049152325 1.152978619675205 0 +119 1.847021380324795 1.631077950847673 0 +120 1.631077950847674 1.847021380324795 0 +121 1.427830058747408 1.651050646667107 0 +122 1.651050646667106 1.427830058747406 0 +123 1.344526687099 1.848147888571915 0 +124 1.848147888571915 1.344526687098998 0 +125 1.649968721817871 1.152131359244414 0 +126 1.152131359244414 1.64996872181787 0 +127 1.534380217637535 1.347146478604053 0 +128 1.347146478604053 1.534380217637533 0 +129 1.889352113273405 1.734479791948566 0 +130 1.110647886726593 1.265520208051434 0 +131 1.734479791948567 1.889352113273407 0 +132 1.265520208051433 1.110647886726593 0 +133 1.10927749111436 1.471326065082704 0 +134 1.528673934917296 1.890722508885641 0 +135 1.890722508885641 1.528673934917296 0 +136 1.471326065082703 1.109277491114359 0 +137 1.744164663543895 1.588463747885493 0 +138 1.588463747885494 1.744164663543894 0 +139 1.255440469871814 1.41331218277418 0 +140 1.41331218277418 1.255440469871815 0 +141 1.266977077412053 1.602231586726453 0 +142 1.602231586726453 1.266977077412054 0 +143 1.39561085927885 1.739168704737854 0 +144 1.739168704737852 1.395610859278847 0 +145 1.747669830000197 1.117793013166533 0 +146 1.117793013166535 1.747669830000198 0 +147 1.256902409288907 1.893259276413316 0 +148 1.893259276413317 1.256902409288907 0 +149 1.531794640854234 1.665223951361509 0 +150 1.665223951361508 1.531794640854234 0 +151 1.901758424048178 1.095293993477543 0 +152 1.095293993477545 1.901758424048178 0 +153 1.57221084811197 1.09566284678802 0 +154 1.904337153211979 1.427789151888031 0 +155 1.427789151888031 1.904337153211979 0 +156 1.095662846788021 1.572210848111967 0 +157 1.798078527818303 1.262616518734937 0 +158 1.262616518734936 1.798078527818304 0 +159 1.202921353097178 1.738564609322744 0 +160 1.738564609322743 1.202921353097178 0 +161 1.558047623785612 1.432895772141976 0 +162 1.432895772141978 1.558047623785612 0 +163 1.450021433783818 1.331171109818391 0 +164 1.33117110981839 1.450021433783816 0 +165 1.344872135603886 1.632644584772506 0 +166 1.632644584772506 1.344872135603887 0 +167 1.178860045828999 1.922403263204486 0 +168 1.922403263204488 1.178860045828998 0 +169 1.82160731126862 1.078212944557763 0 +170 1.078212944557763 1.82160731126862 0 +171 1.921972054863461 1.806981562894044 0 +172 1.806981562894045 1.921972054863463 0 +173 1.078027945136536 1.193018437105953 0 +174 1.193018437105952 1.078027945136537 0 +175 1.816144564361443 1.70408394395173 0 +176 1.704083943951732 1.816144564361444 0 +177 1.183855435638558 1.295916056048269 0 +178 1.295916056048268 1.183855435638558 0 +179 1.92095817130086 1.661175261113056 0 +180 1.07904182869914 1.338824738886945 0 +181 1.338824738886943 1.079041828699139 0 +182 1.661175261113058 1.92095817130086 0 +183 1.669324160777568 1.077838824192506 0 +184 1.077838824192508 1.669324160777567 0 +185 1.922320603372591 1.329986432673865 0 +186 1.329986432673867 1.922320603372589 0 +187 1.439284657412375 1.180653459139306 0 +188 1.815461348980987 1.560748569108928 0 +189 1.180653459139305 1.439284657412375 0 +190 1.560748569108929 1.815461348980985 0 +191 1.468186635345498 1.713193473236357 0 +192 1.71319347323636 1.468186635345496 0 +193 1.576046801042956 1.173000304366177 0 +194 1.423385603130812 1.823354278486964 0 +195 1.173000304366173 1.576046801042952 0 +196 1.823354278486964 1.423385603130812 0 +197 1.494843118518314 1.600520529428728 0 +198 1.600520529428725 1.494843118518312 0 +199 1.661275880337027 1.758625191567975 0 +200 1.758625191567976 1.661275880337028 0 +201 1.241374808432028 1.338724119662973 0 +202 1.338724119662972 1.241374808432026 0 +203 1.589681062443401 1.928819013676577 0 +204 1.928819013676576 1.589681062443401 0 +205 1.071180986323421 1.410318937556601 0 +206 1.410318937556599 1.071180986323423 0 +207 1.663799194525358 1.222357823091753 0 +208 1.222357823091753 1.663799194525361 0 +209 1.282779661788947 1.506035911227879 0 +210 1.506035911227883 1.282779661788948 0 +211 1.331106439824016 1.779463743393781 0 +212 1.779463743393781 1.331106439824014 0 +213 1.387430681930641 1.491712132176408 0 +214 1.491712132176414 1.387430681930643 0 +215 1.679534152670789 1.592976981641047 0 +216 1.592976981641045 1.679534152670789 0 +217 1.737521293460635 1.531953915970157 0 +218 1.531953915970159 1.737521293460633 0 +219 1.482480801017248 1.93222260500507 0 +220 1.51751919898275 1.06777739499493 0 +221 1.93222260500507 1.48248080101725 0 +222 1.06777739499493 1.517519198982751 0 +223 1.521560816990244 1.152704646165468 0 +224 1.850670053686097 1.487989331307737 0 +225 1.487989331307736 1.850670053686097 0 +226 1.152704646165467 1.521560816990241 0 +227 1.383939017789249 1.68314115744443 0 +228 1.683141157444428 1.383939017789249 0 +229 1.383593072003243 1.588607872197797 0 +230 1.588607872197797 1.383593072003243 0 +231 1.398224294654927 1.302052802765113 0 +232 1.302052802765113 1.398224294654927 0 +233 1.621734476217419 1.052059792376795 0 +234 1.378167037132788 1.947962982989646 0 +235 1.947962982989647 1.378167037132786 0 +236 1.052059792376796 1.621734476217417 0 +237 1.058510668753954 1.941489331246046 0 +238 1.941489331246047 1.058510668753954 0 +239 1.94717104962464 1.86008692657443 0 +240 1.860086926574428 1.947171049624639 0 +241 1.052828950375359 1.13991307342557 0 +242 1.139913073425569 1.052828950375358 0 +243 1.955960086598579 1.725576335326328 0 +244 1.044039913401418 1.274423664673669 0 +245 1.274423664673669 1.044039913401419 0 +246 1.725576335326332 1.955960086598582 0 +247 1.854313891738678 1.286684577932552 0 +248 1.286684577932552 1.854313891738678 0 +249 1.152840574410203 1.226469770417378 0 +250 1.773530229582618 1.847159425589794 0 +251 1.226469770417378 1.152840574410204 0 +252 1.847159425589794 1.773530229582619 0 +253 1.774916415941419 1.163716109846591 0 +254 1.163716109846592 1.77491641594142 0 +255 1.850390915933954 1.227047159223457 0 +256 1.227047159223457 1.850390915933954 0 +257 1.726050897194427 1.057603970886109 0 +258 1.274260991338987 1.943996049751356 0 +259 1.05760397088611 1.726050897194427 0 +260 1.943996049751358 1.274260991338987 0 +261 1.316519256191738 1.587129200289769 0 +262 1.587129200289769 1.316519256191739 0 +263 1.463156122960439 1.257885718764508 0 +264 1.257885718764506 1.46315612296044 0 +265 1.143993314916913 1.70445639124501 0 +266 1.704456391245009 1.143993314916911 0 +267 1.682918119422459 1.867846603808816 0 +268 1.867846603808812 1.682918119422458 0 +269 1.132153396191188 1.317081880577542 0 +270 1.31708188057754 1.132153396191188 0 +271 1.791923333687461 1.615024284851209 0 +272 1.211271239038293 1.381474077500972 0 +273 1.615024284851209 1.79192333368746 0 +274 1.381474077500972 1.211271239038293 0 +275 1.126522679120447 1.418158703134009 0 +276 1.581841296865991 1.873477320879554 0 +277 1.873477320879551 1.581841296865989 0 +278 1.418158703134008 1.126522679120446 0 +279 1.130596982688577 1.948313149839124 0 +280 1.047627326801329 1.866278966607692 0 +281 1.948313149839125 1.130596982688574 0 +282 1.866278966607691 1.047627326801329 0 +283 1.250831636230523 1.54580269342243 0 +284 1.545802693422433 1.250831636230524 0 +285 1.252300344186508 1.747321089275218 0 +286 1.747321089275215 1.25230034418651 0 +287 1.292989157390618 1.654828053934523 0 +288 1.654828053934521 1.292989157390619 0 +289 1.641084080272031 1.704710073598156 0 +290 1.704710073598156 1.641084080272033 0 +291 1.789852752494263 1.210819078576449 0 +292 1.210819078576449 1.789852752494265 0 +293 1.440382996630123 1.777736287780227 0 +294 1.777736287780228 1.440382996630121 0 +295 1.532917556203867 1.949566529957859 0 +296 1.949566529957859 1.532917556203866 0 +297 1.467082443796133 1.050433470042141 0 +298 1.05043347004214 1.467082443796134 0 +299 1.718758120831279 1.34261142250016 0 +300 1.34261142250016 1.718758120831282 0 +301 1.851165582159955 1.12274110485503 0 +302 1.122741104855031 1.851165582159954 0 +303 1.383880351066585 1.790795968309861 0 +304 1.790795968309861 1.383880351066584 0 +305 1.957078395644843 1.21597525319468 0 +306 1.784053829407411 1.043234600023431 0 +307 1.043234600023432 1.784053829407412 0 +308 1.215975253194681 1.957078395644841 0 +309 1.611520715429325 1.541790909509032 0 +310 1.541790909509032 1.611520715429328 0 +311 1.210899972333734 1.612486719203788 0 +312 1.612486719203787 1.21089997233374 0 +313 1.892336492868602 1.365414493225321 0 +314 1.632977848588321 1.107578726096997 0 +315 1.107578726096998 1.63297784858832 0 +316 1.365414493225322 1.8923364928686 0 +317 1.432726994044488 1.498310540962406 0 +318 1.498310540962409 1.432726994044489 0 +319 1.482415352578478 1.647381776239105 0 +320 1.647381776239104 1.482415352578477 0 +321 1.482099238168262 1.555629421956253 0 +322 1.555629421956248 1.482099238168258 0 +323 1.346488212676809 1.287747404317843 0 +324 1.287747404317843 1.34648821267681 0 +325 1.223440791704046 1.039364437246172 0 +326 1.960635562753827 1.776559208295951 0 +327 1.039364437246171 1.223440791704047 0 +328 1.776559208295952 1.960635562753828 0 +329 1.564977937623382 1.042264984860064 0 +330 1.435010337775454 1.957737726493082 0 +331 1.957737726493083 1.435010337775452 0 +332 1.042264984860066 1.56497793762338 0 +333 1.876060901468873 1.188596230367579 0 +334 1.18859623036758 1.876060901468873 0 +335 1.873851640366673 1.814067670518256 0 +336 1.814067670518257 1.873851640366676 0 +337 1.185932329481741 1.126148359633325 0 +338 1.126148359633324 1.185932329481741 0 +339 1.70650908689785 1.771193919220809 0 +340 1.77119391922081 1.70650908689785 0 +341 1.22880608077919 1.293490913102149 0 +342 1.293490913102149 1.22880608077919 0 +343 1.954982065583369 1.626916237655493 0 +344 1.045017934416629 1.373083762344508 0 +345 1.626916237655493 1.95498206558337 0 +346 1.373083762344508 1.045017934416629 0 +347 1.453368401366818 1.602644966351887 0 +348 1.602644966351886 1.453368401366815 0 +349 1.447709298902069 1.376019868176358 0 +350 1.376019868176356 1.447709298902065 0 +351 1.388906003544225 1.540965743652531 0 +352 1.540965743652532 1.388906003544225 0 +353 1.123048267714089 1.801540056926871 0 +354 1.80154005692687 1.123048267714088 0 +355 1.80728270068331 1.748171110356596 0 +356 1.192717299316683 1.251828889643401 0 +357 1.251828889643402 1.192717299316683 0 +358 1.748171110356598 1.80728270068331 0 +359 1.892978254917567 1.62404558628459 0 +360 1.624045586284592 1.892978254917566 0 +361 1.375954413715411 1.107021745082432 0 +362 1.107021745082436 1.375954413715407 0 +363 1.28876346446621 1.440827597536062 0 +364 1.440827597536062 1.288763464466211 0 +365 1.427033671358538 1.702723243840398 0 +366 1.702723243840396 1.427033671358536 0 +367 1.928917318253906 1.968735771625016 0 +368 1.071082681746093 1.031264228374981 0 +369 1.968735771625017 1.928917318253907 0 +370 1.031264228374982 1.071082681746093 0 +371 1.32528239691539 1.49824431751533 0 +372 1.498244317515335 1.325282396915392 0 +373 1.659103356681088 1.811538323624262 0 +374 1.188461676375737 1.340896643318912 0 +375 1.340896643318911 1.188461676375738 0 +376 1.811538323624263 1.659103356681089 0 +377 1.145569289979724 1.890301064054547 0 +378 1.890301064054546 1.145569289979722 0 +379 1.889796493972676 1.474480139090939 0 +380 1.474480139090938 1.889796493972676 0 +381 1.108659153470157 1.525773024679662 0 +382 1.525773024679661 1.108659153470157 0 +383 1.328804298634536 1.677636101492689 0 +384 1.677636101492687 1.328804298634538 0 +385 1.390639753204705 1.632250043167761 0 +386 1.63225004316776 1.390639753204704 0 +387 1.472566735639878 1.15427156552039 0 +388 1.527433264360125 1.845728434479609 0 +389 1.845728434479614 1.527433264360125 0 +390 1.154271565520387 1.472566735639874 0 +391 1.771154399540027 1.555094474022461 0 +392 1.555094474022465 1.771154399540027 0 +393 1.303190069487368 1.756927936776172 0 +394 1.75692793677617 1.303190069487368 0 +395 1.693980293886066 1.245334199080409 0 +396 1.245334199080409 1.693980293886068 0 +397 1.296956719087691 1.547009073778456 0 +398 1.54700907377846 1.296956719087688 0 +399 1.040807436706695 1.32319904926385 0 +400 1.959192563293303 1.67680095073615 0 +401 1.676800950736152 1.959192563293305 0 +402 1.323199049263849 1.040807436706695 0 +403 1.705322533766711 1.558267899748807 0 +404 1.558267899748806 1.705322533766708 0 +405 1.964855892163695 1.819683242050706 0 +406 1.180316757949292 1.035144107836304 0 +407 1.035144107836305 1.180316757949293 0 +408 1.819683242050706 1.964855892163694 0 +409 1.883729107123577 1.323043973814478 0 +410 1.32304397381448 1.883729107123576 0 +411 1.90032308690199 1.84795136992838 0 +412 1.152048630071618 1.099676913098009 0 +413 1.099676913098009 1.152048630071618 0 +414 1.84795136992838 1.90032308690199 0 +415 1.636557301413557 1.582192119233065 0 +416 1.582192119233064 1.636557301413558 0 +417 1.690368611584417 1.188535331956359 0 +418 1.18853533195636 1.690368611584419 0 +419 1.175022505295306 1.965008989764131 0 +420 1.824840696963958 1.034828040643561 0 +421 1.965008989764133 1.175022505295306 0 +422 1.034828040643561 1.824840696963961 0 +423 1.203511361452224 1.543357925355743 0 +424 1.543357925355747 1.203511361452227 0 +425 1.108777963959539 1.223244178359285 0 +426 1.223244178359285 1.10877796395954 0 +427 1.891222036040458 1.776755821640712 0 +428 1.776755821640713 1.89122203604046 0 +429 1.152401618561438 1.268314603885517 0 +430 1.268314603885518 1.152401618561437 0 +431 1.847598381438564 1.731685396114476 0 +432 1.731685396114478 1.847598381438568 0 +433 1.115210256043632 1.672788178480366 0 +434 1.672788178480366 1.11521025604363 0 +435 1.674194850427491 1.035779752040723 0 +436 1.325540412619945 1.96486228021083 0 +437 1.035779752040722 1.674194850427492 0 +438 1.964862280210831 1.325540412619944 0 +439 1.857921811274563 1.444140760829791 0 +440 1.444140760829791 1.857921811274564 0 +441 1.136280332156326 1.556341077653394 0 +442 1.556341077653399 1.136280332156332 0 +443 1.40999470344446 1.352208600838271 0 +444 1.352208600838271 1.409994703444459 0 +445 1.309819578068145 1.808335231884333 0 +446 1.808335231884333 1.309819578068144 0 +447 1.084875043988303 1.297372391674491 0 +448 1.297372391674492 1.084875043988305 0 +449 1.915124956011692 1.702627608325507 0 +450 1.702627608325511 1.915124956011696 0 +451 1.96728574125623 1.574999999999999 0 +452 1.575000000000001 1.967285741256231 0 +453 1.424999999999999 1.03271425874377 0 +454 1.032714258743767 1.424999999999999 0 +455 1.698878894912874 1.502517494069071 0 +456 1.502517494069072 1.698878894912873 0 +457 1.846902227097242 1.391116810896704 0 +458 1.607173062007453 1.151428605119598 0 +459 1.391116810896704 1.846902227097242 0 +460 1.151428605119598 1.607173062007452 0 +461 1.9112556142141 1.217730482000099 0 +462 1.217730482000098 1.911255614214093 0 +463 1.777750208680053 1.085243579877084 0 +464 1.085243579877085 1.777750208680054 0 +465 1.484869921367756 1.745949034468182 0 +466 1.745949034468183 1.484869921367753 0 +467 1.429529776034876 1.216884455408708 0 +468 1.216884455408707 1.429529776034876 0 +469 1.473120591081134 1.207732356441173 0 +470 1.207732356441171 1.473120591081135 0 +471 1.466073088742609 1.816450260706398 0 +472 1.816450260706398 1.466073088742608 0 +473 1.704057891539513 1.102487875841138 0 +474 1.102487875841139 1.704057891539513 0 +475 1.904540875113385 1.292794954346816 0 +476 1.292794954346815 1.904540875113382 0 +477 1.728276094446232 1.685517009002038 0 +478 1.685517009002037 1.728276094446232 0 +479 1.314482990997958 1.271723905553764 0 +480 1.271723905553763 1.314482990997958 0 +481 1.80728348784079 1.51988612304525 0 +482 1.519886123045253 1.80728348784079 0 +483 1.173411664763546 1.401998711510771 0 +484 1.401998711510771 1.173411664763546 0 +485 1.831970845968198 1.59717302541845 0 +486 1.59717302541845 1.831970845968198 0 +487 1.569493741514673 1.517597172853245 0 +488 1.51759717285325 1.569493741514683 0 +489 1.680954855266541 1.45984090164441 0 +490 1.459840901644411 1.680954855266542 0 +491 1.092077241209449 1.961652779457519 0 +492 1.907922758790555 1.038347220542478 0 +493 1.038347220542481 1.907922758790553 0 +494 1.961652779457521 1.092077241209446 0 +495 1.252848370831241 1.643328768354839 0 +496 1.643328768354836 1.252848370831241 0 +497 1.644457955401133 1.187364320287396 0 +498 1.187364320287396 1.644457955401135 0 +499 1.901808597998297 1.948249613731501 0 +500 1.051750386268496 1.098191402001701 0 +501 1.098191402001701 1.051750386268495 0 +502 1.948249613731502 1.901808597998299 0 +503 1.923774946108954 1.754719853260077 0 +504 1.24528014673992 1.076225053891043 0 +505 1.076225053891038 1.245280146739918 0 +506 1.75471985326008 1.923774946108961 0 +507 1.411649301933745 1.459960187511624 0 +508 1.459960187511626 1.411649301933746 0 +509 1.226945553051417 1.579306513802602 0 +510 1.579306513802602 1.226945553051422 0 +511 1.25740570291479 1.373535048971603 0 +512 1.373535048971602 1.25740570291479 0 +513 1.743398707548301 1.625589672142274 0 +514 1.625589672142273 1.7433987075483 0 +515 1.472688691578712 1.517753198370208 0 +516 1.517753198370206 1.472688691578714 0 +517 1.086078480510835 1.859477342771219 0 +518 1.859477342771219 1.086078480510835 0 +519 1.753960051115772 1.357324464105859 0 +520 1.357324464105865 1.753960051115777 0 +521 1.028642603999476 1.971357396000524 0 +522 1.968950280659451 1.031049719340549 0 +523 1.519915916080453 1.032095169979427 0 +524 1.967905372291202 1.480081738999314 0 +525 1.032095169979427 1.519915916080453 0 +526 1.480081738999314 1.967905372291202 0 +527 1.500809145802982 1.248314946748099 0 +528 1.248314946748097 1.500809145802985 0 +529 1.34868062903739 1.811421996956815 0 +530 1.811421996956813 1.348680629037391 0 +531 1.558278462607639 1.910646343349908 0 +532 1.910646343349907 1.558278462607638 0 +533 1.441721537392361 1.089353656650092 0 +534 1.089353656650092 1.441721537392362 0 +535 1.443216816796153 1.735754148812604 0 +536 1.735754148812604 1.44321681679615 0 +537 1.619693631669934 1.656168437295585 0 +538 1.656168437295583 1.619693631669937 0 +539 1.484413480003034 1.353410107088967 0 +540 1.353410107088966 1.48441348000303 0 +541 1.472633581888515 1.298472000336959 0 +542 1.298472000336958 1.472633581888514 0 +543 1.187169333441722 1.504826802514988 0 +544 1.50482680251499 1.187169333441724 0 +545 1.351351781039616 1.575259151101544 0 +546 1.575259151101543 1.351351781039617 0 +547 1.620739058144107 1.308759843763949 0 +548 1.308759843763947 1.620739058144109 0 +549 1.41766541149283 1.606520230434032 0 +550 1.606520230434032 1.417665411492829 0 +551 1.630816238376309 1.512133655151921 0 +552 1.512133655151918 1.630816238376305 0 +553 1.516051966046563 1.76949450116318 0 +554 1.769494501163181 1.516051966046561 0 +555 1.73598231203687 1.163291120387096 0 +556 1.163291120387099 1.735982312036869 0 +557 1.713952319979456 1.602699825048679 0 +558 1.602699825048677 1.713952319979454 0 +$EndNodes +$Elements +1063 +1 1 2 2 5 1 81 +2 1 2 2 5 81 82 +3 1 2 2 5 82 83 +4 1 2 2 5 83 84 +5 1 2 2 5 84 85 +6 1 2 2 5 85 86 +7 1 2 2 5 86 87 +8 1 2 2 5 87 88 +9 1 2 2 5 88 89 +10 1 2 2 5 89 90 +11 1 2 2 5 90 91 +12 1 2 2 5 91 92 +13 1 2 2 5 92 93 +14 1 2 2 5 93 94 +15 1 2 2 5 94 95 +16 1 2 2 5 95 96 +17 1 2 2 5 96 97 +18 1 2 2 5 97 98 +19 1 2 2 5 98 99 +20 1 2 2 5 99 100 +21 1 2 2 5 100 101 +22 1 2 2 5 101 102 +23 1 2 2 5 102 103 +24 1 2 2 5 103 104 +25 1 2 2 5 104 105 +26 1 2 2 5 105 106 +27 1 2 2 5 106 107 +28 1 2 2 5 107 108 +29 1 2 2 5 108 3 +30 2 2 1 1 213 128 540 +31 2 2 1 1 127 214 539 +32 2 2 1 1 44 240 499 +33 2 2 1 1 241 79 500 +34 2 2 1 1 6 242 501 +35 2 2 1 1 239 41 502 +36 2 2 1 1 44 45 240 +37 2 2 1 1 40 41 239 +38 2 2 1 1 78 79 241 +39 2 2 1 1 6 7 242 +40 2 2 1 1 128 213 351 +41 2 2 1 1 214 127 352 +42 2 2 1 1 158 248 256 +43 2 2 1 1 247 157 255 +44 2 2 1 1 248 147 256 +45 2 2 1 1 148 247 255 +46 2 2 1 1 68 236 437 +47 2 2 1 1 233 17 435 +48 2 2 1 1 234 55 436 +49 2 2 1 1 30 235 438 +50 2 2 1 1 297 14 523 +51 2 2 1 1 33 296 524 +52 2 2 1 1 295 52 526 +53 2 2 1 1 71 298 525 +54 2 2 1 1 68 69 236 +55 2 2 1 1 54 55 234 +56 2 2 1 1 16 17 233 +57 2 2 1 1 30 31 235 +58 2 2 1 1 223 193 424 +59 2 2 1 1 195 226 423 +60 2 2 1 1 188 271 391 +61 2 2 1 1 273 190 392 +62 2 2 1 1 296 221 524 +63 2 2 1 1 220 297 523 +64 2 2 1 1 219 295 526 +65 2 2 1 1 298 222 525 +66 2 2 1 1 223 424 544 +67 2 2 1 1 423 226 543 +68 2 2 1 1 159 285 292 +69 2 2 1 1 286 160 291 +70 2 2 1 1 148 255 461 +71 2 2 1 1 256 147 462 +72 2 2 1 1 147 258 462 +73 2 2 1 1 260 148 461 +74 2 2 1 1 255 333 461 +75 2 2 1 1 334 256 462 +76 2 2 1 1 258 308 462 +77 2 2 1 1 305 260 461 +78 2 2 1 1 90 231 443 +79 2 2 1 1 232 90 444 +80 2 2 1 1 219 134 295 +81 2 2 1 1 135 221 296 +82 2 2 1 1 133 222 298 +83 2 2 1 1 220 136 297 +84 2 2 1 1 150 215 415 +85 2 2 1 1 216 149 416 +86 2 2 1 1 227 121 365 +87 2 2 1 1 122 228 366 +88 2 2 1 1 96 97 310 +89 2 2 1 1 97 96 309 +90 2 2 1 1 259 146 464 +91 2 2 1 1 145 257 463 +92 2 2 1 1 271 137 391 +93 2 2 1 1 138 273 392 +94 2 2 1 1 184 236 315 +95 2 2 1 1 234 186 316 +96 2 2 1 1 185 235 313 +97 2 2 1 1 233 183 314 +98 2 2 1 1 307 259 464 +99 2 2 1 1 257 306 463 +100 2 2 1 1 115 301 378 +101 2 2 1 1 302 116 377 +102 2 2 1 1 333 115 378 +103 2 2 1 1 116 334 377 +104 2 2 1 1 96 310 488 +105 2 2 1 1 309 96 487 +106 2 2 1 1 285 158 292 +107 2 2 1 1 157 286 291 +108 2 2 1 1 91 90 443 +109 2 2 1 1 90 91 444 +110 2 2 1 1 440 155 459 +111 2 2 1 1 156 441 460 +112 2 2 1 1 442 153 458 +113 2 2 1 1 154 439 457 +114 2 2 1 1 51 52 295 +115 2 2 1 1 13 14 297 +116 2 2 1 1 33 34 296 +117 2 2 1 1 71 72 298 +118 2 2 1 1 261 141 397 +119 2 2 1 1 142 262 398 +120 2 2 1 1 155 316 459 +121 2 2 1 1 315 156 460 +122 2 2 1 1 153 314 458 +123 2 2 1 1 313 154 457 +124 2 2 1 1 341 201 374 +125 2 2 1 1 202 342 375 +126 2 2 1 1 340 200 376 +127 2 2 1 1 199 339 373 +128 2 2 1 1 141 283 397 +129 2 2 1 1 284 142 398 +130 2 2 1 1 113 423 543 +131 2 2 1 1 424 114 544 +132 2 2 1 1 270 181 361 +133 2 2 1 1 179 268 359 +134 2 2 1 1 267 182 360 +135 2 2 1 1 180 269 362 +136 2 2 1 1 236 184 437 +137 2 2 1 1 186 234 436 +138 2 2 1 1 235 185 438 +139 2 2 1 1 183 233 435 +140 2 2 1 1 160 395 417 +141 2 2 1 1 396 159 418 +142 2 2 1 1 215 150 403 +143 2 2 1 1 149 216 404 +144 2 2 1 1 187 278 387 +145 2 2 1 1 190 276 388 +146 2 2 1 1 275 189 390 +147 2 2 1 1 277 188 389 +148 2 2 1 1 254 159 292 +149 2 2 1 1 160 253 291 +150 2 2 1 1 118 270 361 +151 2 2 1 1 268 119 359 +152 2 2 1 1 120 267 360 +153 2 2 1 1 269 117 362 +154 2 2 1 1 90 232 324 +155 2 2 1 1 231 90 323 +156 2 2 1 1 153 220 329 +157 2 2 1 1 221 154 331 +158 2 2 1 1 222 156 332 +159 2 2 1 1 155 219 330 +160 2 2 1 1 316 186 410 +161 2 2 1 1 185 313 409 +162 2 2 1 1 138 218 404 +163 2 2 1 1 217 137 403 +164 2 2 1 1 177 341 374 +165 2 2 1 1 342 178 375 +166 2 2 1 1 175 340 376 +167 2 2 1 1 339 176 373 +168 2 2 1 1 287 109 396 +169 2 2 1 1 110 288 395 +170 2 2 1 1 128 371 540 +171 2 2 1 1 372 127 539 +172 2 2 1 1 144 228 299 +173 2 2 1 1 227 143 300 +174 2 2 1 1 158 256 292 +175 2 2 1 1 255 157 291 +176 2 2 1 1 115 255 291 +177 2 2 1 1 256 116 292 +178 2 2 1 1 287 396 495 +179 2 2 1 1 395 288 496 +180 2 2 1 1 140 231 512 +181 2 2 1 1 232 139 511 +182 2 2 1 1 365 121 490 +183 2 2 1 1 122 366 489 +184 2 2 1 1 134 219 380 +185 2 2 1 1 222 133 381 +186 2 2 1 1 136 220 382 +187 2 2 1 1 221 135 379 +188 2 2 1 1 146 254 353 +189 2 2 1 1 253 145 354 +190 2 2 1 1 54 234 330 +191 2 2 1 1 236 69 332 +192 2 2 1 1 235 31 331 +193 2 2 1 1 16 233 329 +194 2 2 1 1 205 344 362 +195 2 2 1 1 345 203 360 +196 2 2 1 1 204 343 359 +197 2 2 1 1 346 206 361 +198 2 2 1 1 103 104 250 +199 2 2 1 1 85 86 249 +200 2 2 1 1 104 103 252 +201 2 2 1 1 86 85 251 +202 2 2 1 1 37 38 243 +203 2 2 1 1 9 10 245 +204 2 2 1 1 47 48 246 +205 2 2 1 1 75 76 244 +206 2 2 1 1 137 217 391 +207 2 2 1 1 218 138 392 +208 2 2 1 1 162 229 351 +209 2 2 1 1 230 161 352 +210 2 2 1 1 294 144 304 +211 2 2 1 1 143 293 303 +212 2 2 1 1 176 267 373 +213 2 2 1 1 269 177 374 +214 2 2 1 1 178 270 375 +215 2 2 1 1 268 175 376 +216 2 2 1 1 196 294 304 +217 2 2 1 1 293 194 303 +218 2 2 1 1 253 115 291 +219 2 2 1 1 116 254 292 +220 2 2 1 1 184 315 433 +221 2 2 1 1 314 183 434 +222 2 2 1 1 267 120 373 +223 2 2 1 1 270 118 375 +224 2 2 1 1 119 268 376 +225 2 2 1 1 117 269 374 +226 2 2 1 1 344 180 362 +227 2 2 1 1 182 345 360 +228 2 2 1 1 181 346 361 +229 2 2 1 1 343 179 359 +230 2 2 1 1 166 230 546 +231 2 2 1 1 229 165 545 +232 2 2 1 1 238 151 492 +233 2 2 1 1 151 238 494 +234 2 2 1 1 237 152 491 +235 2 2 1 1 152 237 493 +236 2 2 1 1 1 81 80 +237 2 2 1 1 5 81 1 +238 2 2 1 1 42 3 108 +239 2 2 1 1 108 3 43 +240 2 2 1 1 258 57 308 +241 2 2 1 1 66 259 307 +242 2 2 1 1 28 260 305 +243 2 2 1 1 257 19 306 +244 2 2 1 1 165 227 383 +245 2 2 1 1 228 166 384 +246 2 2 1 1 123 248 445 +247 2 2 1 1 247 124 446 +248 2 2 1 1 154 221 379 +249 2 2 1 1 220 153 382 +250 2 2 1 1 219 155 380 +251 2 2 1 1 156 222 381 +252 2 2 1 1 225 134 380 +253 2 2 1 1 223 136 382 +254 2 2 1 1 133 226 381 +255 2 2 1 1 135 224 379 +256 2 2 1 1 56 57 258 +257 2 2 1 1 18 19 257 +258 2 2 1 1 28 29 260 +259 2 2 1 1 66 67 259 +260 2 2 1 1 59 60 279 +261 2 2 1 1 21 22 282 +262 2 2 1 1 25 26 281 +263 2 2 1 1 63 64 280 +264 2 2 1 1 278 136 387 +265 2 2 1 1 276 134 388 +266 2 2 1 1 133 275 390 +267 2 2 1 1 135 277 389 +268 2 2 1 1 225 471 482 +269 2 2 1 1 472 224 481 +270 2 2 1 1 271 188 485 +271 2 2 1 1 190 273 486 +272 2 2 1 1 153 233 314 +273 2 2 1 1 235 154 313 +274 2 2 1 1 155 234 316 +275 2 2 1 1 236 156 315 +276 2 2 1 1 163 231 364 +277 2 2 1 1 232 164 363 +278 2 2 1 1 395 207 417 +279 2 2 1 1 208 396 418 +280 2 2 1 1 139 232 363 +281 2 2 1 1 231 140 364 +282 2 2 1 1 471 111 482 +283 2 2 1 1 112 472 481 +284 2 2 1 1 393 300 520 +285 2 2 1 1 299 394 519 +286 2 2 1 1 151 282 492 +287 2 2 1 1 281 151 494 +288 2 2 1 1 152 279 491 +289 2 2 1 1 280 152 493 +290 2 2 1 1 231 323 512 +291 2 2 1 1 324 232 511 +292 2 2 1 1 160 286 395 +293 2 2 1 1 285 159 396 +294 2 2 1 1 193 223 442 +295 2 2 1 1 226 195 441 +296 2 2 1 1 92 91 443 +297 2 2 1 1 91 92 444 +298 2 2 1 1 94 93 318 +299 2 2 1 1 93 94 317 +300 2 2 1 1 121 227 385 +301 2 2 1 1 228 122 386 +302 2 2 1 1 227 165 385 +303 2 2 1 1 166 228 386 +304 2 2 1 1 99 100 289 +305 2 2 1 1 100 99 290 +306 2 2 1 1 92 350 444 +307 2 2 1 1 349 92 443 +308 2 2 1 1 262 166 546 +309 2 2 1 1 165 261 545 +310 2 2 1 1 281 168 378 +311 2 2 1 1 167 279 377 +312 2 2 1 1 151 281 378 +313 2 2 1 1 279 152 377 +314 2 2 1 1 127 372 398 +315 2 2 1 1 371 128 397 +316 2 2 1 1 87 341 356 +317 2 2 1 1 342 87 357 +318 2 2 1 1 339 102 358 +319 2 2 1 1 102 340 355 +320 2 2 1 1 136 223 387 +321 2 2 1 1 134 225 388 +322 2 2 1 1 224 135 389 +323 2 2 1 1 226 133 390 +324 2 2 1 1 198 322 348 +325 2 2 1 1 321 197 347 +326 2 2 1 1 322 161 348 +327 2 2 1 1 162 321 347 +328 2 2 1 1 93 317 507 +329 2 2 1 1 318 93 508 +330 2 2 1 1 149 310 416 +331 2 2 1 1 309 150 415 +332 2 2 1 1 312 193 458 +333 2 2 1 1 195 311 460 +334 2 2 1 1 115 253 354 +335 2 2 1 1 254 116 353 +336 2 2 1 1 233 153 329 +337 2 2 1 1 154 235 331 +338 2 2 1 1 234 155 330 +339 2 2 1 1 156 236 332 +340 2 2 1 1 276 190 486 +341 2 2 1 1 189 275 483 +342 2 2 1 1 278 187 484 +343 2 2 1 1 188 277 485 +344 2 2 1 1 228 144 366 +345 2 2 1 1 143 227 365 +346 2 2 1 1 247 148 475 +347 2 2 1 1 147 248 476 +348 2 2 1 1 312 458 497 +349 2 2 1 1 460 311 498 +350 2 2 1 1 159 254 556 +351 2 2 1 1 253 160 555 +352 2 2 1 1 211 393 520 +353 2 2 1 1 394 212 519 +354 2 2 1 1 394 157 446 +355 2 2 1 1 158 393 445 +356 2 2 1 1 123 445 529 +357 2 2 1 1 446 124 530 +358 2 2 1 1 320 198 348 +359 2 2 1 1 197 319 347 +360 2 2 1 1 351 229 545 +361 2 2 1 1 230 352 546 +362 2 2 1 1 122 320 348 +363 2 2 1 1 319 121 347 +364 2 2 1 1 286 157 394 +365 2 2 1 1 158 285 393 +366 2 2 1 1 313 124 409 +367 2 2 1 1 123 316 410 +368 2 2 1 1 90 89 323 +369 2 2 1 1 89 90 324 +370 2 2 1 1 111 293 465 +371 2 2 1 1 294 112 466 +372 2 2 1 1 229 162 549 +373 2 2 1 1 161 230 550 +374 2 2 1 1 109 300 393 +375 2 2 1 1 299 110 394 +376 2 2 1 1 455 320 489 +377 2 2 1 1 319 456 490 +378 2 2 1 1 61 237 491 +379 2 2 1 1 237 62 493 +380 2 2 1 1 238 24 494 +381 2 2 1 1 23 238 492 +382 2 2 1 1 62 237 521 +383 2 2 1 1 238 23 522 +384 2 2 1 1 237 61 521 +385 2 2 1 1 24 238 522 +386 2 2 1 1 243 38 326 +387 2 2 1 1 9 245 325 +388 2 2 1 1 47 246 328 +389 2 2 1 1 244 76 327 +390 2 2 1 1 146 353 464 +391 2 2 1 1 354 145 463 +392 2 2 1 1 230 166 386 +393 2 2 1 1 165 229 385 +394 2 2 1 1 193 312 510 +395 2 2 1 1 311 195 509 +396 2 2 1 1 424 193 510 +397 2 2 1 1 195 423 509 +398 2 2 1 1 372 210 398 +399 2 2 1 1 209 371 397 +400 2 2 1 1 57 58 308 +401 2 2 1 1 65 66 307 +402 2 2 1 1 19 20 306 +403 2 2 1 1 27 28 305 +404 2 2 1 1 231 163 443 +405 2 2 1 1 164 232 444 +406 2 2 1 1 69 70 332 +407 2 2 1 1 53 54 330 +408 2 2 1 1 31 32 331 +409 2 2 1 1 15 16 329 +410 2 2 1 1 75 244 399 +411 2 2 1 1 246 48 401 +412 2 2 1 1 37 243 400 +413 2 2 1 1 245 10 402 +414 2 2 1 1 254 146 556 +415 2 2 1 1 145 253 555 +416 2 2 1 1 250 104 336 +417 2 2 1 1 104 252 335 +418 2 2 1 1 85 249 338 +419 2 2 1 1 251 85 337 +420 2 2 1 1 128 351 545 +421 2 2 1 1 352 127 546 +422 2 2 1 1 106 240 414 +423 2 2 1 1 239 106 411 +424 2 2 1 1 241 83 413 +425 2 2 1 1 83 242 412 +426 2 2 1 1 255 115 333 +427 2 2 1 1 116 256 334 +428 2 2 1 1 201 272 374 +429 2 2 1 1 274 202 375 +430 2 2 1 1 200 271 376 +431 2 2 1 1 273 199 373 +432 2 2 1 1 290 99 538 +433 2 2 1 1 99 289 537 +434 2 2 1 1 467 274 484 +435 2 2 1 1 272 468 483 +436 2 2 1 1 209 264 542 +437 2 2 1 1 263 210 541 +438 2 2 1 1 301 115 354 +439 2 2 1 1 116 302 353 +440 2 2 1 1 465 293 535 +441 2 2 1 1 294 466 536 +442 2 2 1 1 146 265 556 +443 2 2 1 1 266 145 555 +444 2 2 1 1 227 300 383 +445 2 2 1 1 299 228 384 +446 2 2 1 1 262 127 398 +447 2 2 1 1 128 261 397 +448 2 2 1 1 341 177 356 +449 2 2 1 1 178 342 357 +450 2 2 1 1 176 339 358 +451 2 2 1 1 340 175 355 +452 2 2 1 1 252 103 355 +453 2 2 1 1 86 251 357 +454 2 2 1 1 249 86 356 +455 2 2 1 1 103 250 358 +456 2 2 1 1 457 304 530 +457 2 2 1 1 303 459 529 +458 2 2 1 1 242 83 501 +459 2 2 1 1 83 241 500 +460 2 2 1 1 106 239 502 +461 2 2 1 1 240 106 499 +462 2 2 1 1 279 60 491 +463 2 2 1 1 282 22 492 +464 2 2 1 1 63 280 493 +465 2 2 1 1 25 281 494 +466 2 2 1 1 258 186 436 +467 2 2 1 1 184 259 437 +468 2 2 1 1 185 260 438 +469 2 2 1 1 257 183 435 +470 2 2 1 1 56 258 436 +471 2 2 1 1 18 257 435 +472 2 2 1 1 260 29 438 +473 2 2 1 1 259 67 437 +474 2 2 1 1 265 126 418 +475 2 2 1 1 125 266 417 +476 2 2 1 1 177 269 429 +477 2 2 1 1 267 176 432 +478 2 2 1 1 270 178 430 +479 2 2 1 1 175 268 431 +480 2 2 1 1 150 320 455 +481 2 2 1 1 319 149 456 +482 2 2 1 1 269 130 429 +483 2 2 1 1 132 270 430 +484 2 2 1 1 268 129 431 +485 2 2 1 1 131 267 432 +486 2 2 1 1 215 290 538 +487 2 2 1 1 289 216 537 +488 2 2 1 1 261 128 545 +489 2 2 1 1 127 262 546 +490 2 2 1 1 315 126 433 +491 2 2 1 1 125 314 434 +492 2 2 1 1 40 239 405 +493 2 2 1 1 240 45 408 +494 2 2 1 1 78 241 407 +495 2 2 1 1 242 7 406 +496 2 2 1 1 275 205 362 +497 2 2 1 1 277 204 359 +498 2 2 1 1 203 276 360 +499 2 2 1 1 206 278 361 +500 2 2 1 1 187 467 484 +501 2 2 1 1 468 189 483 +502 2 2 1 1 166 262 547 +503 2 2 1 1 261 165 548 +504 2 2 1 1 317 213 507 +505 2 2 1 1 214 318 508 +506 2 2 1 1 245 402 448 +507 2 2 1 1 399 244 447 +508 2 2 1 1 246 401 450 +509 2 2 1 1 400 243 449 +510 2 2 1 1 171 239 411 +511 2 2 1 1 242 174 412 +512 2 2 1 1 173 241 413 +513 2 2 1 1 240 172 414 +514 2 2 1 1 264 139 363 +515 2 2 1 1 140 263 364 +516 2 2 1 1 239 171 405 +517 2 2 1 1 174 242 406 +518 2 2 1 1 241 173 407 +519 2 2 1 1 172 240 408 +520 2 2 1 1 196 304 457 +521 2 2 1 1 303 194 459 +522 2 2 1 1 269 180 447 +523 2 2 1 1 181 270 448 +524 2 2 1 1 268 179 449 +525 2 2 1 1 182 267 450 +526 2 2 1 1 140 274 467 +527 2 2 1 1 272 139 468 +528 2 2 1 1 318 214 352 +529 2 2 1 1 213 317 351 +530 2 2 1 1 317 162 351 +531 2 2 1 1 161 318 352 +532 2 2 1 1 124 247 409 +533 2 2 1 1 248 123 410 +534 2 2 1 1 145 266 473 +535 2 2 1 1 265 146 474 +536 2 2 1 1 249 130 425 +537 2 2 1 1 131 250 428 +538 2 2 1 1 132 251 426 +539 2 2 1 1 252 129 427 +540 2 2 1 1 409 247 475 +541 2 2 1 1 248 410 476 +542 2 2 1 1 248 158 445 +543 2 2 1 1 157 247 446 +544 2 2 1 1 119 277 359 +545 2 2 1 1 278 118 361 +546 2 2 1 1 276 120 360 +547 2 2 1 1 117 275 362 +548 2 2 1 1 212 394 446 +549 2 2 1 1 393 211 445 +550 2 2 1 1 198 309 487 +551 2 2 1 1 310 197 488 +552 2 2 1 1 141 261 548 +553 2 2 1 1 262 142 547 +554 2 2 1 1 310 97 416 +555 2 2 1 1 97 309 415 +556 2 2 1 1 102 101 340 +557 2 2 1 1 101 102 339 +558 2 2 1 1 87 88 341 +559 2 2 1 1 88 87 342 +560 2 2 1 1 388 225 482 +561 2 2 1 1 224 389 481 +562 2 2 1 1 210 263 527 +563 2 2 1 1 264 209 528 +564 2 2 1 1 130 249 429 +565 2 2 1 1 129 252 431 +566 2 2 1 1 250 131 432 +567 2 2 1 1 251 132 430 +568 2 2 1 1 168 281 421 +569 2 2 1 1 282 169 420 +570 2 2 1 1 170 280 422 +571 2 2 1 1 279 167 419 +572 2 2 1 1 130 269 447 +573 2 2 1 1 270 132 448 +574 2 2 1 1 129 268 449 +575 2 2 1 1 267 131 450 +576 2 2 1 1 21 282 420 +577 2 2 1 1 281 26 421 +578 2 2 1 1 280 64 422 +579 2 2 1 1 59 279 419 +580 2 2 1 1 35 343 451 +581 2 2 1 1 345 50 452 +582 2 2 1 1 346 12 453 +583 2 2 1 1 73 344 454 +584 2 2 1 1 60 61 491 +585 2 2 1 1 62 63 493 +586 2 2 1 1 24 25 494 +587 2 2 1 1 22 23 492 +588 2 2 1 1 118 274 375 +589 2 2 1 1 120 273 373 +590 2 2 1 1 272 117 374 +591 2 2 1 1 271 119 376 +592 2 2 1 1 289 100 478 +593 2 2 1 1 100 290 477 +594 2 2 1 1 146 259 474 +595 2 2 1 1 257 145 473 +596 2 2 1 1 38 39 326 +597 2 2 1 1 8 9 325 +598 2 2 1 1 46 47 328 +599 2 2 1 1 76 77 327 +600 2 2 1 1 259 184 474 +601 2 2 1 1 183 257 473 +602 2 2 1 1 124 457 530 +603 2 2 1 1 459 123 529 +604 2 2 1 1 35 36 343 +605 2 2 1 1 49 50 345 +606 2 2 1 1 73 74 344 +607 2 2 1 1 11 12 346 +608 2 2 1 1 199 289 478 +609 2 2 1 1 290 200 477 +610 2 2 1 1 118 278 484 +611 2 2 1 1 275 117 483 +612 2 2 1 1 120 276 486 +613 2 2 1 1 277 119 485 +614 2 2 1 1 135 296 532 +615 2 2 1 1 297 136 533 +616 2 2 1 1 133 298 534 +617 2 2 1 1 295 134 531 +618 2 2 1 1 318 161 516 +619 2 2 1 1 162 317 515 +620 2 2 1 1 202 274 512 +621 2 2 1 1 272 201 511 +622 2 2 1 1 192 455 489 +623 2 2 1 1 456 191 490 +624 2 2 1 1 148 260 475 +625 2 2 1 1 258 147 476 +626 2 2 1 1 264 363 542 +627 2 2 1 1 364 263 541 +628 2 2 1 1 186 258 476 +629 2 2 1 1 260 185 475 +630 2 2 1 1 274 140 512 +631 2 2 1 1 139 272 511 +632 2 2 1 1 109 287 383 +633 2 2 1 1 288 110 384 +634 2 2 1 1 271 200 513 +635 2 2 1 1 199 273 514 +636 2 2 1 1 142 288 547 +637 2 2 1 1 287 141 548 +638 2 2 1 1 137 271 513 +639 2 2 1 1 273 138 514 +640 2 2 1 1 104 105 336 +641 2 2 1 1 105 104 335 +642 2 2 1 1 84 85 338 +643 2 2 1 1 85 84 337 +644 2 2 1 1 287 165 383 +645 2 2 1 1 166 288 384 +646 2 2 1 1 285 109 393 +647 2 2 1 1 110 286 394 +648 2 2 1 1 207 395 496 +649 2 2 1 1 396 208 495 +650 2 2 1 1 161 322 516 +651 2 2 1 1 321 162 515 +652 2 2 1 1 169 301 354 +653 2 2 1 1 302 170 353 +654 2 2 1 1 109 285 396 +655 2 2 1 1 286 110 395 +656 2 2 1 1 283 209 397 +657 2 2 1 1 210 284 398 +658 2 2 1 1 311 141 495 +659 2 2 1 1 142 312 496 +660 2 2 1 1 138 404 558 +661 2 2 1 1 403 137 557 +662 2 2 1 1 284 210 527 +663 2 2 1 1 209 283 528 +664 2 2 1 1 208 311 495 +665 2 2 1 1 312 207 496 +666 2 2 1 1 263 140 467 +667 2 2 1 1 139 264 468 +668 2 2 1 1 41 42 369 +669 2 2 1 1 43 44 367 +670 2 2 1 1 79 80 370 +671 2 2 1 1 5 6 368 +672 2 2 1 1 114 284 527 +673 2 2 1 1 283 113 528 +674 2 2 1 1 97 98 416 +675 2 2 1 1 98 97 415 +676 2 2 1 1 126 265 433 +677 2 2 1 1 266 125 434 +678 2 2 1 1 86 87 356 +679 2 2 1 1 87 86 357 +680 2 2 1 1 103 102 355 +681 2 2 1 1 102 103 358 +682 2 2 1 1 165 287 548 +683 2 2 1 1 288 166 547 +684 2 2 1 1 136 278 533 +685 2 2 1 1 277 135 532 +686 2 2 1 1 134 276 531 +687 2 2 1 1 275 133 534 +688 2 2 1 1 179 400 449 +689 2 2 1 1 401 182 450 +690 2 2 1 1 402 181 448 +691 2 2 1 1 180 399 447 +692 2 2 1 1 290 215 557 +693 2 2 1 1 216 289 558 +694 2 2 1 1 204 277 532 +695 2 2 1 1 205 275 534 +696 2 2 1 1 278 206 533 +697 2 2 1 1 276 203 531 +698 2 2 1 1 343 204 451 +699 2 2 1 1 203 345 452 +700 2 2 1 1 344 205 454 +701 2 2 1 1 206 346 453 +702 2 2 1 1 280 170 517 +703 2 2 1 1 169 282 518 +704 2 2 1 1 169 306 420 +705 2 2 1 1 305 168 421 +706 2 2 1 1 307 170 422 +707 2 2 1 1 167 308 419 +708 2 2 1 1 65 307 422 +709 2 2 1 1 306 20 420 +710 2 2 1 1 27 305 421 +711 2 2 1 1 308 58 419 +712 2 2 1 1 284 114 424 +713 2 2 1 1 113 283 423 +714 2 2 1 1 300 109 383 +715 2 2 1 1 110 299 384 +716 2 2 1 1 293 143 535 +717 2 2 1 1 144 294 536 +718 2 2 1 1 152 302 377 +719 2 2 1 1 301 151 378 +720 2 2 1 1 283 141 509 +721 2 2 1 1 142 284 510 +722 2 2 1 1 152 280 517 +723 2 2 1 1 282 151 518 +724 2 2 1 1 321 95 488 +725 2 2 1 1 95 322 487 +726 2 2 1 1 117 272 483 +727 2 2 1 1 274 118 484 +728 2 2 1 1 273 120 486 +729 2 2 1 1 119 271 485 +730 2 2 1 1 293 111 471 +731 2 2 1 1 112 294 472 +732 2 2 1 1 194 293 471 +733 2 2 1 1 294 196 472 +734 2 2 1 1 95 96 488 +735 2 2 1 1 96 95 487 +736 2 2 1 1 105 106 414 +737 2 2 1 1 106 105 411 +738 2 2 1 1 83 84 413 +739 2 2 1 1 84 83 412 +740 2 2 1 1 107 108 367 +741 2 2 1 1 108 107 369 +742 2 2 1 1 82 81 368 +743 2 2 1 1 81 82 370 +744 2 2 1 1 326 171 503 +745 2 2 1 1 174 325 504 +746 2 2 1 1 327 173 505 +747 2 2 1 1 172 328 506 +748 2 2 1 1 289 199 514 +749 2 2 1 1 200 290 513 +750 2 2 1 1 141 287 495 +751 2 2 1 1 288 142 496 +752 2 2 1 1 298 205 534 +753 2 2 1 1 296 204 532 +754 2 2 1 1 206 297 533 +755 2 2 1 1 203 295 531 +756 2 2 1 1 326 39 405 +757 2 2 1 1 8 325 406 +758 2 2 1 1 46 328 408 +759 2 2 1 1 327 77 407 +760 2 2 1 1 194 440 459 +761 2 2 1 1 441 195 460 +762 2 2 1 1 439 196 457 +763 2 2 1 1 193 442 458 +764 2 2 1 1 171 326 405 +765 2 2 1 1 325 174 406 +766 2 2 1 1 173 327 407 +767 2 2 1 1 328 172 408 +768 2 2 1 1 348 161 550 +769 2 2 1 1 162 347 549 +770 2 2 1 1 367 44 499 +771 2 2 1 1 41 369 502 +772 2 2 1 1 79 370 500 +773 2 2 1 1 368 6 501 +774 2 2 1 1 130 447 505 +775 2 2 1 1 448 132 504 +776 2 2 1 1 129 449 503 +777 2 2 1 1 450 131 506 +778 2 2 1 1 304 144 519 +779 2 2 1 1 143 303 520 +780 2 2 1 1 83 82 501 +781 2 2 1 1 82 83 500 +782 2 2 1 1 106 107 499 +783 2 2 1 1 107 106 502 +784 2 2 1 1 212 304 519 +785 2 2 1 1 303 211 520 +786 2 2 1 1 205 298 454 +787 2 2 1 1 295 203 452 +788 2 2 1 1 297 206 453 +789 2 2 1 1 204 296 451 +790 2 2 1 1 51 295 452 +791 2 2 1 1 298 72 454 +792 2 2 1 1 13 297 453 +793 2 2 1 1 296 34 451 +794 2 2 1 1 187 387 469 +795 2 2 1 1 390 189 470 +796 2 2 1 1 447 244 505 +797 2 2 1 1 245 448 504 +798 2 2 1 1 449 243 503 +799 2 2 1 1 246 450 506 +800 2 2 1 1 329 220 523 +801 2 2 1 1 221 331 524 +802 2 2 1 1 330 219 526 +803 2 2 1 1 222 332 525 +804 2 2 1 1 168 333 378 +805 2 2 1 1 334 167 377 +806 2 2 1 1 45 46 408 +807 2 2 1 1 77 78 407 +808 2 2 1 1 7 8 406 +809 2 2 1 1 39 40 405 +810 2 2 1 1 160 417 555 +811 2 2 1 1 418 159 556 +812 2 2 1 1 417 266 555 +813 2 2 1 1 265 418 556 +814 2 2 1 1 105 335 411 +815 2 2 1 1 336 105 414 +816 2 2 1 1 337 84 412 +817 2 2 1 1 84 338 413 +818 2 2 1 1 171 427 503 +819 2 2 1 1 426 174 504 +820 2 2 1 1 173 425 505 +821 2 2 1 1 428 172 506 +822 2 2 1 1 108 43 367 +823 2 2 1 1 42 108 369 +824 2 2 1 1 81 5 368 +825 2 2 1 1 80 81 370 +826 2 2 1 1 302 152 517 +827 2 2 1 1 151 301 518 +828 2 2 1 1 144 299 519 +829 2 2 1 1 300 143 520 +830 2 2 1 1 306 169 463 +831 2 2 1 1 170 307 464 +832 2 2 1 1 168 305 461 +833 2 2 1 1 308 167 462 +834 2 2 1 1 345 182 401 +835 2 2 1 1 346 181 402 +836 2 2 1 1 180 344 399 +837 2 2 1 1 179 343 400 +838 2 2 1 1 418 126 498 +839 2 2 1 1 125 417 497 +840 2 2 1 1 36 37 400 +841 2 2 1 1 74 75 399 +842 2 2 1 1 10 11 402 +843 2 2 1 1 48 49 401 +844 2 2 1 1 126 460 498 +845 2 2 1 1 458 125 497 +846 2 2 1 1 49 345 401 +847 2 2 1 1 343 36 400 +848 2 2 1 1 344 74 399 +849 2 2 1 1 11 346 402 +850 2 2 1 1 309 198 551 +851 2 2 1 1 197 310 552 +852 2 2 1 1 170 302 517 +853 2 2 1 1 301 169 518 +854 2 2 1 1 266 434 473 +855 2 2 1 1 433 265 474 +856 2 2 1 1 163 349 443 +857 2 2 1 1 350 164 444 +858 2 2 1 1 469 387 544 +859 2 2 1 1 390 470 543 +860 2 2 1 1 304 212 530 +861 2 2 1 1 211 303 529 +862 2 2 1 1 124 313 457 +863 2 2 1 1 314 125 458 +864 2 2 1 1 126 315 460 +865 2 2 1 1 316 123 459 +866 2 2 1 1 440 194 471 +867 2 2 1 1 196 439 472 +868 2 2 1 1 439 224 472 +869 2 2 1 1 225 440 471 +870 2 2 1 1 150 309 551 +871 2 2 1 1 310 149 552 +872 2 2 1 1 250 336 428 +873 2 2 1 1 338 249 425 +874 2 2 1 1 251 337 426 +875 2 2 1 1 335 252 427 +876 2 2 1 1 141 311 509 +877 2 2 1 1 312 142 510 +878 2 2 1 1 311 208 498 +879 2 2 1 1 207 312 497 +880 2 2 1 1 174 337 412 +881 2 2 1 1 335 171 411 +882 2 2 1 1 338 173 413 +883 2 2 1 1 172 336 414 +884 2 2 1 1 333 168 461 +885 2 2 1 1 167 334 462 +886 2 2 1 1 58 59 419 +887 2 2 1 1 64 65 422 +888 2 2 1 1 20 21 420 +889 2 2 1 1 26 27 421 +890 2 2 1 1 143 365 535 +891 2 2 1 1 366 144 536 +892 2 2 1 1 340 101 477 +893 2 2 1 1 101 339 478 +894 2 2 1 1 88 342 479 +895 2 2 1 1 341 88 480 +896 2 2 1 1 200 340 477 +897 2 2 1 1 339 199 478 +898 2 2 1 1 342 202 479 +899 2 2 1 1 201 341 480 +900 2 2 1 1 403 150 455 +901 2 2 1 1 149 404 456 +902 2 2 1 1 404 218 456 +903 2 2 1 1 217 403 455 +904 2 2 1 1 95 321 515 +905 2 2 1 1 322 95 516 +906 2 2 1 1 55 56 436 +907 2 2 1 1 17 18 435 +908 2 2 1 1 29 30 438 +909 2 2 1 1 67 68 437 +910 2 2 1 1 89 324 480 +911 2 2 1 1 323 89 479 +912 2 2 1 1 202 323 479 +913 2 2 1 1 324 201 480 +914 2 2 1 1 94 95 515 +915 2 2 1 1 95 94 516 +916 2 2 1 1 23 2 522 +917 2 2 1 1 61 4 521 +918 2 2 1 1 2 24 522 +919 2 2 1 1 4 62 521 +920 2 2 1 1 15 329 523 +921 2 2 1 1 331 32 524 +922 2 2 1 1 53 330 526 +923 2 2 1 1 332 70 525 +924 2 2 1 1 188 391 481 +925 2 2 1 1 392 190 482 +926 2 2 1 1 215 403 557 +927 2 2 1 1 404 216 558 +928 2 2 1 1 12 13 453 +929 2 2 1 1 34 35 451 +930 2 2 1 1 50 51 452 +931 2 2 1 1 72 73 454 +932 2 2 1 1 317 94 515 +933 2 2 1 1 94 318 516 +934 2 2 1 1 171 335 427 +935 2 2 1 1 337 174 426 +936 2 2 1 1 173 338 425 +937 2 2 1 1 336 172 428 +938 2 2 1 1 358 250 432 +939 2 2 1 1 252 355 431 +940 2 2 1 1 357 251 430 +941 2 2 1 1 249 356 429 +942 2 2 1 1 176 358 432 +943 2 2 1 1 355 175 431 +944 2 2 1 1 356 177 429 +945 2 2 1 1 178 357 430 +946 2 2 1 1 121 319 490 +947 2 2 1 1 320 122 489 +948 2 2 1 1 93 92 508 +949 2 2 1 1 92 93 507 +950 2 2 1 1 470 264 528 +951 2 2 1 1 263 469 527 +952 2 2 1 1 366 192 489 +953 2 2 1 1 191 365 490 +954 2 2 1 1 113 470 528 +955 2 2 1 1 469 114 527 +956 2 2 1 1 456 218 465 +957 2 2 1 1 217 455 466 +958 2 2 1 1 320 150 551 +959 2 2 1 1 149 319 552 +960 2 2 1 1 191 456 465 +961 2 2 1 1 455 192 466 +962 2 2 1 1 322 198 487 +963 2 2 1 1 197 321 488 +964 2 2 1 1 198 320 551 +965 2 2 1 1 319 197 552 +966 2 2 1 1 101 100 477 +967 2 2 1 1 88 89 480 +968 2 2 1 1 100 101 478 +969 2 2 1 1 89 88 479 +970 2 2 1 1 223 382 442 +971 2 2 1 1 381 226 441 +972 2 2 1 1 379 224 439 +973 2 2 1 1 225 380 440 +974 2 2 1 1 382 153 442 +975 2 2 1 1 154 379 439 +976 2 2 1 1 156 381 441 +977 2 2 1 1 380 155 440 +978 2 2 1 1 122 348 550 +979 2 2 1 1 347 121 549 +980 2 2 1 1 325 245 504 +981 2 2 1 1 328 246 506 +982 2 2 1 1 244 327 505 +983 2 2 1 1 243 326 503 +984 2 2 1 1 323 202 512 +985 2 2 1 1 201 324 511 +986 2 2 1 1 92 349 508 +987 2 2 1 1 350 92 507 +988 2 2 1 1 52 53 526 +989 2 2 1 1 70 71 525 +990 2 2 1 1 14 15 523 +991 2 2 1 1 32 33 524 +992 2 2 1 1 169 354 463 +993 2 2 1 1 353 170 464 +994 2 2 1 1 98 99 537 +995 2 2 1 1 99 98 538 +996 2 2 1 1 466 192 536 +997 2 2 1 1 191 465 535 +998 2 2 1 1 208 418 498 +999 2 2 1 1 417 207 497 +1000 2 2 1 1 121 385 549 +1001 2 2 1 1 386 122 550 +1002 2 2 1 1 481 391 554 +1003 2 2 1 1 392 482 553 +1004 2 2 1 1 214 349 539 +1005 2 2 1 1 350 213 540 +1006 2 2 1 1 349 163 539 +1007 2 2 1 1 164 350 540 +1008 2 2 1 1 213 350 507 +1009 2 2 1 1 349 214 508 +1010 2 2 1 1 212 446 530 +1011 2 2 1 1 445 211 529 +1012 2 2 1 1 389 188 481 +1013 2 2 1 1 190 388 482 +1014 2 2 1 1 98 415 538 +1015 2 2 1 1 416 98 537 +1016 2 2 1 1 163 364 541 +1017 2 2 1 1 363 164 542 +1018 2 2 1 1 415 215 538 +1019 2 2 1 1 216 416 537 +1020 2 2 1 1 468 264 470 +1021 2 2 1 1 263 467 469 +1022 2 2 1 1 189 468 470 +1023 2 2 1 1 467 187 469 +1024 2 2 1 1 410 186 476 +1025 2 2 1 1 185 409 475 +1026 2 2 1 1 192 366 536 +1027 2 2 1 1 365 191 535 +1028 2 2 1 1 372 163 541 +1029 2 2 1 1 164 371 542 +1030 2 2 1 1 434 183 473 +1031 2 2 1 1 184 433 474 +1032 2 2 1 1 163 372 539 +1033 2 2 1 1 371 164 540 +1034 2 2 1 1 210 372 541 +1035 2 2 1 1 371 209 542 +1036 2 2 1 1 425 130 505 +1037 2 2 1 1 132 426 504 +1038 2 2 1 1 131 428 506 +1039 2 2 1 1 427 129 503 +1040 2 2 1 1 226 390 543 +1041 2 2 1 1 387 223 544 +1042 2 2 1 1 385 229 549 +1043 2 2 1 1 230 386 550 +1044 2 2 1 1 423 283 509 +1045 2 2 1 1 284 424 510 +1046 2 2 1 1 218 392 553 +1047 2 2 1 1 391 217 554 +1048 2 2 1 1 107 367 499 +1049 2 2 1 1 82 368 501 +1050 2 2 1 1 370 82 500 +1051 2 2 1 1 369 107 502 +1052 2 2 1 1 465 218 553 +1053 2 2 1 1 217 466 554 +1054 2 2 1 1 111 465 553 +1055 2 2 1 1 466 112 554 +1056 2 2 1 1 289 514 558 +1057 2 2 1 1 513 290 557 +1058 2 2 1 1 114 469 544 +1059 2 2 1 1 470 113 543 +1060 2 2 1 1 514 138 558 +1061 2 2 1 1 137 513 557 +1062 2 2 1 1 482 111 553 +1063 2 2 1 1 112 481 554 +$EndElements diff --git a/src/grids/mmdg5_8e-2.msh b/src/grids/mmdg5_8e-2.msh new file mode 100644 index 0000000000000000000000000000000000000000..e9179092d064f4b626f95901a32432a991e4e24c --- /dev/null +++ b/src/grids/mmdg5_8e-2.msh @@ -0,0 +1,732 @@ +$MeshFormat +2.2 0 8 +$EndMeshFormat +$Nodes +253 +1 1 1 0 +2 2 1 0 +3 2 2 0 +4 1 2 0 +5 1.076923076923077 1 0 +6 1.153846153846154 1 0 +7 1.23076923076923 1 0 +8 1.307692307692307 1 0 +9 1.384615384615384 1 0 +10 1.461538461538461 1 0 +11 1.538461538461538 1 0 +12 1.615384615384615 1 0 +13 1.692307692307692 1 0 +14 1.769230769230769 1 0 +15 1.846153846153846 1 0 +16 1.923076923076923 1 0 +17 2 1.076923076923077 0 +18 2 1.153846153846154 0 +19 2 1.23076923076923 0 +20 2 1.307692307692307 0 +21 2 1.384615384615384 0 +22 2 1.461538461538461 0 +23 2 1.538461538461538 0 +24 2 1.615384615384615 0 +25 2 1.692307692307692 0 +26 2 1.769230769230769 0 +27 2 1.846153846153846 0 +28 2 1.923076923076923 0 +29 1.923076923076923 2 0 +30 1.846153846153846 2 0 +31 1.76923076923077 2 0 +32 1.692307692307693 2 0 +33 1.615384615384616 2 0 +34 1.538461538461539 2 0 +35 1.461538461538462 2 0 +36 1.384615384615385 2 0 +37 1.307692307692308 2 0 +38 1.230769230769231 2 0 +39 1.153846153846154 2 0 +40 1.076923076923077 2 0 +41 1 1.923076923076923 0 +42 1 1.846153846153846 0 +43 1 1.76923076923077 0 +44 1 1.692307692307693 0 +45 1 1.615384615384616 0 +46 1 1.538461538461539 0 +47 1 1.461538461538462 0 +48 1 1.384615384615385 0 +49 1 1.307692307692308 0 +50 1 1.230769230769231 0 +51 1 1.153846153846154 0 +52 1 1.076923076923077 0 +53 1.055555555555555 1.055555555555555 0 +54 1.111111111111111 1.111111111111111 0 +55 1.166666666666666 1.166666666666666 0 +56 1.222222222222222 1.222222222222222 0 +57 1.277777777777777 1.277777777777777 0 +58 1.333333333333333 1.333333333333333 0 +59 1.388888888888888 1.388888888888888 0 +60 1.444444444444444 1.444444444444444 0 +61 1.5 1.5 0 +62 1.555555555555555 1.555555555555555 0 +63 1.611111111111111 1.611111111111111 0 +64 1.666666666666667 1.666666666666667 0 +65 1.722222222222222 1.722222222222222 0 +66 1.777777777777778 1.777777777777778 0 +67 1.833333333333333 1.833333333333333 0 +68 1.888888888888889 1.888888888888889 0 +69 1.944444444444444 1.944444444444444 0 +70 1.292735042735043 1.707264957264957 0 +71 1.707264957264957 1.292735042735043 0 +72 1.784037132601749 1.477004077657108 0 +73 1.215962867398251 1.522995922342892 0 +74 1.509080253312596 1.210921374546512 0 +75 1.490919746687404 1.789078625453489 0 +76 1.157094605601142 1.842575858398543 0 +77 1.842905394398859 1.157424141601458 0 +78 1.647955721656381 1.145624217732997 0 +79 1.843552075828363 1.352861733798269 0 +80 1.156447924171637 1.647138266201729 0 +81 1.352044278343618 1.854375782267003 0 +82 1.628360347171247 1.855457832670777 0 +83 1.855457832670777 1.628360347171247 0 +84 1.144542167329223 1.371639652828754 0 +85 1.371639652828753 1.144542167329223 0 +86 1.640889549620963 1.445222819474063 0 +87 1.359110450379039 1.554777180525937 0 +88 1.542280086952971 1.342687701922401 0 +89 1.457719913047029 1.657312298077602 0 +90 1.264686469823076 1.106050847713209 0 +91 1.10605084771321 1.264686469823077 0 +92 1.735313530176923 1.893949152286791 0 +93 1.893949152286791 1.735313530176922 0 +94 1.889918574771269 1.447497704616807 0 +95 1.110081425228731 1.552502295383191 0 +96 1.447497704616809 1.88991857477127 0 +97 1.55250229538319 1.11008142522873 0 +98 1.590378949319125 1.744909338559649 0 +99 1.409621050680876 1.255090661440351 0 +100 1.745981648046965 1.591497416834326 0 +101 1.254018351953037 1.408502583165676 0 +102 1.240533918361546 1.605681472424173 0 +103 1.610581881761596 1.246679000678086 0 +104 1.759466081638456 1.394318527575826 0 +105 1.389418118238403 1.753320999321916 0 +106 1.741544480719785 1.100231242248597 0 +107 1.258455519280216 1.899768757751404 0 +108 1.100661785599341 1.743538594990269 0 +109 1.899338214400659 1.25646140500973 0 +110 1.895777248266989 1.10427767439973 0 +111 1.104222751733011 1.89572232560027 0 +112 1.900755103145281 1.548497428435631 0 +113 1.54849742843563 1.90075510314528 0 +114 1.099244896854717 1.451502571564369 0 +115 1.451502571564369 1.099244896854719 0 +116 1.667140155582831 1.525906925251218 0 +117 1.332859844417172 1.474093074748781 0 +118 1.738847285916229 1.196742009829521 0 +119 1.261152714083772 1.803257990170479 0 +120 1.801650031022268 1.260904241868706 0 +121 1.198349968977732 1.739095758131293 0 +122 1.368612742054867 1.656830851793847 0 +123 1.631387257945133 1.343169148206154 0 +124 1.442028256032895 1.567192701072768 0 +125 1.557971743967108 1.432807298927234 0 +126 1.546170184125098 1.670528446152427 0 +127 1.453829815874902 1.329471553847575 0 +128 1.192556203740086 1.07975963554001 0 +129 1.079759635540011 1.192556203740086 0 +130 1.920240364459988 1.807443796259913 0 +131 1.807443796259913 1.920240364459988 0 +132 1.335028777208945 1.071701590747538 0 +133 1.071701590747538 1.335028777208946 0 +134 1.664971222791054 1.928298409252462 0 +135 1.928298409252462 1.664971222791054 0 +136 1.303795906834999 1.201338625923711 0 +137 1.201338625923711 1.303795906834999 0 +138 1.696204093165 1.798661374076288 0 +139 1.798661374076288 1.696204093165 0 +140 1.810519230109831 1.075860981249611 0 +141 1.075717774081228 1.81084663576033 0 +142 1.189480769890169 1.92413901875039 0 +143 1.924282225918772 1.189153364239668 0 +144 1.667271005647754 1.076570716988083 0 +145 1.077636898236009 1.667254308686309 0 +146 1.332728994352246 1.923429283011917 0 +147 1.922363101763991 1.332745691313689 0 +148 1.809449134198363 1.557087325894348 0 +149 1.558707229609365 1.809518315649139 0 +150 1.441292770390636 1.190481684350861 0 +151 1.190550865801637 1.442912674105652 0 +152 1.287397779691152 1.541538881919513 0 +153 1.712602220308849 1.458461118080487 0 +154 1.584316597547853 1.176543803693791 0 +155 1.415683402452147 1.823456196306209 0 +156 1.167729419705464 1.575952341096263 0 +157 1.832270580294537 1.424047658903735 0 +158 1.402056861108611 1.503136839232101 0 +159 1.597943138891394 1.496863160767897 0 +160 1.049285260441335 1.117998420235197 0 +161 1.117998420235197 1.049285260441335 0 +162 1.950714739558665 1.882001579764803 0 +163 1.882001579764803 1.950714739558664 0 +164 1.546488303604623 1.269885974351029 0 +165 1.453511696395375 1.730114025648972 0 +166 1.512843587935723 1.609938442438252 0 +167 1.487156412064275 1.390061557561752 0 +168 1.230442321291473 1.678868609247887 0 +169 1.769557678708527 1.321131390752112 0 +170 1.93348726929183 1.487174123171693 0 +171 1.487174123171694 1.933487269291831 0 +172 1.066512730708169 1.512825876828306 0 +173 1.512825876828306 1.066512730708169 0 +174 1.67491032891115 1.218691868470552 0 +175 1.325089671088849 1.781308131529449 0 +176 1.679557567570835 1.595697507169844 0 +177 1.320442432429167 1.404302492830155 0 +178 1.265931963945192 1.350442837539444 0 +179 1.348752924244887 1.265144591899373 0 +180 1.651247075755112 1.734855408100626 0 +181 1.734068036054809 1.649557162460555 0 +182 1.68854854223913 1.378873421102834 0 +183 1.311451457760871 1.621126578897167 0 +184 1.480689170290468 1.266477918749209 0 +185 1.519310829709533 1.733522081250792 0 +186 1.059969938072509 1.587522032361673 0 +187 1.587523046511419 1.059905179736265 0 +188 1.940030061927489 1.412477967638325 0 +189 1.41247695348858 1.940094820263734 0 +190 1.855088820555016 1.495370539403854 0 +191 1.498310484351775 1.853279762836327 0 +192 1.501689515648227 1.146720237163676 0 +193 1.144911179444985 1.504629460596146 0 +194 1.736870093611234 1.52845922983359 0 +195 1.263129906388763 1.471540770166412 0 +196 1.390799887040673 1.314551107582015 0 +197 1.609200112959327 1.685448892417985 0 +198 1.594128409232116 1.942163728191685 0 +199 1.40587159076788 1.057836271808317 0 +200 1.942163728191687 1.594128409232115 0 +201 1.057836271808316 1.405871590767883 0 +202 1.841706168072586 1.767439903813117 0 +203 1.767439903813117 1.841706168072585 0 +204 1.158293831927412 1.23256009618688 0 +205 1.232560096186879 1.158293831927413 0 +206 1.048916891783945 1.268106800067489 0 +207 1.268106800067489 1.048916891783944 0 +208 1.73189319993251 1.951083108216055 0 +209 1.951083108216056 1.73189319993251 0 +210 1.943336116381163 1.056663883618839 0 +211 1.056663883618839 1.943336116381163 0 +212 1.730769230769231 1.045943707292523 0 +213 1.954056292707478 1.269230769230769 0 +214 1.26923076923077 1.954056292707478 0 +215 1.045943707292522 1.73076923076923 0 +216 1.884396019814875 1.044482907084298 0 +217 1.044472592119955 1.88441020572615 0 +218 1.955527407880046 1.115589794273848 0 +219 1.115603980185125 1.955517092915703 0 +220 1.397580896177419 1.605672880351046 0 +221 1.602419103822584 1.394327119648958 0 +222 1.226842513276214 1.84445915344006 0 +223 1.773157486723788 1.155540846559942 0 +224 1.707281051595848 1.148900150304949 0 +225 1.292718948404152 1.851099849695052 0 +226 1.617075249664581 1.554191158637219 0 +227 1.382924750335423 1.44580884136278 0 +228 1.607913733349319 1.113745068675973 0 +229 1.113919253067287 1.608059587365035 0 +230 1.39208626665068 1.886254931324027 0 +231 1.886080746932711 1.391940412634963 0 +232 1.308031416720957 1.140794529799419 0 +233 1.140794529799419 1.308031416720959 0 +234 1.691968583279043 1.85920547020058 0 +235 1.859205470200581 1.691968583279041 0 +236 1.154349104529174 1.774167076324334 0 +237 1.845650895470826 1.225832923675664 0 +238 1.645458100049992 1.282762799951112 0 +239 1.354541899950008 1.717237200048889 0 +240 1.852960335633813 1.290145656873144 0 +241 1.147039664366187 1.709854343126855 0 +242 1.747132107189374 1.250239542652416 0 +243 1.252867892810626 1.749760457347583 0 +244 1.595239126062863 1.297036925021756 0 +245 1.404760873937136 1.702963074978245 0 +246 1.374998881081373 1.210565750778825 0 +247 1.790643203401739 1.624921322453663 0 +248 1.209356796598262 1.375078677546336 0 +249 1.625001118918627 1.789434249221176 0 +250 1.209075088101206 1.790025346050695 0 +251 1.790924911898794 1.209974653949303 0 +252 1.501020579521085 1.55756934217292 0 +253 1.498979420478917 1.442430657827082 0 +$EndNodes +$Elements +470 +1 1 2 2 5 1 53 +2 1 2 2 5 53 54 +3 1 2 2 5 54 55 +4 1 2 2 5 55 56 +5 1 2 2 5 56 57 +6 1 2 2 5 57 58 +7 1 2 2 5 58 59 +8 1 2 2 5 59 60 +9 1 2 2 5 60 61 +10 1 2 2 5 61 62 +11 1 2 2 5 62 63 +12 1 2 2 5 63 64 +13 1 2 2 5 64 65 +14 1 2 2 5 65 66 +15 1 2 2 5 66 67 +16 1 2 2 5 67 68 +17 1 2 2 5 68 69 +18 1 2 2 5 69 3 +19 2 2 1 1 77 110 143 +20 2 2 1 1 110 77 140 +21 2 2 1 1 76 111 141 +22 2 2 1 1 111 76 142 +23 2 2 1 1 104 72 153 +24 2 2 1 1 102 73 152 +25 2 2 1 1 72 104 157 +26 2 2 1 1 73 102 156 +27 2 2 1 1 80 102 168 +28 2 2 1 1 79 104 169 +29 2 2 1 1 102 80 156 +30 2 2 1 1 104 79 157 +31 2 2 1 1 83 112 200 +32 2 2 1 1 113 82 198 +33 2 2 1 1 84 114 201 +34 2 2 1 1 115 85 199 +35 2 2 1 1 74 154 164 +36 2 2 1 1 75 155 165 +37 2 2 1 1 154 103 164 +38 2 2 1 1 155 105 165 +39 2 2 1 1 144 106 224 +40 2 2 1 1 146 107 225 +41 2 2 1 1 89 126 185 +42 2 2 1 1 88 127 184 +43 2 2 1 1 165 89 185 +44 2 2 1 1 164 88 184 +45 2 2 1 1 61 125 159 +46 2 2 1 1 61 124 158 +47 2 2 1 1 109 147 240 +48 2 2 1 1 108 145 241 +49 2 2 1 1 135 83 200 +50 2 2 1 1 82 134 198 +51 2 2 1 1 133 84 201 +52 2 2 1 1 85 132 199 +53 2 2 1 1 141 43 215 +54 2 2 1 1 38 142 214 +55 2 2 1 1 14 140 212 +56 2 2 1 1 143 19 213 +57 2 2 1 1 71 169 182 +58 2 2 1 1 70 168 183 +59 2 2 1 1 106 144 212 +60 2 2 1 1 145 108 215 +61 2 2 1 1 147 109 213 +62 2 2 1 1 107 146 214 +63 2 2 1 1 112 83 148 +64 2 2 1 1 85 115 150 +65 2 2 1 1 82 113 149 +66 2 2 1 1 114 84 151 +67 2 2 1 1 62 61 159 +68 2 2 1 1 60 61 158 +69 2 2 1 1 34 113 198 +70 2 2 1 1 112 23 200 +71 2 2 1 1 10 115 199 +72 2 2 1 1 114 47 201 +73 2 2 1 1 113 34 171 +74 2 2 1 1 23 112 170 +75 2 2 1 1 115 10 173 +76 2 2 1 1 47 114 172 +77 2 2 1 1 63 126 166 +78 2 2 1 1 59 127 167 +79 2 2 1 1 117 87 152 +80 2 2 1 1 116 86 153 +81 2 2 1 1 5 53 1 +82 2 2 1 1 28 3 69 +83 2 2 1 1 1 53 52 +84 2 2 1 1 69 3 29 +85 2 2 1 1 87 117 158 +86 2 2 1 1 86 116 159 +87 2 2 1 1 140 77 223 +88 2 2 1 1 142 76 222 +89 2 2 1 1 169 104 182 +90 2 2 1 1 168 102 183 +91 2 2 1 1 108 141 215 +92 2 2 1 1 140 106 212 +93 2 2 1 1 142 107 214 +94 2 2 1 1 109 143 213 +95 2 2 1 1 42 43 141 +96 2 2 1 1 14 15 140 +97 2 2 1 1 38 39 142 +98 2 2 1 1 18 19 143 +99 2 2 1 1 77 143 237 +100 2 2 1 1 76 141 236 +101 2 2 1 1 67 131 203 +102 2 2 1 1 130 67 202 +103 2 2 1 1 55 128 205 +104 2 2 1 1 129 55 204 +105 2 2 1 1 99 150 184 +106 2 2 1 1 98 149 185 +107 2 2 1 1 62 63 166 +108 2 2 1 1 60 59 167 +109 2 2 1 1 78 144 224 +110 2 2 1 1 81 146 225 +111 2 2 1 1 67 68 131 +112 2 2 1 1 68 67 130 +113 2 2 1 1 54 55 129 +114 2 2 1 1 55 54 128 +115 2 2 1 1 102 152 183 +116 2 2 1 1 104 153 182 +117 2 2 1 1 123 71 182 +118 2 2 1 1 122 70 183 +119 2 2 1 1 150 74 184 +120 2 2 1 1 149 75 185 +121 2 2 1 1 6 7 128 +122 2 2 1 1 50 51 129 +123 2 2 1 1 26 27 130 +124 2 2 1 1 30 31 131 +125 2 2 1 1 80 168 241 +126 2 2 1 1 79 169 240 +127 2 2 1 1 131 92 203 +128 2 2 1 1 93 130 202 +129 2 2 1 1 128 90 205 +130 2 2 1 1 91 129 204 +131 2 2 1 1 151 101 195 +132 2 2 1 1 148 100 194 +133 2 2 1 1 145 80 241 +134 2 2 1 1 147 79 240 +135 2 2 1 1 154 78 174 +136 2 2 1 1 155 81 175 +137 2 2 1 1 106 140 223 +138 2 2 1 1 107 142 222 +139 2 2 1 1 103 154 174 +140 2 2 1 1 105 155 175 +141 2 2 1 1 64 176 181 +142 2 2 1 1 58 177 178 +143 2 2 1 1 71 123 238 +144 2 2 1 1 70 122 239 +145 2 2 1 1 128 7 207 +146 2 2 1 1 50 129 206 +147 2 2 1 1 26 130 209 +148 2 2 1 1 131 31 208 +149 2 2 1 1 152 87 183 +150 2 2 1 1 153 86 182 +151 2 2 1 1 48 49 133 +152 2 2 1 1 24 25 135 +153 2 2 1 1 8 9 132 +154 2 2 1 1 32 33 134 +155 2 2 1 1 169 120 240 +156 2 2 1 1 168 121 241 +157 2 2 1 1 141 108 236 +158 2 2 1 1 143 109 237 +159 2 2 1 1 72 157 190 +160 2 2 1 1 155 75 191 +161 2 2 1 1 73 156 193 +162 2 2 1 1 154 74 192 +163 2 2 1 1 72 148 194 +164 2 2 1 1 73 151 195 +165 2 2 1 1 88 125 167 +166 2 2 1 1 89 124 166 +167 2 2 1 1 124 87 158 +168 2 2 1 1 125 86 159 +169 2 2 1 1 70 175 243 +170 2 2 1 1 71 174 242 +171 2 2 1 1 93 135 209 +172 2 2 1 1 91 133 206 +173 2 2 1 1 134 92 208 +174 2 2 1 1 132 90 207 +175 2 2 1 1 56 57 137 +176 2 2 1 1 57 56 136 +177 2 2 1 1 66 65 139 +178 2 2 1 1 65 66 138 +179 2 2 1 1 88 123 221 +180 2 2 1 1 89 122 220 +181 2 2 1 1 177 101 178 +182 2 2 1 1 176 100 181 +183 2 2 1 1 90 128 207 +184 2 2 1 1 92 131 208 +185 2 2 1 1 129 91 206 +186 2 2 1 1 130 93 209 +187 2 2 1 1 36 37 146 +188 2 2 1 1 20 21 147 +189 2 2 1 1 12 13 144 +190 2 2 1 1 44 45 145 +191 2 2 1 1 157 94 190 +192 2 2 1 1 96 155 191 +193 2 2 1 1 156 95 193 +194 2 2 1 1 97 154 192 +195 2 2 1 1 135 93 235 +196 2 2 1 1 92 134 234 +197 2 2 1 1 133 91 233 +198 2 2 1 1 90 132 232 +199 2 2 1 1 123 88 244 +200 2 2 1 1 122 89 245 +201 2 2 1 1 177 117 195 +202 2 2 1 1 176 116 194 +203 2 2 1 1 125 61 253 +204 2 2 1 1 124 61 252 +205 2 2 1 1 22 170 188 +206 2 2 1 1 46 172 186 +207 2 2 1 1 171 35 189 +208 2 2 1 1 173 11 187 +209 2 2 1 1 172 114 193 +210 2 2 1 1 115 173 192 +211 2 2 1 1 170 112 190 +212 2 2 1 1 113 171 191 +213 2 2 1 1 101 177 195 +214 2 2 1 1 100 176 194 +215 2 2 1 1 95 172 193 +216 2 2 1 1 173 97 192 +217 2 2 1 1 94 170 190 +218 2 2 1 1 171 96 191 +219 2 2 1 1 125 88 221 +220 2 2 1 1 124 89 220 +221 2 2 1 1 54 129 160 +222 2 2 1 1 128 54 161 +223 2 2 1 1 131 68 163 +224 2 2 1 1 68 130 162 +225 2 2 1 1 87 124 220 +226 2 2 1 1 86 125 221 +227 2 2 1 1 127 88 167 +228 2 2 1 1 126 89 166 +229 2 2 1 1 6 128 161 +230 2 2 1 1 129 51 160 +231 2 2 1 1 130 27 162 +232 2 2 1 1 30 131 163 +233 2 2 1 1 54 53 161 +234 2 2 1 1 53 54 160 +235 2 2 1 1 68 69 163 +236 2 2 1 1 69 68 162 +237 2 2 1 1 135 25 209 +238 2 2 1 1 133 49 206 +239 2 2 1 1 8 132 207 +240 2 2 1 1 32 134 208 +241 2 2 1 1 175 119 243 +242 2 2 1 1 174 118 242 +243 2 2 1 1 126 63 197 +244 2 2 1 1 127 59 196 +245 2 2 1 1 175 70 239 +246 2 2 1 1 174 71 238 +247 2 2 1 1 127 99 184 +248 2 2 1 1 126 98 185 +249 2 2 1 1 204 137 233 +250 2 2 1 1 136 205 232 +251 2 2 1 1 202 139 235 +252 2 2 1 1 138 203 234 +253 2 2 1 1 56 137 204 +254 2 2 1 1 136 56 205 +255 2 2 1 1 66 139 202 +256 2 2 1 1 138 66 203 +257 2 2 1 1 99 127 196 +258 2 2 1 1 98 126 197 +259 2 2 1 1 29 30 163 +260 2 2 1 1 27 28 162 +261 2 2 1 1 51 52 160 +262 2 2 1 1 5 6 161 +263 2 2 1 1 144 13 212 +264 2 2 1 1 146 37 214 +265 2 2 1 1 20 147 213 +266 2 2 1 1 44 145 215 +267 2 2 1 1 64 65 180 +268 2 2 1 1 57 58 178 +269 2 2 1 1 58 57 179 +270 2 2 1 1 65 64 181 +271 2 2 1 1 170 94 188 +272 2 2 1 1 172 95 186 +273 2 2 1 1 97 173 187 +274 2 2 1 1 96 171 189 +275 2 2 1 1 134 82 234 +276 2 2 1 1 132 85 232 +277 2 2 1 1 84 133 233 +278 2 2 1 1 83 135 235 +279 2 2 1 1 9 10 199 +280 2 2 1 1 47 48 201 +281 2 2 1 1 33 34 198 +282 2 2 1 1 23 24 200 +283 2 2 1 1 21 22 188 +284 2 2 1 1 35 36 189 +285 2 2 1 1 45 46 186 +286 2 2 1 1 11 12 187 +287 2 2 1 1 64 180 197 +288 2 2 1 1 58 179 196 +289 2 2 1 1 122 183 220 +290 2 2 1 1 123 182 221 +291 2 2 1 1 16 2 210 +292 2 2 1 1 4 41 211 +293 2 2 1 1 40 4 211 +294 2 2 1 1 2 17 210 +295 2 2 1 1 183 87 220 +296 2 2 1 1 182 86 221 +297 2 2 1 1 101 151 248 +298 2 2 1 1 150 99 246 +299 2 2 1 1 149 98 249 +300 2 2 1 1 100 148 247 +301 2 2 1 1 34 35 171 +302 2 2 1 1 22 23 170 +303 2 2 1 1 10 11 173 +304 2 2 1 1 46 47 172 +305 2 2 1 1 28 69 162 +306 2 2 1 1 69 29 163 +307 2 2 1 1 52 53 160 +308 2 2 1 1 53 5 161 +309 2 2 1 1 57 136 179 +310 2 2 1 1 137 57 178 +311 2 2 1 1 65 138 180 +312 2 2 1 1 139 65 181 +313 2 2 1 1 138 82 249 +314 2 2 1 1 136 85 246 +315 2 2 1 1 83 139 247 +316 2 2 1 1 84 137 248 +317 2 2 1 1 132 9 199 +318 2 2 1 1 48 133 201 +319 2 2 1 1 24 135 200 +320 2 2 1 1 134 33 198 +321 2 2 1 1 205 90 232 +322 2 2 1 1 91 204 233 +323 2 2 1 1 93 202 235 +324 2 2 1 1 203 92 234 +325 2 2 1 1 64 63 176 +326 2 2 1 1 58 59 177 +327 2 2 1 1 63 64 197 +328 2 2 1 1 59 58 196 +329 2 2 1 1 142 39 219 +330 2 2 1 1 140 15 216 +331 2 2 1 1 18 143 218 +332 2 2 1 1 42 141 217 +333 2 2 1 1 110 140 216 +334 2 2 1 1 143 110 218 +335 2 2 1 1 111 142 219 +336 2 2 1 1 141 111 217 +337 2 2 1 1 167 125 253 +338 2 2 1 1 166 124 252 +339 2 2 1 1 36 146 189 +340 2 2 1 1 147 21 188 +341 2 2 1 1 12 144 187 +342 2 2 1 1 145 45 186 +343 2 2 1 1 179 99 196 +344 2 2 1 1 180 98 197 +345 2 2 1 1 112 148 190 +346 2 2 1 1 150 115 192 +347 2 2 1 1 114 151 193 +348 2 2 1 1 149 113 191 +349 2 2 1 1 151 84 248 +350 2 2 1 1 85 150 246 +351 2 2 1 1 82 149 249 +352 2 2 1 1 148 83 247 +353 2 2 1 1 82 138 234 +354 2 2 1 1 85 136 232 +355 2 2 1 1 137 84 233 +356 2 2 1 1 139 83 235 +357 2 2 1 1 175 81 225 +358 2 2 1 1 174 78 224 +359 2 2 1 1 119 175 225 +360 2 2 1 1 118 174 224 +361 2 2 1 1 117 152 195 +362 2 2 1 1 116 153 194 +363 2 2 1 1 144 78 228 +364 2 2 1 1 80 145 229 +365 2 2 1 1 146 81 230 +366 2 2 1 1 79 147 231 +367 2 2 1 1 75 165 185 +368 2 2 1 1 74 164 184 +369 2 2 1 1 148 72 190 +370 2 2 1 1 75 149 191 +371 2 2 1 1 74 150 192 +372 2 2 1 1 151 73 193 +373 2 2 1 1 117 177 227 +374 2 2 1 1 116 176 226 +375 2 2 1 1 177 59 227 +376 2 2 1 1 176 63 226 +377 2 2 1 1 103 174 238 +378 2 2 1 1 105 175 239 +379 2 2 1 1 153 72 194 +380 2 2 1 1 152 73 195 +381 2 2 1 1 66 67 203 +382 2 2 1 1 67 66 202 +383 2 2 1 1 55 56 204 +384 2 2 1 1 56 55 205 +385 2 2 1 1 31 32 208 +386 2 2 1 1 49 50 206 +387 2 2 1 1 7 8 207 +388 2 2 1 1 25 26 209 +389 2 2 1 1 13 14 212 +390 2 2 1 1 19 20 213 +391 2 2 1 1 43 44 215 +392 2 2 1 1 37 38 214 +393 2 2 1 1 156 80 229 +394 2 2 1 1 81 155 230 +395 2 2 1 1 78 154 228 +396 2 2 1 1 157 79 231 +397 2 2 1 1 95 156 229 +398 2 2 1 1 94 157 231 +399 2 2 1 1 154 97 228 +400 2 2 1 1 155 96 230 +401 2 2 1 1 41 42 217 +402 2 2 1 1 15 16 216 +403 2 2 1 1 17 18 218 +404 2 2 1 1 39 40 219 +405 2 2 1 1 63 62 226 +406 2 2 1 1 59 60 227 +407 2 2 1 1 61 60 253 +408 2 2 1 1 61 62 252 +409 2 2 1 1 179 136 246 +410 2 2 1 1 180 138 249 +411 2 2 1 1 139 181 247 +412 2 2 1 1 137 178 248 +413 2 2 1 1 89 165 245 +414 2 2 1 1 88 164 244 +415 2 2 1 1 120 169 242 +416 2 2 1 1 121 168 243 +417 2 2 1 1 169 71 242 +418 2 2 1 1 168 70 243 +419 2 2 1 1 158 117 227 +420 2 2 1 1 159 116 226 +421 2 2 1 1 62 159 226 +422 2 2 1 1 60 158 227 +423 2 2 1 1 164 103 244 +424 2 2 1 1 165 105 245 +425 2 2 1 1 189 146 230 +426 2 2 1 1 187 144 228 +427 2 2 1 1 147 188 231 +428 2 2 1 1 145 186 229 +429 2 2 1 1 62 166 252 +430 2 2 1 1 60 167 253 +431 2 2 1 1 178 101 248 +432 2 2 1 1 181 100 247 +433 2 2 1 1 99 179 246 +434 2 2 1 1 98 180 249 +435 2 2 1 1 186 95 229 +436 2 2 1 1 188 94 231 +437 2 2 1 1 97 187 228 +438 2 2 1 1 96 189 230 +439 2 2 1 1 237 109 240 +440 2 2 1 1 236 108 241 +441 2 2 1 1 120 237 240 +442 2 2 1 1 121 236 241 +443 2 2 1 1 211 41 217 +444 2 2 1 1 16 210 216 +445 2 2 1 1 40 211 219 +446 2 2 1 1 210 17 218 +447 2 2 1 1 210 110 216 +448 2 2 1 1 111 211 217 +449 2 2 1 1 110 210 218 +450 2 2 1 1 211 111 219 +451 2 2 1 1 106 223 224 +452 2 2 1 1 107 222 225 +453 2 2 1 1 222 119 225 +454 2 2 1 1 223 118 224 +455 2 2 1 1 120 242 251 +456 2 2 1 1 121 243 250 +457 2 2 1 1 243 119 250 +458 2 2 1 1 242 118 251 +459 2 2 1 1 222 76 250 +460 2 2 1 1 223 77 251 +461 2 2 1 1 119 222 250 +462 2 2 1 1 118 223 251 +463 2 2 1 1 76 236 250 +464 2 2 1 1 77 237 251 +465 2 2 1 1 103 238 244 +466 2 2 1 1 105 239 245 +467 2 2 1 1 237 120 251 +468 2 2 1 1 236 121 250 +469 2 2 1 1 238 123 244 +470 2 2 1 1 239 122 245 +$EndElements diff --git a/src/mmdg.cc b/src/mmdg.cc index 20d7d6de0762ab0e1d9ae12edd65b5b65976439d..894c3044684293d6a78f6c7bdf46c0387cd4122e 100644 --- a/src/mmdg.cc +++ b/src/mmdg.cc @@ -24,6 +24,7 @@ #include <dune/mmdg/problems/mmdgproblem2.hh> #include <dune/mmdg/problems/mmdgproblem3.hh> #include <dune/mmdg/problems/mmdgproblem4.hh> +#include <dune/mmdg/problems/mmdgproblem5.hh> int main(int argc, char** argv) { @@ -82,6 +83,11 @@ int main(int argc, char** argv) problem = new MMDGProblem4<Coordinate>(aperture); gridType = "mmdg2_5e-2"; } + else if (pt["problem"] == "mmdg5") + { + problem = new MMDGProblem5<Coordinate>(aperture); + gridType = "mmdg5_5e-2"; + } else { DUNE_THROW(Dune::Exception, diff --git a/src/mmdgAnalysis.py b/src/mmdgAnalysis.py index c66dd514341f1a88a1b75593be2f175655f24a4b..89c1f6e9214ab090e843489e82ed89c44bdd929f 100644 --- a/src/mmdgAnalysis.py +++ b/src/mmdgAnalysis.py @@ -63,6 +63,7 @@ numberOfMeans = 1 # names of the grid files sorted by refinement gridfiles = ["mmdg2_1", "mmdg2_5e-1", "mmdg2_2e-1", "mmdg2_8e-2", "mmdg2_5e-2"] # gridfiles = ["mmdg3_1", "mmdg3_5e-1", "mmdg3_2e-1", "mmdg3_8e-2", "mmdg3_5e-2"] +# gridfiles = ["mmdg5_1", "mmdg5_5e-1", "mmdg5_2e-1", "mmdg5_8e-2", "mmdg5_5e-2"] dim = 2 diff --git a/src/parameterMMDG.ini b/src/parameterMMDG.ini index a6283fc3114613f0a0e3f66adc1fcb6338c5992e..d214f65a5b92af753727c8050b00ed0444bc76c1 100644 --- a/src/parameterMMDG.ini +++ b/src/parameterMMDG.ini @@ -1,5 +1,5 @@ mu = 1000 xi = 0.75 +d = 1 problem = mmdg2 gridfile = mmdg2_5e-2 -d = 1