From ce573fe94c2962b7e5430cb111ae80316a4bfc9d Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Maximilian=20H=C3=B6rl?=
<maximilian.hoerl@mathematik.uni-stuttgart.de>
Date: Tue, 17 Mar 2020 14:50:36 +0100
Subject: [PATCH] add vtk output to mmdg
---
dune/mmdg/dg.hh | 19 +-
dune/mmdg/mmdg.hh | 83 +-
dune/mmdg/problems/coupleddgproblem.hh | 12 +-
src/grids/sphere2d.msh | 1744 ++++++++++++++++++++++++
src/mmdg.cc | 4 +-
5 files changed, 1844 insertions(+), 18 deletions(-)
create mode 100644 src/grids/sphere2d.msh
diff --git a/dune/mmdg/dg.hh b/dune/mmdg/dg.hh
index 83676e7..b82886f 100644
--- a/dune/mmdg/dg.hh
+++ b/dune/mmdg/dg.hh
@@ -38,7 +38,7 @@ public:
static constexpr auto edge = Dune::GeometryTypes::simplex(dim-1);
//constructor
- DG (const Grid& grid, const GridView& gridView, const Mapper mapper, const Problem& problem) :
+ DG (const Grid& grid, const GridView& gridView, const Mapper& mapper, const Problem& problem) :
//NOTE: why does this not work?
//grid_(grid), gridView_(grid.leafGridView()), mapper_(*(new Mapper(gridView_, Dune::mcmgElementLayout()))), problem_(problem), dof_((1 + dim) * gridView_.size(0))
grid_(grid), gridView_(gridView), mapper_(mapper), problem_(problem),
@@ -72,7 +72,7 @@ public:
const int dof_;
protected:
- DG (const Grid& grid, const GridView& gridView, const Mapper mapper,
+ DG (const Grid& grid, const GridView& gridView, const Mapper& mapper,
const Problem& problem, const int dof) :
grid_(grid), gridView_(gridView), mapper_(mapper), problem_(problem),
dof_(dof)
@@ -451,13 +451,13 @@ protected:
}
//writes the solution to a vtk file
- void writeVTKoutput () const
+ void writeVTKoutput (const int bulkDOF) const
{
//storage for pressure data, for each element we store the
//pressure at the corners of the element, the VTKFunction will
//create the output using a linear interpolation for each element
- Dune::DynamicVector<Scalar> pressure(dof_, 0.0);
- Dune::DynamicVector<Scalar> exactPressure(dof_, 0.0);
+ Dune::DynamicVector<Scalar> pressure(bulkDOF, 0.0);
+ Dune::DynamicVector<Scalar> exactPressure(bulkDOF, 0.0);
for (const auto& elem : elements(gridView_))
{
@@ -502,10 +502,17 @@ protected:
vtkWriter.pwrite("pressureData", "", "data");
}
-
std::shared_ptr<Matrix> A; //stiffness matrix
Vector b; //load vector
Vector d; //solution vector
+
+private:
+ //writes the solution to a vtk file
+ void writeVTKoutput () const
+ {
+ writeVTKoutput(dof_);
+ }
+
};
#endif
diff --git a/dune/mmdg/mmdg.hh b/dune/mmdg/mmdg.hh
index 5892fa8..69a48fe 100644
--- a/dune/mmdg/mmdg.hh
+++ b/dune/mmdg/mmdg.hh
@@ -9,6 +9,7 @@ class MMDG : public DG<Grid, GridView, Mapper, Problem>
{
public:
static constexpr int dim = GridView::dimension;
+ static constexpr auto iSimplex = Dune::GeometryTypes::simplex(dim-1);
using Base = DG<Grid, GridView, Mapper, Problem>;
using Scalar = typename Base::Scalar; //NOTE using Base::Scalar
@@ -18,7 +19,9 @@ public:
using InterfaceBasis = std::array<Coordinate, dim-1>;
using QRuleInterface = typename Base::QRuleEdge;
using QRulesInterface = typename Base::QRulesEdge;
- static constexpr auto iSimplex = Dune::GeometryTypes::simplex(dim-1);
+ using VTKFunction = Dune::VTKFunction<InterfaceGridView>;
+ using P1Function = Dune::NonconformingP1VTKFunction<InterfaceGridView,
+ Dune::DynamicVector<Scalar>>;
//constructor
MMDG (const Grid& grid, const GridView& gridView, const Mapper& mapper,
@@ -32,6 +35,13 @@ public:
{
assembleSLE(mu, xi);
Base::A->compress();
+
+ Dune::InverseOperatorResult result;
+ Dune::UMFPack<Matrix> solver(*Base::A);
+ solver.apply(Base::d, Base::b, result);
+
+ //write solution to a vtk file
+ writeVTKoutput();
}
const InterfaceGridView& iGridView_;
@@ -96,8 +106,8 @@ private:
//get basis of the interface coordinate system for evaluation of
//basis function phi_iElem,i with i=1,...,dim-1
- const InterfaceBasis iFrame = interfaceFrame(
- Base::grid_.asIntersection(iElem).centerUnitOuterNormal());
+ const InterfaceBasis iFrame =
+ interfaceFrame(intersct.centerUnitOuterNormal());
//in the total system of linear equations (SLE) Ad = b,
//the index iElemIdxSLE refers to the basis function (0, phi_iElem,0)
@@ -558,7 +568,7 @@ private:
//3d: tau_1 ~ [ normal[1] + normal[2], tau_2 ~ [ -normal[1],
// -normal[0], normal[0] + normal[2]
// -normal[0] ], -normal[1] ]
- InterfaceBasis interfaceFrame(const Coordinate& normal)
+ InterfaceBasis interfaceFrame(const Coordinate& normal) const
{
InterfaceBasis iBasis;
@@ -587,6 +597,71 @@ private:
return iBasis;
}
+ void writeVTKoutput () const
+ {
+ //degrees of freedom of bulk and interface grids
+ const int bulkDOF = (dim + 1) * Base::gridView_.size(0);
+ const int interfaceDOF = dim * iGridView_.size(0);
+
+ //=== write bulk solution to file ===
+ Base::writeVTKoutput(bulkDOF);
+
+ //=== write interface solution to file ===
+
+ //storage for pressure data, for each interface element we store the
+ //pressure at the corners of the interface element, the VTKFunction will
+ //create the output using a linear interpolation for each interface element
+ Dune::DynamicVector<Scalar> iPressure(interfaceDOF, 0.0);
+ Dune::DynamicVector<Scalar> exactIPressure(interfaceDOF, 0.0);
+
+ for (const auto& iElem : elements(iGridView_))
+ {
+ const int iElemIdxSLE = bulkDOF + dim * iMapper_.index(iElem);
+ const auto& iGeo = iElem.geometry();
+
+ //get basis of the interface coordinate system for evaluation of
+ //basis function phi_iElem,i with i=1,...,dim-1
+ const InterfaceBasis iFrame = interfaceFrame(
+ Base::grid_.asIntersection(iElem).centerUnitOuterNormal());
+
+ for (int k = 0; k < iGeo.corners(); k++)
+ {
+ if (Base::problem_.hasExactSolution())
+ {
+ exactIPressure[iElemIdxSLE + k] =
+ Base::problem_.exactInterfaceSolution(iGeo.corner(k));
+ }
+
+ //contribution of the basis function
+ // phi_iElem,0 (x) = indicator(iElem);
+ //at the kth corner of iElem
+ iPressure[iElemIdxSLE + k] = Base::d[iElemIdxSLE];
+
+ for (int i = 0; i < dim - 1; i++)
+ {
+ //contribution of the basis function
+ // phi_iElem,i (x) = (x * tau[i]) * indicator(iElem);
+ //at the kth corner of iElem
+ iPressure[iElemIdxSLE + k] +=
+ Base::d[iElemIdxSLE + i + 1] * (iGeo.corner(k) * iFrame[i]);
+ }
+ }
+ }
+
+ Dune::VTKWriter<InterfaceGridView>
+ vtkWriter(iGridView_, Dune::VTK::nonconforming);
+ vtkWriter.addVertexData( std::shared_ptr<const VTKFunction>(
+ new P1Function(iGridView_, iPressure, "interfacePressure")));
+
+ if (Base::problem_.hasExactSolution())
+ {
+ vtkWriter.addVertexData( std::shared_ptr<const VTKFunction>(
+ new P1Function(iGridView_, exactIPressure, "exactInterfacePressure")));
+ }
+
+ vtkWriter.pwrite("interfacePressureData", "", "data");
+ }
+
};
#endif
diff --git a/dune/mmdg/problems/coupleddgproblem.hh b/dune/mmdg/problems/coupleddgproblem.hh
index 6e23c89..e6f8d6e 100644
--- a/dune/mmdg/problems/coupleddgproblem.hh
+++ b/dune/mmdg/problems/coupleddgproblem.hh
@@ -14,6 +14,12 @@ class CoupledDGProblem : public DGProblem<Vector, Scalar>
using Base = DGProblem<Vector, Scalar>;
using Matrix = typename Base::Matrix;
+ //the exact solution on the interface at position pos
+ virtual Scalar exactInterfaceSolution (const Vector& pos) const
+ {
+ return Scalar(0.0);
+ }
+
//interface source term at position pos
virtual Scalar qInterface (const Vector& pos) const
{
@@ -26,12 +32,6 @@ class CoupledDGProblem : public DGProblem<Vector, Scalar>
return Scalar(1.0);
}
- //(tangential) gradient of the aperture at position pos
- virtual Vector gradAperture (const Vector& pos) const
- {
- return Vector(0.0);
- }
-
//tangential permeability tensor of the interface at position pos
virtual Matrix Kparallel (const Vector& pos) const
{
diff --git a/src/grids/sphere2d.msh b/src/grids/sphere2d.msh
new file mode 100644
index 0000000..f24bcc9
--- /dev/null
+++ b/src/grids/sphere2d.msh
@@ -0,0 +1,1744 @@
+$MeshFormat
+2.2 0 8
+$EndMeshFormat
+$Nodes
+599
+1 0 0 0
+2 1 0 0
+3 1 1 0
+4 0 1 0
+5 0.5 0.9 0
+6 0.35 0.75 0
+7 0.5 0.6 0
+8 0.65 0.75 0
+9 0.04999999999989965 0 0
+10 0.09999999999981467 0 0
+11 0.1499999999997036 0 0
+12 0.1999999999995579 0 0
+13 0.2499999999994121 0 0
+14 0.2999999999992664 0 0
+15 0.3499999999991206 0 0
+16 0.3999999999989749 0 0
+17 0.4499999999988292 0 0
+18 0.4999999999986943 0 0
+19 0.5499999999988151 0 0
+20 0.5999999999989468 0 0
+21 0.6499999999990784 0 0
+22 0.69999999999921 0 0
+23 0.7499999999993417 0 0
+24 0.7999999999994734 0 0
+25 0.8499999999996051 0 0
+26 0.8999999999997368 0 0
+27 0.9499999999998684 0 0
+28 1 0.04999999999989965 0
+29 1 0.09999999999981467 0
+30 1 0.1499999999997036 0
+31 1 0.1999999999995579 0
+32 1 0.2499999999994121 0
+33 1 0.2999999999992664 0
+34 1 0.3499999999991206 0
+35 1 0.3999999999989749 0
+36 1 0.4499999999988292 0
+37 1 0.4999999999986943 0
+38 1 0.5499999999988151 0
+39 1 0.5999999999989468 0
+40 1 0.6499999999990784 0
+41 1 0.69999999999921 0
+42 1 0.7499999999993417 0
+43 1 0.7999999999994734 0
+44 1 0.8499999999996051 0
+45 1 0.8999999999997368 0
+46 1 0.9499999999998684 0
+47 0.9499999999997918 1 0
+48 0.8999999999995836 1 0
+49 0.8499999999996529 1 0
+50 0.7999999999999998 1 0
+51 0.7500000000003466 1 0
+52 0.7000000000006934 1 0
+53 0.6500000000010401 1 0
+54 0.6000000000013869 1 0
+55 0.5500000000017335 1 0
+56 0.5000000000020587 1 0
+57 0.4500000000018723 1 0
+58 0.4000000000016644 1 0
+59 0.3500000000014564 1 0
+60 0.3000000000012483 1 0
+61 0.2500000000010403 1 0
+62 0.2000000000008322 1 0
+63 0.1500000000006241 1 0
+64 0.100000000000416 1 0
+65 0.05000000000020799 1 0
+66 0 0.9499999999997918 0
+67 0 0.8999999999995836 0
+68 0 0.8499999999996529 0
+69 0 0.7999999999999998 0
+70 0 0.7500000000003466 0
+71 0 0.7000000000006934 0
+72 0 0.6500000000010401 0
+73 0 0.6000000000013869 0
+74 0 0.5500000000017335 0
+75 0 0.5000000000020587 0
+76 0 0.4500000000018723 0
+77 0 0.4000000000016644 0
+78 0 0.3500000000014564 0
+79 0 0.3000000000012483 0
+80 0 0.2500000000010403 0
+81 0 0.2000000000008322 0
+82 0 0.1500000000006241 0
+83 0 0.100000000000416 0
+84 0 0.05000000000020799 0
+85 0.453647450780622 0.8926584774237589 0
+86 0.4118322119679112 0.8713525490194939 0
+87 0.3786474507039843 0.838167787651489 0
+88 0.3573415225082845 0.7963525490102292 0
+89 0.3573415225779883 0.7036474507752443 0
+90 0.3786474509488204 0.661832212011523 0
+91 0.4118322123654675 0.6286474506916646 0
+92 0.4536474509691241 0.607341522514993 0
+93 0.5463525492164328 0.6073415225752841 0
+94 0.5881677879880485 0.6286474509485089 0
+95 0.6213525492671707 0.6618322123088091 0
+96 0.6426584774848487 0.7036474509686369 0
+97 0.6426584774247159 0.7963525492164328 0
+98 0.6213525490514911 0.8381677879880485 0
+99 0.5881677876911909 0.8713525492671707 0
+100 0.5463525490313631 0.8926584774848487 0
+101 0.3189945241507618 0.3226544027841542 0
+102 0.6809451753627315 0.3216309813694689 0
+103 0.7805334726857065 0.5133890875186577 0
+104 0.2183677555839037 0.5129602409185885 0
+105 0.5016439734449822 0.2081018505106127 0
+106 0.4974562453527041 0.4144703900691039 0
+107 0.1863308474258644 0.1863308474259515 0
+108 0.8136691525801667 0.1863308474197672 0
+109 0.1783722942414739 0.8012703432443191 0
+110 0.8213015477609455 0.6991108343266341 0
+111 0.1632611783956815 0.3574406982371006 0
+112 0.837240224019091 0.3550135516200195 0
+113 0.3472452761822939 0.1595219365168092 0
+114 0.6524343282655077 0.1605429510991677 0
+115 0.6222892879088412 0.4792881469340461 0
+116 0.4999999999999952 0.7500000000000011 0
+117 0.3735377092606734 0.4869360140506486 0
+118 0.7898318142287654 0.8467603203808851 0
+119 0.1478820549245306 0.6456000630487561 0
+120 0.8745223014858544 0.5883422466248326 0
+121 0.7324633546320717 0.6207404191435332 0
+122 0.2671093988320164 0.8815511761092146 0
+123 0.2656817742754761 0.621013759861046 0
+124 0.5634935229448634 0.3108099532405584 0
+125 0.4253346695637267 0.3113921947150382 0
+126 0.8869945694868784 0.4634762130147316 0
+127 0.1165580953314957 0.4621773835834631 0
+128 0.8862559425024215 0.7927644037853313 0
+129 0.2679036265989492 0.4296960361924566 0
+130 0.7233959523143316 0.421127840485136 0
+131 0.4470235716015805 0.11309562190905 0
+132 0.1188668157297 0.8829432820388647 0
+133 0.1073347617383978 0.2635640234119433 0
+134 0.8926652382547051 0.2635640234288273 0
+135 0.1027723608102328 0.1117912262095486 0
+136 0.8972276391917119 0.1117912262093265 0
+137 0.5594137712199742 0.1072546418762209 0
+138 0.2475927791169805 0.7195969072975297 0
+139 0.5265702225016204 0.5059358183417253 0
+140 0.2546951476327677 0.1087192772254437 0
+141 0.7452637437460657 0.1089852422020708 0
+142 0.8756350681059081 0.8964377270299716 0
+143 0.2269965668654961 0.2775294316963477 0
+144 0.771490431353876 0.2763963684049586 0
+145 0.7024280031570532 0.900707575645413 0
+146 0.09866512435377567 0.737978102258491 0
+147 0.09544738889438467 0.5659290206401522 0
+148 0.3845739835541037 0.3904076301375724 0
+149 0.7461691414887199 0.7599620234098473 0
+150 0.6572371102154313 0.5673057753427317 0
+151 0.6169921530277097 0.3765321598896578 0
+152 0.9114521970759928 0.6742816671079419 0
+153 0.3414436824425633 0.5724743506861438 0
+154 0.6351720850975711 0.2451236672873143 0
+155 0.4495750388406308 0.5193419111988838 0
+156 0.3593991471830846 0.244207268115211 0
+157 0.2747393284570373 0.2083291992076511 0
+158 0.7232884812288243 0.2091370455115789 0
+159 0.3517712746283457 0.9242834039054917 0
+160 0.1853252696370091 0.9135709213536087 0
+161 0.2680464157665426 0.7958961097422692 0
+162 0.08115114570433962 0.3891507218777193 0
+163 0.9152100140347055 0.3891680421001668 0
+164 0.3672863835869281 0.08012508261734277 0
+165 0.1973572211544855 0.4342339915220699 0
+166 0.8040453101058151 0.4347551240039552 0
+167 0.07644291682884963 0.1866238622210362 0
+168 0.9235570831677304 0.1866238622273756 0
+169 0.6358893637253391 0.08165077921021982 0
+170 0.07924524015977584 0.814967045448018 0
+171 0.4239179540535095 0.1972979398289636 0
+172 0.7052869630975451 0.5021422883726643 0
+173 0.5736714898933901 0.1904704214619407 0
+174 0.4249999999999975 0.7618788329849965 0
+175 0.4881211670094411 0.675 0
+176 0.5764838143172589 0.7414018079610822 0
+177 0.5062835352933159 0.8285320053734242 0
+178 0.1919083949384496 0.5859235062107955 0
+179 0.1830209620514912 0.0804073476665709 0
+180 0.8169790379498036 0.08040734766613497 0
+181 0.7783440999795792 0.9231378488578053 0
+182 0.2943256441907511 0.5064828201676498 0
+183 0.8137131248178493 0.6187333045050307 0
+184 0.07198329968605219 0.6637232325325355 0
+185 0.7125553915779134 0.8255538488749383 0
+186 0.2420653812016603 0.3556915009318665 0
+187 0.7162853305035023 0.6910600345791089 0
+188 0.5614089287487484 0.4337624682611698 0
+189 0.9266309636474888 0.5325956485576818 0
+190 0.4326394759116713 0.4495534935171767 0
+191 0.1739139344176251 0.7225584697165623 0
+192 0.7585599693183998 0.3535781065164037 0
+193 0.5869534582661667 0.5470707302448944 0
+194 0.6320633904518222 0.9286056789882188 0
+195 0.8185727349725702 0.7756483731211691 0
+196 0.4997955311277253 0.3288124707775839 0
+197 0.8488538957455855 0.5275465730472561 0
+198 0.9306053690018283 0.8479472552320556 0
+199 0.06347874322572213 0.3217248706963154 0
+200 0.9365207075795818 0.3217252115110646 0
+201 0.9349886966491463 0.749011569482589 0
+202 0.2114171012111613 0.6594092944091384 0
+203 0.5016588869397121 0.06797218363804416 0
+204 0.1539812358975107 0.5195017016706647 0
+205 0.2967853160746573 0.6792234771733605 0
+206 0.4017562661221452 0.9355971982277441 0
+207 0.9268501882076796 0.9257227039290514 0
+208 0.4521720276710546 0.8167865974752775 0
+209 0.5510660968194043 0.6797135476535654 0
+210 0.5707236023072578 0.7996005435106294 0
+211 0.4381373385673064 0.7024990533039687 0
+212 0.06635917244426495 0.5181322413008239 0
+213 0.06611516951872876 0.9346190325398593 0
+214 0.5543716319693202 0.2507432879078373 0
+215 0.5036804263027204 0.1459746499361169 0
+216 0.3073163351524544 0.06513924970157009 0
+217 0.9353922225431665 0.6227473793370543 0
+218 0.7350090277987589 0.5687727272497948 0
+219 0.4006976059931963 0.5671395704540068 0
+220 0.5490886029692275 0.3733419394036243 0
+221 0.4871409253493302 0.2749191543540659 0
+222 0.2121291590007329 0.8638777218942592 0
+223 0.06327788539350787 0.06327788539342601 0
+224 0.9367221146084541 0.06327788539160781 0
+225 0.4408453192569133 0.3782357053890621 0
+226 0.672265162832817 0.6260525092862941 0
+227 0.3175238472732977 0.8460929421150545 0
+228 0.690844559035799 0.07324980066571454 0
+229 0.1485070800631576 0.2969615021959272 0
+230 0.8515889397175606 0.2965266306768867 0
+231 0.6862730010668697 0.2512313359680488 0
+232 0.3409693371981573 0.4446099852805522 0
+233 0.8409599685162412 0.8477843183722567 0
+234 0.3104699612440741 0.2510212862046502 0
+235 0.2617709125442271 0.5619213222720593 0
+236 0.6635652629321784 0.4256567033348759 0
+237 0.2851247162560687 0.9304034211976683 0
+238 0.3750924296143534 0.3235916737233546 0
+239 0.3328410440182527 0.6266250913165181 0
+240 0.8267567819907106 0.9333325436254075 0
+241 0.8731071350283545 0.7339817586611024 0
+242 0.5756785237919199 0.9450786574931633 0
+243 0.6216707706196916 0.3258829750713973 0
+244 0.7660667422145149 0.6709887092082645 0
+245 0.4320968816896103 0.2448530419385232 0
+246 0.4271581750771073 0.05202140637405563 0
+247 0.4844243257236938 0.4618766879791599 0
+248 0.9431495850766271 0.4775380703472147 0
+249 0.3135419049022754 0.372111320949921 0
+250 0.8087666142117923 0.563355196163648 0
+251 0.05898910560640176 0.4422673185029301 0
+252 0.1420165037221835 0.9432689791196578 0
+253 0.05852230625539075 0.6072681929335997 0
+254 0.03989836924342544 0.873505928323029 0
+255 0.161585131926244 0.2391857946243138 0
+256 0.8384148680534601 0.2391857946436191 0
+257 0.1605019015588762 0.1327874198246194 0
+258 0.8394980984457867 0.1327874198210596 0
+259 0.8546982616315805 0.6567907673047465 0
+260 0.5801094089017853 0.05519988950025151 0
+261 0.4688130687149384 0.946906534357117 0
+262 0.1458935810696156 0.4160586787075894 0
+263 0.1332611389930252 0.8357289212824203 0
+264 0.8599803548467403 0.4139820253212571 0
+265 0.05809049950584526 0.1396850039264114 0
+266 0.941909500492071 0.1396850039331882 0
+267 0.6875591900709364 0.3773195842445296 0
+268 0.4948819351903923 0.5440302349095326 0
+269 0.1260030920151941 0.03924500965653527 0
+270 0.8739969079849601 0.03924500965583505 0
+271 0.2837450967389632 0.1534509574488906 0
+272 0.715889371218871 0.1538115739137546 0
+273 0.05689393008678054 0.2392007003479759 0
+274 0.9431060699081403 0.2392007003574637 0
+275 0.2564698821974526 0.4745014878049475 0
+276 0.2272508895731511 0.7745050539575379 0
+277 0.7434751496091593 0.4747776748391724 0
+278 0.230095007397805 0.05527116088276794 0
+279 0.7713718482344473 0.05568096101351351 0
+280 0.3954895358822054 0.1311946474491396 0
+281 0.2975329254004855 0.7517533657618267 0
+282 0.2312484260837477 0.1548659677766349 0
+283 0.7701477087279389 0.153649072736263 0
+284 0.7246615011906069 0.9499865872288887 0
+285 0.05228162893836207 0.7689765800051673 0
+286 0.5234740235903357 0.9474406115558548 0
+287 0.1434928364647098 0.6008169115510076 0
+288 0.9369143619260438 0.7963424734034621 0
+289 0.6056935254245887 0.1265287360215735 0
+290 0.694799648716804 0.7780647687088221 0
+291 0.2195533029098034 0.2322181345580657 0
+292 0.7796613335576197 0.2323431189059976 0
+293 0.1284043423964295 0.7826218269009529 0
+294 0.6672484030990004 0.8649456527811299 0
+295 0.05257044232818209 0.7119689820929865 0
+296 0.5768097512884905 0.5005928579182558 0
+297 0.8177449505079945 0.8887008563606282 0
+298 0.2294882552756616 0.9549993967365497 0
+299 0.1229000152380553 0.3341367260228106 0
+300 0.8784021734240119 0.3314180266104568 0
+301 0.2612852635045506 0.3106588653245886 0
+302 0.1338391216935454 0.1775558540666768 0
+303 0.8661608783066284 0.1775558540597449 0
+304 0.7362799479076265 0.3095610459081599 0
+305 0.3298691987089148 0.2094081037923071 0
+306 0.8405619397184987 0.4754487586667506 0
+307 0.6749006885114874 0.1148165124130752 0
+308 0.6662408589113106 0.210027510078292 0
+309 0.3254241043492987 0.1136919987351632 0
+310 0.7398309408910323 0.86763111489038 0
+311 0.7750351953203204 0.7187944341426856 0
+312 0.1252582898663277 0.6904780465402182 0
+313 0.8756060054718642 0.9573551674591529 0
+314 0.1630281605397375 0.4727664473336795 0
+315 0.9542570184033007 0.4234226026146065 0
+316 0.7770489429100018 0.7966942205010731 0
+317 0.6129423969816983 0.584448125931658 0
+318 0.6229307686840463 0.5286705925933066 0
+319 0.0468667719792867 0.5585293346167686 0
+320 0.6710340972836625 0.9503793447962866 0
+321 0.16351444840473 0.8748630676245268 0
+322 0.5409128530145592 0.5598664917682544 0
+323 0.6180557007651564 0.4282727335387148 0
+324 0.3029177324594935 0.463906074489936 0
+325 0.3677573223375093 0.5326984032651564 0
+326 0.9596026746388204 0.5759470607763799 0
+327 0.386659609547144 0.4391546871797619 0
+328 0.6988567673691469 0.5877809739252305 0
+329 0.7699016756352965 0.398516480583833 0
+330 0.2935585545982071 0.5900449787731803 0
+331 0.1038443414274288 0.6239159461890235 0
+332 0.2303613598020513 0.4004222160535651 0
+333 0.3759913958355357 0.2008753378910623 0
+334 0.2092667149691265 0.1172373126918267 0
+335 0.790924558041748 0.1163157894181978 0
+336 0.7335262162965845 0.2613649833701958 0
+337 0.3780284240725885 0.6089205280021748 0
+338 0.7722316906008789 0.5919001002102559 0
+339 0.204204591054406 0.3193076314496321 0
+340 0.9637164607887825 0.6645399502951359 0
+341 0.6752906645349077 0.6731013751859178 0
+342 0.9584397092973465 0.8803586587207719 0
+343 0.4144145905079837 0.4867929429676162 0
+344 0.4071906266093612 0.8032841959364372 0
+345 0.4477615637324044 0.6520967215558213 0
+346 0.5919241503074778 0.6993248921312549 0
+347 0.550406892613776 0.8420073842713645 0
+348 0.4788496094589973 0.8595498742121117 0
+349 0.5213644246187449 0.6428123010426474 0
+350 0.6083983921963612 0.7713276330303331 0
+351 0.3952673185812053 0.7280033917497746 0
+352 0.6946295406314056 0.4623734776278599 0
+353 0.7964147896075295 0.3204578454089672 0
+354 0.6928239467305521 0.7280123650345414 0
+355 0.628187546919812 0.8864955313482064 0
+356 0.2637903916353117 0.2590618002233394 0
+357 0.4678410189532725 0.7698723693925276 0
+358 0.5190768058978629 0.7157603295484917 0
+359 0.5366629222045098 0.7585063407504753 0
+360 0.4616770325977443 0.7369342997859611 0
+361 0.3990795894476696 0.2732787185143976 0
+362 0.4744916117858474 0.1766977717550853 0
+363 0.2500846528867531 0.8379561173758063 0
+364 0.6021354104131599 0.2750256555632312 0
+365 0.6529666302828054 0.285562216418078 0
+366 0.6208584498738083 0.2003601180756036 0
+367 0.358868909639378 0.8713028563654119 0
+368 0.3220296400713081 0.965096374256715 0
+369 0.0450289126262645 0.3607604138323677 0
+370 0.9549710873749327 0.3607604138653258 0
+371 0.4494834520346518 0.5653640774679654 0
+372 0.1154017348982857 0.2184380366574663 0
+373 0.8845982650969715 0.2184380366670949 0
+374 0.1813660057272621 0.9623678594419632 0
+375 0.05053477354373449 0.2855334669670243 0
+376 0.9494651479974092 0.2855335156609285 0
+377 0.3452707255374174 0.2846309735189408 0
+378 0.03671708261646304 0.1758256035554163 0
+379 0.963282917377656 0.1758256035652568 0
+380 0.580711692084902 0.3534220668924306 0
+381 0.8938335439772058 0.6278666316511751 0
+382 0.2261306160633835 0.6245770892644745 0
+383 0.5317517780571531 0.18017010714328 0
+384 0.5382517311015336 0.04319803631442524 0
+385 0.18547155080015 0.5414054709994462 0
+386 0.3109721755075838 0.8812147116513576 0
+387 0.1108468018487458 0.5093351955856934 0
+388 0.6249999999990127 0.0368995104443293 0
+389 0.5270362527499111 0.2890114230556702 0
+390 0.6985827545497421 0.5420700460840071 0
+391 0.3174534614518505 0.8002951775851398 0
+392 0.3435007446427953 0.03901085329331016 0
+393 0.2806337743862168 0.3508239784579412 0
+394 0.03815739922519534 0.8255254899704985 0
+395 0.4574241457449562 0.4157646110794656 0
+396 0.3764347381353931 0.9648485910158631 0
+397 0.4674128648176603 0.04043219862157266 0
+398 0.651028594780229 0.3580375376812758 0
+399 0.2710147269579462 0.03850981891594153 0
+400 0.8160547675285058 0.8144748372609291 0
+401 0.1898828083564689 0.3924677648585352 0
+402 0.4840327136217024 0.5045899035376922 0
+403 0.3184271287711127 0.713553021289589 0
+404 0.2335349266971445 0.9073606538120773 0
+405 0.5321150042679361 0.4579975750317197 0
+406 0.9116850701940886 0.4304794630957428 0
+407 0.8850661428994187 0.5488096409595046 0
+408 0.2546522150490536 0.6705671661752909 0
+409 0.8120964339104122 0.3924554715241102 0
+410 0.1735160514492026 0.03954570748268497 0
+411 0.8266086811722826 0.03958055443240611 0
+412 0.2177033844180772 0.8148011089179827 0
+413 0.2140459534889999 0.6991120268244346 0
+414 0.03855836599309355 0.4069297364830907 0
+415 0.2108552637392916 0.4726719043460517 0
+416 0.03807980624245771 0.4828928843747823 0
+417 0.03185040451496983 0.6723545193255792 0
+418 0.1173408236332442 0.3787075016549468 0
+419 0.3865628287010288 0.7679666828284335 0
+420 0.4821789212659429 0.6354501090226924 0
+421 0.6138929668611892 0.7331403568182614 0
+422 0.5168771995047183 0.8645153440195986 0
+423 0.09962748576320227 0.4271054605938512 0
+424 0.3316376776348741 0.6708530868159079 0
+425 0.7254198094717239 0.03713715208398233 0
+426 0.8905187830110065 0.8367392614492802 0
+427 0.879055278955432 0.3703673244309111 0
+428 0.1807399821219172 0.7609410195444795 0
+429 0.6632877144010623 0.5011655953647051 0
+430 0.3294846476477217 0.5339859469527758 0
+431 0.8918576915949723 0.5065850912218047 0
+432 0.6308846443884883 0.6183388729246823 0
+433 0.5956168520498442 0.9077605865243321 0
+434 0.1063293600794699 0.1496886732496585 0
+435 0.8936706399207855 0.149688673250139 0
+436 0.7304293753751753 0.3789493549468333 0
+437 0.7923887859133348 0.4742481721121519 0
+438 0.1767916661041777 0.6761089307016153 0
+439 0.8696274884022378 0.697568862851349 0
+440 0.8405575370915923 0.5913446637488688 0
+441 0.5796386579344532 0.3978245567832364 0
+442 0.3853751124708639 0.04198269809718994 0
+443 0.3546475230466961 0.3631741336976982 0
+444 0.09746792111223453 0.9605693208032362 0
+445 0.7394039505379052 0.5251333957571672 0
+446 0.777375849609712 0.9660432626799348 0
+447 0.9665895474867356 0.5258076370020353 0
+448 0.7109263696910029 0.6542635890729319 0
+449 0.4847660365164375 0.1104679975442236 0
+450 0.1833780091978758 0.628775024955162 0
+451 0.4530172289211364 0.4844309878401058 0
+452 0.08894353958045299 0.8546627798510918 0
+453 0.5577651960046823 0.1475884667674519 0
+454 0.4083574267873839 0.08929280826436128 0
+455 0.1399007730583397 0.7440597624318674 0
+456 0.5886343710279837 0.2297807499736143 0
+457 0.7506789193613408 0.8279186190688776 0
+458 0.5364880554007446 0.4080136013586033 0
+459 0.4695925342774275 0.2320679337304671 0
+460 0.6250000000012135 0.9682588137628259 0
+461 0.9138928550856427 0.579567879061424 0
+462 0.4338905703499352 0.1544234824313795 0
+463 0.9540885663612039 0.7163323595564562 0
+464 0.04464615467308888 0.09074374766342681 0
+465 0.955353845327473 0.09074374766311784 0
+466 0.5032038283620132 0.7879332745942292 0
+467 0.4278277089335074 0.9631296728486094 0
+468 0.7395149931396 0.9072882462612831 0
+469 0.6653780528260819 0.0373375754746653 0
+470 0.7750688261756686 0.8858267004049365 0
+471 0.2567734003242691 0.5259099192617107 0
+472 0.9652142268716739 0.8316968223602572 0
+473 0.6736564052901084 0.8234244082007277 0
+474 0.1434297035809975 0.08842916792419363 0
+475 0.8565702964204095 0.08842916792412682 0
+476 0.4092442864957742 0.5263788199007129 0
+477 0.575065809093381 0.5854748642937199 0
+478 0.8067712456040665 0.6587319812397521 0
+479 0.7782492315738763 0.6322189028613673 0
+480 0.4069270085608257 0.3588586923128377 0
+481 0.1403336027720304 0.5570542547337253 0
+482 0.5157500776285544 0.3727351005722654 0
+483 0.7158038865114336 0.3462471362169249 0
+484 0.331786725709314 0.4952643566508288 0
+485 0.8281381749050026 0.735634863709433 0
+486 0.9625760449027807 0.9625760449027999 0
+487 0.8568761422096569 0.7684865009985571 0
+488 0.5125627948588626 0.2451228721333351 0
+489 0.5952765173179219 0.09265851165206641 0
+490 0.7380694708452544 0.7208034650606205 0
+491 0.08964908396208576 0.7761358886531573 0
+492 0.66019228818227 0.9062267567118509 0
+493 0.1899778659372953 0.2704304225676201 0
+494 0.8098606070430447 0.2704005515172035 0
+495 0.5374249729195943 0.3372459303803913 0
+496 0.4555417927657185 0.3291218831712886 0
+497 0.7355682586351789 0.7970729957346112 0
+498 0.04030148585790171 0.960287862177903 0
+499 0.4166890859861498 0.4081255394616878 0
+500 0.9070465576574881 0.7215574852985747 0
+501 0.1754812924386141 0.8387741038691691 0
+502 0.6908823940641997 0.1845666086614767 0
+503 0.788992837919323 0.7573467829768417 0
+504 0.7671889511670084 0.5525101013799046 0
+505 0.915717896104047 0.8806524239225173 0
+506 0.2336706046120287 0.4417925865415113 0
+507 0.02972586850406834 0.5243795319986964 0
+508 0.03985705480607428 0.03894801433812083 0
+509 0.96014294519457 0.03894801433742239 0
+510 0.3980969936418819 0.2321024612576316 0
+511 0.5565468856024529 0.5268919713531723 0
+512 0.9643449430747179 0.762201524860582 0
+513 0.03038191512837742 0.324287627704582 0
+514 0.9696180208258252 0.324287667452228 0
+515 0.2584880578049887 0.753672445543504 0
+516 0.03715874684545626 0.9169940794964864 0
+517 0.9111391282344532 0.3486390746416438 0
+518 0.08255670524725733 0.03631144004051522 0
+519 0.9174432947531812 0.03631144003965196 0
+520 0.08849659720937329 0.3493075797449283 0
+521 0.9235760340830108 0.9697786262318169 0
+522 0.9695731884815614 0.9237314815104457 0
+523 0.607700302689504 0.164103575997529 0
+524 0.4560628269566469 0.07888036939188454 0
+525 0.2956300706452174 0.6432331172024737 0
+526 0.2741904758143821 0.3865770219354968 0
+527 0.8180257967188529 0.5081292117368242 0
+528 0.7025590660962384 0.2885706333473892 0
+529 0.03310751022537751 0.6379295298279174 0
+530 0.7674178084270689 0.4394306280361689 0
+531 0.08062591178552957 0.4738287221558386 0
+532 0.442428875330346 0.8551913848218944 0
+533 0.5571292620462622 0.6451076684028514 0
+534 0.6048923315972 0.8071292620461874 0
+535 0.3956040537695126 0.6925172798097742 0
+536 0.8990152624279947 0.7581686941320099 0
+537 0.4325751928810223 0.9274359877404146 0
+538 0.2932457120113867 0.5487564584370643 0
+539 0.1863407933121187 0.5038611530536861 0
+540 0.1436599711899129 0.905909199317847 0
+541 0.5859901685374369 0.4597727176952695 0
+542 0.2822740097288993 0.7115578587885797 0
+543 0.8515397231234003 0.8059829491645871 0
+544 0.09695097059521789 0.3058827842120544 0
+545 0.9031720420490227 0.305427821538674 0
+546 0.3620232047225836 0.1161728771032239 0
+547 0.9728617755964831 0.6249999999990127 0
+548 0.3027470965174698 0.1865780603515497 0
+549 0.8520228034260743 0.5632002997237465 0
+550 0.2898659448168561 0.1005990892932667 0
+551 0.1961555405127517 0.1486247306413643 0
+552 0.8045336721190954 0.148937218707601 0
+553 0.226536034345496 0.5536959767022478 0
+554 0.7967854315953429 0.361336700378885 0
+555 0.7050699236367185 0.8607487986687201 0
+556 0.2384398876454908 0.1923262439244666 0
+557 0.7618006503195552 0.1941868937850041 0
+558 0.5437726732086153 0.2173982148551033 0
+559 0.6551705316447196 0.4598190425309201 0
+560 0.5852934889063478 0.6629473186059692 0
+561 0.5870614764120536 0.8352215230928042 0
+562 0.4160385149224072 0.8369997903855178 0
+563 0.4129551185506674 0.6635221362219785 0
+564 0.2999621732144231 0.2856054656111346 0
+565 0.3953069464606959 0.1686626688234709 0
+566 0.217005266087277 0.738397653814008 0
+567 0.7310494322793162 0.07228378747237291 0
+568 0.5538878033903838 0.7174861347338046 0
+569 0.5334561561561746 0.8033159097000246 0
+570 0.4822350848385833 0.7124753842711691 0
+571 0.2444528984859503 0.5908225707626107 0
+572 0.535810055103262 0.07923194399908105 0
+573 0.07698090601083714 0.9009940884824853 0
+574 0.3982403975286148 0.9028075058241124 0
+575 0.348504861040586 0.4033624727121226 0
+576 0.6488509972522262 0.3929331869422485 0
+577 0.03401078590641864 0.2103300332252601 0
+578 0.9659892140907053 0.2103300332298132 0
+579 0.2676622922407929 0.0721386839953696 0
+580 0.08721872242227677 0.7032151726457829 0
+581 0.4216161225992181 0.5961629782558722 0
+582 0.4686273629688805 0.1405831218583108 0
+583 0.4447977321822471 0.2776054877372967 0
+584 0.09538386437663225 0.07339038223312919 0
+585 0.904616135624607 0.07339038223085059 0
+586 0.2758485868524275 0.9728772371413709 0
+587 0.6412073410895367 0.1271296457272502 0
+588 0.2870148820303137 0.8371476176036567 0
+589 0.346294612797947 0.8283168062512931 0
+590 0.2036635418869092 0.3600469172533407 0
+591 0.8254835343950109 0.9729374286521076 0
+592 0.7099379584235467 0.1016466023858731 0
+593 0.4778021853127619 0.3731900268431282 0
+594 0.1075178812139328 0.9213839837169916 0
+595 0.8567499905195095 0.6175190826458412 0
+596 0.3092938734749872 0.4137749393501131 0
+597 0.1492892089859849 0.2053776331936021 0
+598 0.8507107910093066 0.2053776331975565 0
+599 0.5235735920373428 0.1117318466784937 0
+$EndNodes
+$Elements
+1136
+1 1 2 2 5 5 85
+2 1 2 2 5 85 86
+3 1 2 2 5 86 87
+4 1 2 2 5 87 88
+5 1 2 2 5 88 6
+6 1 2 2 6 6 89
+7 1 2 2 6 89 90
+8 1 2 2 6 90 91
+9 1 2 2 6 91 92
+10 1 2 2 6 92 7
+11 1 2 2 7 7 93
+12 1 2 2 7 93 94
+13 1 2 2 7 94 95
+14 1 2 2 7 95 96
+15 1 2 2 7 96 8
+16 1 2 2 8 8 97
+17 1 2 2 8 97 98
+18 1 2 2 8 98 99
+19 1 2 2 8 99 100
+20 1 2 2 8 100 5
+21 2 2 1 1 158 231 308
+22 2 2 1 1 177 208 466
+23 2 2 1 1 119 312 331
+24 2 2 1 1 158 308 502
+25 2 2 1 1 100 242 286
+26 2 2 1 1 273 133 375
+27 2 2 1 1 134 274 376
+28 2 2 1 1 229 133 255
+29 2 2 1 1 134 230 256
+30 2 2 1 1 221 196 496
+31 2 2 1 1 408 138 413
+32 2 2 1 1 375 133 544
+33 2 2 1 1 134 376 545
+34 2 2 1 1 135 265 464
+35 2 2 1 1 266 136 465
+36 2 2 1 1 135 464 584
+37 2 2 1 1 465 136 585
+38 2 2 1 1 152 217 340
+39 2 2 1 1 166 264 306
+40 2 2 1 1 130 236 267
+41 2 2 1 1 232 117 484
+42 2 2 1 1 203 18 384
+43 2 2 1 1 264 126 306
+44 2 2 1 1 312 184 331
+45 2 2 1 1 124 214 364
+46 2 2 1 1 247 106 405
+47 2 2 1 1 278 140 334
+48 2 2 1 1 141 279 335
+49 2 2 1 1 253 73 319
+50 2 2 1 1 216 14 392
+51 2 2 1 1 150 226 432
+52 2 2 1 1 247 190 395
+53 2 2 1 1 179 278 334
+54 2 2 1 1 279 180 335
+55 2 2 1 1 106 247 395
+56 2 2 1 1 128 288 426
+57 2 2 1 1 202 408 413
+58 2 2 1 1 18 203 397
+59 2 2 1 1 142 207 313
+60 2 2 1 1 221 496 583
+61 2 2 1 1 216 164 309
+62 2 2 1 1 285 70 295
+63 2 2 1 1 196 221 389
+64 2 2 1 1 233 142 297
+65 2 2 1 1 390 150 429
+66 2 2 1 1 324 232 484
+67 2 2 1 1 288 198 426
+68 2 2 1 1 133 229 544
+69 2 2 1 1 230 134 545
+70 2 2 1 1 208 357 466
+71 2 2 1 1 240 181 297
+72 2 2 1 1 234 157 305
+73 2 2 1 1 255 107 291
+74 2 2 1 1 108 256 292
+75 2 2 1 1 271 157 556
+76 2 2 1 1 158 272 557
+77 2 2 1 1 206 159 574
+78 2 2 1 1 150 318 429
+79 2 2 1 1 405 106 458
+80 2 2 1 1 385 178 481
+81 2 2 1 1 317 193 318
+82 2 2 1 1 282 271 556
+83 2 2 1 1 272 283 557
+84 2 2 1 1 242 55 286
+85 2 2 1 1 150 317 318
+86 2 2 1 1 227 367 386
+87 2 2 1 1 178 287 481
+88 2 2 1 1 181 240 446
+89 2 2 1 1 122 237 404
+90 2 2 1 1 221 459 488
+91 2 2 1 1 142 240 297
+92 2 2 1 1 35 315 370
+93 2 2 1 1 142 233 426
+94 2 2 1 1 187 244 490
+95 2 2 1 1 159 206 396
+96 2 2 1 1 121 244 448
+97 2 2 1 1 56 261 286
+98 2 2 1 1 207 142 505
+99 2 2 1 1 265 82 464
+100 2 2 1 1 30 266 465
+101 2 2 1 1 147 212 387
+102 2 2 1 1 95 226 341
+103 2 2 1 1 96 341 354
+104 2 2 1 1 314 204 387
+105 2 2 1 1 153 219 337
+106 2 2 1 1 211 175 570
+107 2 2 1 1 359 358 568
+108 2 2 1 1 19 20 260
+109 2 2 1 1 174 208 344
+110 2 2 1 1 175 211 345
+111 2 2 1 1 208 177 348
+112 2 2 1 1 209 175 349
+113 2 2 1 1 210 176 350
+114 2 2 1 1 211 174 351
+115 2 2 1 1 208 174 357
+116 2 2 1 1 175 209 358
+117 2 2 1 1 176 210 359
+118 2 2 1 1 174 211 360
+119 2 2 1 1 116 358 359
+120 2 2 1 1 229 111 299
+121 2 2 1 1 112 230 300
+122 2 2 1 1 127 314 387
+123 2 2 1 1 252 63 444
+124 2 2 1 1 296 405 541
+125 2 2 1 1 296 139 405
+126 2 2 1 1 82 83 464
+127 2 2 1 1 29 30 465
+128 2 2 1 1 212 147 319
+129 2 2 1 1 244 187 448
+130 2 2 1 1 459 105 488
+131 2 2 1 1 54 55 242
+132 2 2 1 1 12 278 410
+133 2 2 1 1 279 24 411
+134 2 2 1 1 118 233 297
+135 2 2 1 1 262 165 314
+136 2 2 1 1 16 17 246
+137 2 2 1 1 149 290 354
+138 2 2 1 1 301 186 339
+139 2 2 1 1 444 65 498
+140 2 2 1 1 192 304 353
+141 2 2 1 1 127 262 314
+142 2 2 1 1 145 284 320
+143 2 2 1 1 70 71 295
+144 2 2 1 1 191 312 438
+145 2 2 1 1 263 170 293
+146 2 2 1 1 69 285 394
+147 2 2 1 1 12 13 278
+148 2 2 1 1 23 24 279
+149 2 2 1 1 268 139 322
+150 2 2 1 1 143 301 339
+151 2 2 1 1 229 299 544
+152 2 2 1 1 300 230 545
+153 2 2 1 1 214 124 389
+154 2 2 1 1 304 144 353
+155 2 2 1 1 5 100 286
+156 2 2 1 1 264 163 406
+157 2 2 1 1 111 229 339
+158 2 2 1 1 230 112 353
+159 2 2 1 1 73 74 319
+160 2 2 1 1 169 228 307
+161 2 2 1 1 152 463 500
+162 2 2 1 1 218 121 328
+163 2 2 1 1 459 221 583
+164 2 2 1 1 231 154 308
+165 2 2 1 1 36 37 248
+166 2 2 1 1 156 234 305
+167 2 2 1 1 226 95 432
+168 2 2 1 1 164 216 392
+169 2 2 1 1 313 49 591
+170 2 2 1 1 367 227 589
+171 2 2 1 1 446 240 591
+172 2 2 1 1 222 122 404
+173 2 2 1 1 322 139 511
+174 2 2 1 1 244 311 490
+175 2 2 1 1 317 150 432
+176 2 2 1 1 56 57 261
+177 2 2 1 1 217 152 381
+178 2 2 1 1 109 293 428
+179 2 2 1 1 140 271 282
+180 2 2 1 1 272 141 283
+181 2 2 1 1 197 306 431
+182 2 2 1 1 312 191 455
+183 2 2 1 1 7 268 322
+184 2 2 1 1 121 218 338
+185 2 2 1 1 153 239 330
+186 2 2 1 1 126 248 431
+187 2 2 1 1 14 216 399
+188 2 2 1 1 69 70 285
+189 2 2 1 1 163 264 427
+190 2 2 1 1 219 153 325
+191 2 2 1 1 63 64 444
+192 2 2 1 1 248 189 431
+193 2 2 1 1 138 408 542
+194 2 2 1 1 251 76 414
+195 2 2 1 1 160 222 404
+196 2 2 1 1 152 340 463
+197 2 2 1 1 149 354 490
+198 2 2 1 1 159 237 386
+199 2 2 1 1 245 459 583
+200 2 2 1 1 237 122 386
+201 2 2 1 1 261 5 286
+202 2 2 1 1 222 160 321
+203 2 2 1 1 244 121 479
+204 2 2 1 1 253 147 331
+205 2 2 1 1 240 313 591
+206 2 2 1 1 428 293 455
+207 2 2 1 1 44 342 472
+208 2 2 1 1 306 126 431
+209 2 2 1 1 90 239 337
+210 2 2 1 1 135 257 434
+211 2 2 1 1 258 136 435
+212 2 2 1 1 240 142 313
+213 2 2 1 1 85 5 261
+214 2 2 1 1 67 68 254
+215 2 2 1 1 154 231 365
+216 2 2 1 1 238 125 480
+217 2 2 1 1 226 150 328
+218 2 2 1 1 122 222 363
+219 2 2 1 1 121 226 328
+220 2 2 1 1 109 263 293
+221 2 2 1 1 268 7 371
+222 2 2 1 1 101 238 443
+223 2 2 1 1 396 206 467
+224 2 2 1 1 237 298 404
+225 2 2 1 1 133 273 372
+226 2 2 1 1 274 134 373
+227 2 2 1 1 117 232 327
+228 2 2 1 1 273 167 372
+229 2 2 1 1 168 274 373
+230 2 2 1 1 142 426 505
+231 2 2 1 1 10 11 269
+232 2 2 1 1 25 26 270
+233 2 2 1 1 28 465 509
+234 2 2 1 1 464 84 508
+235 2 2 1 1 341 187 354
+236 2 2 1 1 244 110 311
+237 2 2 1 1 97 8 290
+238 2 2 1 1 228 169 469
+239 2 2 1 1 238 101 377
+240 2 2 1 1 231 158 336
+241 2 2 1 1 58 396 467
+242 2 2 1 1 226 121 448
+243 2 2 1 1 242 100 433
+244 2 2 1 1 480 125 496
+245 2 2 1 1 276 161 412
+246 2 2 1 1 64 65 444
+247 2 2 1 1 361 245 583
+248 2 2 1 1 257 107 302
+249 2 2 1 1 108 258 303
+250 2 2 1 1 105 362 383
+251 2 2 1 1 55 56 286
+252 2 2 1 1 145 320 492
+253 2 2 1 1 362 215 383
+254 2 2 1 1 14 15 392
+255 2 2 1 1 213 444 498
+256 2 2 1 1 315 163 370
+257 2 2 1 1 278 179 410
+258 2 2 1 1 180 279 411
+259 2 2 1 1 257 302 434
+260 2 2 1 1 303 258 435
+261 2 2 1 1 147 287 331
+262 2 2 1 1 193 296 318
+263 2 2 1 1 36 248 315
+264 2 2 1 1 61 62 298
+265 2 2 1 1 285 170 394
+266 2 2 1 1 107 255 597
+267 2 2 1 1 256 108 598
+268 2 2 1 1 272 114 307
+269 2 2 1 1 19 260 384
+270 2 2 1 1 227 386 588
+271 2 2 1 1 234 156 377
+272 2 2 1 1 249 101 443
+273 2 2 1 1 196 389 495
+274 2 2 1 1 233 118 400
+275 2 2 1 1 296 115 318
+276 2 2 1 1 132 263 321
+277 2 2 1 1 364 214 456
+278 2 2 1 1 125 361 583
+279 2 2 1 1 146 285 295
+280 2 2 1 1 48 49 313
+281 2 2 1 1 22 425 469
+282 2 2 1 1 157 234 356
+283 2 2 1 1 386 122 588
+284 2 2 1 1 425 228 469
+285 2 2 1 1 236 130 352
+286 2 2 1 1 330 239 525
+287 2 2 1 1 76 77 414
+288 2 2 1 1 245 361 510
+289 2 2 1 1 301 143 356
+290 2 2 1 1 44 45 342
+291 2 2 1 1 298 62 374
+292 2 2 1 1 242 194 460
+293 2 2 1 1 255 291 493
+294 2 2 1 1 292 256 494
+295 2 2 1 1 93 7 322
+296 2 2 1 1 311 485 503
+297 2 2 1 1 237 159 368
+298 2 2 1 1 250 103 527
+299 2 2 1 1 340 217 547
+300 2 2 1 1 239 153 337
+301 2 2 1 1 113 271 309
+302 2 2 1 1 250 183 338
+303 2 2 1 1 371 219 476
+304 2 2 1 1 125 238 361
+305 2 2 1 1 391 227 588
+306 2 2 1 1 325 117 476
+307 2 2 1 1 147 253 319
+308 2 2 1 1 248 126 406
+309 2 2 1 1 197 250 527
+310 2 2 1 1 464 223 584
+311 2 2 1 1 224 465 585
+312 2 2 1 1 117 343 476
+313 2 2 1 1 480 225 499
+314 2 2 1 1 124 243 380
+315 2 2 1 1 269 410 474
+316 2 2 1 1 411 270 475
+317 2 2 1 1 144 304 336
+318 2 2 1 1 275 129 324
+319 2 2 1 1 243 102 398
+320 2 2 1 1 51 52 284
+321 2 2 1 1 248 37 447
+322 2 2 1 1 410 179 474
+323 2 2 1 1 180 411 475
+324 2 2 1 1 288 128 536
+325 2 2 1 1 311 110 485
+326 2 2 1 1 54 242 460
+327 2 2 1 1 194 242 433
+328 2 2 1 1 225 480 496
+329 2 2 1 1 287 119 331
+330 2 2 1 1 110 244 478
+331 2 2 1 1 189 248 447
+332 2 2 1 1 323 236 559
+333 2 2 1 1 243 124 364
+334 2 2 1 1 239 90 424
+335 2 2 1 1 102 243 365
+336 2 2 1 1 7 92 371
+337 2 2 1 1 161 363 412
+338 2 2 1 1 126 264 406
+339 2 2 1 1 243 151 380
+340 2 2 1 1 165 262 401
+341 2 2 1 1 465 224 509
+342 2 2 1 1 223 464 508
+343 2 2 1 1 227 391 589
+344 2 2 1 1 28 29 465
+345 2 2 1 1 83 84 464
+346 2 2 1 1 236 323 576
+347 2 2 1 1 342 198 472
+348 2 2 1 1 254 394 452
+349 2 2 1 1 155 371 476
+350 2 2 1 1 151 243 398
+351 2 2 1 1 262 111 401
+352 2 2 1 1 305 157 548
+353 2 2 1 1 367 159 386
+354 2 2 1 1 16 246 442
+355 2 2 1 1 101 249 393
+356 2 2 1 1 184 253 331
+357 2 2 1 1 107 257 551
+358 2 2 1 1 258 108 552
+359 2 2 1 1 312 119 438
+360 2 2 1 1 86 85 574
+361 2 2 1 1 110 439 485
+362 2 2 1 1 115 323 559
+363 2 2 1 1 103 250 504
+364 2 2 1 1 254 68 394
+365 2 2 1 1 147 387 481
+366 2 2 1 1 123 330 525
+367 2 2 1 1 77 78 369
+368 2 2 1 1 34 35 370
+369 2 2 1 1 85 537 574
+370 2 2 1 1 255 133 372
+371 2 2 1 1 134 256 373
+372 2 2 1 1 277 172 352
+373 2 2 1 1 257 179 334
+374 2 2 1 1 180 258 335
+375 2 2 1 1 182 275 324
+376 2 2 1 1 264 166 409
+377 2 2 1 1 97 290 473
+378 2 2 1 1 130 277 352
+379 2 2 1 1 202 438 450
+380 2 2 1 1 301 356 564
+381 2 2 1 1 201 288 536
+382 2 2 1 1 247 139 402
+383 2 2 1 1 496 125 583
+384 2 2 1 1 112 264 409
+385 2 2 1 1 8 96 354
+386 2 2 1 1 87 367 589
+387 2 2 1 1 110 259 439
+388 2 2 1 1 274 32 376
+389 2 2 1 1 80 273 375
+390 2 2 1 1 394 170 452
+391 2 2 1 1 96 95 341
+392 2 2 1 1 323 151 576
+393 2 2 1 1 63 252 374
+394 2 2 1 1 269 11 410
+395 2 2 1 1 25 270 411
+396 2 2 1 1 284 52 320
+397 2 2 1 1 294 98 473
+398 2 2 1 1 246 17 397
+399 2 2 1 1 245 171 459
+400 2 2 1 1 148 480 499
+401 2 2 1 1 173 456 558
+402 2 2 1 1 171 245 510
+403 2 2 1 1 85 261 537
+404 2 2 1 1 439 241 485
+405 2 2 1 1 361 156 510
+406 2 2 1 1 281 6 391
+407 2 2 1 1 102 267 398
+408 2 2 1 1 67 254 516
+409 2 2 1 1 139 247 405
+410 2 2 1 1 293 170 491
+411 2 2 1 1 252 160 374
+412 2 2 1 1 190 247 451
+413 2 2 1 1 160 298 374
+414 2 2 1 1 251 162 423
+415 2 2 1 1 91 90 337
+416 2 2 1 1 146 312 455
+417 2 2 1 1 297 181 470
+418 2 2 1 1 152 259 381
+419 2 2 1 1 259 110 478
+420 2 2 1 1 183 250 440
+421 2 2 1 1 183 259 478
+422 2 2 1 1 157 271 548
+423 2 2 1 1 160 252 540
+424 2 2 1 1 76 251 416
+425 2 2 1 1 162 251 414
+426 2 2 1 1 77 369 414
+427 2 2 1 1 153 330 538
+428 2 2 1 1 169 260 388
+429 2 2 1 1 18 19 384
+430 2 2 1 1 253 184 529
+431 2 2 1 1 430 153 538
+432 2 2 1 1 140 282 334
+433 2 2 1 1 283 141 335
+434 2 2 1 1 250 197 549
+435 2 2 1 1 337 219 581
+436 2 2 1 1 264 112 427
+437 2 2 1 1 155 268 371
+438 2 2 1 1 405 188 541
+439 2 2 1 1 327 190 343
+440 2 2 1 1 260 20 388
+441 2 2 1 1 139 268 402
+442 2 2 1 1 117 327 343
+443 2 2 1 1 229 255 493
+444 2 2 1 1 256 230 494
+445 2 2 1 1 336 304 528
+446 2 2 1 1 164 442 454
+447 2 2 1 1 271 140 550
+448 2 2 1 1 217 381 461
+449 2 2 1 1 73 253 529
+450 2 2 1 1 259 183 595
+451 2 2 1 1 148 327 575
+452 2 2 1 1 167 265 434
+453 2 2 1 1 266 168 435
+454 2 2 1 1 265 135 434
+455 2 2 1 1 136 266 435
+456 2 2 1 1 35 36 315
+457 2 2 1 1 98 97 473
+458 2 2 1 1 10 269 518
+459 2 2 1 1 270 26 519
+460 2 2 1 1 442 246 454
+461 2 2 1 1 260 137 572
+462 2 2 1 1 284 181 446
+463 2 2 1 1 257 135 474
+464 2 2 1 1 136 258 475
+465 2 2 1 1 179 257 474
+466 2 2 1 1 258 180 475
+467 2 2 1 1 51 284 446
+468 2 2 1 1 265 167 378
+469 2 2 1 1 168 266 379
+470 2 2 1 1 6 281 403
+471 2 2 1 1 82 265 378
+472 2 2 1 1 266 30 379
+473 2 2 1 1 137 260 489
+474 2 2 1 1 326 39 547
+475 2 2 1 1 271 113 548
+476 2 2 1 1 117 325 484
+477 2 2 1 1 365 231 528
+478 2 2 1 1 167 273 577
+479 2 2 1 1 274 168 578
+480 2 2 1 1 263 132 452
+481 2 2 1 1 292 144 336
+482 2 2 1 1 275 182 471
+483 2 2 1 1 170 263 452
+484 2 2 1 1 137 289 453
+485 2 2 1 1 114 308 366
+486 2 2 1 1 158 292 336
+487 2 2 1 1 38 39 326
+488 2 2 1 1 324 129 596
+489 2 2 1 1 114 272 502
+490 2 2 1 1 104 275 471
+491 2 2 1 1 263 109 501
+492 2 2 1 1 267 236 576
+493 2 2 1 1 381 120 461
+494 2 2 1 1 141 272 592
+495 2 2 1 1 143 339 493
+496 2 2 1 1 100 99 433
+497 2 2 1 1 273 80 577
+498 2 2 1 1 32 274 578
+499 2 2 1 1 339 229 493
+500 2 2 1 1 260 169 489
+501 2 2 1 1 463 201 500
+502 2 2 1 1 268 155 402
+503 2 2 1 1 315 248 406
+504 2 2 1 1 98 294 355
+505 2 2 1 1 290 8 354
+506 2 2 1 1 272 158 502
+507 2 2 1 1 172 277 445
+508 2 2 1 1 259 152 439
+509 2 2 1 1 325 430 484
+510 2 2 1 1 111 262 418
+511 2 2 1 1 277 103 445
+512 2 2 1 1 262 127 423
+513 2 2 1 1 87 88 344
+514 2 2 1 1 91 92 345
+515 2 2 1 1 99 100 347
+516 2 2 1 1 95 96 346
+517 2 2 1 1 6 89 351
+518 2 2 1 1 5 85 348
+519 2 2 1 1 7 93 349
+520 2 2 1 1 8 97 350
+521 2 2 1 1 130 267 436
+522 2 2 1 1 352 172 429
+523 2 2 1 1 291 143 493
+524 2 2 1 1 144 292 494
+525 2 2 1 1 498 66 516
+526 2 2 1 1 107 282 556
+527 2 2 1 1 283 108 557
+528 2 2 1 1 267 102 483
+529 2 2 1 1 508 9 518
+530 2 2 1 1 27 509 519
+531 2 2 1 1 327 232 575
+532 2 2 1 1 438 119 450
+533 2 2 1 1 124 380 495
+534 2 2 1 1 199 375 544
+535 2 2 1 1 376 200 545
+536 2 2 1 1 353 144 494
+537 2 2 1 1 305 113 333
+538 2 2 1 1 217 326 547
+539 2 2 1 1 230 353 494
+540 2 2 1 1 159 367 574
+541 2 2 1 1 301 101 393
+542 2 2 1 1 143 291 356
+543 2 2 1 1 156 305 333
+544 2 2 1 1 291 157 356
+545 2 2 1 1 387 204 481
+546 2 2 1 1 361 238 377
+547 2 2 1 1 225 395 499
+548 2 2 1 1 52 53 320
+549 2 2 1 1 102 365 528
+550 2 2 1 1 287 178 450
+551 2 2 1 1 161 391 588
+552 2 2 1 1 456 214 558
+553 2 2 1 1 176 359 568
+554 2 2 1 1 360 211 570
+555 2 2 1 1 119 287 450
+556 2 2 1 1 232 324 596
+557 2 2 1 1 43 288 512
+558 2 2 1 1 129 275 506
+559 2 2 1 1 161 276 515
+560 2 2 1 1 275 104 415
+561 2 2 1 1 309 164 546
+562 2 2 1 1 154 364 456
+563 2 2 1 1 290 185 473
+564 2 2 1 1 140 278 579
+565 2 2 1 1 161 281 391
+566 2 2 1 1 408 205 542
+567 2 2 1 1 277 130 530
+568 2 2 1 1 278 13 399
+569 2 2 1 1 6 88 391
+570 2 2 1 1 352 429 559
+571 2 2 1 1 109 276 412
+572 2 2 1 1 218 390 445
+573 2 2 1 1 353 112 554
+574 2 2 1 1 284 145 468
+575 2 2 1 1 380 220 495
+576 2 2 1 1 276 109 428
+577 2 2 1 1 299 111 418
+578 2 2 1 1 291 107 556
+579 2 2 1 1 108 292 557
+580 2 2 1 1 15 16 442
+581 2 2 1 1 32 33 376
+582 2 2 1 1 79 80 375
+583 2 2 1 1 453 289 523
+584 2 2 1 1 304 192 483
+585 2 2 1 1 279 141 567
+586 2 2 1 1 261 57 467
+587 2 2 1 1 113 280 565
+588 2 2 1 1 280 113 546
+589 2 2 1 1 390 172 445
+590 2 2 1 1 181 284 468
+591 2 2 1 1 314 165 415
+592 2 2 1 1 23 279 425
+593 2 2 1 1 103 277 437
+594 2 2 1 1 281 161 515
+595 2 2 1 1 111 339 590
+596 2 2 1 1 320 194 492
+597 2 2 1 1 91 337 581
+598 2 2 1 1 309 271 550
+599 2 2 1 1 280 131 462
+600 2 2 1 1 172 390 429
+601 2 2 1 1 131 280 454
+602 2 2 1 1 282 107 551
+603 2 2 1 1 108 283 552
+604 2 2 1 1 66 67 516
+605 2 2 1 1 308 154 366
+606 2 2 1 1 167 302 372
+607 2 2 1 1 303 168 373
+608 2 2 1 1 99 98 355
+609 2 2 1 1 17 18 397
+610 2 2 1 1 454 246 524
+611 2 2 1 1 382 123 408
+612 2 2 1 1 9 10 518
+613 2 2 1 1 26 27 519
+614 2 2 1 1 113 305 548
+615 2 2 1 1 40 41 340
+616 2 2 1 1 194 320 460
+617 2 2 1 1 202 382 408
+618 2 2 1 1 288 43 472
+619 2 2 1 1 104 471 553
+620 2 2 1 1 368 159 396
+621 2 2 1 1 294 145 492
+622 2 2 1 1 526 249 596
+623 2 2 1 1 285 146 491
+624 2 2 1 1 383 215 453
+625 2 2 1 1 186 332 590
+626 2 2 1 1 219 371 581
+627 2 2 1 1 330 123 571
+628 2 2 1 1 170 285 491
+629 2 2 1 1 342 207 505
+630 2 2 1 1 272 307 592
+631 2 2 1 1 251 423 531
+632 2 2 1 1 235 330 571
+633 2 2 1 1 426 198 505
+634 2 2 1 1 453 215 599
+635 2 2 1 1 118 297 470
+636 2 2 1 1 131 454 524
+637 2 2 1 1 115 296 541
+638 2 2 1 1 89 6 403
+639 2 2 1 1 357 174 360
+640 2 2 1 1 116 357 360
+641 2 2 1 1 320 53 460
+642 2 2 1 1 166 306 437
+643 2 2 1 1 62 63 374
+644 2 2 1 1 3 47 486
+645 2 2 1 1 46 3 486
+646 2 2 1 1 392 15 442
+647 2 2 1 1 288 201 512
+648 2 2 1 1 289 169 587
+649 2 2 1 1 298 160 404
+650 2 2 1 1 308 114 502
+651 2 2 1 1 198 288 472
+652 2 2 1 1 298 237 586
+653 2 2 1 1 329 192 554
+654 2 2 1 1 186 301 393
+655 2 2 1 1 132 321 540
+656 2 2 1 1 59 60 368
+657 2 2 1 1 389 124 495
+658 2 2 1 1 169 289 489
+659 2 2 1 1 423 127 531
+660 2 2 1 1 290 149 497
+661 2 2 1 1 295 71 417
+662 2 2 1 1 287 147 481
+663 2 2 1 1 185 290 497
+664 2 2 1 1 61 298 586
+665 2 2 1 1 87 86 367
+666 2 2 1 1 332 401 590
+667 2 2 1 1 184 295 417
+668 2 2 1 1 58 59 396
+669 2 2 1 1 234 377 564
+670 2 2 1 1 289 137 489
+671 2 2 1 1 42 43 512
+672 2 2 1 1 41 42 463
+673 2 2 1 1 65 4 498
+674 2 2 1 1 4 66 498
+675 2 2 1 1 382 202 450
+676 2 2 1 1 146 295 580
+677 2 2 1 1 188 405 458
+678 2 2 1 1 344 88 419
+679 2 2 1 1 345 92 420
+680 2 2 1 1 347 100 422
+681 2 2 1 1 346 96 421
+682 2 2 1 1 6 351 419
+683 2 2 1 1 5 348 422
+684 2 2 1 1 7 349 420
+685 2 2 1 1 8 350 421
+686 2 2 1 1 21 22 469
+687 2 2 1 1 112 300 427
+688 2 2 1 1 317 94 477
+689 2 2 1 1 368 60 586
+690 2 2 1 1 164 454 546
+691 2 2 1 1 157 291 556
+692 2 2 1 1 292 158 557
+693 2 2 1 1 454 280 546
+694 2 2 1 1 307 114 587
+695 2 2 1 1 293 146 455
+696 2 2 1 1 237 368 586
+697 2 2 1 1 114 523 587
+698 2 2 1 1 485 195 503
+699 2 2 1 1 169 388 469
+700 2 2 1 1 114 366 523
+701 2 2 1 1 146 293 491
+702 2 2 1 1 281 515 542
+703 2 2 1 1 296 193 511
+704 2 2 1 1 139 296 511
+705 2 2 1 1 118 316 400
+706 2 2 1 1 243 364 365
+707 2 2 1 1 223 508 518
+708 2 2 1 1 509 224 519
+709 2 2 1 1 84 1 508
+710 2 2 1 1 1 9 508
+711 2 2 1 1 27 2 509
+712 2 2 1 1 2 28 509
+713 2 2 1 1 295 184 580
+714 2 2 1 1 364 154 365
+715 2 2 1 1 302 167 434
+716 2 2 1 1 168 303 435
+717 2 2 1 1 80 81 577
+718 2 2 1 1 31 32 578
+719 2 2 1 1 94 93 477
+720 2 2 1 1 218 328 390
+721 2 2 1 1 409 329 554
+722 2 2 1 1 204 385 481
+723 2 2 1 1 523 289 587
+724 2 2 1 1 294 473 555
+725 2 2 1 1 341 226 448
+726 2 2 1 1 323 188 441
+727 2 2 1 1 327 148 499
+728 2 2 1 1 471 235 553
+729 2 2 1 1 94 317 432
+730 2 2 1 1 30 31 379
+731 2 2 1 1 81 82 378
+732 2 2 1 1 302 107 597
+733 2 2 1 1 108 303 598
+734 2 2 1 1 356 234 564
+735 2 2 1 1 473 185 555
+736 2 2 1 1 515 138 542
+737 2 2 1 1 377 101 564
+738 2 2 1 1 57 58 467
+739 2 2 1 1 213 498 516
+740 2 2 1 1 151 323 441
+741 2 2 1 1 40 340 547
+742 2 2 1 1 391 88 589
+743 2 2 1 1 122 363 588
+744 2 2 1 1 102 304 483
+745 2 2 1 1 163 315 406
+746 2 2 1 1 145 294 555
+747 2 2 1 1 89 403 424
+748 2 2 1 1 113 309 546
+749 2 2 1 1 72 73 529
+750 2 2 1 1 316 195 400
+751 2 2 1 1 418 262 423
+752 2 2 1 1 191 428 455
+753 2 2 1 1 169 307 587
+754 2 2 1 1 168 379 578
+755 2 2 1 1 378 167 577
+756 2 2 1 1 101 301 564
+757 2 2 1 1 20 21 388
+758 2 2 1 1 393 249 526
+759 2 2 1 1 314 415 539
+760 2 2 1 1 403 205 424
+761 2 2 1 1 13 14 399
+762 2 2 1 1 162 418 423
+763 2 2 1 1 323 115 541
+764 2 2 1 1 326 189 447
+765 2 2 1 1 90 89 424
+766 2 2 1 1 304 102 528
+767 2 2 1 1 95 94 432
+768 2 2 1 1 75 76 416
+769 2 2 1 1 38 326 447
+770 2 2 1 1 68 69 394
+771 2 2 1 1 22 23 425
+772 2 2 1 1 463 42 512
+773 2 2 1 1 189 326 461
+774 2 2 1 1 321 160 540
+775 2 2 1 1 328 150 390
+776 2 2 1 1 338 183 479
+777 2 2 1 1 192 436 483
+778 2 2 1 1 316 149 503
+779 2 2 1 1 148 443 480
+780 2 2 1 1 329 130 436
+781 2 2 1 1 355 294 492
+782 2 2 1 1 443 238 480
+783 2 2 1 1 313 207 521
+784 2 2 1 1 192 329 436
+785 2 2 1 1 449 203 599
+786 2 2 1 1 48 313 521
+787 2 2 1 1 203 572 599
+788 2 2 1 1 195 316 503
+789 2 2 1 1 310 185 457
+790 2 2 1 1 354 187 490
+791 2 2 1 1 156 361 377
+792 2 2 1 1 231 336 528
+793 2 2 1 1 130 329 530
+794 2 2 1 1 149 316 497
+795 2 2 1 1 99 355 433
+796 2 2 1 1 326 217 461
+797 2 2 1 1 71 72 417
+798 2 2 1 1 173 453 523
+799 2 2 1 1 43 44 472
+800 2 2 1 1 118 310 457
+801 2 2 1 1 255 372 597
+802 2 2 1 1 373 256 598
+803 2 2 1 1 93 322 477
+804 2 2 1 1 129 526 596
+805 2 2 1 1 11 12 410
+806 2 2 1 1 24 25 411
+807 2 2 1 1 306 197 527
+808 2 2 1 1 318 115 429
+809 2 2 1 1 183 478 479
+810 2 2 1 1 332 165 401
+811 2 2 1 1 145 310 468
+812 2 2 1 1 384 260 572
+813 2 2 1 1 478 244 479
+814 2 2 1 1 322 193 477
+815 2 2 1 1 183 440 595
+816 2 2 1 1 310 118 470
+817 2 2 1 1 311 149 490
+818 2 2 1 1 307 228 592
+819 2 2 1 1 88 6 419
+820 2 2 1 1 92 7 420
+821 2 2 1 1 96 8 421
+822 2 2 1 1 100 5 422
+823 2 2 1 1 166 329 409
+824 2 2 1 1 332 129 506
+825 2 2 1 1 355 194 433
+826 2 2 1 1 149 311 503
+827 2 2 1 1 325 153 430
+828 2 2 1 1 185 310 555
+829 2 2 1 1 316 118 457
+830 2 2 1 1 371 92 581
+831 2 2 1 1 517 300 545
+832 2 2 1 1 299 520 544
+833 2 2 1 1 184 312 580
+834 2 2 1 1 310 145 555
+835 2 2 1 1 441 188 458
+836 2 2 1 1 383 173 558
+837 2 2 1 1 333 113 565
+838 2 2 1 1 299 418 520
+839 2 2 1 1 362 171 462
+840 2 2 1 1 200 517 545
+841 2 2 1 1 520 199 544
+842 2 2 1 1 173 366 456
+843 2 2 1 1 171 362 459
+844 2 2 1 1 220 441 458
+845 2 2 1 1 321 263 501
+846 2 2 1 1 216 309 550
+847 2 2 1 1 400 195 543
+848 2 2 1 1 437 306 527
+849 2 2 1 1 50 51 446
+850 2 2 1 1 193 317 477
+851 2 2 1 1 37 38 447
+852 2 2 1 1 174 344 419
+853 2 2 1 1 175 345 420
+854 2 2 1 1 177 347 422
+855 2 2 1 1 176 346 421
+856 2 2 1 1 351 174 419
+857 2 2 1 1 348 177 422
+858 2 2 1 1 349 175 420
+859 2 2 1 1 350 176 421
+860 2 2 1 1 312 146 580
+861 2 2 1 1 116 359 466
+862 2 2 1 1 347 177 569
+863 2 2 1 1 346 176 568
+864 2 2 1 1 209 346 568
+865 2 2 1 1 210 347 569
+866 2 2 1 1 333 171 510
+867 2 2 1 1 412 222 501
+868 2 2 1 1 204 314 539
+869 2 2 1 1 59 368 396
+870 2 2 1 1 239 424 525
+871 2 2 1 1 424 205 525
+872 2 2 1 1 109 412 501
+873 2 2 1 1 220 380 441
+874 2 2 1 1 53 54 460
+875 2 2 1 1 380 151 441
+876 2 2 1 1 537 206 574
+877 2 2 1 1 496 196 593
+878 2 2 1 1 184 417 529
+879 2 2 1 1 358 116 570
+880 2 2 1 1 116 360 570
+881 2 2 1 1 382 178 571
+882 2 2 1 1 155 343 451
+883 2 2 1 1 363 222 412
+884 2 2 1 1 195 487 543
+885 2 2 1 1 219 325 476
+886 2 2 1 1 182 324 484
+887 2 2 1 1 222 321 501
+888 2 2 1 1 50 446 591
+889 2 2 1 1 343 190 451
+890 2 2 1 1 212 319 507
+891 2 2 1 1 369 162 414
+892 2 2 1 1 319 74 507
+893 2 2 1 1 385 104 553
+894 2 2 1 1 85 86 532
+895 2 2 1 1 93 94 533
+896 2 2 1 1 97 98 534
+897 2 2 1 1 89 90 535
+898 2 2 1 1 338 218 504
+899 2 2 1 1 407 197 431
+900 2 2 1 1 407 189 461
+901 2 2 1 1 178 553 571
+902 2 2 1 1 190 327 499
+903 2 2 1 1 187 341 448
+904 2 2 1 1 250 338 504
+905 2 2 1 1 193 322 511
+906 2 2 1 1 74 75 507
+907 2 2 1 1 428 191 566
+908 2 2 1 1 343 155 476
+909 2 2 1 1 191 413 566
+910 2 2 1 1 78 79 513
+911 2 2 1 1 33 34 514
+912 2 2 1 1 188 323 541
+913 2 2 1 1 370 163 517
+914 2 2 1 1 162 369 520
+915 2 2 1 1 332 186 526
+916 2 2 1 1 127 387 531
+917 2 2 1 1 340 41 463
+918 2 2 1 1 47 48 521
+919 2 2 1 1 45 46 522
+920 2 2 1 1 200 370 517
+921 2 2 1 1 369 199 520
+922 2 2 1 1 121 338 479
+923 2 2 1 1 452 132 573
+924 2 2 1 1 387 212 531
+925 2 2 1 1 165 332 506
+926 2 2 1 1 164 392 442
+927 2 2 1 1 254 452 573
+928 2 2 1 1 39 40 547
+929 2 2 1 1 90 91 563
+930 2 2 1 1 86 87 562
+931 2 2 1 1 98 99 561
+932 2 2 1 1 94 95 560
+933 2 2 1 1 60 61 586
+934 2 2 1 1 362 105 459
+935 2 2 1 1 357 116 466
+936 2 2 1 1 49 50 591
+937 2 2 1 1 366 154 456
+938 2 2 1 1 334 282 551
+939 2 2 1 1 283 335 552
+940 2 2 1 1 375 199 513
+941 2 2 1 1 200 376 514
+942 2 2 1 1 92 91 581
+943 2 2 1 1 79 375 513
+944 2 2 1 1 376 33 514
+945 2 2 1 1 416 212 507
+946 2 2 1 1 88 87 589
+947 2 2 1 1 348 85 532
+948 2 2 1 1 349 93 533
+949 2 2 1 1 350 97 534
+950 2 2 1 1 351 89 535
+951 2 2 1 1 257 334 551
+952 2 2 1 1 335 258 552
+953 2 2 1 1 75 416 507
+954 2 2 1 1 426 233 543
+955 2 2 1 1 128 426 543
+956 2 2 1 1 103 437 527
+957 2 2 1 1 156 333 510
+958 2 2 1 1 344 208 562
+959 2 2 1 1 346 209 560
+960 2 2 1 1 347 210 561
+961 2 2 1 1 345 211 563
+962 2 2 1 1 129 332 526
+963 2 2 1 1 329 166 530
+964 2 2 1 1 202 413 438
+965 2 2 1 1 198 342 505
+966 2 2 1 1 137 453 599
+967 2 2 1 1 372 302 597
+968 2 2 1 1 303 373 598
+969 2 2 1 1 207 342 522
+970 2 2 1 1 218 445 504
+971 2 2 1 1 330 235 538
+972 2 2 1 1 342 45 522
+973 2 2 1 1 474 135 584
+974 2 2 1 1 136 475 585
+975 2 2 1 1 269 474 584
+976 2 2 1 1 475 270 585
+977 2 2 1 1 445 103 504
+978 2 2 1 1 246 397 524
+979 2 2 1 1 201 463 512
+980 2 2 1 1 189 407 431
+981 2 2 1 1 415 104 539
+982 2 2 1 1 178 382 450
+983 2 2 1 1 173 383 453
+984 2 2 1 1 194 355 492
+985 2 2 1 1 440 250 549
+986 2 2 1 1 247 402 451
+987 2 2 1 1 171 333 565
+988 2 2 1 1 413 191 438
+989 2 2 1 1 208 348 532
+990 2 2 1 1 209 349 533
+991 2 2 1 1 210 350 534
+992 2 2 1 1 211 351 535
+993 2 2 1 1 402 155 451
+994 2 2 1 1 339 186 590
+995 2 2 1 1 138 515 566
+996 2 2 1 1 91 345 563
+997 2 2 1 1 87 344 562
+998 2 2 1 1 99 347 561
+999 2 2 1 1 95 346 560
+1000 2 2 1 1 515 276 566
+1001 2 2 1 1 388 21 469
+1002 2 2 1 1 104 385 539
+1003 2 2 1 1 192 353 554
+1004 2 2 1 1 370 200 514
+1005 2 2 1 1 199 369 513
+1006 2 2 1 1 106 395 593
+1007 2 2 1 1 236 352 559
+1008 2 2 1 1 418 162 520
+1009 2 2 1 1 34 370 514
+1010 2 2 1 1 369 78 513
+1011 2 2 1 1 395 225 593
+1012 2 2 1 1 385 204 539
+1013 2 2 1 1 120 407 461
+1014 2 2 1 1 389 221 488
+1015 2 2 1 1 366 173 523
+1016 2 2 1 1 196 482 593
+1017 2 2 1 1 214 389 488
+1018 2 2 1 1 205 408 525
+1019 2 2 1 1 403 281 542
+1020 2 2 1 1 458 106 482
+1021 2 2 1 1 186 393 526
+1022 2 2 1 1 175 358 570
+1023 2 2 1 1 358 209 568
+1024 2 2 1 1 359 210 569
+1025 2 2 1 1 316 457 497
+1026 2 2 1 1 220 458 482
+1027 2 2 1 1 408 123 525
+1028 2 2 1 1 362 462 582
+1029 2 2 1 1 220 482 495
+1030 2 2 1 1 482 196 495
+1031 2 2 1 1 225 496 593
+1032 2 2 1 1 182 430 538
+1033 2 2 1 1 235 471 538
+1034 2 2 1 1 471 182 538
+1035 2 2 1 1 363 161 588
+1036 2 2 1 1 215 362 582
+1037 2 2 1 1 487 241 536
+1038 2 2 1 1 468 310 470
+1039 2 2 1 1 367 86 574
+1040 2 2 1 1 128 487 536
+1041 2 2 1 1 395 190 499
+1042 2 2 1 1 444 213 594
+1043 2 2 1 1 398 267 576
+1044 2 2 1 1 466 359 569
+1045 2 2 1 1 181 468 470
+1046 2 2 1 1 212 416 531
+1047 2 2 1 1 241 500 536
+1048 2 2 1 1 416 251 531
+1049 2 2 1 1 401 111 590
+1050 2 2 1 1 500 201 536
+1051 2 2 1 1 430 182 484
+1052 2 2 1 1 415 165 506
+1053 2 2 1 1 436 267 483
+1054 2 2 1 1 399 216 579
+1055 2 2 1 1 429 115 559
+1056 2 2 1 1 397 203 524
+1057 2 2 1 1 462 171 565
+1058 2 2 1 1 112 409 554
+1059 2 2 1 1 381 259 595
+1060 2 2 1 1 178 385 553
+1061 2 2 1 1 280 462 565
+1062 2 2 1 1 437 277 530
+1063 2 2 1 1 120 381 595
+1064 2 2 1 1 439 152 500
+1065 2 2 1 1 105 383 558
+1066 2 2 1 1 275 415 506
+1067 2 2 1 1 249 575 596
+1068 2 2 1 1 195 485 487
+1069 2 2 1 1 123 382 571
+1070 2 2 1 1 241 439 500
+1071 2 2 1 1 166 437 530
+1072 2 2 1 1 81 378 577
+1073 2 2 1 1 379 31 578
+1074 2 2 1 1 203 384 572
+1075 2 2 1 1 485 241 487
+1076 2 2 1 1 213 573 594
+1077 2 2 1 1 278 399 579
+1078 2 2 1 1 407 120 549
+1079 2 2 1 1 249 443 575
+1080 2 2 1 1 457 185 497
+1081 2 2 1 1 425 279 567
+1082 2 2 1 1 205 403 542
+1083 2 2 1 1 233 400 543
+1084 2 2 1 1 427 300 517
+1085 2 2 1 1 163 427 517
+1086 2 2 1 1 197 407 549
+1087 2 2 1 1 151 398 576
+1088 2 2 1 1 261 467 537
+1089 2 2 1 1 467 206 537
+1090 2 2 1 1 449 131 524
+1091 2 2 1 1 417 72 529
+1092 2 2 1 1 203 449 524
+1093 2 2 1 1 553 235 571
+1094 2 2 1 1 449 215 582
+1095 2 2 1 1 131 449 582
+1096 2 2 1 1 120 440 549
+1097 2 2 1 1 276 428 566
+1098 2 2 1 1 413 138 566
+1099 2 2 1 1 228 425 567
+1100 2 2 1 1 575 232 596
+1101 2 2 1 1 209 533 560
+1102 2 2 1 1 210 534 561
+1103 2 2 1 1 208 532 562
+1104 2 2 1 1 211 535 563
+1105 2 2 1 1 207 486 521
+1106 2 2 1 1 486 207 522
+1107 2 2 1 1 516 254 573
+1108 2 2 1 1 486 47 521
+1109 2 2 1 1 46 486 522
+1110 2 2 1 1 213 516 573
+1111 2 2 1 1 443 148 575
+1112 2 2 1 1 132 540 594
+1113 2 2 1 1 177 466 569
+1114 2 2 1 1 540 252 594
+1115 2 2 1 1 462 131 582
+1116 2 2 1 1 487 128 543
+1117 2 2 1 1 252 444 594
+1118 2 2 1 1 216 550 579
+1119 2 2 1 1 440 120 595
+1120 2 2 1 1 214 488 558
+1121 2 2 1 1 215 449 599
+1122 2 2 1 1 488 105 558
+1123 2 2 1 1 518 269 584
+1124 2 2 1 1 270 519 585
+1125 2 2 1 1 223 518 584
+1126 2 2 1 1 519 224 585
+1127 2 2 1 1 572 137 599
+1128 2 2 1 1 482 106 593
+1129 2 2 1 1 228 567 592
+1130 2 2 1 1 532 86 562
+1131 2 2 1 1 533 94 560
+1132 2 2 1 1 534 98 561
+1133 2 2 1 1 535 90 563
+1134 2 2 1 1 573 132 594
+1135 2 2 1 1 550 140 579
+1136 2 2 1 1 567 141 592
+$EndElements
diff --git a/src/mmdg.cc b/src/mmdg.cc
index 32140f3..4d567e2 100644
--- a/src/mmdg.cc
+++ b/src/mmdg.cc
@@ -33,7 +33,7 @@ int main(int argc, char** argv)
static constexpr int dim = GRIDDIM;
using Grid = Dune::MovingMesh<dim>;
- using GridFactory = Dune::DGFGridFactory<Grid>;
+ using GridFactory = Dune::GmshGridFactory<Grid, /*implicit=*/false>;
using GridView = typename Grid::LeafGridView;
using Mapper = typename Dune::MultipleCodimMultipleGeomTypeMapper<GridView>;
using IGrid = typename Grid::InterfaceGrid;
@@ -53,7 +53,7 @@ int main(int argc, char** argv)
const double xi = std::stod(pt["xi"]); //coupling parameter
//create a grid from .dgf file
- GridFactory gridFactory( "grids/cube" + std::to_string(dim) + "d.dgf" );
+ GridFactory gridFactory( "grids/sphere" + std::to_string(dim) + "d.msh" );
const Grid& grid = *gridFactory.grid();
const GridView& gridView = grid.leafGridView();
const IGrid& iGrid = grid.interfaceGrid();
--
GitLab