diff --git a/NonLinearSolver/field/nlsField.cpp b/NonLinearSolver/field/nlsField.cpp
index cafbf9e3a1e6083bed34bb0753c0a1c79ec0535e..d14f1de4d154e5339db9af61c995a45a2edf0103 100644
--- a/NonLinearSolver/field/nlsField.cpp
+++ b/NonLinearSolver/field/nlsField.cpp
@@ -90,6 +90,13 @@ void nlsField::createFile(){
}
}
+void nlsField::resetFile(){
+ if(_vBuildView->size() > 0)
+ {
+ this->createFile();
+ }
+};
+
nlsField::nlsField(const std::string &fnn, const uint32_t fms, const int ncomp,
const std::vector<dataBuildView> &dbview_) : numfile(0), fmaxsize(fms),
totelem(0), numcomp(ncomp), view(false),
diff --git a/NonLinearSolver/field/nlsField.h b/NonLinearSolver/field/nlsField.h
index 29cdaadc43d173d7d800d8cc38606ec3dd14e7d9..87d1ba3d47aa7d0fde81aad3d9e4b6a1003ed907 100644
--- a/NonLinearSolver/field/nlsField.h
+++ b/NonLinearSolver/field/nlsField.h
@@ -93,6 +93,7 @@ class nlsField{
virtual void onelabView(const std::string &meshFileName, const std::string &fname,
const std::vector<partDomain*> &vdom,const std::string &valuename,
const int comp,const double time, const int nstep) const;
+ void resetFile();
#endif // SWIG
};
#endif // nlsField_H_
diff --git a/NonLinearSolver/field/unknownField.cpp b/NonLinearSolver/field/unknownField.cpp
index c6083e9354fc94ed85a4bbaf7df7805b81e3b8d2..c9b2f58affb20d660346c5627c3d47c46cca3c3c 100644
--- a/NonLinearSolver/field/unknownField.cpp
+++ b/NonLinearSolver/field/unknownField.cpp
@@ -289,3 +289,11 @@ unknownField::restart()
nodesField::restart();
return;
}
+
+void
+unknownField::openArchivingFiles(){
+ for(int i=0;i<varch.size();i++)
+ {
+ varch[i].openFile();
+ }
+};
diff --git a/NonLinearSolver/field/unknownField.h b/NonLinearSolver/field/unknownField.h
index 7cb9a9e3082d3b3edefa64e095a86769a24fe9bf..242fbe4057cb5bde896339310ab1f1364e1cfda0 100644
--- a/NonLinearSolver/field/unknownField.h
+++ b/NonLinearSolver/field/unknownField.h
@@ -42,13 +42,19 @@ class unknownField : public nodesField{
archiveNode(Dof R, const int pnum, int n, int comp,nonLinearBoundaryCondition::whichCondition wv=nonLinearBoundaryCondition::position,
int nstep_=1,bool isc=false) : physnum(pnum),nodenum(n), D(R),_comp(comp), wc(wv), nstep(nstep_), iscontact(isc),fp(NULL)
{
- // open File
+ openFile();
+ }
+ archiveNode(const archiveNode &source) : physnum(source.physnum), nodenum(source.nodenum), D(source.D), _comp(source._comp),
+ nstep(source.nstep), wc(source.wc), fp(source.fp), iscontact(source.iscontact){}
+ ~archiveNode(){};
+ void openFile(){
+ // open File
std::ostringstream oss;
oss << nodenum;
std::string s = oss.str();
// component of displacement
oss.str("");
- oss << comp;
+ oss << _comp;
std::string s2 = oss.str();
#if defined(HAVE_MPI)
if(Msg::GetCommSize() != 1){
@@ -70,10 +76,7 @@ class unknownField : public nodesField{
break;
}
fp = fopen(fname.c_str(),"a");
- }
- archiveNode(const archiveNode &source) : physnum(source.physnum), nodenum(source.nodenum), D(source.D), _comp(source._comp),
- nstep(source.nstep), wc(source.wc), fp(source.fp), iscontact(source.iscontact){}
- ~archiveNode(){};
+ };
void closeFile(){fclose(fp);}
};
protected:
@@ -99,5 +102,6 @@ class unknownField : public nodesField{
virtual void archive(const double time,const int step=1);
virtual void valuesForOnelab(fullVector<double> &nodalValues) const;
virtual void restart();
+ void openArchivingFiles();
};
#endif // _UNKNOWNFIELD_H_
diff --git a/NonLinearSolver/internalPoints/ipField.cpp b/NonLinearSolver/internalPoints/ipField.cpp
index 6758b6ac168dc3bd95f2cc71f37cb586b653da90..d2af88d95eeffa587733be50f0312912b12a6040 100644
--- a/NonLinearSolver/internalPoints/ipField.cpp
+++ b/NonLinearSolver/internalPoints/ipField.cpp
@@ -1068,16 +1068,6 @@ double IPField::getIPcompAndVolume(const partDomain *ef,MElement *ele, const IPS
}
}
-void IPField::resetUnknownField()
-{
- // reset archiving file if needed
- if(_vBuildView->size() > 0)
- {
-// FILE *fp = fopen(this->fileName.c_str(),"w"); WHY ??
- this->createFile();
- }
-}
-
void IPField::setFactorOnArchivingFiles(const int fact)
{
for(int i=0;i<viparch.size();i++)
diff --git a/NonLinearSolver/internalPoints/ipField.h b/NonLinearSolver/internalPoints/ipField.h
index 495705646a9e2bad3a4b234948c7d9454643f225..20f610d56e1f041173df437f8f4935fae3441617 100644
--- a/NonLinearSolver/internalPoints/ipField.h
+++ b/NonLinearSolver/internalPoints/ipField.h
@@ -216,8 +216,6 @@ class IPField : public elementsField {
// initial broken
void initialBroken(GModel* pModel, std::vector<int> &vnumphys);
void initialBroken(MElement *iele, materialLaw *mlaw);
- // used in case of switch between scheme
- void resetUnknownField();
void setFactorOnArchivingFiles(const int fact);
// close all files
diff --git a/NonLinearSolver/nlsolver/nonLinearMechSolver.cpp b/NonLinearSolver/nlsolver/nonLinearMechSolver.cpp
index 6616d7ceadb5c5b627ab41f7e3b0fbdaa7414e89..2b402a1ff3c51280d89db67a6ccaf852b8723e24 100644
--- a/NonLinearSolver/nlsolver/nonLinearMechSolver.cpp
+++ b/NonLinearSolver/nlsolver/nonLinearMechSolver.cpp
@@ -890,19 +890,20 @@ void nonLinearMechSolver::init2()
}
#endif // HAVE_MPI
- /* fields creation */
- if(_ufield) delete _ufield;
- _ufield = new unknownField(this,3,anoded,unknownView);
- // this->restart(ufield); // CHANGE THIS
-
- if (OneUnknownView.size() > 0){
- if (_oneUnknownField) delete _oneUnknownField;
- _oneUnknownField = new OneUnknownField(this,OneUnknownView);
- }
bool stiffcomputation=true;
if(whatScheme==Explicit) stiffcomputation=false;
if(!_previousInit){
+ /* fields creation */
+ if(_ufield) delete _ufield;
+ _ufield = new unknownField(this,3,anoded,unknownView);
+ // this->restart(ufield); // CHANGE THIS
+
+ if (OneUnknownView.size() > 0){
+ if (_oneUnknownField) delete _oneUnknownField;
+ _oneUnknownField = new OneUnknownField(this,OneUnknownView);
+ }
+
Msg::Info("Creation of IPVariable");
_ipf = new IPField(this,vaip,ipView); // Field for GaussPoint
Msg::Info("IPVariable are created");
@@ -924,18 +925,44 @@ void nonLinearMechSolver::init2()
}
}
-
+
+ if(_energField) delete _energField;
+ _energField = new energeticField(this, energyView,_energyComputation,_fractureEnergyComputation,getScheme());
}
else{ // as nextStep is performed at the end of previous scheme the last step is in previous (and not in current) !!
- _ipf->resetUnknownField();
+ // all files are closed from previous scheme
+ // archived files need to be reopenned
+ _ipf->resetFile();
_ipf->copy(IPStateBase::previous,IPStateBase::initial);
_ipf->copy(IPStateBase::previous,IPStateBase::current); // same value to avois initialization problem ??
_ipf->openArchivingFiles();
+
+ if (_unresetUnknowns){
+ _ufield->resetFile();
+ _ufield->openArchivingFiles();
+ if (_energField){
+ _energField->resetFile();
+ }
+ if (_oneUnknownField){
+ _oneUnknownField->resetFile();
+ }
+ }
+ else{
+ // recreate fields as unknowns are changed
+ if(_ufield) delete _ufield;
+ _ufield = new unknownField(this,3,anoded,unknownView);
+ // this->restart(ufield); // CHANGE THIS
+
+ if (OneUnknownView.size() > 0){
+ if (_oneUnknownField) delete _oneUnknownField;
+ _oneUnknownField = new OneUnknownField(this,OneUnknownView);
+ }
+
+ if(_energField) delete _energField;
+ _energField = new energeticField(this, energyView,_energyComputation,_fractureEnergyComputation,getScheme());
+ }
}
- if(_energField) delete _energField;
- _energField = new energeticField(this, energyView,_energyComputation,_fractureEnergyComputation,getScheme());
-
/* terms initialization */
this->initTerms(_ufield,_ipf);
_previousInit = true;
@@ -974,6 +1001,7 @@ void nonLinearMechSolver::endOfScheme(unknownField *ufield,energeticField *efiel
if (_oneUnknownField!= NULL){
_oneUnknownField->forceView();
_oneUnknownField->archive(curtime,step);
+ _oneUnknownField->closeFile();
}
_ipf->forceView();
diff --git a/dG3D/benchmarks/CMakeLists.txt b/dG3D/benchmarks/CMakeLists.txt
index 480d1b239d7261153d2a38a1baf249aa50fe651c..97bc8a904b89cb48648595499a6ca5c5dc919d20 100644
--- a/dG3D/benchmarks/CMakeLists.txt
+++ b/dG3D/benchmarks/CMakeLists.txt
@@ -112,3 +112,4 @@ add_subdirectory(multiscaleEnhanceStrain)
add_subdirectory(nonLinearMixedBC_2D_DirectionFollowing)
add_subdirectory(interpolationPBC_2DShifted)
add_subdirectory(DG_PRI6_resetBC)
+add_subdirectory(DG_PRI6_DiriToNeum)
diff --git a/dG3D/benchmarks/DG_PRI6_DiriToNeum/CMakeLists.txt b/dG3D/benchmarks/DG_PRI6_DiriToNeum/CMakeLists.txt
new file mode 100644
index 0000000000000000000000000000000000000000..b2f3e302ec34f718b47d65153ccec9ba3fcf16bc
--- /dev/null
+++ b/dG3D/benchmarks/DG_PRI6_DiriToNeum/CMakeLists.txt
@@ -0,0 +1,12 @@
+# test file
+
+set(PYFILE twoHole.py)
+
+set(FILES2DELETE
+ *.csv
+ disp*
+ stress*
+ previousScheme*
+)
+
+add_cm3python_test(${PYFILE} "${FILES2DELETE}")
diff --git a/dG3D/benchmarks/DG_PRI6_DiriToNeum/twoHole.geo b/dG3D/benchmarks/DG_PRI6_DiriToNeum/twoHole.geo
new file mode 100644
index 0000000000000000000000000000000000000000..9b42aefe356d6eff249f6f74ccae56f8f5f0ae66
--- /dev/null
+++ b/dG3D/benchmarks/DG_PRI6_DiriToNeum/twoHole.geo
@@ -0,0 +1,72 @@
+
+mm=1.0e-3;
+n=1;
+L=10*mm;
+sl1=0.4*L/n;
+Lx = 2*L;
+Ly = 4*L;
+Point(1)={0,0,0,sl1};
+Point(2)={Lx,0,0,sl1};
+Point(3)={Lx,Ly,0,sl1};
+Point(4)={0,Ly,0,sl1};
+Line(1)={1,2};
+Line(2)={2,3};
+Line(3)={3,4};
+Line(4)={4,1};
+x0=0.5*Lx/n; y0=0.5*Ly/n; r=0.2*L/n;
+
+x=0.3*Lx;
+y=0.4*Ly;
+r = 0.2*L;
+sl2=0.7*sl1;
+p1=newp; Point(p1)={x-r,y,0,sl2};
+p2=newp; Point(p2)={x,y+r,0,sl2};
+p3=newp; Point(p3)={x+r,y,0,sl2};
+p4=newp; Point(p4)={x,y-r,0,sl2};
+pc=newp; Point(pc)={x,y,0,sl2};
+c1 = newreg; Circle(c1) = {p1,pc,p2};
+c2 = newreg; Circle(c2) = {p2,pc,p3};
+c3 = newreg; Circle(c3) = {p3,pc,p4};
+c4 = newreg; Circle(c4) = {p4,pc,p1};
+l[1]=newreg; Line Loop(l[1]) = {c1,c2,c3,c4};
+
+
+x=0.7*Lx;
+y=0.6*Ly;
+r = 0.2*L;
+sl2=0.4*sl1;
+p1=newp; Point(p1)={x-r,y,0,sl2};
+p2=newp; Point(p2)={x,y+r,0,sl2};
+p3=newp; Point(p3)={x+r,y,0,sl2};
+p4=newp; Point(p4)={x,y-r,0,sl2};
+pc=newp; Point(pc)={x,y,0,sl2};
+c1 = newreg; Circle(c1) = {p1,pc,p2};
+c2 = newreg; Circle(c2) = {p2,pc,p3};
+c3 = newreg; Circle(c3) = {p3,pc,p4};
+c4 = newreg; Circle(c4) = {p4,pc,p1};
+l[2]=newreg; Line Loop(l[2]) = {c1,c2,c3,c4};
+
+
+l[0]=newreg;
+Line Loop(l[0])={1,2,3,4};
+Plane Surface(11)={l[]};
+//Recombine Surface{11};
+
+Physical Surface(16) = {11};
+Physical Line(17) = {1};
+Physical Line(18) = {3};
+Physical Point(19) = {4};
+Physical Point(20) = {1};
+//+
+//Recombine Surface {11};
+//+
+Extrude {0, 0, 1e-3} {
+ Surface{11}; Layers{2}; Recombine;
+}
+
+//+
+Physical Volume(83) = {1};
+//+
+Physical Surface(84) = {37};
+//+
+Physical Surface(85) = {45};
diff --git a/dG3D/benchmarks/DG_PRI6_DiriToNeum/twoHole.msh b/dG3D/benchmarks/DG_PRI6_DiriToNeum/twoHole.msh
new file mode 100644
index 0000000000000000000000000000000000000000..88eb048998fa937733cb23c3cce59073082c8ddc
--- /dev/null
+++ b/dG3D/benchmarks/DG_PRI6_DiriToNeum/twoHole.msh
@@ -0,0 +1,1313 @@
+$MeshFormat
+2.2 0 8
+$EndMeshFormat
+$Nodes
+468
+1 0 0 0
+2 0.02 0 0
+3 0.02 0.04 0
+4 0 0.04 0
+5 0.004 0.016 0
+6 0.006 0.018 0
+7 0.008 0.016 0
+8 0.006 0.014 0
+9 0.012 0.024 0
+10 0.014 0.026 0
+11 0.016 0.024 0
+12 0.014 0.022 0
+13 0 0 0.001
+14 0.02 0 0.001
+15 0.02 0.04 0.001
+16 0 0.04 0.001
+17 0.004 0.016 0.001
+18 0.006 0.014 0.001
+19 0.008 0.016 0.001
+20 0.006 0.018 0.001
+21 0.012 0.024 0.001
+22 0.014 0.022 0.001
+23 0.016 0.024 0.001
+24 0.014 0.026 0.001
+25 0.003999999999991045 0 0
+26 0.007999999999979423 0 0
+27 0.01199999999997901 0 0
+28 0.01599999999998917 0 0
+29 0.02 0.003999999999992779 0
+30 0.02 0.007999999999982091 0
+31 0.02 0.0119999999999705 0
+32 0.02 0.01599999999995885 0
+33 0.02 0.01999999999994786 0
+34 0.02 0.02399999999995802 0
+35 0.02 0.02799999999996886 0
+36 0.02 0.03199999999997834 0
+37 0.02 0.03599999999998998 0
+38 0.01600000000001043 0.04 0
+39 0.01200000000002099 0.04 0
+40 0.008000000000021266 0.04 0
+41 0.004000000000010967 0.04 0
+42 0 0.0360000000000088 0
+43 0 0.03200000000002086 0
+44 0 0.02800000000003047 0
+45 0 0.02400000000004199 0
+46 0 0.02000000000005283 0
+47 0 0.01600000000004253 0
+48 0 0.01200000000003156 0
+49 0 0.008000000000021934 0
+50 0 0.004000000000010155 0
+51 0.004585786437623847 0.01741421356237004 0
+52 0.007414213562370057 0.01741421356237614 0
+53 0.007414213562374714 0.01458578643762852 0
+54 0.004585786437629543 0.01458578643762427 0
+55 0.01258578643762438 0.02541421356237057 0
+56 0.01541421356237008 0.02541421356237611 0
+57 0.01541421356237474 0.02258578643762854 0
+58 0.01258578643763034 0.02258578643762347 0
+59 0.003999999999991045 0 0.001
+60 0.007999999999979423 0 0.001
+61 0.01199999999997901 0 0.001
+62 0.01599999999998917 0 0.001
+63 0.02 0.003999999999992779 0.001
+64 0.02 0.007999999999982091 0.001
+65 0.02 0.0119999999999705 0.001
+66 0.02 0.01599999999995885 0.001
+67 0.02 0.01999999999994786 0.001
+68 0.02 0.02399999999995802 0.001
+69 0.02 0.02799999999996886 0.001
+70 0.02 0.03199999999997834 0.001
+71 0.02 0.03599999999998998 0.001
+72 0.01600000000001043 0.04 0.001
+73 0.01200000000002099 0.04 0.001
+74 0.008000000000021266 0.04 0.001
+75 0.004000000000010967 0.04 0.001
+76 0 0.0360000000000088 0.001
+77 0 0.03200000000002086 0.001
+78 0 0.02800000000003047 0.001
+79 0 0.02400000000004199 0.001
+80 0 0.02000000000005283 0.001
+81 0 0.01600000000004253 0.001
+82 0 0.01200000000003156 0.001
+83 0 0.008000000000021934 0.001
+84 0 0.004000000000010155 0.001
+85 0.004585786437629543 0.01458578643762427 0.001
+86 0.007414213562374714 0.01458578643762852 0.001
+87 0.007414213562370057 0.01741421356237614 0.001
+88 0.004585786437623847 0.01741421356237004 0.001
+89 0.01258578643763034 0.02258578643762347 0.001
+90 0.01541421356237474 0.02258578643762854 0.001
+91 0.01541421356237008 0.02541421356237611 0.001
+92 0.01258578643762438 0.02541421356237057 0.001
+93 0 0 0.0005
+94 0.02 0 0.0005
+95 0.02 0.04 0.0005
+96 0 0.04 0.0005
+97 0.004 0.016 0.0005
+98 0.006 0.014 0.0005
+99 0.008 0.016 0.0005
+100 0.006 0.018 0.0005
+101 0.012 0.024 0.0005
+102 0.014 0.022 0.0005
+103 0.016 0.024 0.0005
+104 0.014 0.026 0.0005
+105 0.01069491464816358 0.03275706766290423 0
+106 0.009305085351833412 0.007242932337095425 0
+107 0.01195336568324535 0.01850333395213949 0
+108 0.01169692766308795 0.01273881552073794 0
+109 0.008268255376061504 0.02156696680671535 0
+110 0.008303072336911341 0.0272611844792605 0
+111 0.01187364030974172 0.02913348640269208 0
+112 0.008126359690256567 0.01086651359730744 0
+113 0.003890646163783284 0.02210218278604862 0
+114 0.01610129463824834 0.01788473837026788 0
+115 0.00527471461890471 0.0312856638364231 0
+116 0.01472528538109378 0.008714336163572449 0
+117 0.01623634084102215 0.03000849594336865 0
+118 0.003763659158977146 0.009991504056632143 0
+119 0.01443926882698468 0.03459017605072595 0
+120 0.005560731173012773 0.005409823949275548 0
+121 0.0107927702276238 0.02079995987352916 0
+122 0.009311828935983524 0.02433847876298567 0
+123 0.010650184926707 0.01597034408016556 0
+124 0.009488676316585601 0.01873394905844635 0
+125 0.007122237134825811 0.0361465517015278 0
+126 0.01287776286517033 0.003853448298470206 0
+127 0.008170792723072185 0.02983066727881701 0
+128 0.01182920727692591 0.01016933272118098 0
+129 0.006054294952975383 0.02377714633014847 0
+130 0.01401598402785485 0.01623087328408666 0
+131 0.007294908290742031 0.01983327956720814 0
+132 0.003027281658933691 0.02681472674903151 0
+133 0.01697271834106652 0.0131852732509674 0
+134 0.01069523874921455 0.02713050105997452 0
+135 0.009304761250784256 0.01286949894002464 0
+136 0.01275290459629124 0.02014826212529715 0
+137 0.01086648419919771 0.02253493272183302 0
+138 0.008075124555504954 0.003752623987231804 0
+139 0.01192487544449155 0.0362473760127683 0
+140 0.004476527263781469 0.01967798854167316 0
+141 0.01563081723820658 0.02008827489546555 0
+142 0.01340073119304642 0.03170462382804348 0
+143 0.006599268806951707 0.00829537617195707 0
+144 0.01261348915701591 0.02734733328150114 0
+145 0.00738651084298308 0.01265266671849831 0
+146 0.002370572056322577 0.01441224443541275 0
+147 0.01762942794367732 0.02558775556458744 0
+148 0.01790285673843247 0.0223785889447873 0
+149 0.002097143261567581 0.01762141105521092 0
+150 0.009240952099717 0.01705617258380825 0
+151 0.01382198503784449 0.01859814134267056 0
+152 0.00625056340314818 0.02138407985407196 0
+153 0.00387890972637342 0.03353716237557244 0
+154 0.01612109027362418 0.006462837624426198 0
+155 0.01088597704851552 0.02528987054004267 0
+156 0.009114022951483685 0.01471012945995643 0
+157 0.01543332795189227 0.02746035978079723 0
+158 0.01396135682223692 0.02895194262264064 0
+159 0.004566672048107451 0.01253964021920341 0
+160 0.006038643177761786 0.01104805737735952 0
+161 0.005504168720476045 0.02848787898487717 0
+162 0.01449583127952336 0.01151212101512062 0
+163 0.007738762426519678 0.03238279071413182 0
+164 0.01226123757347748 0.007617209285866475 0
+165 0.009183244138059184 0.02006818767405194 0
+166 0.01705248754588464 0.009979094675673207 0
+167 0.002947512454114644 0.03002090532432592 0
+168 0.01627553539450034 0.01547401350357579 0
+169 0.003736908562072098 0.02452497465053428 0
+170 0.01210320004370288 0.02141430992964029 0
+171 0.01457530738257165 0.03732467934410703 0
+172 0.005424692617424976 0.002675320655894944 0
+173 0.01688464871030633 0.03342212420093356 0
+174 0.003115351289692559 0.006577875799067628 0
+175 0.0107562083351965 0.01938171418718951 0
+176 0.008019077390322949 0.01859656757431411 0
+177 0.009706092248953281 0.03497450834836441 0
+178 0.01029390775104271 0.005025491651635149 0
+179 0.01041939478605073 0.03065899224442401 0
+180 0.009580605213946953 0.009341007755575252 0
+181 0.005956735004202665 0.01958735532173931 0
+182 0.007863373369673484 0.02328580803017795 0
+183 0.01213662663032722 0.01671419196982106 0
+184 0.009892427780998105 0.02880296629303363 0
+185 0.01010757221900033 0.01119703370696525 0
+186 0.009309166871142308 0.02282005571001185 0
+187 0.0106150246213663 0.02392672936345621 0
+188 0.01237424154086746 0.03390544608957453 0
+189 0.007625758459129293 0.006094553910425497 0
+190 0.01168498587988745 0.02632369803138028 0
+191 0.008315014120111607 0.01367630196861893 0
+192 0.005728816107385676 0.02596630605182034 0
+193 0.0142849710488482 0.01403486190956941 0
+194 0.004600229372246293 0.0371367428154218 0
+195 0.01539977062775254 0.002863257184577836 0
+196 0.01424508695900599 0.02065138529069584 0
+197 0.01070433566529644 0.01772661763859503 0
+198 0.01031585238516193 0.00227338267817154 0
+199 0.009684147614835145 0.0377266173218274 0
+200 0.009388657496714088 0.02579039199186671 0
+201 0.01060374527582381 0.01427137257676252 0
+202 0.01698320415331218 0.03688949659929274 0
+203 0.003016795846686892 0.003110503400708045 0
+204 0.0176937946876085 0.02816647553920276 0
+205 0.002306205312391421 0.01183352446079793 0
+206 0.01771383789648296 0.01965170741782519 0
+207 0.002297981782566726 0.02031890979568184 0
+208 0.01402312665197501 0.02741087956810794 0
+209 0.005976873348024324 0.01258912043189244 0
+210 0.0122314065954415 0.01499340989019052 0
+211 0.007775007199940584 0.02506988594104327 0
+212 0.009683984162416901 0.02155802055722827 0
+213 0.003364359686384905 0.01875813073873399 0
+214 0.0166517312501388 0.02127130301017035 0
+215 0.003999999999991045 0 0.0005
+216 0.007999999999979423 0 0.0005
+217 0.01199999999997901 0 0.0005
+218 0.01599999999998917 0 0.0005
+219 0.02 0.003999999999992779 0.0005
+220 0.02 0.007999999999982091 0.0005
+221 0.02 0.0119999999999705 0.0005
+222 0.02 0.01599999999995885 0.0005
+223 0.02 0.01999999999994786 0.0005
+224 0.02 0.02399999999995802 0.0005
+225 0.02 0.02799999999996886 0.0005
+226 0.02 0.03199999999997834 0.0005
+227 0.02 0.03599999999998998 0.0005
+228 0.01600000000001043 0.04 0.0005
+229 0.01200000000002099 0.04 0.0005
+230 0.008000000000021266 0.04 0.0005
+231 0.004000000000010967 0.04 0.0005
+232 0 0.0360000000000088 0.0005
+233 0 0.03200000000002086 0.0005
+234 0 0.02800000000003047 0.0005
+235 0 0.02400000000004199 0.0005
+236 0 0.02000000000005283 0.0005
+237 0 0.01600000000004253 0.0005
+238 0 0.01200000000003156 0.0005
+239 0 0.008000000000021934 0.0005
+240 0 0.004000000000010155 0.0005
+241 0.004585786437629543 0.01458578643762427 0.0005
+242 0.007414213562374714 0.01458578643762852 0.0005
+243 0.007414213562370057 0.01741421356237614 0.0005
+244 0.004585786437623847 0.01741421356237004 0.0005
+245 0.01258578643763034 0.02258578643762347 0.0005
+246 0.01541421356237474 0.02258578643762854 0.0005
+247 0.01541421356237008 0.02541421356237611 0.0005
+248 0.01258578643762438 0.02541421356237057 0.0005
+249 0.01069491464816358 0.03275706766290423 0.001
+250 0.009305085351833412 0.007242932337095425 0.001
+251 0.01195336568324535 0.01850333395213949 0.001
+252 0.01169692766308795 0.01273881552073794 0.001
+253 0.008268255376061504 0.02156696680671535 0.001
+254 0.008303072336911341 0.0272611844792605 0.001
+255 0.01187364030974172 0.02913348640269208 0.001
+256 0.008126359690256567 0.01086651359730744 0.001
+257 0.003890646163783284 0.02210218278604862 0.001
+258 0.01610129463824834 0.01788473837026788 0.001
+259 0.00527471461890471 0.0312856638364231 0.001
+260 0.01472528538109378 0.008714336163572449 0.001
+261 0.01623634084102215 0.03000849594336865 0.001
+262 0.003763659158977146 0.009991504056632143 0.001
+263 0.01443926882698468 0.03459017605072595 0.001
+264 0.005560731173012773 0.005409823949275548 0.001
+265 0.0107927702276238 0.02079995987352916 0.001
+266 0.009311828935983524 0.02433847876298567 0.001
+267 0.010650184926707 0.01597034408016556 0.001
+268 0.009488676316585601 0.01873394905844635 0.001
+269 0.007122237134825811 0.0361465517015278 0.001
+270 0.01287776286517033 0.003853448298470206 0.001
+271 0.008170792723072185 0.02983066727881701 0.001
+272 0.01182920727692591 0.01016933272118098 0.001
+273 0.006054294952975383 0.02377714633014847 0.001
+274 0.01401598402785485 0.01623087328408666 0.001
+275 0.007294908290742031 0.01983327956720814 0.001
+276 0.003027281658933691 0.02681472674903151 0.001
+277 0.01697271834106652 0.0131852732509674 0.001
+278 0.01069523874921455 0.02713050105997452 0.001
+279 0.009304761250784256 0.01286949894002464 0.001
+280 0.01275290459629124 0.02014826212529715 0.001
+281 0.01086648419919771 0.02253493272183302 0.001
+282 0.008075124555504954 0.003752623987231804 0.001
+283 0.01192487544449155 0.0362473760127683 0.001
+284 0.004476527263781469 0.01967798854167316 0.001
+285 0.01563081723820658 0.02008827489546555 0.001
+286 0.01340073119304642 0.03170462382804348 0.001
+287 0.006599268806951707 0.00829537617195707 0.001
+288 0.01261348915701591 0.02734733328150114 0.001
+289 0.00738651084298308 0.01265266671849831 0.001
+290 0.002370572056322577 0.01441224443541275 0.001
+291 0.01762942794367732 0.02558775556458744 0.001
+292 0.01790285673843247 0.0223785889447873 0.001
+293 0.002097143261567581 0.01762141105521092 0.001
+294 0.009240952099717 0.01705617258380825 0.001
+295 0.01382198503784449 0.01859814134267056 0.001
+296 0.00625056340314818 0.02138407985407196 0.001
+297 0.00387890972637342 0.03353716237557244 0.001
+298 0.01612109027362418 0.006462837624426198 0.001
+299 0.01088597704851552 0.02528987054004267 0.001
+300 0.009114022951483685 0.01471012945995643 0.001
+301 0.01543332795189227 0.02746035978079723 0.001
+302 0.01396135682223692 0.02895194262264064 0.001
+303 0.004566672048107451 0.01253964021920341 0.001
+304 0.006038643177761786 0.01104805737735952 0.001
+305 0.005504168720476045 0.02848787898487717 0.001
+306 0.01449583127952336 0.01151212101512062 0.001
+307 0.007738762426519678 0.03238279071413182 0.001
+308 0.01226123757347748 0.007617209285866475 0.001
+309 0.009183244138059184 0.02006818767405194 0.001
+310 0.01705248754588464 0.009979094675673207 0.001
+311 0.002947512454114644 0.03002090532432592 0.001
+312 0.01627553539450034 0.01547401350357579 0.001
+313 0.003736908562072098 0.02452497465053428 0.001
+314 0.01210320004370288 0.02141430992964029 0.001
+315 0.01457530738257165 0.03732467934410703 0.001
+316 0.005424692617424976 0.002675320655894944 0.001
+317 0.01688464871030633 0.03342212420093356 0.001
+318 0.003115351289692559 0.006577875799067628 0.001
+319 0.0107562083351965 0.01938171418718951 0.001
+320 0.008019077390322949 0.01859656757431411 0.001
+321 0.009706092248953281 0.03497450834836441 0.001
+322 0.01029390775104271 0.005025491651635149 0.001
+323 0.01041939478605073 0.03065899224442401 0.001
+324 0.009580605213946953 0.009341007755575252 0.001
+325 0.005956735004202665 0.01958735532173931 0.001
+326 0.007863373369673484 0.02328580803017795 0.001
+327 0.01213662663032722 0.01671419196982106 0.001
+328 0.009892427780998105 0.02880296629303363 0.001
+329 0.01010757221900033 0.01119703370696525 0.001
+330 0.009309166871142308 0.02282005571001185 0.001
+331 0.0106150246213663 0.02392672936345621 0.001
+332 0.01237424154086746 0.03390544608957453 0.001
+333 0.007625758459129293 0.006094553910425497 0.001
+334 0.01168498587988745 0.02632369803138028 0.001
+335 0.008315014120111607 0.01367630196861893 0.001
+336 0.005728816107385676 0.02596630605182034 0.001
+337 0.0142849710488482 0.01403486190956941 0.001
+338 0.004600229372246293 0.0371367428154218 0.001
+339 0.01539977062775254 0.002863257184577836 0.001
+340 0.01424508695900599 0.02065138529069584 0.001
+341 0.01070433566529644 0.01772661763859503 0.001
+342 0.01031585238516193 0.00227338267817154 0.001
+343 0.009684147614835145 0.0377266173218274 0.001
+344 0.009388657496714088 0.02579039199186671 0.001
+345 0.01060374527582381 0.01427137257676252 0.001
+346 0.01698320415331218 0.03688949659929274 0.001
+347 0.003016795846686892 0.003110503400708045 0.001
+348 0.0176937946876085 0.02816647553920276 0.001
+349 0.002306205312391421 0.01183352446079793 0.001
+350 0.01771383789648296 0.01965170741782519 0.001
+351 0.002297981782566726 0.02031890979568184 0.001
+352 0.01402312665197501 0.02741087956810794 0.001
+353 0.005976873348024324 0.01258912043189244 0.001
+354 0.0122314065954415 0.01499340989019052 0.001
+355 0.007775007199940584 0.02506988594104327 0.001
+356 0.009683984162416901 0.02155802055722827 0.001
+357 0.003364359686384905 0.01875813073873399 0.001
+358 0.0166517312501388 0.02127130301017035 0.001
+359 0.01069491464816358 0.03275706766290423 0.0005
+360 0.009305085351833412 0.007242932337095425 0.0005
+361 0.01195336568324535 0.01850333395213949 0.0005
+362 0.01169692766308795 0.01273881552073794 0.0005
+363 0.008268255376061504 0.02156696680671535 0.0005
+364 0.008303072336911341 0.0272611844792605 0.0005
+365 0.01187364030974172 0.02913348640269208 0.0005
+366 0.008126359690256567 0.01086651359730744 0.0005
+367 0.003890646163783284 0.02210218278604862 0.0005
+368 0.01610129463824834 0.01788473837026788 0.0005
+369 0.00527471461890471 0.0312856638364231 0.0005
+370 0.01472528538109378 0.008714336163572449 0.0005
+371 0.01623634084102215 0.03000849594336865 0.0005
+372 0.003763659158977146 0.009991504056632143 0.0005
+373 0.01443926882698468 0.03459017605072595 0.0005
+374 0.005560731173012773 0.005409823949275548 0.0005
+375 0.0107927702276238 0.02079995987352916 0.0005
+376 0.009311828935983524 0.02433847876298567 0.0005
+377 0.010650184926707 0.01597034408016556 0.0005
+378 0.009488676316585601 0.01873394905844635 0.0005
+379 0.007122237134825811 0.0361465517015278 0.0005
+380 0.01287776286517033 0.003853448298470206 0.0005
+381 0.008170792723072185 0.02983066727881701 0.0005
+382 0.01182920727692591 0.01016933272118098 0.0005
+383 0.006054294952975383 0.02377714633014847 0.0005
+384 0.01401598402785485 0.01623087328408666 0.0005
+385 0.007294908290742031 0.01983327956720814 0.0005
+386 0.003027281658933691 0.02681472674903151 0.0005
+387 0.01697271834106652 0.0131852732509674 0.0005
+388 0.01069523874921455 0.02713050105997452 0.0005
+389 0.009304761250784256 0.01286949894002464 0.0005
+390 0.01275290459629124 0.02014826212529715 0.0005
+391 0.01086648419919771 0.02253493272183302 0.0005
+392 0.008075124555504954 0.003752623987231804 0.0005
+393 0.01192487544449155 0.0362473760127683 0.0005
+394 0.004476527263781469 0.01967798854167316 0.0005
+395 0.01563081723820658 0.02008827489546555 0.0005
+396 0.01340073119304642 0.03170462382804348 0.0005
+397 0.006599268806951707 0.00829537617195707 0.0005
+398 0.01261348915701591 0.02734733328150114 0.0005
+399 0.00738651084298308 0.01265266671849831 0.0005
+400 0.002370572056322577 0.01441224443541275 0.0005
+401 0.01762942794367732 0.02558775556458744 0.0005
+402 0.01790285673843247 0.0223785889447873 0.0005
+403 0.002097143261567581 0.01762141105521092 0.0005
+404 0.009240952099717 0.01705617258380825 0.0005
+405 0.01382198503784449 0.01859814134267056 0.0005
+406 0.00625056340314818 0.02138407985407196 0.0005
+407 0.00387890972637342 0.03353716237557244 0.0005
+408 0.01612109027362418 0.006462837624426198 0.0005
+409 0.01088597704851552 0.02528987054004267 0.0005
+410 0.009114022951483685 0.01471012945995643 0.0005
+411 0.01543332795189227 0.02746035978079723 0.0005
+412 0.01396135682223692 0.02895194262264064 0.0005
+413 0.004566672048107451 0.01253964021920341 0.0005
+414 0.006038643177761786 0.01104805737735952 0.0005
+415 0.005504168720476045 0.02848787898487717 0.0005
+416 0.01449583127952336 0.01151212101512062 0.0005
+417 0.007738762426519678 0.03238279071413182 0.0005
+418 0.01226123757347748 0.007617209285866475 0.0005
+419 0.009183244138059184 0.02006818767405194 0.0005
+420 0.01705248754588464 0.009979094675673207 0.0005
+421 0.002947512454114644 0.03002090532432592 0.0005
+422 0.01627553539450034 0.01547401350357579 0.0005
+423 0.003736908562072098 0.02452497465053428 0.0005
+424 0.01210320004370288 0.02141430992964029 0.0005
+425 0.01457530738257165 0.03732467934410703 0.0005
+426 0.005424692617424976 0.002675320655894944 0.0005
+427 0.01688464871030633 0.03342212420093356 0.0005
+428 0.003115351289692559 0.006577875799067628 0.0005
+429 0.0107562083351965 0.01938171418718951 0.0005
+430 0.008019077390322949 0.01859656757431411 0.0005
+431 0.009706092248953281 0.03497450834836441 0.0005
+432 0.01029390775104271 0.005025491651635149 0.0005
+433 0.01041939478605073 0.03065899224442401 0.0005
+434 0.009580605213946953 0.009341007755575252 0.0005
+435 0.005956735004202665 0.01958735532173931 0.0005
+436 0.007863373369673484 0.02328580803017795 0.0005
+437 0.01213662663032722 0.01671419196982106 0.0005
+438 0.009892427780998105 0.02880296629303363 0.0005
+439 0.01010757221900033 0.01119703370696525 0.0005
+440 0.009309166871142308 0.02282005571001185 0.0005
+441 0.0106150246213663 0.02392672936345621 0.0005
+442 0.01237424154086746 0.03390544608957453 0.0005
+443 0.007625758459129293 0.006094553910425497 0.0005
+444 0.01168498587988745 0.02632369803138028 0.0005
+445 0.008315014120111607 0.01367630196861893 0.0005
+446 0.005728816107385676 0.02596630605182034 0.0005
+447 0.0142849710488482 0.01403486190956941 0.0005
+448 0.004600229372246293 0.0371367428154218 0.0005
+449 0.01539977062775254 0.002863257184577836 0.0005
+450 0.01424508695900599 0.02065138529069584 0.0005
+451 0.01070433566529644 0.01772661763859503 0.0005
+452 0.01031585238516193 0.00227338267817154 0.0005
+453 0.009684147614835145 0.0377266173218274 0.0005
+454 0.009388657496714088 0.02579039199186671 0.0005
+455 0.01060374527582381 0.01427137257676252 0.0005
+456 0.01698320415331218 0.03688949659929274 0.0005
+457 0.003016795846686892 0.003110503400708045 0.0005
+458 0.0176937946876085 0.02816647553920276 0.0005
+459 0.002306205312391421 0.01183352446079793 0.0005
+460 0.01771383789648296 0.01965170741782519 0.0005
+461 0.002297981782566726 0.02031890979568184 0.0005
+462 0.01402312665197501 0.02741087956810794 0.0005
+463 0.005976873348024324 0.01258912043189244 0.0005
+464 0.0122314065954415 0.01499340989019052 0.0005
+465 0.007775007199940584 0.02506988594104327 0.0005
+466 0.009683984162416901 0.02155802055722827 0.0005
+467 0.003364359686384905 0.01875813073873399 0.0005
+468 0.0166517312501388 0.02127130301017035 0.0005
+$EndNodes
+$Elements
+836
+1 15 2 20 1 1
+2 15 2 19 4 4
+3 1 2 17 1 1 25
+4 1 2 17 1 25 26
+5 1 2 17 1 26 27
+6 1 2 17 1 27 28
+7 1 2 17 1 28 2
+8 1 2 18 3 3 38
+9 1 2 18 3 38 39
+10 1 2 18 3 39 40
+11 1 2 18 3 40 41
+12 1 2 18 3 41 4
+13 2 2 16 11 12 196 57
+14 2 2 16 11 153 115 163
+15 2 2 16 11 154 116 164
+16 2 2 16 11 118 205 49
+17 2 2 16 11 117 204 36
+18 2 2 16 11 41 42 194
+19 2 2 16 11 28 29 195
+20 2 2 16 11 4 42 41
+21 2 2 16 11 2 29 28
+22 2 2 16 11 32 168 133
+23 2 2 16 11 45 169 132
+24 2 2 16 11 6 176 131
+25 2 2 16 11 126 154 164
+26 2 2 16 11 125 153 163
+27 2 2 16 11 56 147 157
+28 2 2 16 11 54 146 159
+29 2 2 16 11 121 175 136
+30 2 2 16 11 7 156 150
+31 2 2 16 11 12 170 136
+32 2 2 16 11 107 136 175
+33 2 2 16 11 52 150 176
+34 2 2 16 11 11 148 147
+35 2 2 16 11 5 149 146
+36 2 2 16 11 110 134 184
+37 2 2 16 11 108 135 185
+38 2 2 16 11 6 140 51
+39 2 2 16 11 5 146 54
+40 2 2 16 11 11 147 56
+41 2 2 16 11 5 51 149
+42 2 2 16 11 11 57 148
+43 2 2 16 11 6 52 176
+44 2 2 16 11 110 184 127
+45 2 2 16 11 108 185 128
+46 2 2 16 11 49 174 118
+47 2 2 16 11 36 173 117
+48 2 2 16 11 32 206 114
+49 2 2 16 11 45 207 113
+50 2 2 16 11 7 150 52
+51 2 2 16 11 111 134 144
+52 2 2 16 11 112 135 145
+53 2 2 16 11 131 176 165
+54 2 2 16 11 116 154 166
+55 2 2 16 11 115 153 167
+56 2 2 16 11 111 142 179
+57 2 2 16 11 112 143 180
+58 2 2 16 11 49 205 48
+59 2 2 16 11 36 204 35
+60 2 2 16 11 31 32 133
+61 2 2 16 11 44 45 132
+62 2 2 16 11 45 113 169
+63 2 2 16 11 32 114 168
+64 2 2 16 11 124 165 176
+65 2 2 16 11 117 142 158
+66 2 2 16 11 118 143 160
+67 2 2 16 11 121 136 170
+68 2 2 16 11 12 58 170
+69 2 2 16 11 12 136 196
+70 2 2 16 11 109 131 165
+71 2 2 16 11 109 152 131
+72 2 2 16 11 123 150 156
+73 2 2 16 11 111 158 142
+74 2 2 16 11 112 160 143
+75 2 2 16 11 6 131 181
+76 2 2 16 11 10 144 55
+77 2 2 16 11 8 145 53
+78 2 2 16 11 124 176 150
+79 2 2 16 11 50 203 174
+80 2 2 16 11 37 202 173
+81 2 2 16 11 196 141 57
+82 2 2 16 11 117 173 142
+83 2 2 16 11 118 174 143
+84 2 2 16 11 30 166 154
+85 2 2 16 11 43 167 153
+86 2 2 16 11 107 151 136
+87 2 2 16 11 9 137 58
+88 2 2 16 11 123 197 150
+89 2 2 16 11 141 196 151
+90 2 2 16 11 110 127 161
+91 2 2 16 11 108 128 162
+92 2 2 16 11 135 191 145
+93 2 2 16 11 134 190 144
+94 2 2 16 11 40 194 125
+95 2 2 16 11 27 195 126
+96 2 2 16 11 49 50 174
+97 2 2 16 11 36 37 173
+98 2 2 16 11 1 25 203
+99 2 2 16 11 3 38 202
+100 2 2 16 11 123 183 197
+101 2 2 16 11 136 151 196
+102 2 2 16 11 58 137 170
+103 2 2 16 11 105 179 142
+104 2 2 16 11 106 180 143
+105 2 2 16 11 6 181 140
+106 2 2 16 11 10 56 157
+107 2 2 16 11 8 54 159
+108 2 2 16 11 116 162 128
+109 2 2 16 11 115 161 127
+110 2 2 16 11 110 200 134
+111 2 2 16 11 108 201 135
+112 2 2 16 11 40 125 199
+113 2 2 16 11 27 126 198
+114 2 2 16 11 111 184 134
+115 2 2 16 11 112 185 135
+116 2 2 16 11 29 154 195
+117 2 2 16 11 42 153 194
+118 2 2 16 11 116 128 164
+119 2 2 16 11 115 127 163
+120 2 2 16 11 1 203 50
+121 2 2 16 11 3 202 37
+122 2 2 16 11 121 170 137
+123 2 2 16 11 109 182 152
+124 2 2 16 11 119 142 173
+125 2 2 16 11 120 143 174
+126 2 2 16 11 127 179 163
+127 2 2 16 11 128 180 164
+128 2 2 16 11 57 214 148
+129 2 2 16 11 51 213 149
+130 2 2 16 11 140 181 152
+131 2 2 16 11 29 30 154
+132 2 2 16 11 42 43 153
+133 2 2 16 11 131 152 181
+134 2 2 16 11 133 162 166
+135 2 2 16 11 132 161 167
+136 2 2 16 11 105 163 179
+137 2 2 16 11 106 164 180
+138 2 2 16 11 119 139 188
+139 2 2 16 11 120 138 189
+140 2 2 16 11 40 41 194
+141 2 2 16 11 27 28 195
+142 2 2 16 11 9 55 155
+143 2 2 16 11 7 53 156
+144 2 2 16 11 32 33 206
+145 2 2 16 11 45 46 207
+146 2 2 16 11 107 197 183
+147 2 2 16 11 114 151 130
+148 2 2 16 11 116 166 162
+149 2 2 16 11 115 167 161
+150 2 2 16 11 119 173 202
+151 2 2 16 11 120 174 203
+152 2 2 16 11 113 152 129
+153 2 2 16 11 107 183 151
+154 2 2 16 11 114 141 151
+155 2 2 16 11 26 138 172
+156 2 2 16 11 39 139 171
+157 2 2 16 11 147 204 157
+158 2 2 16 11 146 205 159
+159 2 2 16 11 119 171 139
+160 2 2 16 11 120 172 138
+161 2 2 16 11 55 144 190
+162 2 2 16 11 53 145 191
+163 2 2 16 11 31 133 166
+164 2 2 16 11 44 132 167
+165 2 2 16 11 9 187 137
+166 2 2 16 11 46 47 149
+167 2 2 16 11 33 34 148
+168 2 2 16 11 34 35 147
+169 2 2 16 11 47 48 146
+170 2 2 16 11 51 140 213
+171 2 2 16 11 57 141 214
+172 2 2 16 11 111 144 158
+173 2 2 16 11 112 145 160
+174 2 2 16 11 105 177 163
+175 2 2 16 11 106 178 164
+176 2 2 16 11 113 129 169
+177 2 2 16 11 114 130 168
+178 2 2 16 11 108 162 193
+179 2 2 16 11 110 161 192
+180 2 2 16 11 25 26 172
+181 2 2 16 11 38 39 171
+182 2 2 16 11 105 142 188
+183 2 2 16 11 106 143 189
+184 2 2 16 11 124 175 165
+185 2 2 16 11 113 140 152
+186 2 2 16 11 121 137 212
+187 2 2 16 11 10 208 144
+188 2 2 16 11 8 209 145
+189 2 2 16 11 121 165 175
+190 2 2 16 11 124 150 197
+191 2 2 16 11 125 163 177
+192 2 2 16 11 126 164 178
+193 2 2 16 11 129 152 182
+194 2 2 16 11 43 44 167
+195 2 2 16 11 30 31 166
+196 2 2 16 11 119 188 142
+197 2 2 16 11 120 189 143
+198 2 2 16 11 39 199 139
+199 2 2 16 11 26 198 138
+200 2 2 16 11 105 188 177
+201 2 2 16 11 106 189 178
+202 2 2 16 11 25 172 203
+203 2 2 16 11 38 171 202
+204 2 2 16 11 10 157 208
+205 2 2 16 11 8 159 209
+206 2 2 16 11 39 40 199
+207 2 2 16 11 26 27 198
+208 2 2 16 11 137 187 186
+209 2 2 16 11 133 193 162
+210 2 2 16 11 132 192 161
+211 2 2 16 11 122 187 155
+212 2 2 16 11 34 147 148
+213 2 2 16 11 47 146 149
+214 2 2 16 11 125 194 153
+215 2 2 16 11 126 195 154
+216 2 2 16 11 125 177 199
+217 2 2 16 11 126 178 198
+218 2 2 16 11 130 151 183
+219 2 2 16 11 9 155 187
+220 2 2 16 11 55 190 155
+221 2 2 16 11 53 191 156
+222 2 2 16 11 122 186 187
+223 2 2 16 11 122 155 200
+224 2 2 16 11 123 156 201
+225 2 2 16 11 144 208 158
+226 2 2 16 11 145 209 160
+227 2 2 16 11 139 177 188
+228 2 2 16 11 138 178 189
+229 2 2 16 11 130 193 168
+230 2 2 16 11 129 192 169
+231 2 2 16 11 111 179 184
+232 2 2 16 11 112 180 185
+233 2 2 16 11 113 207 140
+234 2 2 16 11 114 206 141
+235 2 2 16 11 129 211 192
+236 2 2 16 11 130 210 193
+237 2 2 16 11 46 149 207
+238 2 2 16 11 33 148 206
+239 2 2 16 11 117 158 157
+240 2 2 16 11 118 160 159
+241 2 2 16 11 109 165 212
+242 2 2 16 11 149 213 207
+243 2 2 16 11 148 214 206
+244 2 2 16 11 110 192 211
+245 2 2 16 11 108 193 210
+246 2 2 16 11 109 186 182
+247 2 2 16 11 121 212 165
+248 2 2 16 11 124 197 175
+249 2 2 16 11 123 201 210
+250 2 2 16 11 122 200 211
+251 2 2 16 11 35 204 147
+252 2 2 16 11 48 205 146
+253 2 2 16 11 119 202 171
+254 2 2 16 11 120 203 172
+255 2 2 16 11 157 158 208
+256 2 2 16 11 159 160 209
+257 2 2 16 11 140 207 213
+258 2 2 16 11 141 206 214
+259 2 2 16 11 108 210 201
+260 2 2 16 11 110 211 200
+261 2 2 16 11 134 200 155
+262 2 2 16 11 135 201 156
+263 2 2 16 11 133 168 193
+264 2 2 16 11 132 169 192
+265 2 2 16 11 134 155 190
+266 2 2 16 11 135 156 191
+267 2 2 16 11 122 182 186
+268 2 2 16 11 137 186 212
+269 2 2 16 11 138 198 178
+270 2 2 16 11 139 199 177
+271 2 2 16 11 107 175 197
+272 2 2 16 11 109 212 186
+273 2 2 16 11 122 211 182
+274 2 2 16 11 123 210 183
+275 2 2 16 11 127 184 179
+276 2 2 16 11 128 185 180
+277 2 2 16 11 129 182 211
+278 2 2 16 11 130 183 210
+279 2 2 16 11 117 157 204
+280 2 2 16 11 118 159 205
+281 3 2 84 37 1 25 215 93
+282 3 2 84 37 93 215 59 13
+283 3 2 84 37 25 26 216 215
+284 3 2 84 37 215 216 60 59
+285 3 2 84 37 26 27 217 216
+286 3 2 84 37 216 217 61 60
+287 3 2 84 37 27 28 218 217
+288 3 2 84 37 217 218 62 61
+289 3 2 84 37 28 2 94 218
+290 3 2 84 37 218 94 14 62
+291 3 2 85 45 3 38 228 95
+292 3 2 85 45 95 228 72 15
+293 3 2 85 45 38 39 229 228
+294 3 2 85 45 228 229 73 72
+295 3 2 85 45 39 40 230 229
+296 3 2 85 45 229 230 74 73
+297 3 2 85 45 40 41 231 230
+298 3 2 85 45 230 231 75 74
+299 3 2 85 45 41 4 96 231
+300 3 2 85 45 231 96 16 75
+301 6 2 83 1 57 12 196 246 102 450
+302 6 2 83 1 246 102 450 90 22 340
+303 6 2 83 1 163 153 115 417 407 369
+304 6 2 83 1 417 407 369 307 297 259
+305 6 2 83 1 164 154 116 418 408 370
+306 6 2 83 1 418 408 370 308 298 260
+307 6 2 83 1 49 118 205 239 372 459
+308 6 2 83 1 239 372 459 83 262 349
+309 6 2 83 1 36 117 204 226 371 458
+310 6 2 83 1 226 371 458 70 261 348
+311 6 2 83 1 194 41 42 448 231 232
+312 6 2 83 1 448 231 232 338 75 76
+313 6 2 83 1 195 28 29 449 218 219
+314 6 2 83 1 449 218 219 339 62 63
+315 6 2 83 1 41 4 42 231 96 232
+316 6 2 83 1 231 96 232 75 16 76
+317 6 2 83 1 28 2 29 218 94 219
+318 6 2 83 1 218 94 219 62 14 63
+319 6 2 83 1 133 32 168 387 222 422
+320 6 2 83 1 387 222 422 277 66 312
+321 6 2 83 1 132 45 169 386 235 423
+322 6 2 83 1 386 235 423 276 79 313
+323 6 2 83 1 131 6 176 385 100 430
+324 6 2 83 1 385 100 430 275 20 320
+325 6 2 83 1 164 126 154 418 380 408
+326 6 2 83 1 418 380 408 308 270 298
+327 6 2 83 1 163 125 153 417 379 407
+328 6 2 83 1 417 379 407 307 269 297
+329 6 2 83 1 157 56 147 411 247 401
+330 6 2 83 1 411 247 401 301 91 291
+331 6 2 83 1 159 54 146 413 241 400
+332 6 2 83 1 413 241 400 303 85 290
+333 6 2 83 1 136 121 175 390 375 429
+334 6 2 83 1 390 375 429 280 265 319
+335 6 2 83 1 150 7 156 404 99 410
+336 6 2 83 1 404 99 410 294 19 300
+337 6 2 83 1 136 12 170 390 102 424
+338 6 2 83 1 390 102 424 280 22 314
+339 6 2 83 1 175 107 136 429 361 390
+340 6 2 83 1 429 361 390 319 251 280
+341 6 2 83 1 176 52 150 430 243 404
+342 6 2 83 1 430 243 404 320 87 294
+343 6 2 83 1 147 11 148 401 103 402
+344 6 2 83 1 401 103 402 291 23 292
+345 6 2 83 1 146 5 149 400 97 403
+346 6 2 83 1 400 97 403 290 17 293
+347 6 2 83 1 184 110 134 438 364 388
+348 6 2 83 1 438 364 388 328 254 278
+349 6 2 83 1 185 108 135 439 362 389
+350 6 2 83 1 439 362 389 329 252 279
+351 6 2 83 1 51 6 140 244 100 394
+352 6 2 83 1 244 100 394 88 20 284
+353 6 2 83 1 54 5 146 241 97 400
+354 6 2 83 1 241 97 400 85 17 290
+355 6 2 83 1 56 11 147 247 103 401
+356 6 2 83 1 247 103 401 91 23 291
+357 6 2 83 1 149 5 51 403 97 244
+358 6 2 83 1 403 97 244 293 17 88
+359 6 2 83 1 148 11 57 402 103 246
+360 6 2 83 1 402 103 246 292 23 90
+361 6 2 83 1 176 6 52 430 100 243
+362 6 2 83 1 430 100 243 320 20 87
+363 6 2 83 1 127 110 184 381 364 438
+364 6 2 83 1 381 364 438 271 254 328
+365 6 2 83 1 128 108 185 382 362 439
+366 6 2 83 1 382 362 439 272 252 329
+367 6 2 83 1 118 49 174 372 239 428
+368 6 2 83 1 372 239 428 262 83 318
+369 6 2 83 1 117 36 173 371 226 427
+370 6 2 83 1 371 226 427 261 70 317
+371 6 2 83 1 114 32 206 368 222 460
+372 6 2 83 1 368 222 460 258 66 350
+373 6 2 83 1 113 45 207 367 235 461
+374 6 2 83 1 367 235 461 257 79 351
+375 6 2 83 1 52 7 150 243 99 404
+376 6 2 83 1 243 99 404 87 19 294
+377 6 2 83 1 144 111 134 398 365 388
+378 6 2 83 1 398 365 388 288 255 278
+379 6 2 83 1 145 112 135 399 366 389
+380 6 2 83 1 399 366 389 289 256 279
+381 6 2 83 1 165 131 176 419 385 430
+382 6 2 83 1 419 385 430 309 275 320
+383 6 2 83 1 166 116 154 420 370 408
+384 6 2 83 1 420 370 408 310 260 298
+385 6 2 83 1 167 115 153 421 369 407
+386 6 2 83 1 421 369 407 311 259 297
+387 6 2 83 1 179 111 142 433 365 396
+388 6 2 83 1 433 365 396 323 255 286
+389 6 2 83 1 180 112 143 434 366 397
+390 6 2 83 1 434 366 397 324 256 287
+391 6 2 83 1 48 49 205 238 239 459
+392 6 2 83 1 238 239 459 82 83 349
+393 6 2 83 1 35 36 204 225 226 458
+394 6 2 83 1 225 226 458 69 70 348
+395 6 2 83 1 133 31 32 387 221 222
+396 6 2 83 1 387 221 222 277 65 66
+397 6 2 83 1 132 44 45 386 234 235
+398 6 2 83 1 386 234 235 276 78 79
+399 6 2 83 1 169 45 113 423 235 367
+400 6 2 83 1 423 235 367 313 79 257
+401 6 2 83 1 168 32 114 422 222 368
+402 6 2 83 1 422 222 368 312 66 258
+403 6 2 83 1 176 124 165 430 378 419
+404 6 2 83 1 430 378 419 320 268 309
+405 6 2 83 1 158 117 142 412 371 396
+406 6 2 83 1 412 371 396 302 261 286
+407 6 2 83 1 160 118 143 414 372 397
+408 6 2 83 1 414 372 397 304 262 287
+409 6 2 83 1 170 121 136 424 375 390
+410 6 2 83 1 424 375 390 314 265 280
+411 6 2 83 1 170 12 58 424 102 245
+412 6 2 83 1 424 102 245 314 22 89
+413 6 2 83 1 196 12 136 450 102 390
+414 6 2 83 1 450 102 390 340 22 280
+415 6 2 83 1 165 109 131 419 363 385
+416 6 2 83 1 419 363 385 309 253 275
+417 6 2 83 1 131 109 152 385 363 406
+418 6 2 83 1 385 363 406 275 253 296
+419 6 2 83 1 156 123 150 410 377 404
+420 6 2 83 1 410 377 404 300 267 294
+421 6 2 83 1 142 111 158 396 365 412
+422 6 2 83 1 396 365 412 286 255 302
+423 6 2 83 1 143 112 160 397 366 414
+424 6 2 83 1 397 366 414 287 256 304
+425 6 2 83 1 181 6 131 435 100 385
+426 6 2 83 1 435 100 385 325 20 275
+427 6 2 83 1 55 10 144 248 104 398
+428 6 2 83 1 248 104 398 92 24 288
+429 6 2 83 1 53 8 145 242 98 399
+430 6 2 83 1 242 98 399 86 18 289
+431 6 2 83 1 150 124 176 404 378 430
+432 6 2 83 1 404 378 430 294 268 320
+433 6 2 83 1 174 50 203 428 240 457
+434 6 2 83 1 428 240 457 318 84 347
+435 6 2 83 1 173 37 202 427 227 456
+436 6 2 83 1 427 227 456 317 71 346
+437 6 2 83 1 57 196 141 246 450 395
+438 6 2 83 1 246 450 395 90 340 285
+439 6 2 83 1 142 117 173 396 371 427
+440 6 2 83 1 396 371 427 286 261 317
+441 6 2 83 1 143 118 174 397 372 428
+442 6 2 83 1 397 372 428 287 262 318
+443 6 2 83 1 154 30 166 408 220 420
+444 6 2 83 1 408 220 420 298 64 310
+445 6 2 83 1 153 43 167 407 233 421
+446 6 2 83 1 407 233 421 297 77 311
+447 6 2 83 1 136 107 151 390 361 405
+448 6 2 83 1 390 361 405 280 251 295
+449 6 2 83 1 58 9 137 245 101 391
+450 6 2 83 1 245 101 391 89 21 281
+451 6 2 83 1 150 123 197 404 377 451
+452 6 2 83 1 404 377 451 294 267 341
+453 6 2 83 1 151 141 196 405 395 450
+454 6 2 83 1 405 395 450 295 285 340
+455 6 2 83 1 161 110 127 415 364 381
+456 6 2 83 1 415 364 381 305 254 271
+457 6 2 83 1 162 108 128 416 362 382
+458 6 2 83 1 416 362 382 306 252 272
+459 6 2 83 1 145 135 191 399 389 445
+460 6 2 83 1 399 389 445 289 279 335
+461 6 2 83 1 144 134 190 398 388 444
+462 6 2 83 1 398 388 444 288 278 334
+463 6 2 83 1 125 40 194 379 230 448
+464 6 2 83 1 379 230 448 269 74 338
+465 6 2 83 1 126 27 195 380 217 449
+466 6 2 83 1 380 217 449 270 61 339
+467 6 2 83 1 174 49 50 428 239 240
+468 6 2 83 1 428 239 240 318 83 84
+469 6 2 83 1 173 36 37 427 226 227
+470 6 2 83 1 427 226 227 317 70 71
+471 6 2 83 1 203 1 25 457 93 215
+472 6 2 83 1 457 93 215 347 13 59
+473 6 2 83 1 202 3 38 456 95 228
+474 6 2 83 1 456 95 228 346 15 72
+475 6 2 83 1 197 123 183 451 377 437
+476 6 2 83 1 451 377 437 341 267 327
+477 6 2 83 1 196 136 151 450 390 405
+478 6 2 83 1 450 390 405 340 280 295
+479 6 2 83 1 170 58 137 424 245 391
+480 6 2 83 1 424 245 391 314 89 281
+481 6 2 83 1 142 105 179 396 359 433
+482 6 2 83 1 396 359 433 286 249 323
+483 6 2 83 1 143 106 180 397 360 434
+484 6 2 83 1 397 360 434 287 250 324
+485 6 2 83 1 140 6 181 394 100 435
+486 6 2 83 1 394 100 435 284 20 325
+487 6 2 83 1 157 10 56 411 104 247
+488 6 2 83 1 411 104 247 301 24 91
+489 6 2 83 1 159 8 54 413 98 241
+490 6 2 83 1 413 98 241 303 18 85
+491 6 2 83 1 128 116 162 382 370 416
+492 6 2 83 1 382 370 416 272 260 306
+493 6 2 83 1 127 115 161 381 369 415
+494 6 2 83 1 381 369 415 271 259 305
+495 6 2 83 1 134 110 200 388 364 454
+496 6 2 83 1 388 364 454 278 254 344
+497 6 2 83 1 135 108 201 389 362 455
+498 6 2 83 1 389 362 455 279 252 345
+499 6 2 83 1 199 40 125 453 230 379
+500 6 2 83 1 453 230 379 343 74 269
+501 6 2 83 1 198 27 126 452 217 380
+502 6 2 83 1 452 217 380 342 61 270
+503 6 2 83 1 134 111 184 388 365 438
+504 6 2 83 1 388 365 438 278 255 328
+505 6 2 83 1 135 112 185 389 366 439
+506 6 2 83 1 389 366 439 279 256 329
+507 6 2 83 1 195 29 154 449 219 408
+508 6 2 83 1 449 219 408 339 63 298
+509 6 2 83 1 194 42 153 448 232 407
+510 6 2 83 1 448 232 407 338 76 297
+511 6 2 83 1 164 116 128 418 370 382
+512 6 2 83 1 418 370 382 308 260 272
+513 6 2 83 1 163 115 127 417 369 381
+514 6 2 83 1 417 369 381 307 259 271
+515 6 2 83 1 50 1 203 240 93 457
+516 6 2 83 1 240 93 457 84 13 347
+517 6 2 83 1 37 3 202 227 95 456
+518 6 2 83 1 227 95 456 71 15 346
+519 6 2 83 1 137 121 170 391 375 424
+520 6 2 83 1 391 375 424 281 265 314
+521 6 2 83 1 152 109 182 406 363 436
+522 6 2 83 1 406 363 436 296 253 326
+523 6 2 83 1 173 119 142 427 373 396
+524 6 2 83 1 427 373 396 317 263 286
+525 6 2 83 1 174 120 143 428 374 397
+526 6 2 83 1 428 374 397 318 264 287
+527 6 2 83 1 163 127 179 417 381 433
+528 6 2 83 1 417 381 433 307 271 323
+529 6 2 83 1 164 128 180 418 382 434
+530 6 2 83 1 418 382 434 308 272 324
+531 6 2 83 1 148 57 214 402 246 468
+532 6 2 83 1 402 246 468 292 90 358
+533 6 2 83 1 149 51 213 403 244 467
+534 6 2 83 1 403 244 467 293 88 357
+535 6 2 83 1 152 140 181 406 394 435
+536 6 2 83 1 406 394 435 296 284 325
+537 6 2 83 1 154 29 30 408 219 220
+538 6 2 83 1 408 219 220 298 63 64
+539 6 2 83 1 153 42 43 407 232 233
+540 6 2 83 1 407 232 233 297 76 77
+541 6 2 83 1 181 131 152 435 385 406
+542 6 2 83 1 435 385 406 325 275 296
+543 6 2 83 1 166 133 162 420 387 416
+544 6 2 83 1 420 387 416 310 277 306
+545 6 2 83 1 167 132 161 421 386 415
+546 6 2 83 1 421 386 415 311 276 305
+547 6 2 83 1 179 105 163 433 359 417
+548 6 2 83 1 433 359 417 323 249 307
+549 6 2 83 1 180 106 164 434 360 418
+550 6 2 83 1 434 360 418 324 250 308
+551 6 2 83 1 188 119 139 442 373 393
+552 6 2 83 1 442 373 393 332 263 283
+553 6 2 83 1 189 120 138 443 374 392
+554 6 2 83 1 443 374 392 333 264 282
+555 6 2 83 1 194 40 41 448 230 231
+556 6 2 83 1 448 230 231 338 74 75
+557 6 2 83 1 195 27 28 449 217 218
+558 6 2 83 1 449 217 218 339 61 62
+559 6 2 83 1 155 9 55 409 101 248
+560 6 2 83 1 409 101 248 299 21 92
+561 6 2 83 1 156 7 53 410 99 242
+562 6 2 83 1 410 99 242 300 19 86
+563 6 2 83 1 206 32 33 460 222 223
+564 6 2 83 1 460 222 223 350 66 67
+565 6 2 83 1 207 45 46 461 235 236
+566 6 2 83 1 461 235 236 351 79 80
+567 6 2 83 1 183 107 197 437 361 451
+568 6 2 83 1 437 361 451 327 251 341
+569 6 2 83 1 130 114 151 384 368 405
+570 6 2 83 1 384 368 405 274 258 295
+571 6 2 83 1 162 116 166 416 370 420
+572 6 2 83 1 416 370 420 306 260 310
+573 6 2 83 1 161 115 167 415 369 421
+574 6 2 83 1 415 369 421 305 259 311
+575 6 2 83 1 202 119 173 456 373 427
+576 6 2 83 1 456 373 427 346 263 317
+577 6 2 83 1 203 120 174 457 374 428
+578 6 2 83 1 457 374 428 347 264 318
+579 6 2 83 1 129 113 152 383 367 406
+580 6 2 83 1 383 367 406 273 257 296
+581 6 2 83 1 151 107 183 405 361 437
+582 6 2 83 1 405 361 437 295 251 327
+583 6 2 83 1 151 114 141 405 368 395
+584 6 2 83 1 405 368 395 295 258 285
+585 6 2 83 1 172 26 138 426 216 392
+586 6 2 83 1 426 216 392 316 60 282
+587 6 2 83 1 171 39 139 425 229 393
+588 6 2 83 1 425 229 393 315 73 283
+589 6 2 83 1 157 147 204 411 401 458
+590 6 2 83 1 411 401 458 301 291 348
+591 6 2 83 1 159 146 205 413 400 459
+592 6 2 83 1 413 400 459 303 290 349
+593 6 2 83 1 139 119 171 393 373 425
+594 6 2 83 1 393 373 425 283 263 315
+595 6 2 83 1 138 120 172 392 374 426
+596 6 2 83 1 392 374 426 282 264 316
+597 6 2 83 1 190 55 144 444 248 398
+598 6 2 83 1 444 248 398 334 92 288
+599 6 2 83 1 191 53 145 445 242 399
+600 6 2 83 1 445 242 399 335 86 289
+601 6 2 83 1 166 31 133 420 221 387
+602 6 2 83 1 420 221 387 310 65 277
+603 6 2 83 1 167 44 132 421 234 386
+604 6 2 83 1 421 234 386 311 78 276
+605 6 2 83 1 137 9 187 391 101 441
+606 6 2 83 1 391 101 441 281 21 331
+607 6 2 83 1 149 46 47 403 236 237
+608 6 2 83 1 403 236 237 293 80 81
+609 6 2 83 1 148 33 34 402 223 224
+610 6 2 83 1 402 223 224 292 67 68
+611 6 2 83 1 147 34 35 401 224 225
+612 6 2 83 1 401 224 225 291 68 69
+613 6 2 83 1 146 47 48 400 237 238
+614 6 2 83 1 400 237 238 290 81 82
+615 6 2 83 1 213 51 140 467 244 394
+616 6 2 83 1 467 244 394 357 88 284
+617 6 2 83 1 214 57 141 468 246 395
+618 6 2 83 1 468 246 395 358 90 285
+619 6 2 83 1 158 111 144 412 365 398
+620 6 2 83 1 412 365 398 302 255 288
+621 6 2 83 1 160 112 145 414 366 399
+622 6 2 83 1 414 366 399 304 256 289
+623 6 2 83 1 163 105 177 417 359 431
+624 6 2 83 1 417 359 431 307 249 321
+625 6 2 83 1 164 106 178 418 360 432
+626 6 2 83 1 418 360 432 308 250 322
+627 6 2 83 1 169 113 129 423 367 383
+628 6 2 83 1 423 367 383 313 257 273
+629 6 2 83 1 168 114 130 422 368 384
+630 6 2 83 1 422 368 384 312 258 274
+631 6 2 83 1 193 108 162 447 362 416
+632 6 2 83 1 447 362 416 337 252 306
+633 6 2 83 1 192 110 161 446 364 415
+634 6 2 83 1 446 364 415 336 254 305
+635 6 2 83 1 172 25 26 426 215 216
+636 6 2 83 1 426 215 216 316 59 60
+637 6 2 83 1 171 38 39 425 228 229
+638 6 2 83 1 425 228 229 315 72 73
+639 6 2 83 1 188 105 142 442 359 396
+640 6 2 83 1 442 359 396 332 249 286
+641 6 2 83 1 189 106 143 443 360 397
+642 6 2 83 1 443 360 397 333 250 287
+643 6 2 83 1 165 124 175 419 378 429
+644 6 2 83 1 419 378 429 309 268 319
+645 6 2 83 1 152 113 140 406 367 394
+646 6 2 83 1 406 367 394 296 257 284
+647 6 2 83 1 212 121 137 466 375 391
+648 6 2 83 1 466 375 391 356 265 281
+649 6 2 83 1 144 10 208 398 104 462
+650 6 2 83 1 398 104 462 288 24 352
+651 6 2 83 1 145 8 209 399 98 463
+652 6 2 83 1 399 98 463 289 18 353
+653 6 2 83 1 175 121 165 429 375 419
+654 6 2 83 1 429 375 419 319 265 309
+655 6 2 83 1 197 124 150 451 378 404
+656 6 2 83 1 451 378 404 341 268 294
+657 6 2 83 1 177 125 163 431 379 417
+658 6 2 83 1 431 379 417 321 269 307
+659 6 2 83 1 178 126 164 432 380 418
+660 6 2 83 1 432 380 418 322 270 308
+661 6 2 83 1 182 129 152 436 383 406
+662 6 2 83 1 436 383 406 326 273 296
+663 6 2 83 1 167 43 44 421 233 234
+664 6 2 83 1 421 233 234 311 77 78
+665 6 2 83 1 166 30 31 420 220 221
+666 6 2 83 1 420 220 221 310 64 65
+667 6 2 83 1 142 119 188 396 373 442
+668 6 2 83 1 396 373 442 286 263 332
+669 6 2 83 1 143 120 189 397 374 443
+670 6 2 83 1 397 374 443 287 264 333
+671 6 2 83 1 139 39 199 393 229 453
+672 6 2 83 1 393 229 453 283 73 343
+673 6 2 83 1 138 26 198 392 216 452
+674 6 2 83 1 392 216 452 282 60 342
+675 6 2 83 1 177 105 188 431 359 442
+676 6 2 83 1 431 359 442 321 249 332
+677 6 2 83 1 178 106 189 432 360 443
+678 6 2 83 1 432 360 443 322 250 333
+679 6 2 83 1 203 25 172 457 215 426
+680 6 2 83 1 457 215 426 347 59 316
+681 6 2 83 1 202 38 171 456 228 425
+682 6 2 83 1 456 228 425 346 72 315
+683 6 2 83 1 208 10 157 462 104 411
+684 6 2 83 1 462 104 411 352 24 301
+685 6 2 83 1 209 8 159 463 98 413
+686 6 2 83 1 463 98 413 353 18 303
+687 6 2 83 1 199 39 40 453 229 230
+688 6 2 83 1 453 229 230 343 73 74
+689 6 2 83 1 198 26 27 452 216 217
+690 6 2 83 1 452 216 217 342 60 61
+691 6 2 83 1 186 137 187 440 391 441
+692 6 2 83 1 440 391 441 330 281 331
+693 6 2 83 1 162 133 193 416 387 447
+694 6 2 83 1 416 387 447 306 277 337
+695 6 2 83 1 161 132 192 415 386 446
+696 6 2 83 1 415 386 446 305 276 336
+697 6 2 83 1 155 122 187 409 376 441
+698 6 2 83 1 409 376 441 299 266 331
+699 6 2 83 1 148 34 147 402 224 401
+700 6 2 83 1 402 224 401 292 68 291
+701 6 2 83 1 149 47 146 403 237 400
+702 6 2 83 1 403 237 400 293 81 290
+703 6 2 83 1 153 125 194 407 379 448
+704 6 2 83 1 407 379 448 297 269 338
+705 6 2 83 1 154 126 195 408 380 449
+706 6 2 83 1 408 380 449 298 270 339
+707 6 2 83 1 199 125 177 453 379 431
+708 6 2 83 1 453 379 431 343 269 321
+709 6 2 83 1 198 126 178 452 380 432
+710 6 2 83 1 452 380 432 342 270 322
+711 6 2 83 1 183 130 151 437 384 405
+712 6 2 83 1 437 384 405 327 274 295
+713 6 2 83 1 187 9 155 441 101 409
+714 6 2 83 1 441 101 409 331 21 299
+715 6 2 83 1 155 55 190 409 248 444
+716 6 2 83 1 409 248 444 299 92 334
+717 6 2 83 1 156 53 191 410 242 445
+718 6 2 83 1 410 242 445 300 86 335
+719 6 2 83 1 187 122 186 441 376 440
+720 6 2 83 1 441 376 440 331 266 330
+721 6 2 83 1 200 122 155 454 376 409
+722 6 2 83 1 454 376 409 344 266 299
+723 6 2 83 1 201 123 156 455 377 410
+724 6 2 83 1 455 377 410 345 267 300
+725 6 2 83 1 158 144 208 412 398 462
+726 6 2 83 1 412 398 462 302 288 352
+727 6 2 83 1 160 145 209 414 399 463
+728 6 2 83 1 414 399 463 304 289 353
+729 6 2 83 1 188 139 177 442 393 431
+730 6 2 83 1 442 393 431 332 283 321
+731 6 2 83 1 189 138 178 443 392 432
+732 6 2 83 1 443 392 432 333 282 322
+733 6 2 83 1 168 130 193 422 384 447
+734 6 2 83 1 422 384 447 312 274 337
+735 6 2 83 1 169 129 192 423 383 446
+736 6 2 83 1 423 383 446 313 273 336
+737 6 2 83 1 184 111 179 438 365 433
+738 6 2 83 1 438 365 433 328 255 323
+739 6 2 83 1 185 112 180 439 366 434
+740 6 2 83 1 439 366 434 329 256 324
+741 6 2 83 1 140 113 207 394 367 461
+742 6 2 83 1 394 367 461 284 257 351
+743 6 2 83 1 141 114 206 395 368 460
+744 6 2 83 1 395 368 460 285 258 350
+745 6 2 83 1 192 129 211 446 383 465
+746 6 2 83 1 446 383 465 336 273 355
+747 6 2 83 1 193 130 210 447 384 464
+748 6 2 83 1 447 384 464 337 274 354
+749 6 2 83 1 207 46 149 461 236 403
+750 6 2 83 1 461 236 403 351 80 293
+751 6 2 83 1 206 33 148 460 223 402
+752 6 2 83 1 460 223 402 350 67 292
+753 6 2 83 1 157 117 158 411 371 412
+754 6 2 83 1 411 371 412 301 261 302
+755 6 2 83 1 159 118 160 413 372 414
+756 6 2 83 1 413 372 414 303 262 304
+757 6 2 83 1 212 109 165 466 363 419
+758 6 2 83 1 466 363 419 356 253 309
+759 6 2 83 1 207 149 213 461 403 467
+760 6 2 83 1 461 403 467 351 293 357
+761 6 2 83 1 206 148 214 460 402 468
+762 6 2 83 1 460 402 468 350 292 358
+763 6 2 83 1 211 110 192 465 364 446
+764 6 2 83 1 465 364 446 355 254 336
+765 6 2 83 1 210 108 193 464 362 447
+766 6 2 83 1 464 362 447 354 252 337
+767 6 2 83 1 182 109 186 436 363 440
+768 6 2 83 1 436 363 440 326 253 330
+769 6 2 83 1 165 121 212 419 375 466
+770 6 2 83 1 419 375 466 309 265 356
+771 6 2 83 1 175 124 197 429 378 451
+772 6 2 83 1 429 378 451 319 268 341
+773 6 2 83 1 210 123 201 464 377 455
+774 6 2 83 1 464 377 455 354 267 345
+775 6 2 83 1 211 122 200 465 376 454
+776 6 2 83 1 465 376 454 355 266 344
+777 6 2 83 1 147 35 204 401 225 458
+778 6 2 83 1 401 225 458 291 69 348
+779 6 2 83 1 146 48 205 400 238 459
+780 6 2 83 1 400 238 459 290 82 349
+781 6 2 83 1 171 119 202 425 373 456
+782 6 2 83 1 425 373 456 315 263 346
+783 6 2 83 1 172 120 203 426 374 457
+784 6 2 83 1 426 374 457 316 264 347
+785 6 2 83 1 208 157 158 462 411 412
+786 6 2 83 1 462 411 412 352 301 302
+787 6 2 83 1 209 159 160 463 413 414
+788 6 2 83 1 463 413 414 353 303 304
+789 6 2 83 1 213 140 207 467 394 461
+790 6 2 83 1 467 394 461 357 284 351
+791 6 2 83 1 214 141 206 468 395 460
+792 6 2 83 1 468 395 460 358 285 350
+793 6 2 83 1 201 108 210 455 362 464
+794 6 2 83 1 455 362 464 345 252 354
+795 6 2 83 1 200 110 211 454 364 465
+796 6 2 83 1 454 364 465 344 254 355
+797 6 2 83 1 155 134 200 409 388 454
+798 6 2 83 1 409 388 454 299 278 344
+799 6 2 83 1 156 135 201 410 389 455
+800 6 2 83 1 410 389 455 300 279 345
+801 6 2 83 1 193 133 168 447 387 422
+802 6 2 83 1 447 387 422 337 277 312
+803 6 2 83 1 192 132 169 446 386 423
+804 6 2 83 1 446 386 423 336 276 313
+805 6 2 83 1 190 134 155 444 388 409
+806 6 2 83 1 444 388 409 334 278 299
+807 6 2 83 1 191 135 156 445 389 410
+808 6 2 83 1 445 389 410 335 279 300
+809 6 2 83 1 186 122 182 440 376 436
+810 6 2 83 1 440 376 436 330 266 326
+811 6 2 83 1 212 137 186 466 391 440
+812 6 2 83 1 466 391 440 356 281 330
+813 6 2 83 1 178 138 198 432 392 452
+814 6 2 83 1 432 392 452 322 282 342
+815 6 2 83 1 177 139 199 431 393 453
+816 6 2 83 1 431 393 453 321 283 343
+817 6 2 83 1 197 107 175 451 361 429
+818 6 2 83 1 451 361 429 341 251 319
+819 6 2 83 1 186 109 212 440 363 466
+820 6 2 83 1 440 363 466 330 253 356
+821 6 2 83 1 182 122 211 436 376 465
+822 6 2 83 1 436 376 465 326 266 355
+823 6 2 83 1 183 123 210 437 377 464
+824 6 2 83 1 437 377 464 327 267 354
+825 6 2 83 1 179 127 184 433 381 438
+826 6 2 83 1 433 381 438 323 271 328
+827 6 2 83 1 180 128 185 434 382 439
+828 6 2 83 1 434 382 439 324 272 329
+829 6 2 83 1 211 129 182 465 383 436
+830 6 2 83 1 465 383 436 355 273 326
+831 6 2 83 1 210 130 183 464 384 437
+832 6 2 83 1 464 384 437 354 274 327
+833 6 2 83 1 204 117 157 458 371 411
+834 6 2 83 1 458 371 411 348 261 301
+835 6 2 83 1 205 118 159 459 372 413
+836 6 2 83 1 459 372 413 349 262 303
+$EndElements
diff --git a/dG3D/benchmarks/DG_PRI6_DiriToNeum/twoHole.py b/dG3D/benchmarks/DG_PRI6_DiriToNeum/twoHole.py
new file mode 100644
index 0000000000000000000000000000000000000000..1099e99bdf564019a37512103f2204a23d8e1268
--- /dev/null
+++ b/dG3D/benchmarks/DG_PRI6_DiriToNeum/twoHole.py
@@ -0,0 +1,105 @@
+#coding-Utf-8-*-
+from gmshpy import *
+from dG3Dpy import*
+
+#script to launch beam problem with a python script
+
+
+lawnum1 = 1 # unique number of law
+rho = 7850
+young = 28.9e9
+nu = 0.3
+sy0 = 150.e6
+h = young/50.
+
+
+# geometry
+meshfile="twoHole.msh" # name of mesh file
+
+# solver
+sol = 2 # Gmm=0 (default) Taucs=1 PETsc=2
+soltype = 1 # StaticLinear=0 (default) StaticNonLinear=1
+nstep = 25 # number of step (used only if soltype=1)
+ftime =1. # Final time (used only if soltype=1)
+tol=1.e-6 # relative tolerance for NR scheme (used only if soltype=1)
+nstepArch=1 # Number of step between 2 archiving (used only if soltype=1)
+
+
+
+harden = LinearExponentialJ2IsotropicHardening(1, sy0, h, 0., 10.)
+law1 = J2LinearDG3DMaterialLaw(lawnum1,rho,young,nu,harden)
+
+
+# creation of ElasticField
+beta1 = 1e2
+fullDG = True;
+averageStrainBased = False
+
+myfield1 = dG3DDomain(1000,83,0,lawnum1,fullDG,3)
+myfield1.stabilityParameters(beta1)
+myfield1.averageStrainBased(averageStrainBased)
+
+
+# creation of Solver
+mysolver = nonLinearMechSolver(1000)
+mysolver.loadModel(meshfile)
+mysolver.addDomain(myfield1)
+mysolver.addMaterialLaw(law1)
+
+mysolver.Scheme(soltype)
+mysolver.Solver(sol)
+mysolver.snlData(nstep,ftime,tol)
+mysolver.stepBetweenArchiving(nstepArch)
+mysolver.options("-ksp_type preonly -pc_type lu -pc_factor_mat_solver_package mumps")
+
+# BC
+#mysolver.displacementBC("Volume",83,2,0.)
+mysolver.displacementBC("Face",84,0,0.)
+mysolver.displacementBC("Face",84,1,0.)
+mysolver.displacementBC("Face",84,2,0.)
+mysolver.displacementBC("Face",85,0,0.)
+mysolver.displacementBC("Face",85,1,5e-4)
+mysolver.displacementBC("Face",85,2,0.)
+
+
+mysolver.internalPointBuildView("svm",IPField.SVM, 1, 1)
+mysolver.internalPointBuildView("sig_xx",IPField.SIG_XX, 1, 1)
+mysolver.internalPointBuildView("sig_yy",IPField.SIG_YY, 1, 1)
+mysolver.internalPointBuildView("sig_zz",IPField.SIG_ZZ, 1, 1)
+mysolver.internalPointBuildView("sig_xy",IPField.SIG_XY, 1, 1)
+mysolver.internalPointBuildView("sig_yz",IPField.SIG_YZ, 1, 1)
+mysolver.internalPointBuildView("sig_xz",IPField.SIG_XZ, 1, 1)
+mysolver.internalPointBuildView("epl",IPField.PLASTICSTRAIN, 1, 1)
+mysolver.internalPointBuildView("Damage",IPField.DAMAGE, 1, 1)
+mysolver.energyBuildView("elastic energy",2)
+mysolver.OneUnknownBuildView("disp0",0)
+
+mysolver.archivingForceOnPhysicalGroup("Face", 84, 1)
+mysolver.archivingNodeDisplacement(19,1,1)
+
+t1 = mysolver.solve()
+
+check = TestCheck()
+check.equal(-2.839032e+03,mysolver.getArchivedForceOnPhysicalGroup("Face", 84, 1),1.e-4)
+
+nstep = 5 # number of step (used only if soltype=1)
+ftime =1. # Final time (used only if soltype=1)
+tol=1.e-6 # relative tolerance for NR scheme (used only if soltype=1)
+mysolver.snlData(nstep,ftime,tol)
+
+lastF = 2.839032e+03/(20e-3*1e-3)
+
+fct = LinearFunctionTime(0,lastF,ftime,0.*lastF)
+mysolver.resetBoundaryConditions(False);
+
+mysolver.displacementBC("Face",84,0,0.)
+mysolver.displacementBC("Face",84,1,0.)
+mysolver.displacementBC("Face",84,2,0.)
+mysolver.displacementBC("Face",85,0,0.)
+mysolver.forceBC("Face",85,1,fct)
+mysolver.displacementBC("Face",85,2,0.)
+
+mysolver.solve()
+
+check = TestCheck()
+check.equal(-1.508328e-06,mysolver.getArchivedForceOnPhysicalGroup("Face", 84, 1),1.e-4)