]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Add new testcase intergrid_map.cc
authorwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 14 Nov 2000 16:46:26 +0000 (16:46 +0000)
committerwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 14 Nov 2000 16:46:26 +0000 (16:46 +0000)
git-svn-id: https://svn.dealii.org/trunk@3488 0785d39b-7218-0410-832d-ea1e28bc413d

tests/deal.II/Makefile.in
tests/deal.II/intergrid_map.cc [new file with mode: 0644]
tests/deal.II/intergrid_map.checked [new file with mode: 0644]
tests/deal.II/mglocal.cc

index b7018eb68dd604385630af0b3d427331971047e0..b36ac06ea8b99f27911bbb3e66447bcd3b865a00 100644 (file)
@@ -284,6 +284,21 @@ block_matrices.testcase: $(block_matrices-o-files) $(libraries)
        @$(CXX) $(LDFLAGS) -g  -o $@ $^ $(LIBS)
 
 
+############################################################
+
+
+intergrid_map-cc-files = intergrid_map.cc 
+
+ifeq ($(debug-mode),on)
+intergrid_map-o-files = $(intergrid_map-cc-files:.cc=.go)
+else
+intergrid_map-o-files = $(intergrid_map-cc-files:.cc=.o)
+endif
+
+intergrid_map.testcase: $(intergrid_map-o-files) $(libraries)
+       @echo =====linking======= $<
+       @$(CXX) $(LDFLAGS) -g  -o $@ $^ $(LIBS) $(lib-deal2-1d.g) $(lib-deal2-3d.g)
+
 
 ############################################################
 # Continue with other targets if needed
diff --git a/tests/deal.II/intergrid_map.cc b/tests/deal.II/intergrid_map.cc
new file mode 100644 (file)
index 0000000..7a7125e
--- /dev/null
@@ -0,0 +1,113 @@
+/* $Id$ */
+/* Author: Wolfgang Bangerth, University of Heidelberg, 2000 */
+/* 
+   Purpose: check some things with the intergrid map
+*/
+
+#include <base/logstream.h>
+#include <grid/tria.h>
+#include <dofs/dof_handler.h>
+#include <grid/grid_generator.h>
+#include <grid/tria_accessor.h>
+#include <grid/tria_iterator.h>
+#include <dofs/dof_accessor.h>
+#include <dofs/dof_handler.h>
+#include <grid/intergrid_map.h>
+#include <fe/fe_lib.lagrange.h>
+
+#include <fstream>
+
+
+template <int dim>
+void check ()
+{
+  deallog << "Checking in " << dim << " space dimensions"
+         << endl
+         << "---------------------------------------" << endl;
+  
+                                  // create two grids
+  Triangulation<dim> tria_1, tria_2;
+  GridGenerator::hyper_cube (tria_1, -1, 1);
+  tria_1.refine_global (5-dim);
+  tria_2.copy_triangulation (tria_1);
+
+  FEQ1<dim> fe_1;
+  FEQ2<dim> fe_2;
+  
+                                  // make several loops to refine the
+                                  // two grids
+  for (unsigned int i=0; i<3; ++i)
+    {
+      deallog << "Refinement step " << i << endl;
+      
+      DoFHandler<dim> dof_1 (tria_1);
+      DoFHandler<dim> dof_2 (tria_2);
+
+      dof_1.distribute_dofs (fe_1);
+      dof_2.distribute_dofs (fe_2);
+
+                                      // create some mapping
+      InterGridMap<DoFHandler,dim> intergrid_map_1;
+      InterGridMap<DoFHandler,dim> intergrid_map_2;
+      intergrid_map_1.make_mapping (dof_1, dof_2);
+      intergrid_map_2.make_mapping (dof_2, dof_1);
+
+                                      // write out the mapping
+      typename DoFHandler<dim>::cell_iterator cell=dof_1.begin(),
+                                             endc=dof_1.end();
+      for (; cell!=endc; ++cell)
+       {
+         deallog << cell
+                 << "->"
+                 << intergrid_map_1[cell]
+                 << "->"
+                 << intergrid_map_2[intergrid_map_1[cell]]
+                 << endl;
+// note that not necessarily intergrid_map_2[intergrid_map_1[cell]] ==
+// cell, since the meshes have different refinement steps.
+       };
+      
+
+
+                                      // now refine grids a little,
+                                      // but differently. this
+                                      // produces quite random grids
+      cell = dof_1.begin();
+      for (unsigned int index=0; cell!=endc; ++cell)
+       if (cell->active())
+         {
+           ++index;
+           if (index % 3 == 0)
+             cell->set_refine_flag ();
+         };
+
+      cell = dof_2.begin();
+      endc = dof_2.end();
+      for (unsigned int index=0; cell!=endc; ++cell)
+       if (cell->active())
+         {
+           ++index;
+           if (index % 3 == 1)
+             cell->set_refine_flag ();
+         };
+      
+      tria_1.execute_coarsening_and_refinement ();
+      tria_2.execute_coarsening_and_refinement ();
+    };
+};
+
+
+
+int main ()
+{
+  ofstream logfile("intergrid_map.output");
+  logfile.precision(4);
+  
+  deallog.attach(logfile);
+  deallog.depth_console(0);
+
+  check<1> ();
+  check<2> ();
+  check<3> ();
+};
+
diff --git a/tests/deal.II/intergrid_map.checked b/tests/deal.II/intergrid_map.checked
new file mode 100644 (file)
index 0000000..26ea424
--- /dev/null
@@ -0,0 +1,2157 @@
+
+DEAL::Checking in 1 space dimensions
+DEAL::---------------------------------------
+DEAL::Refinement step 0
+DEAL::0.0->0.0->0.0
+DEAL::1.0->1.0->1.0
+DEAL::1.1->1.1->1.1
+DEAL::2.0->2.0->2.0
+DEAL::2.1->2.1->2.1
+DEAL::2.2->2.2->2.2
+DEAL::2.3->2.3->2.3
+DEAL::3.0->3.0->3.0
+DEAL::3.1->3.1->3.1
+DEAL::3.2->3.2->3.2
+DEAL::3.3->3.3->3.3
+DEAL::3.4->3.4->3.4
+DEAL::3.5->3.5->3.5
+DEAL::3.6->3.6->3.6
+DEAL::3.7->3.7->3.7
+DEAL::4.0->4.0->4.0
+DEAL::4.1->4.1->4.1
+DEAL::4.2->4.2->4.2
+DEAL::4.3->4.3->4.3
+DEAL::4.4->4.4->4.4
+DEAL::4.5->4.5->4.5
+DEAL::4.6->4.6->4.6
+DEAL::4.7->4.7->4.7
+DEAL::4.8->4.8->4.8
+DEAL::4.9->4.9->4.9
+DEAL::4.10->4.10->4.10
+DEAL::4.11->4.11->4.11
+DEAL::4.12->4.12->4.12
+DEAL::4.13->4.13->4.13
+DEAL::4.14->4.14->4.14
+DEAL::4.15->4.15->4.15
+DEAL::Refinement step 1
+DEAL::0.0->0.0->0.0
+DEAL::1.0->1.0->1.0
+DEAL::1.1->1.1->1.1
+DEAL::2.0->2.0->2.0
+DEAL::2.1->2.1->2.1
+DEAL::2.2->2.2->2.2
+DEAL::2.3->2.3->2.3
+DEAL::3.0->3.0->3.0
+DEAL::3.1->3.1->3.1
+DEAL::3.2->3.2->3.2
+DEAL::3.3->3.3->3.3
+DEAL::3.4->3.4->3.4
+DEAL::3.5->3.5->3.5
+DEAL::3.6->3.6->3.6
+DEAL::3.7->3.7->3.7
+DEAL::4.0->4.0->4.0
+DEAL::4.1->4.1->4.1
+DEAL::4.2->4.2->4.2
+DEAL::4.3->4.3->4.3
+DEAL::4.4->4.4->4.4
+DEAL::4.5->4.5->4.5
+DEAL::4.6->4.6->4.6
+DEAL::4.7->4.7->4.7
+DEAL::4.8->4.8->4.8
+DEAL::4.9->4.9->4.9
+DEAL::4.10->4.10->4.10
+DEAL::4.11->4.11->4.11
+DEAL::4.12->4.12->4.12
+DEAL::4.13->4.13->4.13
+DEAL::4.14->4.14->4.14
+DEAL::4.15->4.15->4.15
+DEAL::5.0->4.2->4.2
+DEAL::5.1->4.2->4.2
+DEAL::5.2->4.5->4.5
+DEAL::5.3->4.5->4.5
+DEAL::5.4->4.8->4.8
+DEAL::5.5->4.8->4.8
+DEAL::5.6->4.11->4.11
+DEAL::5.7->4.11->4.11
+DEAL::5.8->4.14->4.14
+DEAL::5.9->4.14->4.14
+DEAL::Refinement step 2
+DEAL::0.0->0.0->0.0
+DEAL::1.0->1.0->1.0
+DEAL::1.1->1.1->1.1
+DEAL::2.0->2.0->2.0
+DEAL::2.1->2.1->2.1
+DEAL::2.2->2.2->2.2
+DEAL::2.3->2.3->2.3
+DEAL::3.0->3.0->3.0
+DEAL::3.1->3.1->3.1
+DEAL::3.2->3.2->3.2
+DEAL::3.3->3.3->3.3
+DEAL::3.4->3.4->3.4
+DEAL::3.5->3.5->3.5
+DEAL::3.6->3.6->3.6
+DEAL::3.7->3.7->3.7
+DEAL::4.0->4.0->4.0
+DEAL::4.1->4.1->4.1
+DEAL::4.2->4.2->4.2
+DEAL::4.3->4.3->4.3
+DEAL::4.4->4.4->4.4
+DEAL::4.5->4.5->4.5
+DEAL::4.6->4.6->4.6
+DEAL::4.7->4.7->4.7
+DEAL::4.8->4.8->4.8
+DEAL::4.9->4.9->4.9
+DEAL::4.10->4.10->4.10
+DEAL::4.11->4.11->4.11
+DEAL::4.12->4.12->4.12
+DEAL::4.13->4.13->4.13
+DEAL::4.14->4.14->4.14
+DEAL::4.15->4.15->4.15
+DEAL::5.0->4.2->4.2
+DEAL::5.1->4.2->4.2
+DEAL::5.2->5.14->5.2
+DEAL::5.3->5.15->5.3
+DEAL::5.4->4.8->4.8
+DEAL::5.5->4.8->4.8
+DEAL::5.6->4.11->4.11
+DEAL::5.7->4.11->4.11
+DEAL::5.8->5.18->5.8
+DEAL::5.9->5.19->5.9
+DEAL::5.10->5.2->5.10
+DEAL::5.11->5.3->5.11
+DEAL::5.12->4.7->4.7
+DEAL::5.13->4.7->4.7
+DEAL::5.14->5.8->5.14
+DEAL::5.15->5.9->5.15
+DEAL::6.0->4.2->4.2
+DEAL::6.1->4.2->4.2
+DEAL::6.2->5.15->5.3
+DEAL::6.3->5.15->5.3
+DEAL::6.4->4.11->4.11
+DEAL::6.5->4.11->4.11
+DEAL::6.6->5.19->5.9
+DEAL::6.7->5.19->5.9
+DEAL::Checking in 2 space dimensions
+DEAL::---------------------------------------
+DEAL::Refinement step 0
+DEAL::0.0->0.0->0.0
+DEAL::1.0->1.0->1.0
+DEAL::1.1->1.1->1.1
+DEAL::1.2->1.2->1.2
+DEAL::1.3->1.3->1.3
+DEAL::2.0->2.0->2.0
+DEAL::2.1->2.1->2.1
+DEAL::2.2->2.2->2.2
+DEAL::2.3->2.3->2.3
+DEAL::2.4->2.4->2.4
+DEAL::2.5->2.5->2.5
+DEAL::2.6->2.6->2.6
+DEAL::2.7->2.7->2.7
+DEAL::2.8->2.8->2.8
+DEAL::2.9->2.9->2.9
+DEAL::2.10->2.10->2.10
+DEAL::2.11->2.11->2.11
+DEAL::2.12->2.12->2.12
+DEAL::2.13->2.13->2.13
+DEAL::2.14->2.14->2.14
+DEAL::2.15->2.15->2.15
+DEAL::3.0->3.0->3.0
+DEAL::3.1->3.1->3.1
+DEAL::3.2->3.2->3.2
+DEAL::3.3->3.3->3.3
+DEAL::3.4->3.4->3.4
+DEAL::3.5->3.5->3.5
+DEAL::3.6->3.6->3.6
+DEAL::3.7->3.7->3.7
+DEAL::3.8->3.8->3.8
+DEAL::3.9->3.9->3.9
+DEAL::3.10->3.10->3.10
+DEAL::3.11->3.11->3.11
+DEAL::3.12->3.12->3.12
+DEAL::3.13->3.13->3.13
+DEAL::3.14->3.14->3.14
+DEAL::3.15->3.15->3.15
+DEAL::3.16->3.16->3.16
+DEAL::3.17->3.17->3.17
+DEAL::3.18->3.18->3.18
+DEAL::3.19->3.19->3.19
+DEAL::3.20->3.20->3.20
+DEAL::3.21->3.21->3.21
+DEAL::3.22->3.22->3.22
+DEAL::3.23->3.23->3.23
+DEAL::3.24->3.24->3.24
+DEAL::3.25->3.25->3.25
+DEAL::3.26->3.26->3.26
+DEAL::3.27->3.27->3.27
+DEAL::3.28->3.28->3.28
+DEAL::3.29->3.29->3.29
+DEAL::3.30->3.30->3.30
+DEAL::3.31->3.31->3.31
+DEAL::3.32->3.32->3.32
+DEAL::3.33->3.33->3.33
+DEAL::3.34->3.34->3.34
+DEAL::3.35->3.35->3.35
+DEAL::3.36->3.36->3.36
+DEAL::3.37->3.37->3.37
+DEAL::3.38->3.38->3.38
+DEAL::3.39->3.39->3.39
+DEAL::3.40->3.40->3.40
+DEAL::3.41->3.41->3.41
+DEAL::3.42->3.42->3.42
+DEAL::3.43->3.43->3.43
+DEAL::3.44->3.44->3.44
+DEAL::3.45->3.45->3.45
+DEAL::3.46->3.46->3.46
+DEAL::3.47->3.47->3.47
+DEAL::3.48->3.48->3.48
+DEAL::3.49->3.49->3.49
+DEAL::3.50->3.50->3.50
+DEAL::3.51->3.51->3.51
+DEAL::3.52->3.52->3.52
+DEAL::3.53->3.53->3.53
+DEAL::3.54->3.54->3.54
+DEAL::3.55->3.55->3.55
+DEAL::3.56->3.56->3.56
+DEAL::3.57->3.57->3.57
+DEAL::3.58->3.58->3.58
+DEAL::3.59->3.59->3.59
+DEAL::3.60->3.60->3.60
+DEAL::3.61->3.61->3.61
+DEAL::3.62->3.62->3.62
+DEAL::3.63->3.63->3.63
+DEAL::Refinement step 1
+DEAL::0.0->0.0->0.0
+DEAL::1.0->1.0->1.0
+DEAL::1.1->1.1->1.1
+DEAL::1.2->1.2->1.2
+DEAL::1.3->1.3->1.3
+DEAL::2.0->2.0->2.0
+DEAL::2.1->2.1->2.1
+DEAL::2.2->2.2->2.2
+DEAL::2.3->2.3->2.3
+DEAL::2.4->2.4->2.4
+DEAL::2.5->2.5->2.5
+DEAL::2.6->2.6->2.6
+DEAL::2.7->2.7->2.7
+DEAL::2.8->2.8->2.8
+DEAL::2.9->2.9->2.9
+DEAL::2.10->2.10->2.10
+DEAL::2.11->2.11->2.11
+DEAL::2.12->2.12->2.12
+DEAL::2.13->2.13->2.13
+DEAL::2.14->2.14->2.14
+DEAL::2.15->2.15->2.15
+DEAL::3.0->3.0->3.0
+DEAL::3.1->3.1->3.1
+DEAL::3.2->3.2->3.2
+DEAL::3.3->3.3->3.3
+DEAL::3.4->3.4->3.4
+DEAL::3.5->3.5->3.5
+DEAL::3.6->3.6->3.6
+DEAL::3.7->3.7->3.7
+DEAL::3.8->3.8->3.8
+DEAL::3.9->3.9->3.9
+DEAL::3.10->3.10->3.10
+DEAL::3.11->3.11->3.11
+DEAL::3.12->3.12->3.12
+DEAL::3.13->3.13->3.13
+DEAL::3.14->3.14->3.14
+DEAL::3.15->3.15->3.15
+DEAL::3.16->3.16->3.16
+DEAL::3.17->3.17->3.17
+DEAL::3.18->3.18->3.18
+DEAL::3.19->3.19->3.19
+DEAL::3.20->3.20->3.20
+DEAL::3.21->3.21->3.21
+DEAL::3.22->3.22->3.22
+DEAL::3.23->3.23->3.23
+DEAL::3.24->3.24->3.24
+DEAL::3.25->3.25->3.25
+DEAL::3.26->3.26->3.26
+DEAL::3.27->3.27->3.27
+DEAL::3.28->3.28->3.28
+DEAL::3.29->3.29->3.29
+DEAL::3.30->3.30->3.30
+DEAL::3.31->3.31->3.31
+DEAL::3.32->3.32->3.32
+DEAL::3.33->3.33->3.33
+DEAL::3.34->3.34->3.34
+DEAL::3.35->3.35->3.35
+DEAL::3.36->3.36->3.36
+DEAL::3.37->3.37->3.37
+DEAL::3.38->3.38->3.38
+DEAL::3.39->3.39->3.39
+DEAL::3.40->3.40->3.40
+DEAL::3.41->3.41->3.41
+DEAL::3.42->3.42->3.42
+DEAL::3.43->3.43->3.43
+DEAL::3.44->3.44->3.44
+DEAL::3.45->3.45->3.45
+DEAL::3.46->3.46->3.46
+DEAL::3.47->3.47->3.47
+DEAL::3.48->3.48->3.48
+DEAL::3.49->3.49->3.49
+DEAL::3.50->3.50->3.50
+DEAL::3.51->3.51->3.51
+DEAL::3.52->3.52->3.52
+DEAL::3.53->3.53->3.53
+DEAL::3.54->3.54->3.54
+DEAL::3.55->3.55->3.55
+DEAL::3.56->3.56->3.56
+DEAL::3.57->3.57->3.57
+DEAL::3.58->3.58->3.58
+DEAL::3.59->3.59->3.59
+DEAL::3.60->3.60->3.60
+DEAL::3.61->3.61->3.61
+DEAL::3.62->3.62->3.62
+DEAL::3.63->3.63->3.63
+DEAL::4.0->3.2->3.2
+DEAL::4.1->3.2->3.2
+DEAL::4.2->3.2->3.2
+DEAL::4.3->3.2->3.2
+DEAL::4.4->3.5->3.5
+DEAL::4.5->3.5->3.5
+DEAL::4.6->3.5->3.5
+DEAL::4.7->3.5->3.5
+DEAL::4.8->3.8->3.8
+DEAL::4.9->3.8->3.8
+DEAL::4.10->3.8->3.8
+DEAL::4.11->3.8->3.8
+DEAL::4.12->3.11->3.11
+DEAL::4.13->3.11->3.11
+DEAL::4.14->3.11->3.11
+DEAL::4.15->3.11->3.11
+DEAL::4.16->3.14->3.14
+DEAL::4.17->3.14->3.14
+DEAL::4.18->3.14->3.14
+DEAL::4.19->3.14->3.14
+DEAL::4.20->3.17->3.17
+DEAL::4.21->3.17->3.17
+DEAL::4.22->3.17->3.17
+DEAL::4.23->3.17->3.17
+DEAL::4.24->3.20->3.20
+DEAL::4.25->3.20->3.20
+DEAL::4.26->3.20->3.20
+DEAL::4.27->3.20->3.20
+DEAL::4.28->3.23->3.23
+DEAL::4.29->3.23->3.23
+DEAL::4.30->3.23->3.23
+DEAL::4.31->3.23->3.23
+DEAL::4.32->3.26->3.26
+DEAL::4.33->3.26->3.26
+DEAL::4.34->3.26->3.26
+DEAL::4.35->3.26->3.26
+DEAL::4.36->3.29->3.29
+DEAL::4.37->3.29->3.29
+DEAL::4.38->3.29->3.29
+DEAL::4.39->3.29->3.29
+DEAL::4.40->3.32->3.32
+DEAL::4.41->3.32->3.32
+DEAL::4.42->3.32->3.32
+DEAL::4.43->3.32->3.32
+DEAL::4.44->3.35->3.35
+DEAL::4.45->3.35->3.35
+DEAL::4.46->3.35->3.35
+DEAL::4.47->3.35->3.35
+DEAL::4.48->3.38->3.38
+DEAL::4.49->3.38->3.38
+DEAL::4.50->3.38->3.38
+DEAL::4.51->3.38->3.38
+DEAL::4.52->3.41->3.41
+DEAL::4.53->3.41->3.41
+DEAL::4.54->3.41->3.41
+DEAL::4.55->3.41->3.41
+DEAL::4.56->3.44->3.44
+DEAL::4.57->3.44->3.44
+DEAL::4.58->3.44->3.44
+DEAL::4.59->3.44->3.44
+DEAL::4.60->3.47->3.47
+DEAL::4.61->3.47->3.47
+DEAL::4.62->3.47->3.47
+DEAL::4.63->3.47->3.47
+DEAL::4.64->3.50->3.50
+DEAL::4.65->3.50->3.50
+DEAL::4.66->3.50->3.50
+DEAL::4.67->3.50->3.50
+DEAL::4.68->3.53->3.53
+DEAL::4.69->3.53->3.53
+DEAL::4.70->3.53->3.53
+DEAL::4.71->3.53->3.53
+DEAL::4.72->3.56->3.56
+DEAL::4.73->3.56->3.56
+DEAL::4.74->3.56->3.56
+DEAL::4.75->3.56->3.56
+DEAL::4.76->3.59->3.59
+DEAL::4.77->3.59->3.59
+DEAL::4.78->3.59->3.59
+DEAL::4.79->3.59->3.59
+DEAL::4.80->3.62->3.62
+DEAL::4.81->3.62->3.62
+DEAL::4.82->3.62->3.62
+DEAL::4.83->3.62->3.62
+DEAL::Refinement step 2
+DEAL::0.0->0.0->0.0
+DEAL::1.0->1.0->1.0
+DEAL::1.1->1.1->1.1
+DEAL::1.2->1.2->1.2
+DEAL::1.3->1.3->1.3
+DEAL::2.0->2.0->2.0
+DEAL::2.1->2.1->2.1
+DEAL::2.2->2.2->2.2
+DEAL::2.3->2.3->2.3
+DEAL::2.4->2.4->2.4
+DEAL::2.5->2.5->2.5
+DEAL::2.6->2.6->2.6
+DEAL::2.7->2.7->2.7
+DEAL::2.8->2.8->2.8
+DEAL::2.9->2.9->2.9
+DEAL::2.10->2.10->2.10
+DEAL::2.11->2.11->2.11
+DEAL::2.12->2.12->2.12
+DEAL::2.13->2.13->2.13
+DEAL::2.14->2.14->2.14
+DEAL::2.15->2.15->2.15
+DEAL::3.0->3.0->3.0
+DEAL::3.1->3.1->3.1
+DEAL::3.2->3.2->3.2
+DEAL::3.3->3.3->3.3
+DEAL::3.4->3.4->3.4
+DEAL::3.5->3.5->3.5
+DEAL::3.6->3.6->3.6
+DEAL::3.7->3.7->3.7
+DEAL::3.8->3.8->3.8
+DEAL::3.9->3.9->3.9
+DEAL::3.10->3.10->3.10
+DEAL::3.11->3.11->3.11
+DEAL::3.12->3.12->3.12
+DEAL::3.13->3.13->3.13
+DEAL::3.14->3.14->3.14
+DEAL::3.15->3.15->3.15
+DEAL::3.16->3.16->3.16
+DEAL::3.17->3.17->3.17
+DEAL::3.18->3.18->3.18
+DEAL::3.19->3.19->3.19
+DEAL::3.20->3.20->3.20
+DEAL::3.21->3.21->3.21
+DEAL::3.22->3.22->3.22
+DEAL::3.23->3.23->3.23
+DEAL::3.24->3.24->3.24
+DEAL::3.25->3.25->3.25
+DEAL::3.26->3.26->3.26
+DEAL::3.27->3.27->3.27
+DEAL::3.28->3.28->3.28
+DEAL::3.29->3.29->3.29
+DEAL::3.30->3.30->3.30
+DEAL::3.31->3.31->3.31
+DEAL::3.32->3.32->3.32
+DEAL::3.33->3.33->3.33
+DEAL::3.34->3.34->3.34
+DEAL::3.35->3.35->3.35
+DEAL::3.36->3.36->3.36
+DEAL::3.37->3.37->3.37
+DEAL::3.38->3.38->3.38
+DEAL::3.39->3.39->3.39
+DEAL::3.40->3.40->3.40
+DEAL::3.41->3.41->3.41
+DEAL::3.42->3.42->3.42
+DEAL::3.43->3.43->3.43
+DEAL::3.44->3.44->3.44
+DEAL::3.45->3.45->3.45
+DEAL::3.46->3.46->3.46
+DEAL::3.47->3.47->3.47
+DEAL::3.48->3.48->3.48
+DEAL::3.49->3.49->3.49
+DEAL::3.50->3.50->3.50
+DEAL::3.51->3.51->3.51
+DEAL::3.52->3.52->3.52
+DEAL::3.53->3.53->3.53
+DEAL::3.54->3.54->3.54
+DEAL::3.55->3.55->3.55
+DEAL::3.56->3.56->3.56
+DEAL::3.57->3.57->3.57
+DEAL::3.58->3.58->3.58
+DEAL::3.59->3.59->3.59
+DEAL::3.60->3.60->3.60
+DEAL::3.61->3.61->3.61
+DEAL::3.62->3.62->3.62
+DEAL::3.63->3.63->3.63
+DEAL::4.0->4.92->4.0
+DEAL::4.1->4.93->4.1
+DEAL::4.2->4.94->4.2
+DEAL::4.3->4.95->4.3
+DEAL::4.4->4.96->4.4
+DEAL::4.5->4.97->4.5
+DEAL::4.6->4.98->4.6
+DEAL::4.7->4.99->4.7
+DEAL::4.8->4.100->4.8
+DEAL::4.9->4.101->4.9
+DEAL::4.10->4.102->4.10
+DEAL::4.11->4.103->4.11
+DEAL::4.12->3.11->3.11
+DEAL::4.13->3.11->3.11
+DEAL::4.14->3.11->3.11
+DEAL::4.15->3.11->3.11
+DEAL::4.16->4.112->4.16
+DEAL::4.17->4.113->4.17
+DEAL::4.18->4.114->4.18
+DEAL::4.19->4.115->4.19
+DEAL::4.20->4.116->4.20
+DEAL::4.21->4.117->4.21
+DEAL::4.22->4.118->4.22
+DEAL::4.23->4.119->4.23
+DEAL::4.24->3.20->3.20
+DEAL::4.25->3.20->3.20
+DEAL::4.26->3.20->3.20
+DEAL::4.27->3.20->3.20
+DEAL::4.28->4.128->4.28
+DEAL::4.29->4.129->4.29
+DEAL::4.30->4.130->4.30
+DEAL::4.31->4.131->4.31
+DEAL::4.32->3.26->3.26
+DEAL::4.33->3.26->3.26
+DEAL::4.34->3.26->3.26
+DEAL::4.35->3.26->3.26
+DEAL::4.36->4.140->4.36
+DEAL::4.37->4.141->4.37
+DEAL::4.38->4.142->4.38
+DEAL::4.39->4.143->4.39
+DEAL::4.40->4.144->4.40
+DEAL::4.41->4.145->4.41
+DEAL::4.42->4.146->4.42
+DEAL::4.43->4.147->4.43
+DEAL::4.44->3.35->3.35
+DEAL::4.45->3.35->3.35
+DEAL::4.46->3.35->3.35
+DEAL::4.47->3.35->3.35
+DEAL::4.48->4.156->4.48
+DEAL::4.49->4.157->4.49
+DEAL::4.50->4.158->4.50
+DEAL::4.51->4.159->4.51
+DEAL::4.52->4.164->4.52
+DEAL::4.53->4.165->4.53
+DEAL::4.54->4.166->4.54
+DEAL::4.55->4.167->4.55
+DEAL::4.56->4.168->4.56
+DEAL::4.57->4.169->4.57
+DEAL::4.58->4.170->4.58
+DEAL::4.59->4.171->4.59
+DEAL::4.60->3.47->3.47
+DEAL::4.61->3.47->3.47
+DEAL::4.62->3.47->3.47
+DEAL::4.63->3.47->3.47
+DEAL::4.64->4.180->4.64
+DEAL::4.65->4.181->4.65
+DEAL::4.66->4.182->4.66
+DEAL::4.67->4.183->4.67
+DEAL::4.68->4.184->4.68
+DEAL::4.69->4.185->4.69
+DEAL::4.70->4.186->4.70
+DEAL::4.71->4.187->4.71
+DEAL::4.72->3.56->3.56
+DEAL::4.73->3.56->3.56
+DEAL::4.74->3.56->3.56
+DEAL::4.75->3.56->3.56
+DEAL::4.76->4.196->4.76
+DEAL::4.77->4.197->4.77
+DEAL::4.78->4.198->4.78
+DEAL::4.79->4.199->4.79
+DEAL::4.80->3.62->3.62
+DEAL::4.81->3.62->3.62
+DEAL::4.82->3.62->3.62
+DEAL::4.83->3.62->3.62
+DEAL::4.84->4.88->4.84
+DEAL::4.85->4.89->4.85
+DEAL::4.86->4.90->4.86
+DEAL::4.87->4.91->4.87
+DEAL::4.88->4.4->4.88
+DEAL::4.89->4.5->4.89
+DEAL::4.90->4.6->4.90
+DEAL::4.91->4.7->4.91
+DEAL::4.92->3.4->3.4
+DEAL::4.93->3.4->3.4
+DEAL::4.94->3.4->3.4
+DEAL::4.95->3.4->3.4
+DEAL::4.96->4.8->4.96
+DEAL::4.97->4.9->4.97
+DEAL::4.98->4.10->4.98
+DEAL::4.99->4.11->4.99
+DEAL::4.100->3.7->3.7
+DEAL::4.101->3.7->3.7
+DEAL::4.102->3.7->3.7
+DEAL::4.103->3.7->3.7
+DEAL::4.104->4.12->4.104
+DEAL::4.105->4.13->4.105
+DEAL::4.106->4.14->4.106
+DEAL::4.107->4.15->4.107
+DEAL::4.108->4.104->4.108
+DEAL::4.109->4.105->4.109
+DEAL::4.110->4.106->4.110
+DEAL::4.111->4.107->4.111
+DEAL::4.112->4.16->4.112
+DEAL::4.113->4.17->4.113
+DEAL::4.114->4.18->4.114
+DEAL::4.115->4.19->4.115
+DEAL::4.116->4.108->4.116
+DEAL::4.117->4.109->4.117
+DEAL::4.118->4.110->4.118
+DEAL::4.119->4.111->4.119
+DEAL::4.120->4.20->4.120
+DEAL::4.121->4.21->4.121
+DEAL::4.122->4.22->4.122
+DEAL::4.123->4.23->4.123
+DEAL::4.124->3.16->3.16
+DEAL::4.125->3.16->3.16
+DEAL::4.126->3.16->3.16
+DEAL::4.127->3.16->3.16
+DEAL::4.128->4.24->4.128
+DEAL::4.129->4.25->4.129
+DEAL::4.130->4.26->4.130
+DEAL::4.131->4.27->4.131
+DEAL::4.132->4.28->4.132
+DEAL::4.133->4.29->4.133
+DEAL::4.134->4.30->4.134
+DEAL::4.135->4.31->4.135
+DEAL::4.136->4.32->4.136
+DEAL::4.137->4.33->4.137
+DEAL::4.138->4.34->4.138
+DEAL::4.139->4.35->4.139
+DEAL::4.140->4.132->4.140
+DEAL::4.141->4.133->4.141
+DEAL::4.142->4.134->4.142
+DEAL::4.143->4.135->4.143
+DEAL::4.144->4.40->4.144
+DEAL::4.145->4.41->4.145
+DEAL::4.146->4.42->4.146
+DEAL::4.147->4.43->4.147
+DEAL::4.148->3.31->3.31
+DEAL::4.149->3.31->3.31
+DEAL::4.150->3.31->3.31
+DEAL::4.151->3.31->3.31
+DEAL::4.152->4.148->4.152
+DEAL::4.153->4.149->4.153
+DEAL::4.154->4.150->4.154
+DEAL::4.155->4.151->4.155
+DEAL::4.156->4.152->4.156
+DEAL::4.157->4.153->4.157
+DEAL::4.158->4.154->4.158
+DEAL::4.159->4.155->4.159
+DEAL::4.160->4.52->4.160
+DEAL::4.161->4.53->4.161
+DEAL::4.162->4.54->4.162
+DEAL::4.163->4.55->4.163
+DEAL::4.164->4.160->4.164
+DEAL::4.165->4.161->4.165
+DEAL::4.166->4.162->4.166
+DEAL::4.167->4.163->4.167
+DEAL::4.168->4.56->4.168
+DEAL::4.169->4.57->4.169
+DEAL::4.170->4.58->4.170
+DEAL::4.171->4.59->4.171
+DEAL::4.172->3.43->3.43
+DEAL::4.173->3.43->3.43
+DEAL::4.174->3.43->3.43
+DEAL::4.175->3.43->3.43
+DEAL::4.176->4.60->4.176
+DEAL::4.177->4.61->4.177
+DEAL::4.178->4.62->4.178
+DEAL::4.179->4.63->4.179
+DEAL::4.180->4.172->4.180
+DEAL::4.181->4.173->4.181
+DEAL::4.182->4.174->4.182
+DEAL::4.183->4.175->4.183
+DEAL::4.184->4.64->4.184
+DEAL::4.185->4.65->4.185
+DEAL::4.186->4.66->4.186
+DEAL::4.187->4.67->4.187
+DEAL::4.188->4.176->4.188
+DEAL::4.189->4.177->4.189
+DEAL::4.190->4.178->4.190
+DEAL::4.191->4.179->4.191
+DEAL::4.192->4.68->4.192
+DEAL::4.193->4.69->4.193
+DEAL::4.194->4.70->4.194
+DEAL::4.195->4.71->4.195
+DEAL::4.196->3.52->3.52
+DEAL::4.197->3.52->3.52
+DEAL::4.198->3.52->3.52
+DEAL::4.199->3.52->3.52
+DEAL::4.200->4.72->4.200
+DEAL::4.201->4.73->4.201
+DEAL::4.202->4.74->4.202
+DEAL::4.203->4.75->4.203
+DEAL::4.204->4.188->4.204
+DEAL::4.205->4.189->4.205
+DEAL::4.206->4.190->4.206
+DEAL::4.207->4.191->4.207
+DEAL::4.208->4.76->4.208
+DEAL::4.209->4.77->4.209
+DEAL::4.210->4.78->4.210
+DEAL::4.211->4.79->4.211
+DEAL::4.212->4.200->4.212
+DEAL::4.213->4.201->4.213
+DEAL::4.214->4.202->4.214
+DEAL::4.215->4.203->4.215
+DEAL::5.0->4.93->4.1
+DEAL::5.1->4.93->4.1
+DEAL::5.2->4.93->4.1
+DEAL::5.3->4.93->4.1
+DEAL::5.4->4.96->4.4
+DEAL::5.5->4.96->4.4
+DEAL::5.6->4.96->4.4
+DEAL::5.7->4.96->4.4
+DEAL::5.8->4.99->4.7
+DEAL::5.9->4.99->4.7
+DEAL::5.10->4.99->4.7
+DEAL::5.11->4.99->4.7
+DEAL::5.12->4.102->4.10
+DEAL::5.13->4.102->4.10
+DEAL::5.14->4.102->4.10
+DEAL::5.15->4.102->4.10
+DEAL::5.16->3.11->3.11
+DEAL::5.17->3.11->3.11
+DEAL::5.18->3.11->3.11
+DEAL::5.19->3.11->3.11
+DEAL::5.20->4.112->4.16
+DEAL::5.21->4.112->4.16
+DEAL::5.22->4.112->4.16
+DEAL::5.23->4.112->4.16
+DEAL::5.24->4.115->4.19
+DEAL::5.25->4.115->4.19
+DEAL::5.26->4.115->4.19
+DEAL::5.27->4.115->4.19
+DEAL::5.28->4.118->4.22
+DEAL::5.29->4.118->4.22
+DEAL::5.30->4.118->4.22
+DEAL::5.31->4.118->4.22
+DEAL::5.32->3.20->3.20
+DEAL::5.33->3.20->3.20
+DEAL::5.34->3.20->3.20
+DEAL::5.35->3.20->3.20
+DEAL::5.36->4.128->4.28
+DEAL::5.37->4.128->4.28
+DEAL::5.38->4.128->4.28
+DEAL::5.39->4.128->4.28
+DEAL::5.40->4.131->4.31
+DEAL::5.41->4.131->4.31
+DEAL::5.42->4.131->4.31
+DEAL::5.43->4.131->4.31
+DEAL::5.44->3.26->3.26
+DEAL::5.45->3.26->3.26
+DEAL::5.46->3.26->3.26
+DEAL::5.47->3.26->3.26
+DEAL::5.48->4.141->4.37
+DEAL::5.49->4.141->4.37
+DEAL::5.50->4.141->4.37
+DEAL::5.51->4.141->4.37
+DEAL::5.52->4.144->4.40
+DEAL::5.53->4.144->4.40
+DEAL::5.54->4.144->4.40
+DEAL::5.55->4.144->4.40
+DEAL::5.56->4.147->4.43
+DEAL::5.57->4.147->4.43
+DEAL::5.58->4.147->4.43
+DEAL::5.59->4.147->4.43
+DEAL::5.60->3.35->3.35
+DEAL::5.61->3.35->3.35
+DEAL::5.62->3.35->3.35
+DEAL::5.63->3.35->3.35
+DEAL::5.64->4.157->4.49
+DEAL::5.65->4.157->4.49
+DEAL::5.66->4.157->4.49
+DEAL::5.67->4.157->4.49
+DEAL::5.68->4.164->4.52
+DEAL::5.69->4.164->4.52
+DEAL::5.70->4.164->4.52
+DEAL::5.71->4.164->4.52
+DEAL::5.72->4.167->4.55
+DEAL::5.73->4.167->4.55
+DEAL::5.74->4.167->4.55
+DEAL::5.75->4.167->4.55
+DEAL::5.76->4.170->4.58
+DEAL::5.77->4.170->4.58
+DEAL::5.78->4.170->4.58
+DEAL::5.79->4.170->4.58
+DEAL::5.80->3.47->3.47
+DEAL::5.81->3.47->3.47
+DEAL::5.82->3.47->3.47
+DEAL::5.83->3.47->3.47
+DEAL::5.84->4.180->4.64
+DEAL::5.85->4.180->4.64
+DEAL::5.86->4.180->4.64
+DEAL::5.87->4.180->4.64
+DEAL::5.88->4.183->4.67
+DEAL::5.89->4.183->4.67
+DEAL::5.90->4.183->4.67
+DEAL::5.91->4.183->4.67
+DEAL::5.92->4.186->4.70
+DEAL::5.93->4.186->4.70
+DEAL::5.94->4.186->4.70
+DEAL::5.95->4.186->4.70
+DEAL::5.96->3.56->3.56
+DEAL::5.97->3.56->3.56
+DEAL::5.98->3.56->3.56
+DEAL::5.99->3.56->3.56
+DEAL::5.100->4.196->4.76
+DEAL::5.101->4.196->4.76
+DEAL::5.102->4.196->4.76
+DEAL::5.103->4.196->4.76
+DEAL::5.104->4.199->4.79
+DEAL::5.105->4.199->4.79
+DEAL::5.106->4.199->4.79
+DEAL::5.107->4.199->4.79
+DEAL::5.108->3.62->3.62
+DEAL::5.109->3.62->3.62
+DEAL::5.110->3.62->3.62
+DEAL::5.111->3.62->3.62
+DEAL::Checking in 3 space dimensions
+DEAL::---------------------------------------
+DEAL::Refinement step 0
+DEAL::0.0->0.0->0.0
+DEAL::1.0->1.0->1.0
+DEAL::1.1->1.1->1.1
+DEAL::1.2->1.2->1.2
+DEAL::1.3->1.3->1.3
+DEAL::1.4->1.4->1.4
+DEAL::1.5->1.5->1.5
+DEAL::1.6->1.6->1.6
+DEAL::1.7->1.7->1.7
+DEAL::2.0->2.0->2.0
+DEAL::2.1->2.1->2.1
+DEAL::2.2->2.2->2.2
+DEAL::2.3->2.3->2.3
+DEAL::2.4->2.4->2.4
+DEAL::2.5->2.5->2.5
+DEAL::2.6->2.6->2.6
+DEAL::2.7->2.7->2.7
+DEAL::2.8->2.8->2.8
+DEAL::2.9->2.9->2.9
+DEAL::2.10->2.10->2.10
+DEAL::2.11->2.11->2.11
+DEAL::2.12->2.12->2.12
+DEAL::2.13->2.13->2.13
+DEAL::2.14->2.14->2.14
+DEAL::2.15->2.15->2.15
+DEAL::2.16->2.16->2.16
+DEAL::2.17->2.17->2.17
+DEAL::2.18->2.18->2.18
+DEAL::2.19->2.19->2.19
+DEAL::2.20->2.20->2.20
+DEAL::2.21->2.21->2.21
+DEAL::2.22->2.22->2.22
+DEAL::2.23->2.23->2.23
+DEAL::2.24->2.24->2.24
+DEAL::2.25->2.25->2.25
+DEAL::2.26->2.26->2.26
+DEAL::2.27->2.27->2.27
+DEAL::2.28->2.28->2.28
+DEAL::2.29->2.29->2.29
+DEAL::2.30->2.30->2.30
+DEAL::2.31->2.31->2.31
+DEAL::2.32->2.32->2.32
+DEAL::2.33->2.33->2.33
+DEAL::2.34->2.34->2.34
+DEAL::2.35->2.35->2.35
+DEAL::2.36->2.36->2.36
+DEAL::2.37->2.37->2.37
+DEAL::2.38->2.38->2.38
+DEAL::2.39->2.39->2.39
+DEAL::2.40->2.40->2.40
+DEAL::2.41->2.41->2.41
+DEAL::2.42->2.42->2.42
+DEAL::2.43->2.43->2.43
+DEAL::2.44->2.44->2.44
+DEAL::2.45->2.45->2.45
+DEAL::2.46->2.46->2.46
+DEAL::2.47->2.47->2.47
+DEAL::2.48->2.48->2.48
+DEAL::2.49->2.49->2.49
+DEAL::2.50->2.50->2.50
+DEAL::2.51->2.51->2.51
+DEAL::2.52->2.52->2.52
+DEAL::2.53->2.53->2.53
+DEAL::2.54->2.54->2.54
+DEAL::2.55->2.55->2.55
+DEAL::2.56->2.56->2.56
+DEAL::2.57->2.57->2.57
+DEAL::2.58->2.58->2.58
+DEAL::2.59->2.59->2.59
+DEAL::2.60->2.60->2.60
+DEAL::2.61->2.61->2.61
+DEAL::2.62->2.62->2.62
+DEAL::2.63->2.63->2.63
+DEAL::Refinement step 1
+DEAL::0.0->0.0->0.0
+DEAL::1.0->1.0->1.0
+DEAL::1.1->1.1->1.1
+DEAL::1.2->1.2->1.2
+DEAL::1.3->1.3->1.3
+DEAL::1.4->1.4->1.4
+DEAL::1.5->1.5->1.5
+DEAL::1.6->1.6->1.6
+DEAL::1.7->1.7->1.7
+DEAL::2.0->2.0->2.0
+DEAL::2.1->2.1->2.1
+DEAL::2.2->2.2->2.2
+DEAL::2.3->2.3->2.3
+DEAL::2.4->2.4->2.4
+DEAL::2.5->2.5->2.5
+DEAL::2.6->2.6->2.6
+DEAL::2.7->2.7->2.7
+DEAL::2.8->2.8->2.8
+DEAL::2.9->2.9->2.9
+DEAL::2.10->2.10->2.10
+DEAL::2.11->2.11->2.11
+DEAL::2.12->2.12->2.12
+DEAL::2.13->2.13->2.13
+DEAL::2.14->2.14->2.14
+DEAL::2.15->2.15->2.15
+DEAL::2.16->2.16->2.16
+DEAL::2.17->2.17->2.17
+DEAL::2.18->2.18->2.18
+DEAL::2.19->2.19->2.19
+DEAL::2.20->2.20->2.20
+DEAL::2.21->2.21->2.21
+DEAL::2.22->2.22->2.22
+DEAL::2.23->2.23->2.23
+DEAL::2.24->2.24->2.24
+DEAL::2.25->2.25->2.25
+DEAL::2.26->2.26->2.26
+DEAL::2.27->2.27->2.27
+DEAL::2.28->2.28->2.28
+DEAL::2.29->2.29->2.29
+DEAL::2.30->2.30->2.30
+DEAL::2.31->2.31->2.31
+DEAL::2.32->2.32->2.32
+DEAL::2.33->2.33->2.33
+DEAL::2.34->2.34->2.34
+DEAL::2.35->2.35->2.35
+DEAL::2.36->2.36->2.36
+DEAL::2.37->2.37->2.37
+DEAL::2.38->2.38->2.38
+DEAL::2.39->2.39->2.39
+DEAL::2.40->2.40->2.40
+DEAL::2.41->2.41->2.41
+DEAL::2.42->2.42->2.42
+DEAL::2.43->2.43->2.43
+DEAL::2.44->2.44->2.44
+DEAL::2.45->2.45->2.45
+DEAL::2.46->2.46->2.46
+DEAL::2.47->2.47->2.47
+DEAL::2.48->2.48->2.48
+DEAL::2.49->2.49->2.49
+DEAL::2.50->2.50->2.50
+DEAL::2.51->2.51->2.51
+DEAL::2.52->2.52->2.52
+DEAL::2.53->2.53->2.53
+DEAL::2.54->2.54->2.54
+DEAL::2.55->2.55->2.55
+DEAL::2.56->2.56->2.56
+DEAL::2.57->2.57->2.57
+DEAL::2.58->2.58->2.58
+DEAL::2.59->2.59->2.59
+DEAL::2.60->2.60->2.60
+DEAL::2.61->2.61->2.61
+DEAL::2.62->2.62->2.62
+DEAL::2.63->2.63->2.63
+DEAL::3.0->2.2->2.2
+DEAL::3.1->2.2->2.2
+DEAL::3.2->2.2->2.2
+DEAL::3.3->2.2->2.2
+DEAL::3.4->2.2->2.2
+DEAL::3.5->2.2->2.2
+DEAL::3.6->2.2->2.2
+DEAL::3.7->2.2->2.2
+DEAL::3.8->2.5->2.5
+DEAL::3.9->2.5->2.5
+DEAL::3.10->2.5->2.5
+DEAL::3.11->2.5->2.5
+DEAL::3.12->2.5->2.5
+DEAL::3.13->2.5->2.5
+DEAL::3.14->2.5->2.5
+DEAL::3.15->2.5->2.5
+DEAL::3.16->2.8->2.8
+DEAL::3.17->2.8->2.8
+DEAL::3.18->2.8->2.8
+DEAL::3.19->2.8->2.8
+DEAL::3.20->2.8->2.8
+DEAL::3.21->2.8->2.8
+DEAL::3.22->2.8->2.8
+DEAL::3.23->2.8->2.8
+DEAL::3.24->2.11->2.11
+DEAL::3.25->2.11->2.11
+DEAL::3.26->2.11->2.11
+DEAL::3.27->2.11->2.11
+DEAL::3.28->2.11->2.11
+DEAL::3.29->2.11->2.11
+DEAL::3.30->2.11->2.11
+DEAL::3.31->2.11->2.11
+DEAL::3.32->2.14->2.14
+DEAL::3.33->2.14->2.14
+DEAL::3.34->2.14->2.14
+DEAL::3.35->2.14->2.14
+DEAL::3.36->2.14->2.14
+DEAL::3.37->2.14->2.14
+DEAL::3.38->2.14->2.14
+DEAL::3.39->2.14->2.14
+DEAL::3.40->2.17->2.17
+DEAL::3.41->2.17->2.17
+DEAL::3.42->2.17->2.17
+DEAL::3.43->2.17->2.17
+DEAL::3.44->2.17->2.17
+DEAL::3.45->2.17->2.17
+DEAL::3.46->2.17->2.17
+DEAL::3.47->2.17->2.17
+DEAL::3.48->2.20->2.20
+DEAL::3.49->2.20->2.20
+DEAL::3.50->2.20->2.20
+DEAL::3.51->2.20->2.20
+DEAL::3.52->2.20->2.20
+DEAL::3.53->2.20->2.20
+DEAL::3.54->2.20->2.20
+DEAL::3.55->2.20->2.20
+DEAL::3.56->2.23->2.23
+DEAL::3.57->2.23->2.23
+DEAL::3.58->2.23->2.23
+DEAL::3.59->2.23->2.23
+DEAL::3.60->2.23->2.23
+DEAL::3.61->2.23->2.23
+DEAL::3.62->2.23->2.23
+DEAL::3.63->2.23->2.23
+DEAL::3.64->2.26->2.26
+DEAL::3.65->2.26->2.26
+DEAL::3.66->2.26->2.26
+DEAL::3.67->2.26->2.26
+DEAL::3.68->2.26->2.26
+DEAL::3.69->2.26->2.26
+DEAL::3.70->2.26->2.26
+DEAL::3.71->2.26->2.26
+DEAL::3.72->2.29->2.29
+DEAL::3.73->2.29->2.29
+DEAL::3.74->2.29->2.29
+DEAL::3.75->2.29->2.29
+DEAL::3.76->2.29->2.29
+DEAL::3.77->2.29->2.29
+DEAL::3.78->2.29->2.29
+DEAL::3.79->2.29->2.29
+DEAL::3.80->2.32->2.32
+DEAL::3.81->2.32->2.32
+DEAL::3.82->2.32->2.32
+DEAL::3.83->2.32->2.32
+DEAL::3.84->2.32->2.32
+DEAL::3.85->2.32->2.32
+DEAL::3.86->2.32->2.32
+DEAL::3.87->2.32->2.32
+DEAL::3.88->2.35->2.35
+DEAL::3.89->2.35->2.35
+DEAL::3.90->2.35->2.35
+DEAL::3.91->2.35->2.35
+DEAL::3.92->2.35->2.35
+DEAL::3.93->2.35->2.35
+DEAL::3.94->2.35->2.35
+DEAL::3.95->2.35->2.35
+DEAL::3.96->2.38->2.38
+DEAL::3.97->2.38->2.38
+DEAL::3.98->2.38->2.38
+DEAL::3.99->2.38->2.38
+DEAL::3.100->2.38->2.38
+DEAL::3.101->2.38->2.38
+DEAL::3.102->2.38->2.38
+DEAL::3.103->2.38->2.38
+DEAL::3.104->2.41->2.41
+DEAL::3.105->2.41->2.41
+DEAL::3.106->2.41->2.41
+DEAL::3.107->2.41->2.41
+DEAL::3.108->2.41->2.41
+DEAL::3.109->2.41->2.41
+DEAL::3.110->2.41->2.41
+DEAL::3.111->2.41->2.41
+DEAL::3.112->2.44->2.44
+DEAL::3.113->2.44->2.44
+DEAL::3.114->2.44->2.44
+DEAL::3.115->2.44->2.44
+DEAL::3.116->2.44->2.44
+DEAL::3.117->2.44->2.44
+DEAL::3.118->2.44->2.44
+DEAL::3.119->2.44->2.44
+DEAL::3.120->2.47->2.47
+DEAL::3.121->2.47->2.47
+DEAL::3.122->2.47->2.47
+DEAL::3.123->2.47->2.47
+DEAL::3.124->2.47->2.47
+DEAL::3.125->2.47->2.47
+DEAL::3.126->2.47->2.47
+DEAL::3.127->2.47->2.47
+DEAL::3.128->2.50->2.50
+DEAL::3.129->2.50->2.50
+DEAL::3.130->2.50->2.50
+DEAL::3.131->2.50->2.50
+DEAL::3.132->2.50->2.50
+DEAL::3.133->2.50->2.50
+DEAL::3.134->2.50->2.50
+DEAL::3.135->2.50->2.50
+DEAL::3.136->2.53->2.53
+DEAL::3.137->2.53->2.53
+DEAL::3.138->2.53->2.53
+DEAL::3.139->2.53->2.53
+DEAL::3.140->2.53->2.53
+DEAL::3.141->2.53->2.53
+DEAL::3.142->2.53->2.53
+DEAL::3.143->2.53->2.53
+DEAL::3.144->2.56->2.56
+DEAL::3.145->2.56->2.56
+DEAL::3.146->2.56->2.56
+DEAL::3.147->2.56->2.56
+DEAL::3.148->2.56->2.56
+DEAL::3.149->2.56->2.56
+DEAL::3.150->2.56->2.56
+DEAL::3.151->2.56->2.56
+DEAL::3.152->2.59->2.59
+DEAL::3.153->2.59->2.59
+DEAL::3.154->2.59->2.59
+DEAL::3.155->2.59->2.59
+DEAL::3.156->2.59->2.59
+DEAL::3.157->2.59->2.59
+DEAL::3.158->2.59->2.59
+DEAL::3.159->2.59->2.59
+DEAL::3.160->2.62->2.62
+DEAL::3.161->2.62->2.62
+DEAL::3.162->2.62->2.62
+DEAL::3.163->2.62->2.62
+DEAL::3.164->2.62->2.62
+DEAL::3.165->2.62->2.62
+DEAL::3.166->2.62->2.62
+DEAL::3.167->2.62->2.62
+DEAL::Refinement step 2
+DEAL::0.0->0.0->0.0
+DEAL::1.0->1.0->1.0
+DEAL::1.1->1.1->1.1
+DEAL::1.2->1.2->1.2
+DEAL::1.3->1.3->1.3
+DEAL::1.4->1.4->1.4
+DEAL::1.5->1.5->1.5
+DEAL::1.6->1.6->1.6
+DEAL::1.7->1.7->1.7
+DEAL::2.0->2.0->2.0
+DEAL::2.1->2.1->2.1
+DEAL::2.2->2.2->2.2
+DEAL::2.3->2.3->2.3
+DEAL::2.4->2.4->2.4
+DEAL::2.5->2.5->2.5
+DEAL::2.6->2.6->2.6
+DEAL::2.7->2.7->2.7
+DEAL::2.8->2.8->2.8
+DEAL::2.9->2.9->2.9
+DEAL::2.10->2.10->2.10
+DEAL::2.11->2.11->2.11
+DEAL::2.12->2.12->2.12
+DEAL::2.13->2.13->2.13
+DEAL::2.14->2.14->2.14
+DEAL::2.15->2.15->2.15
+DEAL::2.16->2.16->2.16
+DEAL::2.17->2.17->2.17
+DEAL::2.18->2.18->2.18
+DEAL::2.19->2.19->2.19
+DEAL::2.20->2.20->2.20
+DEAL::2.21->2.21->2.21
+DEAL::2.22->2.22->2.22
+DEAL::2.23->2.23->2.23
+DEAL::2.24->2.24->2.24
+DEAL::2.25->2.25->2.25
+DEAL::2.26->2.26->2.26
+DEAL::2.27->2.27->2.27
+DEAL::2.28->2.28->2.28
+DEAL::2.29->2.29->2.29
+DEAL::2.30->2.30->2.30
+DEAL::2.31->2.31->2.31
+DEAL::2.32->2.32->2.32
+DEAL::2.33->2.33->2.33
+DEAL::2.34->2.34->2.34
+DEAL::2.35->2.35->2.35
+DEAL::2.36->2.36->2.36
+DEAL::2.37->2.37->2.37
+DEAL::2.38->2.38->2.38
+DEAL::2.39->2.39->2.39
+DEAL::2.40->2.40->2.40
+DEAL::2.41->2.41->2.41
+DEAL::2.42->2.42->2.42
+DEAL::2.43->2.43->2.43
+DEAL::2.44->2.44->2.44
+DEAL::2.45->2.45->2.45
+DEAL::2.46->2.46->2.46
+DEAL::2.47->2.47->2.47
+DEAL::2.48->2.48->2.48
+DEAL::2.49->2.49->2.49
+DEAL::2.50->2.50->2.50
+DEAL::2.51->2.51->2.51
+DEAL::2.52->2.52->2.52
+DEAL::2.53->2.53->2.53
+DEAL::2.54->2.54->2.54
+DEAL::2.55->2.55->2.55
+DEAL::2.56->2.56->2.56
+DEAL::2.57->2.57->2.57
+DEAL::2.58->2.58->2.58
+DEAL::2.59->2.59->2.59
+DEAL::2.60->2.60->2.60
+DEAL::2.61->2.61->2.61
+DEAL::2.62->2.62->2.62
+DEAL::2.63->2.63->2.63
+DEAL::3.0->3.184->3.0
+DEAL::3.1->3.185->3.1
+DEAL::3.2->3.186->3.2
+DEAL::3.3->3.187->3.3
+DEAL::3.4->3.188->3.4
+DEAL::3.5->3.189->3.5
+DEAL::3.6->3.190->3.6
+DEAL::3.7->3.191->3.7
+DEAL::3.8->3.200->3.8
+DEAL::3.9->3.201->3.9
+DEAL::3.10->3.202->3.10
+DEAL::3.11->3.203->3.11
+DEAL::3.12->3.204->3.12
+DEAL::3.13->3.205->3.13
+DEAL::3.14->3.206->3.14
+DEAL::3.15->3.207->3.15
+DEAL::3.16->3.216->3.16
+DEAL::3.17->3.217->3.17
+DEAL::3.18->3.218->3.18
+DEAL::3.19->3.219->3.19
+DEAL::3.20->3.220->3.20
+DEAL::3.21->3.221->3.21
+DEAL::3.22->3.222->3.22
+DEAL::3.23->3.223->3.23
+DEAL::3.24->3.232->3.24
+DEAL::3.25->3.233->3.25
+DEAL::3.26->3.234->3.26
+DEAL::3.27->3.235->3.27
+DEAL::3.28->3.236->3.28
+DEAL::3.29->3.237->3.29
+DEAL::3.30->3.238->3.30
+DEAL::3.31->3.239->3.31
+DEAL::3.32->3.248->3.32
+DEAL::3.33->3.249->3.33
+DEAL::3.34->3.250->3.34
+DEAL::3.35->3.251->3.35
+DEAL::3.36->3.252->3.36
+DEAL::3.37->3.253->3.37
+DEAL::3.38->3.254->3.38
+DEAL::3.39->3.255->3.39
+DEAL::3.40->3.264->3.40
+DEAL::3.41->3.265->3.41
+DEAL::3.42->3.266->3.42
+DEAL::3.43->3.267->3.43
+DEAL::3.44->3.268->3.44
+DEAL::3.45->3.269->3.45
+DEAL::3.46->3.270->3.46
+DEAL::3.47->3.271->3.47
+DEAL::3.48->3.280->3.48
+DEAL::3.49->3.281->3.49
+DEAL::3.50->3.282->3.50
+DEAL::3.51->3.283->3.51
+DEAL::3.52->3.284->3.52
+DEAL::3.53->3.285->3.53
+DEAL::3.54->3.286->3.54
+DEAL::3.55->3.287->3.55
+DEAL::3.56->3.296->3.56
+DEAL::3.57->3.297->3.57
+DEAL::3.58->3.298->3.58
+DEAL::3.59->3.299->3.59
+DEAL::3.60->3.300->3.60
+DEAL::3.61->3.301->3.61
+DEAL::3.62->3.302->3.62
+DEAL::3.63->3.303->3.63
+DEAL::3.64->3.312->3.64
+DEAL::3.65->3.313->3.65
+DEAL::3.66->3.314->3.66
+DEAL::3.67->3.315->3.67
+DEAL::3.68->3.316->3.68
+DEAL::3.69->3.317->3.69
+DEAL::3.70->3.318->3.70
+DEAL::3.71->3.319->3.71
+DEAL::3.72->3.328->3.72
+DEAL::3.73->3.329->3.73
+DEAL::3.74->3.330->3.74
+DEAL::3.75->3.331->3.75
+DEAL::3.76->3.332->3.76
+DEAL::3.77->3.333->3.77
+DEAL::3.78->3.334->3.78
+DEAL::3.79->3.335->3.79
+DEAL::3.80->3.344->3.80
+DEAL::3.81->3.345->3.81
+DEAL::3.82->3.346->3.82
+DEAL::3.83->3.347->3.83
+DEAL::3.84->3.348->3.84
+DEAL::3.85->3.349->3.85
+DEAL::3.86->3.350->3.86
+DEAL::3.87->3.351->3.87
+DEAL::3.88->3.360->3.88
+DEAL::3.89->3.361->3.89
+DEAL::3.90->3.362->3.90
+DEAL::3.91->3.363->3.91
+DEAL::3.92->3.364->3.92
+DEAL::3.93->3.365->3.93
+DEAL::3.94->3.366->3.94
+DEAL::3.95->3.367->3.95
+DEAL::3.96->3.376->3.96
+DEAL::3.97->3.377->3.97
+DEAL::3.98->3.378->3.98
+DEAL::3.99->3.379->3.99
+DEAL::3.100->3.380->3.100
+DEAL::3.101->3.381->3.101
+DEAL::3.102->3.382->3.102
+DEAL::3.103->3.383->3.103
+DEAL::3.104->3.392->3.104
+DEAL::3.105->3.393->3.105
+DEAL::3.106->3.394->3.106
+DEAL::3.107->3.395->3.107
+DEAL::3.108->3.396->3.108
+DEAL::3.109->3.397->3.109
+DEAL::3.110->3.398->3.110
+DEAL::3.111->3.399->3.111
+DEAL::3.112->3.408->3.112
+DEAL::3.113->3.409->3.113
+DEAL::3.114->3.410->3.114
+DEAL::3.115->3.411->3.115
+DEAL::3.116->3.412->3.116
+DEAL::3.117->3.413->3.117
+DEAL::3.118->3.414->3.118
+DEAL::3.119->3.415->3.119
+DEAL::3.120->3.424->3.120
+DEAL::3.121->3.425->3.121
+DEAL::3.122->3.426->3.122
+DEAL::3.123->3.427->3.123
+DEAL::3.124->3.428->3.124
+DEAL::3.125->3.429->3.125
+DEAL::3.126->3.430->3.126
+DEAL::3.127->3.431->3.127
+DEAL::3.128->3.440->3.128
+DEAL::3.129->3.441->3.129
+DEAL::3.130->3.442->3.130
+DEAL::3.131->3.443->3.131
+DEAL::3.132->3.444->3.132
+DEAL::3.133->3.445->3.133
+DEAL::3.134->3.446->3.134
+DEAL::3.135->3.447->3.135
+DEAL::3.136->3.456->3.136
+DEAL::3.137->3.457->3.137
+DEAL::3.138->3.458->3.138
+DEAL::3.139->3.459->3.139
+DEAL::3.140->3.460->3.140
+DEAL::3.141->3.461->3.141
+DEAL::3.142->3.462->3.142
+DEAL::3.143->3.463->3.143
+DEAL::3.144->3.472->3.144
+DEAL::3.145->3.473->3.145
+DEAL::3.146->3.474->3.146
+DEAL::3.147->3.475->3.147
+DEAL::3.148->3.476->3.148
+DEAL::3.149->3.477->3.149
+DEAL::3.150->3.478->3.150
+DEAL::3.151->3.479->3.151
+DEAL::3.152->3.488->3.152
+DEAL::3.153->3.489->3.153
+DEAL::3.154->3.490->3.154
+DEAL::3.155->3.491->3.155
+DEAL::3.156->3.492->3.156
+DEAL::3.157->3.493->3.157
+DEAL::3.158->3.494->3.158
+DEAL::3.159->3.495->3.159
+DEAL::3.160->3.504->3.160
+DEAL::3.161->3.505->3.161
+DEAL::3.162->3.506->3.162
+DEAL::3.163->3.507->3.163
+DEAL::3.164->3.508->3.164
+DEAL::3.165->3.509->3.165
+DEAL::3.166->3.510->3.166
+DEAL::3.167->3.511->3.167
+DEAL::3.168->3.0->3.168
+DEAL::3.169->3.1->3.169
+DEAL::3.170->3.2->3.170
+DEAL::3.171->3.3->3.171
+DEAL::3.172->3.4->3.172
+DEAL::3.173->3.5->3.173
+DEAL::3.174->3.6->3.174
+DEAL::3.175->3.7->3.175
+DEAL::3.176->3.176->3.176
+DEAL::3.177->3.177->3.177
+DEAL::3.178->3.178->3.178
+DEAL::3.179->3.179->3.179
+DEAL::3.180->3.180->3.180
+DEAL::3.181->3.181->3.181
+DEAL::3.182->3.182->3.182
+DEAL::3.183->3.183->3.183
+DEAL::3.184->3.8->3.184
+DEAL::3.185->3.9->3.185
+DEAL::3.186->3.10->3.186
+DEAL::3.187->3.11->3.187
+DEAL::3.188->3.12->3.188
+DEAL::3.189->3.13->3.189
+DEAL::3.190->3.14->3.190
+DEAL::3.191->3.15->3.191
+DEAL::3.192->3.192->3.192
+DEAL::3.193->3.193->3.193
+DEAL::3.194->3.194->3.194
+DEAL::3.195->3.195->3.195
+DEAL::3.196->3.196->3.196
+DEAL::3.197->3.197->3.197
+DEAL::3.198->3.198->3.198
+DEAL::3.199->3.199->3.199
+DEAL::3.200->3.16->3.200
+DEAL::3.201->3.17->3.201
+DEAL::3.202->3.18->3.202
+DEAL::3.203->3.19->3.203
+DEAL::3.204->3.20->3.204
+DEAL::3.205->3.21->3.205
+DEAL::3.206->3.22->3.206
+DEAL::3.207->3.23->3.207
+DEAL::3.208->3.208->3.208
+DEAL::3.209->3.209->3.209
+DEAL::3.210->3.210->3.210
+DEAL::3.211->3.211->3.211
+DEAL::3.212->3.212->3.212
+DEAL::3.213->3.213->3.213
+DEAL::3.214->3.214->3.214
+DEAL::3.215->3.215->3.215
+DEAL::3.216->3.24->3.216
+DEAL::3.217->3.25->3.217
+DEAL::3.218->3.26->3.218
+DEAL::3.219->3.27->3.219
+DEAL::3.220->3.28->3.220
+DEAL::3.221->3.29->3.221
+DEAL::3.222->3.30->3.222
+DEAL::3.223->3.31->3.223
+DEAL::3.224->3.224->3.224
+DEAL::3.225->3.225->3.225
+DEAL::3.226->3.226->3.226
+DEAL::3.227->3.227->3.227
+DEAL::3.228->3.228->3.228
+DEAL::3.229->3.229->3.229
+DEAL::3.230->3.230->3.230
+DEAL::3.231->3.231->3.231
+DEAL::3.232->3.32->3.232
+DEAL::3.233->3.33->3.233
+DEAL::3.234->3.34->3.234
+DEAL::3.235->3.35->3.235
+DEAL::3.236->3.36->3.236
+DEAL::3.237->3.37->3.237
+DEAL::3.238->3.38->3.238
+DEAL::3.239->3.39->3.239
+DEAL::3.240->3.240->3.240
+DEAL::3.241->3.241->3.241
+DEAL::3.242->3.242->3.242
+DEAL::3.243->3.243->3.243
+DEAL::3.244->3.244->3.244
+DEAL::3.245->3.245->3.245
+DEAL::3.246->3.246->3.246
+DEAL::3.247->3.247->3.247
+DEAL::3.248->3.40->3.248
+DEAL::3.249->3.41->3.249
+DEAL::3.250->3.42->3.250
+DEAL::3.251->3.43->3.251
+DEAL::3.252->3.44->3.252
+DEAL::3.253->3.45->3.253
+DEAL::3.254->3.46->3.254
+DEAL::3.255->3.47->3.255
+DEAL::3.256->3.256->3.256
+DEAL::3.257->3.257->3.257
+DEAL::3.258->3.258->3.258
+DEAL::3.259->3.259->3.259
+DEAL::3.260->3.260->3.260
+DEAL::3.261->3.261->3.261
+DEAL::3.262->3.262->3.262
+DEAL::3.263->3.263->3.263
+DEAL::3.264->3.48->3.264
+DEAL::3.265->3.49->3.265
+DEAL::3.266->3.50->3.266
+DEAL::3.267->3.51->3.267
+DEAL::3.268->3.52->3.268
+DEAL::3.269->3.53->3.269
+DEAL::3.270->3.54->3.270
+DEAL::3.271->3.55->3.271
+DEAL::3.272->3.272->3.272
+DEAL::3.273->3.273->3.273
+DEAL::3.274->3.274->3.274
+DEAL::3.275->3.275->3.275
+DEAL::3.276->3.276->3.276
+DEAL::3.277->3.277->3.277
+DEAL::3.278->3.278->3.278
+DEAL::3.279->3.279->3.279
+DEAL::3.280->3.56->3.280
+DEAL::3.281->3.57->3.281
+DEAL::3.282->3.58->3.282
+DEAL::3.283->3.59->3.283
+DEAL::3.284->3.60->3.284
+DEAL::3.285->3.61->3.285
+DEAL::3.286->3.62->3.286
+DEAL::3.287->3.63->3.287
+DEAL::3.288->3.288->3.288
+DEAL::3.289->3.289->3.289
+DEAL::3.290->3.290->3.290
+DEAL::3.291->3.291->3.291
+DEAL::3.292->3.292->3.292
+DEAL::3.293->3.293->3.293
+DEAL::3.294->3.294->3.294
+DEAL::3.295->3.295->3.295
+DEAL::3.296->3.64->3.296
+DEAL::3.297->3.65->3.297
+DEAL::3.298->3.66->3.298
+DEAL::3.299->3.67->3.299
+DEAL::3.300->3.68->3.300
+DEAL::3.301->3.69->3.301
+DEAL::3.302->3.70->3.302
+DEAL::3.303->3.71->3.303
+DEAL::3.304->3.304->3.304
+DEAL::3.305->3.305->3.305
+DEAL::3.306->3.306->3.306
+DEAL::3.307->3.307->3.307
+DEAL::3.308->3.308->3.308
+DEAL::3.309->3.309->3.309
+DEAL::3.310->3.310->3.310
+DEAL::3.311->3.311->3.311
+DEAL::3.312->3.72->3.312
+DEAL::3.313->3.73->3.313
+DEAL::3.314->3.74->3.314
+DEAL::3.315->3.75->3.315
+DEAL::3.316->3.76->3.316
+DEAL::3.317->3.77->3.317
+DEAL::3.318->3.78->3.318
+DEAL::3.319->3.79->3.319
+DEAL::3.320->3.320->3.320
+DEAL::3.321->3.321->3.321
+DEAL::3.322->3.322->3.322
+DEAL::3.323->3.323->3.323
+DEAL::3.324->3.324->3.324
+DEAL::3.325->3.325->3.325
+DEAL::3.326->3.326->3.326
+DEAL::3.327->3.327->3.327
+DEAL::3.328->3.80->3.328
+DEAL::3.329->3.81->3.329
+DEAL::3.330->3.82->3.330
+DEAL::3.331->3.83->3.331
+DEAL::3.332->3.84->3.332
+DEAL::3.333->3.85->3.333
+DEAL::3.334->3.86->3.334
+DEAL::3.335->3.87->3.335
+DEAL::3.336->3.336->3.336
+DEAL::3.337->3.337->3.337
+DEAL::3.338->3.338->3.338
+DEAL::3.339->3.339->3.339
+DEAL::3.340->3.340->3.340
+DEAL::3.341->3.341->3.341
+DEAL::3.342->3.342->3.342
+DEAL::3.343->3.343->3.343
+DEAL::3.344->3.88->3.344
+DEAL::3.345->3.89->3.345
+DEAL::3.346->3.90->3.346
+DEAL::3.347->3.91->3.347
+DEAL::3.348->3.92->3.348
+DEAL::3.349->3.93->3.349
+DEAL::3.350->3.94->3.350
+DEAL::3.351->3.95->3.351
+DEAL::3.352->3.352->3.352
+DEAL::3.353->3.353->3.353
+DEAL::3.354->3.354->3.354
+DEAL::3.355->3.355->3.355
+DEAL::3.356->3.356->3.356
+DEAL::3.357->3.357->3.357
+DEAL::3.358->3.358->3.358
+DEAL::3.359->3.359->3.359
+DEAL::3.360->3.96->3.360
+DEAL::3.361->3.97->3.361
+DEAL::3.362->3.98->3.362
+DEAL::3.363->3.99->3.363
+DEAL::3.364->3.100->3.364
+DEAL::3.365->3.101->3.365
+DEAL::3.366->3.102->3.366
+DEAL::3.367->3.103->3.367
+DEAL::3.368->3.368->3.368
+DEAL::3.369->3.369->3.369
+DEAL::3.370->3.370->3.370
+DEAL::3.371->3.371->3.371
+DEAL::3.372->3.372->3.372
+DEAL::3.373->3.373->3.373
+DEAL::3.374->3.374->3.374
+DEAL::3.375->3.375->3.375
+DEAL::3.376->3.104->3.376
+DEAL::3.377->3.105->3.377
+DEAL::3.378->3.106->3.378
+DEAL::3.379->3.107->3.379
+DEAL::3.380->3.108->3.380
+DEAL::3.381->3.109->3.381
+DEAL::3.382->3.110->3.382
+DEAL::3.383->3.111->3.383
+DEAL::3.384->3.384->3.384
+DEAL::3.385->3.385->3.385
+DEAL::3.386->3.386->3.386
+DEAL::3.387->3.387->3.387
+DEAL::3.388->3.388->3.388
+DEAL::3.389->3.389->3.389
+DEAL::3.390->3.390->3.390
+DEAL::3.391->3.391->3.391
+DEAL::3.392->3.112->3.392
+DEAL::3.393->3.113->3.393
+DEAL::3.394->3.114->3.394
+DEAL::3.395->3.115->3.395
+DEAL::3.396->3.116->3.396
+DEAL::3.397->3.117->3.397
+DEAL::3.398->3.118->3.398
+DEAL::3.399->3.119->3.399
+DEAL::3.400->3.400->3.400
+DEAL::3.401->3.401->3.401
+DEAL::3.402->3.402->3.402
+DEAL::3.403->3.403->3.403
+DEAL::3.404->3.404->3.404
+DEAL::3.405->3.405->3.405
+DEAL::3.406->3.406->3.406
+DEAL::3.407->3.407->3.407
+DEAL::3.408->3.120->3.408
+DEAL::3.409->3.121->3.409
+DEAL::3.410->3.122->3.410
+DEAL::3.411->3.123->3.411
+DEAL::3.412->3.124->3.412
+DEAL::3.413->3.125->3.413
+DEAL::3.414->3.126->3.414
+DEAL::3.415->3.127->3.415
+DEAL::3.416->3.416->3.416
+DEAL::3.417->3.417->3.417
+DEAL::3.418->3.418->3.418
+DEAL::3.419->3.419->3.419
+DEAL::3.420->3.420->3.420
+DEAL::3.421->3.421->3.421
+DEAL::3.422->3.422->3.422
+DEAL::3.423->3.423->3.423
+DEAL::3.424->3.128->3.424
+DEAL::3.425->3.129->3.425
+DEAL::3.426->3.130->3.426
+DEAL::3.427->3.131->3.427
+DEAL::3.428->3.132->3.428
+DEAL::3.429->3.133->3.429
+DEAL::3.430->3.134->3.430
+DEAL::3.431->3.135->3.431
+DEAL::3.432->3.432->3.432
+DEAL::3.433->3.433->3.433
+DEAL::3.434->3.434->3.434
+DEAL::3.435->3.435->3.435
+DEAL::3.436->3.436->3.436
+DEAL::3.437->3.437->3.437
+DEAL::3.438->3.438->3.438
+DEAL::3.439->3.439->3.439
+DEAL::3.440->3.136->3.440
+DEAL::3.441->3.137->3.441
+DEAL::3.442->3.138->3.442
+DEAL::3.443->3.139->3.443
+DEAL::3.444->3.140->3.444
+DEAL::3.445->3.141->3.445
+DEAL::3.446->3.142->3.446
+DEAL::3.447->3.143->3.447
+DEAL::3.448->3.448->3.448
+DEAL::3.449->3.449->3.449
+DEAL::3.450->3.450->3.450
+DEAL::3.451->3.451->3.451
+DEAL::3.452->3.452->3.452
+DEAL::3.453->3.453->3.453
+DEAL::3.454->3.454->3.454
+DEAL::3.455->3.455->3.455
+DEAL::3.456->3.144->3.456
+DEAL::3.457->3.145->3.457
+DEAL::3.458->3.146->3.458
+DEAL::3.459->3.147->3.459
+DEAL::3.460->3.148->3.460
+DEAL::3.461->3.149->3.461
+DEAL::3.462->3.150->3.462
+DEAL::3.463->3.151->3.463
+DEAL::3.464->3.464->3.464
+DEAL::3.465->3.465->3.465
+DEAL::3.466->3.466->3.466
+DEAL::3.467->3.467->3.467
+DEAL::3.468->3.468->3.468
+DEAL::3.469->3.469->3.469
+DEAL::3.470->3.470->3.470
+DEAL::3.471->3.471->3.471
+DEAL::3.472->3.152->3.472
+DEAL::3.473->3.153->3.473
+DEAL::3.474->3.154->3.474
+DEAL::3.475->3.155->3.475
+DEAL::3.476->3.156->3.476
+DEAL::3.477->3.157->3.477
+DEAL::3.478->3.158->3.478
+DEAL::3.479->3.159->3.479
+DEAL::3.480->3.480->3.480
+DEAL::3.481->3.481->3.481
+DEAL::3.482->3.482->3.482
+DEAL::3.483->3.483->3.483
+DEAL::3.484->3.484->3.484
+DEAL::3.485->3.485->3.485
+DEAL::3.486->3.486->3.486
+DEAL::3.487->3.487->3.487
+DEAL::3.488->3.160->3.488
+DEAL::3.489->3.161->3.489
+DEAL::3.490->3.162->3.490
+DEAL::3.491->3.163->3.491
+DEAL::3.492->3.164->3.492
+DEAL::3.493->3.165->3.493
+DEAL::3.494->3.166->3.494
+DEAL::3.495->3.167->3.495
+DEAL::3.496->3.496->3.496
+DEAL::3.497->3.497->3.497
+DEAL::3.498->3.498->3.498
+DEAL::3.499->3.499->3.499
+DEAL::3.500->3.500->3.500
+DEAL::3.501->3.501->3.501
+DEAL::3.502->3.502->3.502
+DEAL::3.503->3.503->3.503
+DEAL::3.504->3.168->3.504
+DEAL::3.505->3.169->3.505
+DEAL::3.506->3.170->3.506
+DEAL::3.507->3.171->3.507
+DEAL::3.508->3.172->3.508
+DEAL::3.509->3.173->3.509
+DEAL::3.510->3.174->3.510
+DEAL::3.511->3.175->3.511
+DEAL::4.0->3.185->3.1
+DEAL::4.1->3.185->3.1
+DEAL::4.2->3.185->3.1
+DEAL::4.3->3.185->3.1
+DEAL::4.4->3.185->3.1
+DEAL::4.5->3.185->3.1
+DEAL::4.6->3.185->3.1
+DEAL::4.7->3.185->3.1
+DEAL::4.8->3.188->3.4
+DEAL::4.9->3.188->3.4
+DEAL::4.10->3.188->3.4
+DEAL::4.11->3.188->3.4
+DEAL::4.12->3.188->3.4
+DEAL::4.13->3.188->3.4
+DEAL::4.14->3.188->3.4
+DEAL::4.15->3.188->3.4
+DEAL::4.16->3.191->3.7
+DEAL::4.17->3.191->3.7
+DEAL::4.18->3.191->3.7
+DEAL::4.19->3.191->3.7
+DEAL::4.20->3.191->3.7
+DEAL::4.21->3.191->3.7
+DEAL::4.22->3.191->3.7
+DEAL::4.23->3.191->3.7
+DEAL::4.24->3.202->3.10
+DEAL::4.25->3.202->3.10
+DEAL::4.26->3.202->3.10
+DEAL::4.27->3.202->3.10
+DEAL::4.28->3.202->3.10
+DEAL::4.29->3.202->3.10
+DEAL::4.30->3.202->3.10
+DEAL::4.31->3.202->3.10
+DEAL::4.32->3.205->3.13
+DEAL::4.33->3.205->3.13
+DEAL::4.34->3.205->3.13
+DEAL::4.35->3.205->3.13
+DEAL::4.36->3.205->3.13
+DEAL::4.37->3.205->3.13
+DEAL::4.38->3.205->3.13
+DEAL::4.39->3.205->3.13
+DEAL::4.40->3.216->3.16
+DEAL::4.41->3.216->3.16
+DEAL::4.42->3.216->3.16
+DEAL::4.43->3.216->3.16
+DEAL::4.44->3.216->3.16
+DEAL::4.45->3.216->3.16
+DEAL::4.46->3.216->3.16
+DEAL::4.47->3.216->3.16
+DEAL::4.48->3.219->3.19
+DEAL::4.49->3.219->3.19
+DEAL::4.50->3.219->3.19
+DEAL::4.51->3.219->3.19
+DEAL::4.52->3.219->3.19
+DEAL::4.53->3.219->3.19
+DEAL::4.54->3.219->3.19
+DEAL::4.55->3.219->3.19
+DEAL::4.56->3.222->3.22
+DEAL::4.57->3.222->3.22
+DEAL::4.58->3.222->3.22
+DEAL::4.59->3.222->3.22
+DEAL::4.60->3.222->3.22
+DEAL::4.61->3.222->3.22
+DEAL::4.62->3.222->3.22
+DEAL::4.63->3.222->3.22
+DEAL::4.64->3.233->3.25
+DEAL::4.65->3.233->3.25
+DEAL::4.66->3.233->3.25
+DEAL::4.67->3.233->3.25
+DEAL::4.68->3.233->3.25
+DEAL::4.69->3.233->3.25
+DEAL::4.70->3.233->3.25
+DEAL::4.71->3.233->3.25
+DEAL::4.72->3.236->3.28
+DEAL::4.73->3.236->3.28
+DEAL::4.74->3.236->3.28
+DEAL::4.75->3.236->3.28
+DEAL::4.76->3.236->3.28
+DEAL::4.77->3.236->3.28
+DEAL::4.78->3.236->3.28
+DEAL::4.79->3.236->3.28
+DEAL::4.80->3.239->3.31
+DEAL::4.81->3.239->3.31
+DEAL::4.82->3.239->3.31
+DEAL::4.83->3.239->3.31
+DEAL::4.84->3.239->3.31
+DEAL::4.85->3.239->3.31
+DEAL::4.86->3.239->3.31
+DEAL::4.87->3.239->3.31
+DEAL::4.88->3.250->3.34
+DEAL::4.89->3.250->3.34
+DEAL::4.90->3.250->3.34
+DEAL::4.91->3.250->3.34
+DEAL::4.92->3.250->3.34
+DEAL::4.93->3.250->3.34
+DEAL::4.94->3.250->3.34
+DEAL::4.95->3.250->3.34
+DEAL::4.96->3.253->3.37
+DEAL::4.97->3.253->3.37
+DEAL::4.98->3.253->3.37
+DEAL::4.99->3.253->3.37
+DEAL::4.100->3.253->3.37
+DEAL::4.101->3.253->3.37
+DEAL::4.102->3.253->3.37
+DEAL::4.103->3.253->3.37
+DEAL::4.104->3.264->3.40
+DEAL::4.105->3.264->3.40
+DEAL::4.106->3.264->3.40
+DEAL::4.107->3.264->3.40
+DEAL::4.108->3.264->3.40
+DEAL::4.109->3.264->3.40
+DEAL::4.110->3.264->3.40
+DEAL::4.111->3.264->3.40
+DEAL::4.112->3.267->3.43
+DEAL::4.113->3.267->3.43
+DEAL::4.114->3.267->3.43
+DEAL::4.115->3.267->3.43
+DEAL::4.116->3.267->3.43
+DEAL::4.117->3.267->3.43
+DEAL::4.118->3.267->3.43
+DEAL::4.119->3.267->3.43
+DEAL::4.120->3.270->3.46
+DEAL::4.121->3.270->3.46
+DEAL::4.122->3.270->3.46
+DEAL::4.123->3.270->3.46
+DEAL::4.124->3.270->3.46
+DEAL::4.125->3.270->3.46
+DEAL::4.126->3.270->3.46
+DEAL::4.127->3.270->3.46
+DEAL::4.128->3.281->3.49
+DEAL::4.129->3.281->3.49
+DEAL::4.130->3.281->3.49
+DEAL::4.131->3.281->3.49
+DEAL::4.132->3.281->3.49
+DEAL::4.133->3.281->3.49
+DEAL::4.134->3.281->3.49
+DEAL::4.135->3.281->3.49
+DEAL::4.136->3.284->3.52
+DEAL::4.137->3.284->3.52
+DEAL::4.138->3.284->3.52
+DEAL::4.139->3.284->3.52
+DEAL::4.140->3.284->3.52
+DEAL::4.141->3.284->3.52
+DEAL::4.142->3.284->3.52
+DEAL::4.143->3.284->3.52
+DEAL::4.144->3.287->3.55
+DEAL::4.145->3.287->3.55
+DEAL::4.146->3.287->3.55
+DEAL::4.147->3.287->3.55
+DEAL::4.148->3.287->3.55
+DEAL::4.149->3.287->3.55
+DEAL::4.150->3.287->3.55
+DEAL::4.151->3.287->3.55
+DEAL::4.152->3.298->3.58
+DEAL::4.153->3.298->3.58
+DEAL::4.154->3.298->3.58
+DEAL::4.155->3.298->3.58
+DEAL::4.156->3.298->3.58
+DEAL::4.157->3.298->3.58
+DEAL::4.158->3.298->3.58
+DEAL::4.159->3.298->3.58
+DEAL::4.160->3.301->3.61
+DEAL::4.161->3.301->3.61
+DEAL::4.162->3.301->3.61
+DEAL::4.163->3.301->3.61
+DEAL::4.164->3.301->3.61
+DEAL::4.165->3.301->3.61
+DEAL::4.166->3.301->3.61
+DEAL::4.167->3.301->3.61
+DEAL::4.168->3.312->3.64
+DEAL::4.169->3.312->3.64
+DEAL::4.170->3.312->3.64
+DEAL::4.171->3.312->3.64
+DEAL::4.172->3.312->3.64
+DEAL::4.173->3.312->3.64
+DEAL::4.174->3.312->3.64
+DEAL::4.175->3.312->3.64
+DEAL::4.176->3.315->3.67
+DEAL::4.177->3.315->3.67
+DEAL::4.178->3.315->3.67
+DEAL::4.179->3.315->3.67
+DEAL::4.180->3.315->3.67
+DEAL::4.181->3.315->3.67
+DEAL::4.182->3.315->3.67
+DEAL::4.183->3.315->3.67
+DEAL::4.184->3.318->3.70
+DEAL::4.185->3.318->3.70
+DEAL::4.186->3.318->3.70
+DEAL::4.187->3.318->3.70
+DEAL::4.188->3.318->3.70
+DEAL::4.189->3.318->3.70
+DEAL::4.190->3.318->3.70
+DEAL::4.191->3.318->3.70
+DEAL::4.192->3.329->3.73
+DEAL::4.193->3.329->3.73
+DEAL::4.194->3.329->3.73
+DEAL::4.195->3.329->3.73
+DEAL::4.196->3.329->3.73
+DEAL::4.197->3.329->3.73
+DEAL::4.198->3.329->3.73
+DEAL::4.199->3.329->3.73
+DEAL::4.200->3.332->3.76
+DEAL::4.201->3.332->3.76
+DEAL::4.202->3.332->3.76
+DEAL::4.203->3.332->3.76
+DEAL::4.204->3.332->3.76
+DEAL::4.205->3.332->3.76
+DEAL::4.206->3.332->3.76
+DEAL::4.207->3.332->3.76
+DEAL::4.208->3.335->3.79
+DEAL::4.209->3.335->3.79
+DEAL::4.210->3.335->3.79
+DEAL::4.211->3.335->3.79
+DEAL::4.212->3.335->3.79
+DEAL::4.213->3.335->3.79
+DEAL::4.214->3.335->3.79
+DEAL::4.215->3.335->3.79
+DEAL::4.216->3.346->3.82
+DEAL::4.217->3.346->3.82
+DEAL::4.218->3.346->3.82
+DEAL::4.219->3.346->3.82
+DEAL::4.220->3.346->3.82
+DEAL::4.221->3.346->3.82
+DEAL::4.222->3.346->3.82
+DEAL::4.223->3.346->3.82
+DEAL::4.224->3.349->3.85
+DEAL::4.225->3.349->3.85
+DEAL::4.226->3.349->3.85
+DEAL::4.227->3.349->3.85
+DEAL::4.228->3.349->3.85
+DEAL::4.229->3.349->3.85
+DEAL::4.230->3.349->3.85
+DEAL::4.231->3.349->3.85
+DEAL::4.232->3.360->3.88
+DEAL::4.233->3.360->3.88
+DEAL::4.234->3.360->3.88
+DEAL::4.235->3.360->3.88
+DEAL::4.236->3.360->3.88
+DEAL::4.237->3.360->3.88
+DEAL::4.238->3.360->3.88
+DEAL::4.239->3.360->3.88
+DEAL::4.240->3.363->3.91
+DEAL::4.241->3.363->3.91
+DEAL::4.242->3.363->3.91
+DEAL::4.243->3.363->3.91
+DEAL::4.244->3.363->3.91
+DEAL::4.245->3.363->3.91
+DEAL::4.246->3.363->3.91
+DEAL::4.247->3.363->3.91
+DEAL::4.248->3.366->3.94
+DEAL::4.249->3.366->3.94
+DEAL::4.250->3.366->3.94
+DEAL::4.251->3.366->3.94
+DEAL::4.252->3.366->3.94
+DEAL::4.253->3.366->3.94
+DEAL::4.254->3.366->3.94
+DEAL::4.255->3.366->3.94
+DEAL::4.256->3.377->3.97
+DEAL::4.257->3.377->3.97
+DEAL::4.258->3.377->3.97
+DEAL::4.259->3.377->3.97
+DEAL::4.260->3.377->3.97
+DEAL::4.261->3.377->3.97
+DEAL::4.262->3.377->3.97
+DEAL::4.263->3.377->3.97
+DEAL::4.264->3.380->3.100
+DEAL::4.265->3.380->3.100
+DEAL::4.266->3.380->3.100
+DEAL::4.267->3.380->3.100
+DEAL::4.268->3.380->3.100
+DEAL::4.269->3.380->3.100
+DEAL::4.270->3.380->3.100
+DEAL::4.271->3.380->3.100
+DEAL::4.272->3.383->3.103
+DEAL::4.273->3.383->3.103
+DEAL::4.274->3.383->3.103
+DEAL::4.275->3.383->3.103
+DEAL::4.276->3.383->3.103
+DEAL::4.277->3.383->3.103
+DEAL::4.278->3.383->3.103
+DEAL::4.279->3.383->3.103
+DEAL::4.280->3.394->3.106
+DEAL::4.281->3.394->3.106
+DEAL::4.282->3.394->3.106
+DEAL::4.283->3.394->3.106
+DEAL::4.284->3.394->3.106
+DEAL::4.285->3.394->3.106
+DEAL::4.286->3.394->3.106
+DEAL::4.287->3.394->3.106
+DEAL::4.288->3.397->3.109
+DEAL::4.289->3.397->3.109
+DEAL::4.290->3.397->3.109
+DEAL::4.291->3.397->3.109
+DEAL::4.292->3.397->3.109
+DEAL::4.293->3.397->3.109
+DEAL::4.294->3.397->3.109
+DEAL::4.295->3.397->3.109
+DEAL::4.296->3.408->3.112
+DEAL::4.297->3.408->3.112
+DEAL::4.298->3.408->3.112
+DEAL::4.299->3.408->3.112
+DEAL::4.300->3.408->3.112
+DEAL::4.301->3.408->3.112
+DEAL::4.302->3.408->3.112
+DEAL::4.303->3.408->3.112
+DEAL::4.304->3.411->3.115
+DEAL::4.305->3.411->3.115
+DEAL::4.306->3.411->3.115
+DEAL::4.307->3.411->3.115
+DEAL::4.308->3.411->3.115
+DEAL::4.309->3.411->3.115
+DEAL::4.310->3.411->3.115
+DEAL::4.311->3.411->3.115
+DEAL::4.312->3.414->3.118
+DEAL::4.313->3.414->3.118
+DEAL::4.314->3.414->3.118
+DEAL::4.315->3.414->3.118
+DEAL::4.316->3.414->3.118
+DEAL::4.317->3.414->3.118
+DEAL::4.318->3.414->3.118
+DEAL::4.319->3.414->3.118
+DEAL::4.320->3.425->3.121
+DEAL::4.321->3.425->3.121
+DEAL::4.322->3.425->3.121
+DEAL::4.323->3.425->3.121
+DEAL::4.324->3.425->3.121
+DEAL::4.325->3.425->3.121
+DEAL::4.326->3.425->3.121
+DEAL::4.327->3.425->3.121
+DEAL::4.328->3.428->3.124
+DEAL::4.329->3.428->3.124
+DEAL::4.330->3.428->3.124
+DEAL::4.331->3.428->3.124
+DEAL::4.332->3.428->3.124
+DEAL::4.333->3.428->3.124
+DEAL::4.334->3.428->3.124
+DEAL::4.335->3.428->3.124
+DEAL::4.336->3.431->3.127
+DEAL::4.337->3.431->3.127
+DEAL::4.338->3.431->3.127
+DEAL::4.339->3.431->3.127
+DEAL::4.340->3.431->3.127
+DEAL::4.341->3.431->3.127
+DEAL::4.342->3.431->3.127
+DEAL::4.343->3.431->3.127
+DEAL::4.344->3.442->3.130
+DEAL::4.345->3.442->3.130
+DEAL::4.346->3.442->3.130
+DEAL::4.347->3.442->3.130
+DEAL::4.348->3.442->3.130
+DEAL::4.349->3.442->3.130
+DEAL::4.350->3.442->3.130
+DEAL::4.351->3.442->3.130
+DEAL::4.352->3.445->3.133
+DEAL::4.353->3.445->3.133
+DEAL::4.354->3.445->3.133
+DEAL::4.355->3.445->3.133
+DEAL::4.356->3.445->3.133
+DEAL::4.357->3.445->3.133
+DEAL::4.358->3.445->3.133
+DEAL::4.359->3.445->3.133
+DEAL::4.360->3.456->3.136
+DEAL::4.361->3.456->3.136
+DEAL::4.362->3.456->3.136
+DEAL::4.363->3.456->3.136
+DEAL::4.364->3.456->3.136
+DEAL::4.365->3.456->3.136
+DEAL::4.366->3.456->3.136
+DEAL::4.367->3.456->3.136
+DEAL::4.368->3.459->3.139
+DEAL::4.369->3.459->3.139
+DEAL::4.370->3.459->3.139
+DEAL::4.371->3.459->3.139
+DEAL::4.372->3.459->3.139
+DEAL::4.373->3.459->3.139
+DEAL::4.374->3.459->3.139
+DEAL::4.375->3.459->3.139
+DEAL::4.376->3.462->3.142
+DEAL::4.377->3.462->3.142
+DEAL::4.378->3.462->3.142
+DEAL::4.379->3.462->3.142
+DEAL::4.380->3.462->3.142
+DEAL::4.381->3.462->3.142
+DEAL::4.382->3.462->3.142
+DEAL::4.383->3.462->3.142
+DEAL::4.384->3.473->3.145
+DEAL::4.385->3.473->3.145
+DEAL::4.386->3.473->3.145
+DEAL::4.387->3.473->3.145
+DEAL::4.388->3.473->3.145
+DEAL::4.389->3.473->3.145
+DEAL::4.390->3.473->3.145
+DEAL::4.391->3.473->3.145
+DEAL::4.392->3.476->3.148
+DEAL::4.393->3.476->3.148
+DEAL::4.394->3.476->3.148
+DEAL::4.395->3.476->3.148
+DEAL::4.396->3.476->3.148
+DEAL::4.397->3.476->3.148
+DEAL::4.398->3.476->3.148
+DEAL::4.399->3.476->3.148
+DEAL::4.400->3.479->3.151
+DEAL::4.401->3.479->3.151
+DEAL::4.402->3.479->3.151
+DEAL::4.403->3.479->3.151
+DEAL::4.404->3.479->3.151
+DEAL::4.405->3.479->3.151
+DEAL::4.406->3.479->3.151
+DEAL::4.407->3.479->3.151
+DEAL::4.408->3.490->3.154
+DEAL::4.409->3.490->3.154
+DEAL::4.410->3.490->3.154
+DEAL::4.411->3.490->3.154
+DEAL::4.412->3.490->3.154
+DEAL::4.413->3.490->3.154
+DEAL::4.414->3.490->3.154
+DEAL::4.415->3.490->3.154
+DEAL::4.416->3.493->3.157
+DEAL::4.417->3.493->3.157
+DEAL::4.418->3.493->3.157
+DEAL::4.419->3.493->3.157
+DEAL::4.420->3.493->3.157
+DEAL::4.421->3.493->3.157
+DEAL::4.422->3.493->3.157
+DEAL::4.423->3.493->3.157
+DEAL::4.424->3.504->3.160
+DEAL::4.425->3.504->3.160
+DEAL::4.426->3.504->3.160
+DEAL::4.427->3.504->3.160
+DEAL::4.428->3.504->3.160
+DEAL::4.429->3.504->3.160
+DEAL::4.430->3.504->3.160
+DEAL::4.431->3.504->3.160
+DEAL::4.432->3.507->3.163
+DEAL::4.433->3.507->3.163
+DEAL::4.434->3.507->3.163
+DEAL::4.435->3.507->3.163
+DEAL::4.436->3.507->3.163
+DEAL::4.437->3.507->3.163
+DEAL::4.438->3.507->3.163
+DEAL::4.439->3.507->3.163
+DEAL::4.440->3.510->3.166
+DEAL::4.441->3.510->3.166
+DEAL::4.442->3.510->3.166
+DEAL::4.443->3.510->3.166
+DEAL::4.444->3.510->3.166
+DEAL::4.445->3.510->3.166
+DEAL::4.446->3.510->3.166
+DEAL::4.447->3.510->3.166
index 90fff768cac3eedb622c7877859d9d8277cff56d..d750477fdd934cc32fe49a6e7624725f346bf816 100644 (file)
@@ -23,6 +23,7 @@
 #include <lac/precondition.h>
 #include <grid/tria.h>
 #include <grid/tria_iterator.h>
+#include <grid/grid_out.h>
 #include <multigrid/mg_dof_handler.h>
 #include <dofs/dof_constraints.h>
 #include <dofs/dof_accessor.h>
@@ -73,7 +74,7 @@ int main()
   logfile.precision (3);
   deallog.attach(logfile);
 //  deallog.log_execution_time(true);
-  deallog.depth_console(0);
+  deallog.depth_console(10);
 
   Helmholtz equation;
   RHSFunction<2> rhs;
@@ -83,9 +84,9 @@ int main()
   FEQ2<2> fe2;
   FEQ3<2> fe3;
   FEQ4<2> fe4;
-  for (unsigned int degree=1;degree<=3;degree++)
+  for (unsigned int degree=1;degree<=1;degree++)
     {
-      Triangulation<2> tr;
+      Triangulation<2> tr(Triangulation<2>::maximum_smoothing);
       MGDoFHandler<2> mgdof(tr);
       DoFHandler<2>& dof(mgdof);
   
@@ -101,11 +102,38 @@ int main()
        }
 
       tr.refine_global(1);
-      Triangulation<2>::active_cell_iterator cell = tr.begin_active();
-      cell->set_refine_flag();
+
+      Triangulation<2>::active_cell_iterator cell, endc;
+      for (int i=0; i<8; ++i) 
+       {
+         int n_levels = tr.n_levels();
+         cell = tr.begin_active();
+         endc = tr.end();
+         
+         for (; cell!=endc; ++cell) 
+           {
+             double r      = rand()*1.0/RAND_MAX,
+                    weight = 1.*
+                             (cell->level()*cell->level()) /
+                             (n_levels*n_levels);
+             
+             if (r <= 0.5*weight)
+               cell->set_refine_flag ();
+           };
+         
+         tr.execute_coarsening_and_refinement ();
+       };
+//        Triangulation<2>::active_cell_iterator cell = tr.begin_active();
+//        cell->set_refine_flag();
       tr.execute_coarsening_and_refinement();
 
-      tr.refine_global(1);
+//      tr.refine_global(4);
+      if (true)
+       {
+         ofstream x ("T");
+         GridOut().write_eps (tr, x);
+       };
+      
       dof.distribute_dofs(*fe);
       const unsigned int size = dof.n_dofs();
       deallog << "DoFs " << size << endl;
@@ -166,7 +194,7 @@ int main()
          SolverCG<> cgsolver(cgcontrol, cgmem);
          PreconditionIdentity cgprec;
          MGCoarseGridLACIteration<SolverCG<>, SparseMatrix<double>, PreconditionIdentity>
-           coarse(cgsolver, mgA[tr.n_levels()-2], cgprec);
+           coarse(cgsolver, mgA[0], cgprec);
          
          MGSmootherRelaxation<double>
            smoother(mgdof, mgA, &SparseMatrix<double>::template precondition_SSOR<double>,
@@ -177,7 +205,7 @@ int main()
          transfer.build_matrices(mgdof);
 
 
-Multigrid<2> multigrid(mgdof, hanging_nodes, mgstruct, mgA, transfer, tr.n_levels()-2);
+Multigrid<2> multigrid(mgdof, hanging_nodes, mgstruct, mgA, transfer, 0);
          PreconditionMG<Multigrid<2> >
            mgprecondition(multigrid, smoother, smoother, coarse);
 

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.