--- /dev/null
+//---------------------------- rt_11.cc ---------------------------
+// rt_11.cc,v 1.3 2003/06/09 16:00:38 wolf Exp
+// Version:
+//
+// Copyright (C) 2003, 2005, 2006 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//---------------------------- rt_11.cc ---------------------------
+
+// Observe how the values of the shape functions change as we make a
+// cell smaller and smaller. Evaluate the values with FEFaceValues, to
+// make sure the values scale as in rt_10 where we used FEValues
+
+#include "../tests.h"
+#include <base/quadrature_lib.h>
+#include <base/logstream.h>
+#include <lac/vector.h>
+#include <lac/solver_cg.h>
+#include <lac/precondition.h>
+#include <lac/vector_memory.h>
+#include <lac/sparse_matrix.h>
+#include <lac/sparsity_pattern.h>
+#include <grid/tria.h>
+#include <grid/tria_iterator.h>
+#include <dofs/dof_accessor.h>
+#include <grid/grid_generator.h>
+#include <grid/grid_tools.h>
+#include <fe/fe_raviart_thomas.h>
+#include <fe/fe_values.h>
+
+#include <vector>
+#include <fstream>
+#include <string>
+
+#define PRECISION 2
+
+
+std::ofstream logfile ("rt_11/output");
+
+template<int dim>
+void
+test (const unsigned int degree)
+{
+ FE_RaviartThomas<dim> fe_rt(degree);
+
+ deallog << "Degree=" << degree
+ << std::endl;
+
+ for (double h=1; h>1./128; h/=2)
+ {
+ deallog << " h=" << h
+ << std::endl;
+
+ Triangulation<dim> tr;
+ GridGenerator::hyper_cube(tr, 0., h);
+
+ DoFHandler<dim> dof(tr);
+ dof.distribute_dofs(fe_rt);
+
+ QTrapez<dim-1> quadrature;
+
+ FEFaceValues<dim> fe_values (fe_rt, quadrature, update_values);
+ fe_values.reinit (dof.begin_active(), 0);
+ for (unsigned int q=0; q<quadrature.n_quadrature_points; ++q)
+ {
+ deallog << " Quadrature point " << q << ": ";
+ for (unsigned int i=0; i<fe_rt.dofs_per_cell; ++i)
+ {
+ deallog << '[';
+ for (unsigned int c=0; c<fe_rt.n_components(); ++c)
+ deallog << fe_values.shape_value_component(i,q,c) << ' ';
+ deallog << ']';
+ }
+ deallog << std::endl;
+ }
+ }
+}
+
+
+int
+main()
+{
+ logfile.precision (PRECISION);
+ logfile.setf(std::ios::fixed);
+ deallog.attach(logfile);
+ deallog.depth_console(0);
+ deallog.threshold_double(1.e-10);
+
+ for (unsigned int i=0; i<4; ++i)
+ test<2>(i);
+
+ return 0;
+}
+
+
+
--- /dev/null
+
+DEAL::Degree=0
+DEAL:: h=1.00
+DEAL:: Quadrature point 0: [1.00 0 ][0 0 ][0 1.00 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][1.00 0 ][0 1.00 ][0 0 ]
+DEAL:: Quadrature point 2: [1.00 0 ][0 0 ][0 0 ][0 1.00 ]
+DEAL:: Quadrature point 3: [0 0 ][1.00 0 ][0 0 ][0 1.00 ]
+DEAL:: h=0.50
+DEAL:: Quadrature point 0: [2.00 0 ][0 0 ][0 2.00 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][2.00 0 ][0 2.00 ][0 0 ]
+DEAL:: Quadrature point 2: [2.00 0 ][0 0 ][0 0 ][0 2.00 ]
+DEAL:: Quadrature point 3: [0 0 ][2.00 0 ][0 0 ][0 2.00 ]
+DEAL:: h=0.25
+DEAL:: Quadrature point 0: [4.00 0 ][0 0 ][0 4.00 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][4.00 0 ][0 4.00 ][0 0 ]
+DEAL:: Quadrature point 2: [4.00 0 ][0 0 ][0 0 ][0 4.00 ]
+DEAL:: Quadrature point 3: [0 0 ][4.00 0 ][0 0 ][0 4.00 ]
+DEAL:: h=0.12
+DEAL:: Quadrature point 0: [8.00 0 ][0 0 ][0 8.00 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][8.00 0 ][0 8.00 ][0 0 ]
+DEAL:: Quadrature point 2: [8.00 0 ][0 0 ][0 0 ][0 8.00 ]
+DEAL:: Quadrature point 3: [0 0 ][8.00 0 ][0 0 ][0 8.00 ]
+DEAL:: h=0.06
+DEAL:: Quadrature point 0: [16.00 0 ][0 0 ][0 16.00 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][16.00 0 ][0 16.00 ][0 0 ]
+DEAL:: Quadrature point 2: [16.00 0 ][0 0 ][0 0 ][0 16.00 ]
+DEAL:: Quadrature point 3: [0 0 ][16.00 0 ][0 0 ][0 16.00 ]
+DEAL:: h=0.03
+DEAL:: Quadrature point 0: [32.00 0 ][0 0 ][0 32.00 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][32.00 0 ][0 32.00 ][0 0 ]
+DEAL:: Quadrature point 2: [32.00 0 ][0 0 ][0 0 ][0 32.00 ]
+DEAL:: Quadrature point 3: [0 0 ][32.00 0 ][0 0 ][0 32.00 ]
+DEAL:: h=0.02
+DEAL:: Quadrature point 0: [64.00 0 ][0 0 ][0 64.00 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][64.00 0 ][0 64.00 ][0 0 ]
+DEAL:: Quadrature point 2: [64.00 0 ][0 0 ][0 0 ][0 64.00 ]
+DEAL:: Quadrature point 3: [0 0 ][64.00 0 ][0 0 ][0 64.00 ]
+DEAL::Degree=1
+DEAL:: h=1.00
+DEAL:: Quadrature point 0: [1.00 0 ][-1.73 0 ][0 0 ][0 0 ][0 1.00 ][0 -1.73 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][0 0 ][1.00 0 ][-1.73 0 ][0 1.00 ][0 1.73 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 2: [1.00 0 ][1.73 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 1.00 ][0 -1.73 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 3: [0 0 ][0 0 ][1.00 0 ][1.73 0 ][0 0 ][0 0 ][0 1.00 ][0 1.73 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: h=0.50
+DEAL:: Quadrature point 0: [2.00 0 ][-3.46 0 ][0 0 ][0 0 ][0 2.00 ][0 -3.46 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][0 0 ][2.00 0 ][-3.46 0 ][0 2.00 ][0 3.46 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 2: [2.00 0 ][3.46 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 2.00 ][0 -3.46 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 3: [0 0 ][0 0 ][2.00 0 ][3.46 0 ][0 0 ][0 0 ][0 2.00 ][0 3.46 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: h=0.25
+DEAL:: Quadrature point 0: [4.00 0 ][-6.93 0 ][0 0 ][0 0 ][0 4.00 ][0 -6.93 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][0 0 ][4.00 0 ][-6.93 0 ][0 4.00 ][0 6.93 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 2: [4.00 0 ][6.93 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 4.00 ][0 -6.93 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 3: [0 0 ][0 0 ][4.00 0 ][6.93 0 ][0 0 ][0 0 ][0 4.00 ][0 6.93 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: h=0.12
+DEAL:: Quadrature point 0: [8.00 0 ][-13.86 0 ][0 0 ][0 0 ][0 8.00 ][0 -13.86 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][0 0 ][8.00 0 ][-13.86 0 ][0 8.00 ][0 13.86 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 2: [8.00 0 ][13.86 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 8.00 ][0 -13.86 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 3: [0 0 ][0 0 ][8.00 0 ][13.86 0 ][0 0 ][0 0 ][0 8.00 ][0 13.86 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: h=0.06
+DEAL:: Quadrature point 0: [16.00 0 ][-27.71 0 ][0 0 ][0 0 ][0 16.00 ][0 -27.71 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][0 0 ][16.00 0 ][-27.71 0 ][0 16.00 ][0 27.71 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 2: [16.00 0 ][27.71 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 16.00 ][0 -27.71 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 3: [0 0 ][0 0 ][16.00 0 ][27.71 0 ][0 0 ][0 0 ][0 16.00 ][0 27.71 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: h=0.03
+DEAL:: Quadrature point 0: [32.00 0 ][-55.43 0 ][0 0 ][0 0 ][0 32.00 ][0 -55.43 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][0 0 ][32.00 0 ][-55.43 0 ][0 32.00 ][0 55.43 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 2: [32.00 0 ][55.43 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 32.00 ][0 -55.43 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 3: [0 0 ][0 0 ][32.00 0 ][55.43 0 ][0 0 ][0 0 ][0 32.00 ][0 55.43 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: h=0.02
+DEAL:: Quadrature point 0: [64.00 0 ][-110.85 0 ][0 0 ][0 0 ][0 64.00 ][0 -110.85 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][0 0 ][64.00 0 ][-110.85 0 ][0 64.00 ][0 110.85 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 2: [64.00 0 ][110.85 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 64.00 ][0 -110.85 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 3: [0 0 ][0 0 ][64.00 0 ][110.85 0 ][0 0 ][0 0 ][0 64.00 ][0 110.85 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL::Degree=2
+DEAL:: h=1.00
+DEAL:: Quadrature point 0: [1.00 0 ][-1.73 0 ][2.24 0 ][0 0 ][0 0 ][0 0 ][0 1.00 ][0 -1.73 ][0 2.24 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][1.00 0 ][-1.73 0 ][2.24 0 ][0 1.00 ][0 1.73 ][0 2.24 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 2: [1.00 0 ][1.73 0 ][2.24 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 1.00 ][0 -1.73 ][0 2.24 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][1.00 0 ][1.73 0 ][2.24 0 ][0 0 ][0 0 ][0 0 ][0 1.00 ][0 1.73 ][0 2.24 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: h=0.50
+DEAL:: Quadrature point 0: [2.00 0 ][-3.46 0 ][4.47 0 ][0 0 ][0 0 ][0 0 ][0 2.00 ][0 -3.46 ][0 4.47 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][2.00 0 ][-3.46 0 ][4.47 0 ][0 2.00 ][0 3.46 ][0 4.47 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 2: [2.00 0 ][3.46 0 ][4.47 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 2.00 ][0 -3.46 ][0 4.47 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][2.00 0 ][3.46 0 ][4.47 0 ][0 0 ][0 0 ][0 0 ][0 2.00 ][0 3.46 ][0 4.47 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: h=0.25
+DEAL:: Quadrature point 0: [4.00 0 ][-6.93 0 ][8.94 0 ][0 0 ][0 0 ][0 0 ][0 4.00 ][0 -6.93 ][0 8.94 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][4.00 0 ][-6.93 0 ][8.94 0 ][0 4.00 ][0 6.93 ][0 8.94 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 2: [4.00 0 ][6.93 0 ][8.94 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 4.00 ][0 -6.93 ][0 8.94 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][4.00 0 ][6.93 0 ][8.94 0 ][0 0 ][0 0 ][0 0 ][0 4.00 ][0 6.93 ][0 8.94 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: h=0.12
+DEAL:: Quadrature point 0: [8.00 0 ][-13.86 0 ][17.89 0 ][0 0 ][0 0 ][0 0 ][0 8.00 ][0 -13.86 ][0 17.89 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][8.00 0 ][-13.86 0 ][17.89 0 ][0 8.00 ][0 13.86 ][0 17.89 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 2: [8.00 0 ][13.86 0 ][17.89 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 8.00 ][0 -13.86 ][0 17.89 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][8.00 0 ][13.86 0 ][17.89 0 ][0 0 ][0 0 ][0 0 ][0 8.00 ][0 13.86 ][0 17.89 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: h=0.06
+DEAL:: Quadrature point 0: [16.00 0 ][-27.71 0 ][35.78 0 ][0 0 ][0 0 ][0 0 ][0 16.00 ][0 -27.71 ][0 35.78 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][16.00 0 ][-27.71 0 ][35.78 0 ][0 16.00 ][0 27.71 ][0 35.78 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 2: [16.00 0 ][27.71 0 ][35.78 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 16.00 ][0 -27.71 ][0 35.78 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][16.00 0 ][27.71 0 ][35.78 0 ][0 0 ][0 0 ][0 0 ][0 16.00 ][0 27.71 ][0 35.78 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: h=0.03
+DEAL:: Quadrature point 0: [32.00 0 ][-55.43 0 ][71.55 0 ][0 0 ][0 0 ][0 0 ][0 32.00 ][0 -55.43 ][0 71.55 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][32.00 0 ][-55.43 0 ][71.55 0 ][0 32.00 ][0 55.43 ][0 71.55 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 2: [32.00 0 ][55.43 0 ][71.55 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 32.00 ][0 -55.43 ][0 71.55 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][32.00 0 ][55.43 0 ][71.55 0 ][0 0 ][0 0 ][0 0 ][0 32.00 ][0 55.43 ][0 71.55 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: h=0.02
+DEAL:: Quadrature point 0: [64.00 0 ][-110.85 0 ][143.11 0 ][0 0 ][0 0 ][0 0 ][0 64.00 ][0 -110.85 ][0 143.11 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][64.00 0 ][-110.85 0 ][143.11 0 ][0 64.00 ][0 110.85 ][0 143.11 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 2: [64.00 0 ][110.85 0 ][143.11 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 64.00 ][0 -110.85 ][0 143.11 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][64.00 0 ][110.85 0 ][143.11 0 ][0 0 ][0 0 ][0 0 ][0 64.00 ][0 110.85 ][0 143.11 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL::Degree=3
+DEAL:: h=1.00
+DEAL:: Quadrature point 0: [1.00 0 ][-1.73 0 ][2.24 0 ][-2.65 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 1.00 ][0 -1.73 ][0 2.24 ][0 -2.65 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][0 0 ][1.00 0 ][-1.73 0 ][2.24 0 ][-2.65 0 ][0 1.00 ][0 1.73 ][0 2.24 ][0 2.65 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 2: [1.00 0 ][1.73 0 ][2.24 0 ][2.65 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 1.00 ][0 -1.73 ][0 2.24 ][0 -2.65 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][0 0 ][1.00 0 ][1.73 0 ][2.24 0 ][2.65 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 1.00 ][0 1.73 ][0 2.24 ][0 2.65 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: h=0.50
+DEAL:: Quadrature point 0: [2.00 0 ][-3.46 0 ][4.47 0 ][-5.29 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 2.00 ][0 -3.46 ][0 4.47 ][0 -5.29 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][0 0 ][2.00 0 ][-3.46 0 ][4.47 0 ][-5.29 0 ][0 2.00 ][0 3.46 ][0 4.47 ][0 5.29 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 2: [2.00 0 ][3.46 0 ][4.47 0 ][5.29 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 2.00 ][0 -3.46 ][0 4.47 ][0 -5.29 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][0 0 ][2.00 0 ][3.46 0 ][4.47 0 ][5.29 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 2.00 ][0 3.46 ][0 4.47 ][0 5.29 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: h=0.25
+DEAL:: Quadrature point 0: [4.00 0 ][-6.93 0 ][8.94 0 ][-10.58 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 4.00 ][0 -6.93 ][0 8.94 ][0 -10.58 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][0 0 ][4.00 0 ][-6.93 0 ][8.94 0 ][-10.58 0 ][0 4.00 ][0 6.93 ][0 8.94 ][0 10.58 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 2: [4.00 0 ][6.93 0 ][8.94 0 ][10.58 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 4.00 ][0 -6.93 ][0 8.94 ][0 -10.58 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][0 0 ][4.00 0 ][6.93 0 ][8.94 0 ][10.58 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 4.00 ][0 6.93 ][0 8.94 ][0 10.58 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: h=0.12
+DEAL:: Quadrature point 0: [8.00 0 ][-13.86 0 ][17.89 0 ][-21.17 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 8.00 ][0 -13.86 ][0 17.89 ][0 -21.17 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][0 0 ][8.00 0 ][-13.86 0 ][17.89 0 ][-21.17 0 ][0 8.00 ][0 13.86 ][0 17.89 ][0 21.17 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 2: [8.00 0 ][13.86 0 ][17.89 0 ][21.17 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 8.00 ][0 -13.86 ][0 17.89 ][0 -21.17 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][0 0 ][8.00 0 ][13.86 0 ][17.89 0 ][21.17 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 8.00 ][0 13.86 ][0 17.89 ][0 21.17 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: h=0.06
+DEAL:: Quadrature point 0: [16.00 0 ][-27.71 0 ][35.78 0 ][-42.33 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 16.00 ][0 -27.71 ][0 35.78 ][0 -42.33 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][0 0 ][16.00 0 ][-27.71 0 ][35.78 0 ][-42.33 0 ][0 16.00 ][0 27.71 ][0 35.78 ][0 42.33 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 2: [16.00 0 ][27.71 0 ][35.78 0 ][42.33 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 16.00 ][0 -27.71 ][0 35.78 ][0 -42.33 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][0 0 ][16.00 0 ][27.71 0 ][35.78 0 ][42.33 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 16.00 ][0 27.71 ][0 35.78 ][0 42.33 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: h=0.03
+DEAL:: Quadrature point 0: [32.00 0 ][-55.43 0 ][71.55 0 ][-84.66 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 32.00 ][0 -55.43 ][0 71.55 ][0 -84.66 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][0 0 ][32.00 0 ][-55.43 0 ][71.55 0 ][-84.66 0 ][0 32.00 ][0 55.43 ][0 71.55 ][0 84.66 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 2: [32.00 0 ][55.43 0 ][71.55 0 ][84.66 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 32.00 ][0 -55.43 ][0 71.55 ][0 -84.66 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][0 0 ][32.00 0 ][55.43 0 ][71.55 0 ][84.66 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 32.00 ][0 55.43 ][0 71.55 ][0 84.66 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: h=0.02
+DEAL:: Quadrature point 0: [64.00 0 ][-110.85 0 ][143.11 0 ][-169.33 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 64.00 ][0 -110.85 ][0 143.11 ][0 -169.33 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][0 0 ][64.00 0 ][-110.85 0 ][143.11 0 ][-169.33 0 ][0 64.00 ][0 110.85 ][0 143.11 ][0 169.33 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 2: [64.00 0 ][110.85 0 ][143.11 0 ][169.33 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 64.00 ][0 -110.85 ][0 143.11 ][0 -169.33 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][0 0 ][64.00 0 ][110.85 0 ][143.11 0 ][169.33 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 64.00 ][0 110.85 ][0 143.11 ][0 169.33 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]