]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Patch by Denis Davydov: Get the tensor extractors working right.
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Tue, 9 Jul 2013 21:37:20 +0000 (21:37 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Tue, 9 Jul 2013 21:37:20 +0000 (21:37 +0000)
git-svn-id: https://svn.dealii.org/trunk@29959 0785d39b-7218-0410-832d-ea1e28bc413d

14 files changed:
deal.II/include/deal.II/fe/fe_values.h
deal.II/source/fe/fe_values.cc
tests/deal.II/fe_values_view_21_nonsymmetric/cmp/generic
tests/deal.II/fe_values_view_23_nonsymmetric.cc [new file with mode: 0644]
tests/deal.II/fe_values_view_23_nonsymmetric/cmp/generic [new file with mode: 0644]
tests/deal.II/fe_values_view_23_nonsymmetric/fe_values_view_09/cmp/generic [new file with mode: 0644]
tests/deal.II/fe_values_view_24_nonsymmetric.cc [new file with mode: 0644]
tests/deal.II/fe_values_view_24_nonsymmetric/cmp/generic [new file with mode: 0644]
tests/deal.II/fe_values_view_24_nonsymmetric/fe_values_view_09/cmp/generic [new file with mode: 0644]
tests/deal.II/fe_values_view_25_nonsymmetric.cc [new file with mode: 0644]
tests/deal.II/fe_values_view_25_nonsymmetric/cmp/generic [new file with mode: 0644]
tests/deal.II/fe_values_view_25_nonsymmetric/fe_values_view_09/cmp/generic [new file with mode: 0644]
tests/deal.II/project_common.h
tests/deal.II/project_q_01.cc

index b9da03e4df24d23827509b67dd9fd7d0d3001b31..e3edf8f049e5e25d564ac85014bdb0f8c77025e1 100644 (file)
@@ -4651,7 +4651,7 @@ namespace FEValuesViews
         const unsigned int comp =
           shape_function_data[shape_function].single_nonzero_component_index;
         const TableIndices<2> indices = dealii::Tensor<2,spacedim>::unrolled_to_component_indices(comp);
-        return_value[indices] = fe_values.shape_values(snc,q_point);//last index first [jj][ii]
+        return_value[indices] = fe_values.shape_values(snc,q_point);
         return return_value;
       }
     else
@@ -4661,7 +4661,7 @@ namespace FEValuesViews
           if (shape_function_data[shape_function].is_nonzero_shape_function_component[d]) {
             const TableIndices<2> indices = dealii::Tensor<2,spacedim>::unrolled_to_component_indices(d);
             return_value[indices]
-              = fe_values.shape_values(shape_function_data[shape_function].row_index[d],q_point);//last index first [jj][ii]
+              = fe_values.shape_values(shape_function_data[shape_function].row_index[d],q_point);
           }
         return return_value;
       }
index 23b1b8c22617d003c4ac80f4ddde090c9a5088e7..1a54152551fc137abf52a7c8d2ec193e4c9a9d16 100644 (file)
@@ -1173,7 +1173,7 @@ namespace FEValuesViews
 
                   const double *shape_value_ptr = &shape_values(snc,0);
                   for (unsigned int q_point=0; q_point<n_quadrature_points; ++q_point)
-                    values[q_point][indices] += value * *shape_value_ptr++;//last index first [j][i]
+                    values[q_point][indices] += value * *shape_value_ptr++;
                 }
               else
                 for (unsigned int d=0;
@@ -1185,7 +1185,7 @@ namespace FEValuesViews
                       const double *shape_value_ptr =
                         &shape_values(shape_function_data[shape_function].row_index[d],0);
                       for (unsigned int q_point=0; q_point<n_quadrature_points; ++q_point)
-                        values[q_point][indices] += value * *shape_value_ptr++;//last index first [j][i]
+                        values[q_point][indices] += value * *shape_value_ptr++;
                     }
             }
         }
index 80180cbeef65da584a2f61c078b1bc38912f2942..84c5609545e34aac2e5639a0f5c0e14949115d56 100644 (file)
@@ -6,4 +6,4 @@ DEAL::Number of degrees of freedom: 80 = (72 + 8)
 DEAL::dofs_per_cell: 80
 DEAL::fe.dofs_per_cell: 80     quadrature_formula.size(): 1
 DEAL::3.50 11.5 19.5 27.5 35.5 43.5 51.5 59.5 67.5
-DEAL::6.50 6.00 5.00
+DEAL::3.50 3.50 3.50
diff --git a/tests/deal.II/fe_values_view_23_nonsymmetric.cc b/tests/deal.II/fe_values_view_23_nonsymmetric.cc
new file mode 100644 (file)
index 0000000..9adfc38
--- /dev/null
@@ -0,0 +1,110 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007, 2008, 2010, 2011 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.
+//
+//----------------------------------------------------------------------
+
+
+// there was a bug in getting the divergence of shape functions for the
+// Tensor extractors. test that it is fixed by comparing with
+// get_function_divergences
+
+#include "../tests.h"
+#include <deal.II/base/logstream.h>
+#include <deal.II/base/function.h>
+#include <deal.II/base/quadrature_lib.h>
+#include <deal.II/lac/vector.h>
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/dofs/dof_handler.h>
+#include <deal.II/fe/fe_q.h>
+#include <deal.II/fe/fe_system.h>
+#include <deal.II/fe/fe_values.h>
+#include <deal.II/fe/mapping_q1.h>
+
+#include <fstream>
+
+
+
+template<int dim>
+void test (const Triangulation<dim>& tr,
+          const FiniteElement<dim>& fe)
+{
+  deallog << "FE=" << fe.get_name()
+         << std::endl;
+
+  DoFHandler<dim> dof(tr);
+  dof.distribute_dofs(fe);
+
+  Vector<double> fe_function(dof.n_dofs());
+  for (unsigned int i=0; i<dof.n_dofs(); ++i)
+    fe_function(i) = (i+1)*(i+2);
+  
+  const QGauss<dim> quadrature(2);
+  FEValues<dim> fe_values (fe, quadrature,
+                          update_values | update_gradients);
+  fe_values.reinit (dof.begin_active());
+
+                                  // let the FEValues object compute the
+                                  // divergences at quadrature points
+  std::vector<Tensor<1,dim> > divergences (quadrature.size());
+  FEValuesExtractors::Tensor<2> extractor(0);
+  fe_values[extractor]
+    .get_function_divergences (fe_function, divergences);
+
+                                  // now do the same "by hand"
+  std::vector<unsigned int> local_dof_indices (fe.dofs_per_cell);
+  dof.begin_active()->get_dof_indices (local_dof_indices);
+  
+  for (unsigned int q=0; q<quadrature.size(); ++q)
+    {
+      Tensor<1,dim> div_alt;
+      for (unsigned int i=0; i<fe.dofs_per_cell; ++i)
+       div_alt += fe_values[extractor].divergence (i,q) *
+                  fe_function(local_dof_indices[i]);
+      
+      deallog << "q_point=" << q << std::endl
+             << "   method 1: " << divergences[q] << std::endl
+             << "   method 2: " << div_alt << std::endl
+             << std::endl;
+      Assert ((divergences[q] - div_alt).norm() <= divergences[q].norm(),
+             ExcInternalError());
+    }
+}
+
+
+
+template<int dim>
+void test_hyper_sphere()
+{
+  Triangulation<dim> tr;
+  GridGenerator::hyper_ball(tr);
+
+  static const HyperBallBoundary<dim> boundary;
+  tr.set_boundary (0, boundary);
+
+  FESystem<dim> fe (FE_Q<dim>(1),
+                   Tensor<2,dim>::n_independent_components);
+  test(tr, fe);
+}
+
+
+int main()
+{
+  std::ofstream logfile ("fe_values_view_23_nonsymmetric/output");
+  deallog << std::setprecision (3);
+
+  deallog.attach(logfile);
+  deallog.depth_console (0);
+  deallog.threshold_double(1.e-7);
+
+  test_hyper_sphere<2>();
+  test_hyper_sphere<3>();
+}
diff --git a/tests/deal.II/fe_values_view_23_nonsymmetric/cmp/generic b/tests/deal.II/fe_values_view_23_nonsymmetric/cmp/generic
new file mode 100644 (file)
index 0000000..870eda4
--- /dev/null
@@ -0,0 +1,51 @@
+
+DEAL::FE=FESystem<2>[FE_Q<2>(1)^4]
+DEAL::q_point=0
+DEAL::   method 1: 329. 370.
+DEAL::   method 2: 329. 370.
+DEAL::
+DEAL::q_point=1
+DEAL::   method 1: 472. 521.
+DEAL::   method 2: 472. 521.
+DEAL::
+DEAL::q_point=2
+DEAL::   method 1: 354. 397.
+DEAL::   method 2: 354. 397.
+DEAL::
+DEAL::q_point=3
+DEAL::   method 1: 586. 642.
+DEAL::   method 2: 586. 642.
+DEAL::
+DEAL::FE=FESystem<3>[FE_Q<3>(1)^9]
+DEAL::q_point=0
+DEAL::   method 1: 8.02e+03 8.32e+03 8.61e+03
+DEAL::   method 2: 8.02e+03 8.32e+03 8.61e+03
+DEAL::
+DEAL::q_point=1
+DEAL::   method 1: 9.35e+03 9.64e+03 9.94e+03
+DEAL::   method 2: 9.35e+03 9.64e+03 9.94e+03
+DEAL::
+DEAL::q_point=2
+DEAL::   method 1: 1.02e+04 1.05e+04 1.08e+04
+DEAL::   method 2: 1.02e+04 1.05e+04 1.08e+04
+DEAL::
+DEAL::q_point=3
+DEAL::   method 1: 1.16e+04 1.19e+04 1.22e+04
+DEAL::   method 2: 1.16e+04 1.19e+04 1.22e+04
+DEAL::
+DEAL::q_point=4
+DEAL::   method 1: 1.07e+04 1.10e+04 1.13e+04
+DEAL::   method 2: 1.07e+04 1.10e+04 1.13e+04
+DEAL::
+DEAL::q_point=5
+DEAL::   method 1: 1.20e+04 1.23e+04 1.26e+04
+DEAL::   method 2: 1.20e+04 1.23e+04 1.26e+04
+DEAL::
+DEAL::q_point=6
+DEAL::   method 1: 1.29e+04 1.32e+04 1.35e+04
+DEAL::   method 2: 1.29e+04 1.32e+04 1.35e+04
+DEAL::
+DEAL::q_point=7
+DEAL::   method 1: 1.42e+04 1.45e+04 1.48e+04
+DEAL::   method 2: 1.42e+04 1.45e+04 1.48e+04
+DEAL::
diff --git a/tests/deal.II/fe_values_view_23_nonsymmetric/fe_values_view_09/cmp/generic b/tests/deal.II/fe_values_view_23_nonsymmetric/fe_values_view_09/cmp/generic
new file mode 100644 (file)
index 0000000..977609c
--- /dev/null
@@ -0,0 +1,91 @@
+
+DEAL::FE=FESystem<2>[FE_Q<2>(1)-FE_RaviartThomas<2>(1)-FE_Nedelec<2>(1)]
+DEAL::component=0
+DEAL::0.807 4.36
+DEAL::0.807 5.29
+DEAL::1.31 4.07
+DEAL::1.31 5.59
+DEAL::component=1
+DEAL::18.2 352.
+DEAL::-18.2 337.
+DEAL::333. 180.
+DEAL::-295. 296.
+DEAL::component=2
+DEAL::43.2 -31.7
+DEAL::43.2 118.
+DEAL::125. -61.1
+DEAL::125. -40.8
+DEAL::component=3
+DEAL::2.58e-16 5.84
+DEAL::2.58e-16 5.84
+DEAL::3.58e-16 9.52
+DEAL::3.58e-16 9.52
+DEAL::component=4
+DEAL::5.84 -6.75
+DEAL::5.84 6.75
+DEAL::9.52 -11.0
+DEAL::9.52 11.0
+DEAL::FE=FESystem<3>[FE_Q<3>(1)-FE_RaviartThomas<3>(1)-FE_Nedelec<3>(1)]
+DEAL::component=0
+DEAL::2.37 4.73 9.46
+DEAL::2.37 4.73 9.46
+DEAL::2.37 4.73 9.46
+DEAL::2.37 4.73 9.46
+DEAL::2.37 4.73 9.46
+DEAL::2.37 4.73 9.46
+DEAL::2.37 4.73 9.46
+DEAL::2.37 4.73 9.46
+DEAL::component=1
+DEAL::8.18e-13 -275. -138.
+DEAL::-5.71e-13 -275. -138.
+DEAL::-367. -275. 4.71e+03
+DEAL::367. -275. 4.92e+03
+DEAL::-184. 4.57e+03 -138.
+DEAL::184. 4.79e+03 -138.
+DEAL::4.80e+03 4.57e+03 4.71e+03
+DEAL::-4.38e+03 4.79e+03 4.92e+03
+DEAL::component=2
+DEAL::-275. 1.30e-12 -138.
+DEAL::-275. -459. 4.80e+03
+DEAL::-275. -1.40e-12 -138.
+DEAL::-275. 459. 5.02e+03
+DEAL::4.67e+03 -91.8 -138.
+DEAL::4.67e+03 3.52e+03 4.80e+03
+DEAL::4.88e+03 91.8 -138.
+DEAL::4.88e+03 -3.09e+03 5.02e+03
+DEAL::component=3
+DEAL::-275. -138. 1.36e-12
+DEAL::-275. 4.90e+03 -367.
+DEAL::4.76e+03 -138. -184.
+DEAL::4.76e+03 4.90e+03 2.23e+03
+DEAL::-275. -138. -3.64e-12
+DEAL::-275. 5.11e+03 367.
+DEAL::4.97e+03 -138. 184.
+DEAL::4.97e+03 5.11e+03 -1.81e+03
+DEAL::component=4
+DEAL::-7.70e-16 5.60 22.4
+DEAL::-7.70e-16 5.60 22.4
+DEAL::-6.69e-16 5.60 22.4
+DEAL::-6.69e-16 5.60 22.4
+DEAL::1.58e-16 5.60 22.4
+DEAL::1.58e-16 5.60 22.4
+DEAL::9.52e-16 5.60 22.4
+DEAL::9.52e-16 5.60 22.4
+DEAL::component=5
+DEAL::5.60 5.78e-16 22.4
+DEAL::5.60 -1.39e-15 22.4
+DEAL::5.60 5.78e-16 22.4
+DEAL::5.60 -1.39e-15 22.4
+DEAL::5.60 -6.60e-16 22.4
+DEAL::5.60 3.53e-16 22.4
+DEAL::5.60 -6.60e-16 22.4
+DEAL::5.60 3.53e-16 22.4
+DEAL::component=6
+DEAL::5.60 11.2 -1.38e-15
+DEAL::5.60 11.2 1.96e-16
+DEAL::5.60 11.2 1.23e-17
+DEAL::5.60 11.2 1.50e-15
+DEAL::5.60 11.2 -1.38e-15
+DEAL::5.60 11.2 1.96e-16
+DEAL::5.60 11.2 1.23e-17
+DEAL::5.60 11.2 1.50e-15
diff --git a/tests/deal.II/fe_values_view_24_nonsymmetric.cc b/tests/deal.II/fe_values_view_24_nonsymmetric.cc
new file mode 100644 (file)
index 0000000..7d1335f
--- /dev/null
@@ -0,0 +1,104 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007, 2008, 2010, 2011 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.
+//
+//----------------------------------------------------------------------
+
+
+// there was a bug in getting the divergence of shape functions for the
+// Tensor extractors. test that it is fixed by comparing with
+// get_function_divergences
+
+#include "../tests.h"
+#include <deal.II/base/logstream.h>
+#include <deal.II/base/function.h>
+#include <deal.II/base/quadrature_lib.h>
+#include <deal.II/lac/vector.h>
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/dofs/dof_handler.h>
+#include <deal.II/fe/fe_q.h>
+#include <deal.II/fe/fe_system.h>
+#include <deal.II/fe/fe_values.h>
+#include <deal.II/fe/mapping_q1.h>
+
+#include <fstream>
+
+
+
+template<int dim>
+void test (const Triangulation<dim>& tr,
+          const FiniteElement<dim>& fe)
+{
+  deallog << "FE=" << fe.get_name()
+         << std::endl;
+
+  DoFHandler<dim> dof(tr);
+  dof.distribute_dofs(fe);
+
+  Vector<double> fe_function(dof.n_dofs());
+  for (unsigned int i=0; i<dof.n_dofs(); ++i)
+    fe_function(i) = i+1;
+  
+  const QGauss<dim> quadrature(2);
+  FEValues<dim> fe_values (fe, quadrature,
+                          update_values | update_gradients);
+  fe_values.reinit (dof.begin_active());
+
+                                  // let the FEValues object compute the
+                                  // divergences at quadrature points
+  std::vector<Tensor<1,dim> > divergences (quadrature.size());
+  FEValuesExtractors::Tensor<2> extractor(0);
+  fe_values[extractor]
+    .get_function_divergences (fe_function, divergences);
+
+                                  // now do the same "by hand"
+  std::vector<unsigned int> local_dof_indices (fe.dofs_per_cell);
+  dof.begin_active()->get_dof_indices (local_dof_indices);
+  
+  for (unsigned int i=0; i<fe.dofs_per_cell; ++i)
+    {
+      deallog << "i=" << i << std::endl;
+      
+      for (unsigned int q=0; q<quadrature.size(); ++q)
+       deallog << "  q_point=" << q << std::endl
+               << "    div= " << fe_values[extractor].divergence (i,q) << std::endl;
+    }
+}
+
+
+
+template<int dim>
+void test_hyper_sphere()
+{
+  Triangulation<dim> tr;
+  GridGenerator::hyper_ball(tr);
+
+  static const HyperBallBoundary<dim> boundary;
+  tr.set_boundary (0, boundary);
+
+  FESystem<dim> fe (FE_Q<dim>(1),
+                   Tensor<2,dim>::n_independent_components);
+  test(tr, fe);
+}
+
+
+int main()
+{
+  std::ofstream logfile ("fe_values_view_24_nonsymmetric/output");
+  deallog << std::setprecision (3);
+
+  deallog.attach(logfile);
+  deallog.depth_console (0);
+  deallog.threshold_double(1.e-7);
+
+  test_hyper_sphere<2>();
+  test_hyper_sphere<3>();
+}
diff --git a/tests/deal.II/fe_values_view_24_nonsymmetric/cmp/generic b/tests/deal.II/fe_values_view_24_nonsymmetric/cmp/generic
new file mode 100644 (file)
index 0000000..c48eecd
--- /dev/null
@@ -0,0 +1,1371 @@
+
+DEAL::FE=FESystem<2>[FE_Q<2>(1)^4]
+DEAL::i=0
+DEAL::  q_point=0
+DEAL::    div= -0.636 0.00
+DEAL::  q_point=1
+DEAL::    div= -0.636 0.00
+DEAL::  q_point=2
+DEAL::    div= -0.278 0.00
+DEAL::  q_point=3
+DEAL::    div= -0.278 0.00
+DEAL::i=1
+DEAL::  q_point=0
+DEAL::    div= 0.00 -0.636
+DEAL::  q_point=1
+DEAL::    div= 0.00 -0.636
+DEAL::  q_point=2
+DEAL::    div= 0.00 -0.278
+DEAL::  q_point=3
+DEAL::    div= 0.00 -0.278
+DEAL::i=2
+DEAL::  q_point=0
+DEAL::    div= -1.54 0.00
+DEAL::  q_point=1
+DEAL::    div= -0.878 0.00
+DEAL::  q_point=2
+DEAL::    div= -1.74 0.00
+DEAL::  q_point=3
+DEAL::    div= -0.671 0.00
+DEAL::i=3
+DEAL::  q_point=0
+DEAL::    div= 0.00 -1.54
+DEAL::  q_point=1
+DEAL::    div= 0.00 -0.878
+DEAL::  q_point=2
+DEAL::    div= 0.00 -1.74
+DEAL::  q_point=3
+DEAL::    div= 0.00 -0.671
+DEAL::i=4
+DEAL::  q_point=0
+DEAL::    div= 0.636 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.636 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.278 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.278 0.00
+DEAL::i=5
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.636
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.636
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.278
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.278
+DEAL::i=6
+DEAL::  q_point=0
+DEAL::    div= -0.878 0.00
+DEAL::  q_point=1
+DEAL::    div= -1.54 0.00
+DEAL::  q_point=2
+DEAL::    div= -0.671 0.00
+DEAL::  q_point=3
+DEAL::    div= -1.74 0.00
+DEAL::i=7
+DEAL::  q_point=0
+DEAL::    div= 0.00 -0.878
+DEAL::  q_point=1
+DEAL::    div= 0.00 -1.54
+DEAL::  q_point=2
+DEAL::    div= 0.00 -0.671
+DEAL::  q_point=3
+DEAL::    div= 0.00 -1.74
+DEAL::i=8
+DEAL::  q_point=0
+DEAL::    div= -0.171 0.00
+DEAL::  q_point=1
+DEAL::    div= -0.171 0.00
+DEAL::  q_point=2
+DEAL::    div= -1.04 0.00
+DEAL::  q_point=3
+DEAL::    div= -1.04 0.00
+DEAL::i=9
+DEAL::  q_point=0
+DEAL::    div= 0.00 -0.171
+DEAL::  q_point=1
+DEAL::    div= 0.00 -0.171
+DEAL::  q_point=2
+DEAL::    div= 0.00 -1.04
+DEAL::  q_point=3
+DEAL::    div= 0.00 -1.04
+DEAL::i=10
+DEAL::  q_point=0
+DEAL::    div= 2.00 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.412 0.00
+DEAL::  q_point=2
+DEAL::    div= 2.50 0.00
+DEAL::  q_point=3
+DEAL::    div= -0.0883 0.00
+DEAL::i=11
+DEAL::  q_point=0
+DEAL::    div= 0.00 2.00
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.412
+DEAL::  q_point=2
+DEAL::    div= 0.00 2.50
+DEAL::  q_point=3
+DEAL::    div= 0.00 -0.0883
+DEAL::i=12
+DEAL::  q_point=0
+DEAL::    div= 0.171 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.171 0.00
+DEAL::  q_point=2
+DEAL::    div= 1.04 0.00
+DEAL::  q_point=3
+DEAL::    div= 1.04 0.00
+DEAL::i=13
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.171
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.171
+DEAL::  q_point=2
+DEAL::    div= 0.00 1.04
+DEAL::  q_point=3
+DEAL::    div= 0.00 1.04
+DEAL::i=14
+DEAL::  q_point=0
+DEAL::    div= 0.412 0.00
+DEAL::  q_point=1
+DEAL::    div= 2.00 0.00
+DEAL::  q_point=2
+DEAL::    div= -0.0883 0.00
+DEAL::  q_point=3
+DEAL::    div= 2.50 0.00
+DEAL::i=15
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.412
+DEAL::  q_point=1
+DEAL::    div= 0.00 2.00
+DEAL::  q_point=2
+DEAL::    div= 0.00 -0.0883
+DEAL::  q_point=3
+DEAL::    div= 0.00 2.50
+DEAL::FE=FESystem<3>[FE_Q<3>(1)^9]
+DEAL::i=0
+DEAL::  q_point=0
+DEAL::    div= -1.47 0.00 0.00
+DEAL::  q_point=1
+DEAL::    div= -1.47 0.00 0.00
+DEAL::  q_point=2
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=3
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=4
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=5
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=6
+DEAL::    div= -0.106 0.00 0.00
+DEAL::  q_point=7
+DEAL::    div= -0.106 0.00 0.00
+DEAL::i=1
+DEAL::  q_point=0
+DEAL::    div= 0.00 -1.47 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.00 -1.47 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.00 -0.106 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.00 -0.106 0.00
+DEAL::i=2
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.00 -1.47
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.00 -1.47
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.00 -0.106
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.00 -0.106
+DEAL::i=3
+DEAL::  q_point=0
+DEAL::    div= -1.47 0.00 0.00
+DEAL::  q_point=1
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=2
+DEAL::    div= -1.47 0.00 0.00
+DEAL::  q_point=3
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=4
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=5
+DEAL::    div= -0.106 0.00 0.00
+DEAL::  q_point=6
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=7
+DEAL::    div= -0.106 0.00 0.00
+DEAL::i=4
+DEAL::  q_point=0
+DEAL::    div= 0.00 -1.47 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.00 -1.47 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.00 -0.106 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.00 -0.106 0.00
+DEAL::i=5
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.00 -1.47
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.00 -1.47
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.00 -0.106
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.00 -0.106
+DEAL::i=6
+DEAL::  q_point=0
+DEAL::    div= -1.47 0.00 0.00
+DEAL::  q_point=1
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=2
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=3
+DEAL::    div= -0.106 0.00 0.00
+DEAL::  q_point=4
+DEAL::    div= -1.47 0.00 0.00
+DEAL::  q_point=5
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=6
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=7
+DEAL::    div= -0.106 0.00 0.00
+DEAL::i=7
+DEAL::  q_point=0
+DEAL::    div= 0.00 -1.47 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.00 -0.106 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.00 -1.47 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.00 -0.106 0.00
+DEAL::i=8
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.00 -1.47
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.00 -0.106
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.00 -1.47
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.00 -0.106
+DEAL::i=9
+DEAL::  q_point=0
+DEAL::    div= 1.47 0.00 0.00
+DEAL::  q_point=1
+DEAL::    div= 1.47 0.00 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.106 0.00 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.106 0.00 0.00
+DEAL::i=10
+DEAL::  q_point=0
+DEAL::    div= 0.00 1.47 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.00 1.47 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.106 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.106 0.00
+DEAL::i=11
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.00 1.47
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.00 1.47
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.00 0.106
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.00 0.106
+DEAL::i=12
+DEAL::  q_point=0
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=1
+DEAL::    div= -1.47 0.00 0.00
+DEAL::  q_point=2
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=3
+DEAL::    div= -1.47 0.00 0.00
+DEAL::  q_point=4
+DEAL::    div= -0.106 0.00 0.00
+DEAL::  q_point=5
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=6
+DEAL::    div= -0.106 0.00 0.00
+DEAL::  q_point=7
+DEAL::    div= -0.394 0.00 0.00
+DEAL::i=13
+DEAL::  q_point=0
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.00 -1.47 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.00 -1.47 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.00 -0.106 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.00 -0.106 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::i=14
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.00 -1.47
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.00 -1.47
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.00 -0.106
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.00 -0.106
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::i=15
+DEAL::  q_point=0
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=1
+DEAL::    div= -1.47 0.00 0.00
+DEAL::  q_point=2
+DEAL::    div= -0.106 0.00 0.00
+DEAL::  q_point=3
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=4
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=5
+DEAL::    div= -1.47 0.00 0.00
+DEAL::  q_point=6
+DEAL::    div= -0.106 0.00 0.00
+DEAL::  q_point=7
+DEAL::    div= -0.394 0.00 0.00
+DEAL::i=16
+DEAL::  q_point=0
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.00 -1.47 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.00 -0.106 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.00 -1.47 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.00 -0.106 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::i=17
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.00 -1.47
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.00 -0.106
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.00 -1.47
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.00 -0.106
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::i=18
+DEAL::  q_point=0
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=1
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=2
+DEAL::    div= -1.47 0.00 0.00
+DEAL::  q_point=3
+DEAL::    div= -1.47 0.00 0.00
+DEAL::  q_point=4
+DEAL::    div= -0.106 0.00 0.00
+DEAL::  q_point=5
+DEAL::    div= -0.106 0.00 0.00
+DEAL::  q_point=6
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=7
+DEAL::    div= -0.394 0.00 0.00
+DEAL::i=19
+DEAL::  q_point=0
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.00 -1.47 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.00 -1.47 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.00 -0.106 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.00 -0.106 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::i=20
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.00 -1.47
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.00 -1.47
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.00 -0.106
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.00 -0.106
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::i=21
+DEAL::  q_point=0
+DEAL::    div= 1.47 0.00 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=2
+DEAL::    div= 1.47 0.00 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.106 0.00 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.106 0.00 0.00
+DEAL::i=22
+DEAL::  q_point=0
+DEAL::    div= 0.00 1.47 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.00 1.47 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.106 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.106 0.00
+DEAL::i=23
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.00 1.47
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.00 1.47
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.00 0.106
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.00 0.106
+DEAL::i=24
+DEAL::  q_point=0
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=1
+DEAL::    div= -0.106 0.00 0.00
+DEAL::  q_point=2
+DEAL::    div= -1.47 0.00 0.00
+DEAL::  q_point=3
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=4
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=5
+DEAL::    div= -0.106 0.00 0.00
+DEAL::  q_point=6
+DEAL::    div= -1.47 0.00 0.00
+DEAL::  q_point=7
+DEAL::    div= -0.394 0.00 0.00
+DEAL::i=25
+DEAL::  q_point=0
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.00 -0.106 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.00 -1.47 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.00 -0.106 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.00 -1.47 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::i=26
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.00 -0.106
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.00 -1.47
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.00 -0.106
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.00 -1.47
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::i=27
+DEAL::  q_point=0
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=2
+DEAL::    div= 1.47 0.00 0.00
+DEAL::  q_point=3
+DEAL::    div= 1.47 0.00 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.106 0.00 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.106 0.00 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.394 0.00 0.00
+DEAL::i=28
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.00 1.47 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.00 1.47 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.106 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.106 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.394 0.00
+DEAL::i=29
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.00 1.47
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.00 1.47
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.00 0.106
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.00 0.106
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.00 0.394
+DEAL::i=30
+DEAL::  q_point=0
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=1
+DEAL::    div= 1.47 0.00 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=3
+DEAL::    div= 1.47 0.00 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.106 0.00 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.106 0.00 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.394 0.00 0.00
+DEAL::i=31
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.00 1.47 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.00 1.47 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.106 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.106 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.394 0.00
+DEAL::i=32
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.00 1.47
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.00 1.47
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.00 0.106
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.00 0.106
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.00 0.394
+DEAL::i=33
+DEAL::  q_point=0
+DEAL::    div= -0.106 0.00 0.00
+DEAL::  q_point=1
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=2
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=3
+DEAL::    div= -1.47 0.00 0.00
+DEAL::  q_point=4
+DEAL::    div= -0.106 0.00 0.00
+DEAL::  q_point=5
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=6
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=7
+DEAL::    div= -1.47 0.00 0.00
+DEAL::i=34
+DEAL::  q_point=0
+DEAL::    div= 0.00 -0.106 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.00 -1.47 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.00 -0.106 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.00 -1.47 0.00
+DEAL::i=35
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.00 -0.106
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.00 -1.47
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.00 -0.106
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.00 -1.47
+DEAL::i=36
+DEAL::  q_point=0
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=1
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=2
+DEAL::    div= -0.106 0.00 0.00
+DEAL::  q_point=3
+DEAL::    div= -0.106 0.00 0.00
+DEAL::  q_point=4
+DEAL::    div= -1.47 0.00 0.00
+DEAL::  q_point=5
+DEAL::    div= -1.47 0.00 0.00
+DEAL::  q_point=6
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=7
+DEAL::    div= -0.394 0.00 0.00
+DEAL::i=37
+DEAL::  q_point=0
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.00 -0.106 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.00 -0.106 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.00 -1.47 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.00 -1.47 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::i=38
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.00 -0.106
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.00 -0.106
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.00 -1.47
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.00 -1.47
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::i=39
+DEAL::  q_point=0
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=1
+DEAL::    div= -0.106 0.00 0.00
+DEAL::  q_point=2
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=3
+DEAL::    div= -0.106 0.00 0.00
+DEAL::  q_point=4
+DEAL::    div= -1.47 0.00 0.00
+DEAL::  q_point=5
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=6
+DEAL::    div= -1.47 0.00 0.00
+DEAL::  q_point=7
+DEAL::    div= -0.394 0.00 0.00
+DEAL::i=40
+DEAL::  q_point=0
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.00 -0.106 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.00 -0.106 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.00 -1.47 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.00 -1.47 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::i=41
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.00 -0.106
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.00 -0.106
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.00 -1.47
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.00 -1.47
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::i=42
+DEAL::  q_point=0
+DEAL::    div= 1.47 0.00 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.106 0.00 0.00
+DEAL::  q_point=4
+DEAL::    div= 1.47 0.00 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.106 0.00 0.00
+DEAL::i=43
+DEAL::  q_point=0
+DEAL::    div= 0.00 1.47 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.106 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.00 1.47 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.106 0.00
+DEAL::i=44
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.00 1.47
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.00 0.106
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.00 1.47
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.00 0.106
+DEAL::i=45
+DEAL::  q_point=0
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.106 0.00 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.106 0.00 0.00
+DEAL::  q_point=4
+DEAL::    div= 1.47 0.00 0.00
+DEAL::  q_point=5
+DEAL::    div= 1.47 0.00 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.394 0.00 0.00
+DEAL::i=46
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.106 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.106 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.00 1.47 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.00 1.47 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.394 0.00
+DEAL::i=47
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.00 0.106
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.00 0.106
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.00 1.47
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.00 1.47
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.00 0.394
+DEAL::i=48
+DEAL::  q_point=0
+DEAL::    div= -0.106 0.00 0.00
+DEAL::  q_point=1
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=2
+DEAL::    div= -0.106 0.00 0.00
+DEAL::  q_point=3
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=4
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=5
+DEAL::    div= -1.47 0.00 0.00
+DEAL::  q_point=6
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=7
+DEAL::    div= -1.47 0.00 0.00
+DEAL::i=49
+DEAL::  q_point=0
+DEAL::    div= 0.00 -0.106 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.00 -0.106 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.00 -1.47 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.00 -1.47 0.00
+DEAL::i=50
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.00 -0.106
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.00 -0.106
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.00 -1.47
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.00 -1.47
+DEAL::i=51
+DEAL::  q_point=0
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=1
+DEAL::    div= 1.47 0.00 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.106 0.00 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=5
+DEAL::    div= 1.47 0.00 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.106 0.00 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.394 0.00 0.00
+DEAL::i=52
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.00 1.47 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.106 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.00 1.47 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.106 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.394 0.00
+DEAL::i=53
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.00 1.47
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.00 0.106
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.00 1.47
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.00 0.106
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.00 0.394
+DEAL::i=54
+DEAL::  q_point=0
+DEAL::    div= -0.106 0.00 0.00
+DEAL::  q_point=1
+DEAL::    div= -0.106 0.00 0.00
+DEAL::  q_point=2
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=3
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=4
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=5
+DEAL::    div= -0.394 0.00 0.00
+DEAL::  q_point=6
+DEAL::    div= -1.47 0.00 0.00
+DEAL::  q_point=7
+DEAL::    div= -1.47 0.00 0.00
+DEAL::i=55
+DEAL::  q_point=0
+DEAL::    div= 0.00 -0.106 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.00 -0.106 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.00 -0.394 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.00 -1.47 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.00 -1.47 0.00
+DEAL::i=56
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.00 -0.106
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.00 -0.106
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.00 -0.394
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.00 -1.47
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.00 -1.47
+DEAL::i=57
+DEAL::  q_point=0
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.106 0.00 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.106 0.00 0.00
+DEAL::  q_point=4
+DEAL::    div= 1.47 0.00 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=6
+DEAL::    div= 1.47 0.00 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.394 0.00 0.00
+DEAL::i=58
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.106 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.106 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.00 1.47 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.00 1.47 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.394 0.00
+DEAL::i=59
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.00 0.106
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.00 0.106
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.00 1.47
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.00 1.47
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.00 0.394
+DEAL::i=60
+DEAL::  q_point=0
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.106 0.00 0.00
+DEAL::  q_point=2
+DEAL::    div= 1.47 0.00 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.106 0.00 0.00
+DEAL::  q_point=6
+DEAL::    div= 1.47 0.00 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.394 0.00 0.00
+DEAL::i=61
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.106 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.00 1.47 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.106 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.00 1.47 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.394 0.00
+DEAL::i=62
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.00 0.106
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.00 1.47
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.00 0.106
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.00 1.47
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.00 0.394
+DEAL::i=63
+DEAL::  q_point=0
+DEAL::    div= 0.106 0.00 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.106 0.00 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=6
+DEAL::    div= 1.47 0.00 0.00
+DEAL::  q_point=7
+DEAL::    div= 1.47 0.00 0.00
+DEAL::i=64
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.106 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.106 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.00 1.47 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.00 1.47 0.00
+DEAL::i=65
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.00 0.106
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.00 0.106
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.00 1.47
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.00 1.47
+DEAL::i=66
+DEAL::  q_point=0
+DEAL::    div= 0.106 0.00 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.106 0.00 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=5
+DEAL::    div= 1.47 0.00 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=7
+DEAL::    div= 1.47 0.00 0.00
+DEAL::i=67
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.106 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.106 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.00 1.47 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.00 1.47 0.00
+DEAL::i=68
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.00 0.106
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.00 0.106
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.00 1.47
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.00 1.47
+DEAL::i=69
+DEAL::  q_point=0
+DEAL::    div= 0.106 0.00 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=3
+DEAL::    div= 1.47 0.00 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.106 0.00 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.394 0.00 0.00
+DEAL::  q_point=7
+DEAL::    div= 1.47 0.00 0.00
+DEAL::i=70
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.106 0.00
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=3
+DEAL::    div= 0.00 1.47 0.00
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.106 0.00
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.394 0.00
+DEAL::  q_point=7
+DEAL::    div= 0.00 1.47 0.00
+DEAL::i=71
+DEAL::  q_point=0
+DEAL::    div= 0.00 0.00 0.106
+DEAL::  q_point=1
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=2
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=3
+DEAL::    div= 0.00 0.00 1.47
+DEAL::  q_point=4
+DEAL::    div= 0.00 0.00 0.106
+DEAL::  q_point=5
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=6
+DEAL::    div= 0.00 0.00 0.394
+DEAL::  q_point=7
+DEAL::    div= 0.00 0.00 1.47
diff --git a/tests/deal.II/fe_values_view_24_nonsymmetric/fe_values_view_09/cmp/generic b/tests/deal.II/fe_values_view_24_nonsymmetric/fe_values_view_09/cmp/generic
new file mode 100644 (file)
index 0000000..977609c
--- /dev/null
@@ -0,0 +1,91 @@
+
+DEAL::FE=FESystem<2>[FE_Q<2>(1)-FE_RaviartThomas<2>(1)-FE_Nedelec<2>(1)]
+DEAL::component=0
+DEAL::0.807 4.36
+DEAL::0.807 5.29
+DEAL::1.31 4.07
+DEAL::1.31 5.59
+DEAL::component=1
+DEAL::18.2 352.
+DEAL::-18.2 337.
+DEAL::333. 180.
+DEAL::-295. 296.
+DEAL::component=2
+DEAL::43.2 -31.7
+DEAL::43.2 118.
+DEAL::125. -61.1
+DEAL::125. -40.8
+DEAL::component=3
+DEAL::2.58e-16 5.84
+DEAL::2.58e-16 5.84
+DEAL::3.58e-16 9.52
+DEAL::3.58e-16 9.52
+DEAL::component=4
+DEAL::5.84 -6.75
+DEAL::5.84 6.75
+DEAL::9.52 -11.0
+DEAL::9.52 11.0
+DEAL::FE=FESystem<3>[FE_Q<3>(1)-FE_RaviartThomas<3>(1)-FE_Nedelec<3>(1)]
+DEAL::component=0
+DEAL::2.37 4.73 9.46
+DEAL::2.37 4.73 9.46
+DEAL::2.37 4.73 9.46
+DEAL::2.37 4.73 9.46
+DEAL::2.37 4.73 9.46
+DEAL::2.37 4.73 9.46
+DEAL::2.37 4.73 9.46
+DEAL::2.37 4.73 9.46
+DEAL::component=1
+DEAL::8.18e-13 -275. -138.
+DEAL::-5.71e-13 -275. -138.
+DEAL::-367. -275. 4.71e+03
+DEAL::367. -275. 4.92e+03
+DEAL::-184. 4.57e+03 -138.
+DEAL::184. 4.79e+03 -138.
+DEAL::4.80e+03 4.57e+03 4.71e+03
+DEAL::-4.38e+03 4.79e+03 4.92e+03
+DEAL::component=2
+DEAL::-275. 1.30e-12 -138.
+DEAL::-275. -459. 4.80e+03
+DEAL::-275. -1.40e-12 -138.
+DEAL::-275. 459. 5.02e+03
+DEAL::4.67e+03 -91.8 -138.
+DEAL::4.67e+03 3.52e+03 4.80e+03
+DEAL::4.88e+03 91.8 -138.
+DEAL::4.88e+03 -3.09e+03 5.02e+03
+DEAL::component=3
+DEAL::-275. -138. 1.36e-12
+DEAL::-275. 4.90e+03 -367.
+DEAL::4.76e+03 -138. -184.
+DEAL::4.76e+03 4.90e+03 2.23e+03
+DEAL::-275. -138. -3.64e-12
+DEAL::-275. 5.11e+03 367.
+DEAL::4.97e+03 -138. 184.
+DEAL::4.97e+03 5.11e+03 -1.81e+03
+DEAL::component=4
+DEAL::-7.70e-16 5.60 22.4
+DEAL::-7.70e-16 5.60 22.4
+DEAL::-6.69e-16 5.60 22.4
+DEAL::-6.69e-16 5.60 22.4
+DEAL::1.58e-16 5.60 22.4
+DEAL::1.58e-16 5.60 22.4
+DEAL::9.52e-16 5.60 22.4
+DEAL::9.52e-16 5.60 22.4
+DEAL::component=5
+DEAL::5.60 5.78e-16 22.4
+DEAL::5.60 -1.39e-15 22.4
+DEAL::5.60 5.78e-16 22.4
+DEAL::5.60 -1.39e-15 22.4
+DEAL::5.60 -6.60e-16 22.4
+DEAL::5.60 3.53e-16 22.4
+DEAL::5.60 -6.60e-16 22.4
+DEAL::5.60 3.53e-16 22.4
+DEAL::component=6
+DEAL::5.60 11.2 -1.38e-15
+DEAL::5.60 11.2 1.96e-16
+DEAL::5.60 11.2 1.23e-17
+DEAL::5.60 11.2 1.50e-15
+DEAL::5.60 11.2 -1.38e-15
+DEAL::5.60 11.2 1.96e-16
+DEAL::5.60 11.2 1.23e-17
+DEAL::5.60 11.2 1.50e-15
diff --git a/tests/deal.II/fe_values_view_25_nonsymmetric.cc b/tests/deal.II/fe_values_view_25_nonsymmetric.cc
new file mode 100644 (file)
index 0000000..04fb1e2
--- /dev/null
@@ -0,0 +1,102 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 2007, 2008, 2010, 2011 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.
+//
+//----------------------------------------------------------------------
+
+
+// like _24, but for a simpler mesh for which the output has been
+// verified to be correct
+
+#include "../tests.h"
+#include <deal.II/base/logstream.h>
+#include <deal.II/base/function.h>
+#include <deal.II/base/quadrature_lib.h>
+#include <deal.II/lac/vector.h>
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/dofs/dof_handler.h>
+#include <deal.II/fe/fe_q.h>
+#include <deal.II/fe/fe_system.h>
+#include <deal.II/fe/fe_values.h>
+#include <deal.II/fe/mapping_q1.h>
+
+#include <fstream>
+
+
+
+template<int dim>
+void test (const Triangulation<dim>& tr,
+          const FiniteElement<dim>& fe)
+{
+  deallog << "FE=" << fe.get_name()
+         << std::endl;
+
+  DoFHandler<dim> dof(tr);
+  dof.distribute_dofs(fe);
+
+  Vector<double> fe_function(dof.n_dofs());
+  for (unsigned int i=0; i<dof.n_dofs(); ++i)
+    fe_function(i) = i+1;
+
+  const QGauss<dim> quadrature(2);
+  FEValues<dim> fe_values (fe, quadrature,
+                          update_values | update_gradients | update_q_points);
+  fe_values.reinit (dof.begin_active());
+
+                                  // let the FEValues object compute the
+                                  // divergences at quadrature points
+  std::vector<Tensor<1,dim> > divergences (quadrature.size());
+  FEValuesExtractors::Tensor<2> extractor(0);
+  fe_values[extractor]
+    .get_function_divergences (fe_function, divergences);
+
+                                  // now do the same "by hand"
+  std::vector<unsigned int> local_dof_indices (fe.dofs_per_cell);
+  dof.begin_active()->get_dof_indices (local_dof_indices);
+
+  for (unsigned int i=0; i<fe.dofs_per_cell; ++i)
+    {
+      deallog << "i=" << i << std::endl;
+
+      for (unsigned int q=0; q<quadrature.size(); ++q)
+       deallog << "  q_point=" << fe_values.quadrature_point(q) << std::endl
+               << "    value= " << fe_values[extractor].value (i,q) << std::endl
+               << "    div= " << fe_values[extractor].divergence (i,q) << std::endl;
+    }
+}
+
+
+
+template<int dim>
+void test_hyper_cube()
+{
+  Triangulation<dim> tr;
+  GridGenerator::hyper_cube(tr);
+
+  FESystem<dim> fe (FE_Q<dim>(1),
+                   Tensor<2,dim>::n_independent_components);
+  test(tr, fe);
+}
+
+
+int main()
+{
+  std::ofstream logfile ("fe_values_view_25_nonsymmetric/output");
+  deallog << std::setprecision (3);
+
+  deallog.attach(logfile);
+  deallog.depth_console (0);
+  deallog.threshold_double(1.e-7);
+
+  test_hyper_cube<1>();
+  test_hyper_cube<2>();
+  test_hyper_cube<3>();
+}
diff --git a/tests/deal.II/fe_values_view_25_nonsymmetric/cmp/generic b/tests/deal.II/fe_values_view_25_nonsymmetric/cmp/generic
new file mode 100644 (file)
index 0000000..440f957
--- /dev/null
@@ -0,0 +1,2026 @@
+
+DEAL::FE=FESystem<1>[FE_Q<1>(1)]
+DEAL::i=0
+DEAL::  q_point=0.211
+DEAL::    value= 0.789
+DEAL::    div= -1.00
+DEAL::  q_point=0.789
+DEAL::    value= 0.211
+DEAL::    div= -1.00
+DEAL::i=1
+DEAL::  q_point=0.211
+DEAL::    value= 0.211
+DEAL::    div= 1.00
+DEAL::  q_point=0.789
+DEAL::    value= 0.789
+DEAL::    div= 1.00
+DEAL::FE=FESystem<2>[FE_Q<2>(1)^4]
+DEAL::i=0
+DEAL::  q_point=0.211 0.211
+DEAL::    value= 0.622 0.00 0.00 0.00
+DEAL::    div= -0.789 0.00
+DEAL::  q_point=0.789 0.211
+DEAL::    value= 0.167 0.00 0.00 0.00
+DEAL::    div= -0.789 0.00
+DEAL::  q_point=0.211 0.789
+DEAL::    value= 0.167 0.00 0.00 0.00
+DEAL::    div= -0.211 0.00
+DEAL::  q_point=0.789 0.789
+DEAL::    value= 0.0447 0.00 0.00 0.00
+DEAL::    div= -0.211 0.00
+DEAL::i=1
+DEAL::  q_point=0.211 0.211
+DEAL::    value= 0.00 0.622 0.00 0.00
+DEAL::    div= 0.00 -0.789
+DEAL::  q_point=0.789 0.211
+DEAL::    value= 0.00 0.167 0.00 0.00
+DEAL::    div= 0.00 -0.789
+DEAL::  q_point=0.211 0.789
+DEAL::    value= 0.00 0.167 0.00 0.00
+DEAL::    div= 0.00 -0.211
+DEAL::  q_point=0.789 0.789
+DEAL::    value= 0.00 0.0447 0.00 0.00
+DEAL::    div= 0.00 -0.211
+DEAL::i=2
+DEAL::  q_point=0.211 0.211
+DEAL::    value= 0.00 0.00 0.622 0.00
+DEAL::    div= -0.789 0.00
+DEAL::  q_point=0.789 0.211
+DEAL::    value= 0.00 0.00 0.167 0.00
+DEAL::    div= -0.211 0.00
+DEAL::  q_point=0.211 0.789
+DEAL::    value= 0.00 0.00 0.167 0.00
+DEAL::    div= -0.789 0.00
+DEAL::  q_point=0.789 0.789
+DEAL::    value= 0.00 0.00 0.0447 0.00
+DEAL::    div= -0.211 0.00
+DEAL::i=3
+DEAL::  q_point=0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.622
+DEAL::    div= 0.00 -0.789
+DEAL::  q_point=0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.167
+DEAL::    div= 0.00 -0.211
+DEAL::  q_point=0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.167
+DEAL::    div= 0.00 -0.789
+DEAL::  q_point=0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.0447
+DEAL::    div= 0.00 -0.211
+DEAL::i=4
+DEAL::  q_point=0.211 0.211
+DEAL::    value= 0.167 0.00 0.00 0.00
+DEAL::    div= 0.789 0.00
+DEAL::  q_point=0.789 0.211
+DEAL::    value= 0.622 0.00 0.00 0.00
+DEAL::    div= 0.789 0.00
+DEAL::  q_point=0.211 0.789
+DEAL::    value= 0.0447 0.00 0.00 0.00
+DEAL::    div= 0.211 0.00
+DEAL::  q_point=0.789 0.789
+DEAL::    value= 0.167 0.00 0.00 0.00
+DEAL::    div= 0.211 0.00
+DEAL::i=5
+DEAL::  q_point=0.211 0.211
+DEAL::    value= 0.00 0.167 0.00 0.00
+DEAL::    div= 0.00 0.789
+DEAL::  q_point=0.789 0.211
+DEAL::    value= 0.00 0.622 0.00 0.00
+DEAL::    div= 0.00 0.789
+DEAL::  q_point=0.211 0.789
+DEAL::    value= 0.00 0.0447 0.00 0.00
+DEAL::    div= 0.00 0.211
+DEAL::  q_point=0.789 0.789
+DEAL::    value= 0.00 0.167 0.00 0.00
+DEAL::    div= 0.00 0.211
+DEAL::i=6
+DEAL::  q_point=0.211 0.211
+DEAL::    value= 0.00 0.00 0.167 0.00
+DEAL::    div= -0.211 0.00
+DEAL::  q_point=0.789 0.211
+DEAL::    value= 0.00 0.00 0.622 0.00
+DEAL::    div= -0.789 0.00
+DEAL::  q_point=0.211 0.789
+DEAL::    value= 0.00 0.00 0.0447 0.00
+DEAL::    div= -0.211 0.00
+DEAL::  q_point=0.789 0.789
+DEAL::    value= 0.00 0.00 0.167 0.00
+DEAL::    div= -0.789 0.00
+DEAL::i=7
+DEAL::  q_point=0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.167
+DEAL::    div= 0.00 -0.211
+DEAL::  q_point=0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.622
+DEAL::    div= 0.00 -0.789
+DEAL::  q_point=0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.0447
+DEAL::    div= 0.00 -0.211
+DEAL::  q_point=0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.167
+DEAL::    div= 0.00 -0.789
+DEAL::i=8
+DEAL::  q_point=0.211 0.211
+DEAL::    value= 0.167 0.00 0.00 0.00
+DEAL::    div= -0.211 0.00
+DEAL::  q_point=0.789 0.211
+DEAL::    value= 0.0447 0.00 0.00 0.00
+DEAL::    div= -0.211 0.00
+DEAL::  q_point=0.211 0.789
+DEAL::    value= 0.622 0.00 0.00 0.00
+DEAL::    div= -0.789 0.00
+DEAL::  q_point=0.789 0.789
+DEAL::    value= 0.167 0.00 0.00 0.00
+DEAL::    div= -0.789 0.00
+DEAL::i=9
+DEAL::  q_point=0.211 0.211
+DEAL::    value= 0.00 0.167 0.00 0.00
+DEAL::    div= 0.00 -0.211
+DEAL::  q_point=0.789 0.211
+DEAL::    value= 0.00 0.0447 0.00 0.00
+DEAL::    div= 0.00 -0.211
+DEAL::  q_point=0.211 0.789
+DEAL::    value= 0.00 0.622 0.00 0.00
+DEAL::    div= 0.00 -0.789
+DEAL::  q_point=0.789 0.789
+DEAL::    value= 0.00 0.167 0.00 0.00
+DEAL::    div= 0.00 -0.789
+DEAL::i=10
+DEAL::  q_point=0.211 0.211
+DEAL::    value= 0.00 0.00 0.167 0.00
+DEAL::    div= 0.789 0.00
+DEAL::  q_point=0.789 0.211
+DEAL::    value= 0.00 0.00 0.0447 0.00
+DEAL::    div= 0.211 0.00
+DEAL::  q_point=0.211 0.789
+DEAL::    value= 0.00 0.00 0.622 0.00
+DEAL::    div= 0.789 0.00
+DEAL::  q_point=0.789 0.789
+DEAL::    value= 0.00 0.00 0.167 0.00
+DEAL::    div= 0.211 0.00
+DEAL::i=11
+DEAL::  q_point=0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.167
+DEAL::    div= 0.00 0.789
+DEAL::  q_point=0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.0447
+DEAL::    div= 0.00 0.211
+DEAL::  q_point=0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.622
+DEAL::    div= 0.00 0.789
+DEAL::  q_point=0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.167
+DEAL::    div= 0.00 0.211
+DEAL::i=12
+DEAL::  q_point=0.211 0.211
+DEAL::    value= 0.0447 0.00 0.00 0.00
+DEAL::    div= 0.211 0.00
+DEAL::  q_point=0.789 0.211
+DEAL::    value= 0.167 0.00 0.00 0.00
+DEAL::    div= 0.211 0.00
+DEAL::  q_point=0.211 0.789
+DEAL::    value= 0.167 0.00 0.00 0.00
+DEAL::    div= 0.789 0.00
+DEAL::  q_point=0.789 0.789
+DEAL::    value= 0.622 0.00 0.00 0.00
+DEAL::    div= 0.789 0.00
+DEAL::i=13
+DEAL::  q_point=0.211 0.211
+DEAL::    value= 0.00 0.0447 0.00 0.00
+DEAL::    div= 0.00 0.211
+DEAL::  q_point=0.789 0.211
+DEAL::    value= 0.00 0.167 0.00 0.00
+DEAL::    div= 0.00 0.211
+DEAL::  q_point=0.211 0.789
+DEAL::    value= 0.00 0.167 0.00 0.00
+DEAL::    div= 0.00 0.789
+DEAL::  q_point=0.789 0.789
+DEAL::    value= 0.00 0.622 0.00 0.00
+DEAL::    div= 0.00 0.789
+DEAL::i=14
+DEAL::  q_point=0.211 0.211
+DEAL::    value= 0.00 0.00 0.0447 0.00
+DEAL::    div= 0.211 0.00
+DEAL::  q_point=0.789 0.211
+DEAL::    value= 0.00 0.00 0.167 0.00
+DEAL::    div= 0.789 0.00
+DEAL::  q_point=0.211 0.789
+DEAL::    value= 0.00 0.00 0.167 0.00
+DEAL::    div= 0.211 0.00
+DEAL::  q_point=0.789 0.789
+DEAL::    value= 0.00 0.00 0.622 0.00
+DEAL::    div= 0.789 0.00
+DEAL::i=15
+DEAL::  q_point=0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.0447
+DEAL::    div= 0.00 0.211
+DEAL::  q_point=0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.167
+DEAL::    div= 0.00 0.789
+DEAL::  q_point=0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.167
+DEAL::    div= 0.00 0.211
+DEAL::  q_point=0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.622
+DEAL::    div= 0.00 0.789
+DEAL::FE=FESystem<3>[FE_Q<3>(1)^9]
+DEAL::i=0
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.491 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.622 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.622 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.0447 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00944 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.0447 0.00 0.00
+DEAL::i=1
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.491 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.622 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.622 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.0447 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00944 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.0447 0.00
+DEAL::i=2
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.491 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.622
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.622
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.0447
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00944 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.0447
+DEAL::i=3
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.491 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.622 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.622 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.0447 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00944 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.0447 0.00 0.00
+DEAL::i=4
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.491 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.622 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.622 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.0447 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00944 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.0447 0.00
+DEAL::i=5
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.491 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.622
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.622
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.0447
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00944 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.0447
+DEAL::i=6
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.491 0.00 0.00
+DEAL::    div= -0.622 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00
+DEAL::    div= -0.0447 0.00 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00
+DEAL::    div= -0.622 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00944 0.00 0.00
+DEAL::    div= -0.0447 0.00 0.00
+DEAL::i=7
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.491 0.00
+DEAL::    div= 0.00 -0.622 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00
+DEAL::    div= 0.00 -0.0447 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00
+DEAL::    div= 0.00 -0.622 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00944 0.00
+DEAL::    div= 0.00 -0.0447 0.00
+DEAL::i=8
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.491
+DEAL::    div= 0.00 0.00 -0.622
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352
+DEAL::    div= 0.00 0.00 -0.0447
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131
+DEAL::    div= 0.00 0.00 -0.622
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00944
+DEAL::    div= 0.00 0.00 -0.0447
+DEAL::i=9
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.622 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.491 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.622 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00944 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.0447 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.0447 0.00 0.00
+DEAL::i=10
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.622 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.491 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.622 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00944 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.0447 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.0447 0.00
+DEAL::i=11
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.622
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.491 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.622
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00944 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.0447
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.0447
+DEAL::i=12
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.491 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.622 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.622 0.00 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.0447 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00944 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.0447 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::i=13
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.491 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.622 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.622 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.0447 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00944 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.0447 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::i=14
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.491 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.622
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.622
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.0447
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00944 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.0447
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::i=15
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.491 0.00 0.00
+DEAL::    div= -0.622 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00
+DEAL::    div= -0.0447 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00
+DEAL::    div= -0.622 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00944 0.00 0.00
+DEAL::    div= -0.0447 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::i=16
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.491 0.00
+DEAL::    div= 0.00 -0.622 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00
+DEAL::    div= 0.00 -0.0447 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00
+DEAL::    div= 0.00 -0.622 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00944 0.00
+DEAL::    div= 0.00 -0.0447 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::i=17
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.491
+DEAL::    div= 0.00 0.00 -0.622
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352
+DEAL::    div= 0.00 0.00 -0.0447
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131
+DEAL::    div= 0.00 0.00 -0.622
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00944
+DEAL::    div= 0.00 0.00 -0.0447
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::i=18
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.491 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.622 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.622 0.00 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.0447 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00944 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.0447 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::i=19
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.491 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.622 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.622 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.0447 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00944 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.0447 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::i=20
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.491 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.622
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.622
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.0447
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00944 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.0447
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::i=21
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.622 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.491 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.622 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00944 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.0447 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.0447 0.00 0.00
+DEAL::i=22
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.622 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.491 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.622 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00944 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.0447 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.0447 0.00
+DEAL::i=23
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.622
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.491 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.622
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00944 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.0447
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.0447
+DEAL::i=24
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00
+DEAL::    div= -0.0447 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.491 0.00 0.00
+DEAL::    div= -0.622 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00944 0.00 0.00
+DEAL::    div= -0.0447 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00
+DEAL::    div= -0.622 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::i=25
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00
+DEAL::    div= 0.00 -0.0447 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.491 0.00
+DEAL::    div= 0.00 -0.622 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00944 0.00
+DEAL::    div= 0.00 -0.0447 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00
+DEAL::    div= 0.00 -0.622 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::i=26
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352
+DEAL::    div= 0.00 0.00 -0.0447
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.491
+DEAL::    div= 0.00 0.00 -0.622
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00944
+DEAL::    div= 0.00 0.00 -0.0447
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131
+DEAL::    div= 0.00 0.00 -0.622
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::i=27
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.622 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.491 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.622 0.00 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00944 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.0447 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.0447 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::i=28
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.622 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.491 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.622 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00944 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.0447 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.0447 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::i=29
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.622
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.491 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.622
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00944 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.0447
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.0447
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::i=30
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.622 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.491 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.622 0.00 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00944 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.0447 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.0447 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::i=31
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.622 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.491 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.622 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00944 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.0447 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.0447 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::i=32
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.622
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.491 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.622
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00944 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.0447
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.0447
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::i=33
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00
+DEAL::    div= -0.0447 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.491 0.00 0.00
+DEAL::    div= -0.622 0.00 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00944 0.00 0.00
+DEAL::    div= -0.0447 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00
+DEAL::    div= -0.622 0.00 0.00
+DEAL::i=34
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00
+DEAL::    div= 0.00 -0.0447 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.491 0.00
+DEAL::    div= 0.00 -0.622 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00944 0.00
+DEAL::    div= 0.00 -0.0447 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00
+DEAL::    div= 0.00 -0.622 0.00
+DEAL::i=35
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352
+DEAL::    div= 0.00 0.00 -0.0447
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.491
+DEAL::    div= 0.00 0.00 -0.622
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00944
+DEAL::    div= 0.00 0.00 -0.0447
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131
+DEAL::    div= 0.00 0.00 -0.622
+DEAL::i=36
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.0447 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00944 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.0447 0.00 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.491 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.622 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.622 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::i=37
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.0447 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00944 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.0447 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.491 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.622 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.622 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::i=38
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.0447
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00944 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.0447
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.491 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.622
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.622
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::i=39
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.0447 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00944 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.0447 0.00 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.491 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.622 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.622 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::i=40
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.0447 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00944 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.0447 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.491 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.622 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.622 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::i=41
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.0447
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00944 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.0447
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.491 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.622
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.622
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::i=42
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00
+DEAL::    div= 0.622 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00944 0.00 0.00
+DEAL::    div= 0.0447 0.00 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.491 0.00 0.00
+DEAL::    div= 0.622 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00
+DEAL::    div= 0.0447 0.00 0.00
+DEAL::i=43
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00
+DEAL::    div= 0.00 0.622 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00944 0.00
+DEAL::    div= 0.00 0.0447 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.491 0.00
+DEAL::    div= 0.00 0.622 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00
+DEAL::    div= 0.00 0.0447 0.00
+DEAL::i=44
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131
+DEAL::    div= 0.00 0.00 0.622
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00944
+DEAL::    div= 0.00 0.00 0.0447
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.491
+DEAL::    div= 0.00 0.00 0.622
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352
+DEAL::    div= 0.00 0.00 0.0447
+DEAL::i=45
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00944 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.0447 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.0447 0.00 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.622 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.491 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.622 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::i=46
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00944 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.0447 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.0447 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.622 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.491 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.622 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::i=47
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00944 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.0447
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.0447
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.622
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.491 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.622
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::i=48
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.0447 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00944 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.0447 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.491 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.622 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.622 0.00 0.00
+DEAL::i=49
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.0447 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00944 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.0447 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.491 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.622 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.622 0.00
+DEAL::i=50
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.0447
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00944 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.0447
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.491 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.622
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.622
+DEAL::i=51
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00
+DEAL::    div= 0.622 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00944 0.00 0.00
+DEAL::    div= 0.0447 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.491 0.00 0.00
+DEAL::    div= 0.622 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00
+DEAL::    div= 0.0447 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::i=52
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00
+DEAL::    div= 0.00 0.622 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00944 0.00
+DEAL::    div= 0.00 0.0447 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.491 0.00
+DEAL::    div= 0.00 0.622 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00
+DEAL::    div= 0.00 0.0447 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::i=53
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131
+DEAL::    div= 0.00 0.00 0.622
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00944
+DEAL::    div= 0.00 0.00 0.0447
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.491
+DEAL::    div= 0.00 0.00 0.622
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352
+DEAL::    div= 0.00 0.00 0.0447
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131
+DEAL::    div= 0.00 0.00 0.167
+DEAL::i=54
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.0447 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00944 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.0447 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.167 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.491 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.622 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= -0.622 0.00 0.00
+DEAL::i=55
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.0447 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00944 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.0447 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.167 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.491 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.622 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 -0.622 0.00
+DEAL::i=56
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.0447
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00944 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.0447
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.167
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.491 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.622
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 -0.622
+DEAL::i=57
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00944 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.0447 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.0447 0.00 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.622 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.491 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.622 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::i=58
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00944 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.0447 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.0447 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.622 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.491 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.622 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::i=59
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00944 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.0447
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.0447
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.622
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.491 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.622
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::i=60
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00944 0.00 0.00
+DEAL::    div= 0.0447 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00
+DEAL::    div= 0.622 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00
+DEAL::    div= 0.0447 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.491 0.00 0.00
+DEAL::    div= 0.622 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::i=61
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00944 0.00
+DEAL::    div= 0.00 0.0447 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00
+DEAL::    div= 0.00 0.622 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00
+DEAL::    div= 0.00 0.0447 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.491 0.00
+DEAL::    div= 0.00 0.622 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::i=62
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00944
+DEAL::    div= 0.00 0.00 0.0447
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131
+DEAL::    div= 0.00 0.00 0.622
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352
+DEAL::    div= 0.00 0.00 0.0447
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.491
+DEAL::    div= 0.00 0.00 0.622
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131
+DEAL::    div= 0.00 0.00 0.167
+DEAL::i=63
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00944 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.0447 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.0447 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.622 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.491 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.622 0.00 0.00
+DEAL::i=64
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00944 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.0447 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.0447 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.622 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.491 0.00 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.622 0.00
+DEAL::i=65
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00944 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.0447
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.0447
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.622
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.491 0.00 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.622
+DEAL::i=66
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00944 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.0447 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.0447 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.622 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.491 0.00 0.00 0.00 0.00 0.00
+DEAL::    div= 0.622 0.00 0.00
+DEAL::i=67
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00944 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.0447 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.0447 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.622 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.491 0.00 0.00 0.00 0.00
+DEAL::    div= 0.00 0.622 0.00
+DEAL::i=68
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00944 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.0447
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.0447
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.622
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.491 0.00 0.00 0.00
+DEAL::    div= 0.00 0.00 0.622
+DEAL::i=69
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00944 0.00 0.00
+DEAL::    div= 0.0447 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00
+DEAL::    div= 0.622 0.00 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00 0.00
+DEAL::    div= 0.0447 0.00 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00 0.00
+DEAL::    div= 0.167 0.00 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.491 0.00 0.00
+DEAL::    div= 0.622 0.00 0.00
+DEAL::i=70
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00944 0.00
+DEAL::    div= 0.00 0.0447 0.00
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00
+DEAL::    div= 0.00 0.622 0.00
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352 0.00
+DEAL::    div= 0.00 0.0447 0.00
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131 0.00
+DEAL::    div= 0.00 0.167 0.00
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.491 0.00
+DEAL::    div= 0.00 0.622 0.00
+DEAL::i=71
+DEAL::  q_point=0.211 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00944
+DEAL::    div= 0.00 0.00 0.0447
+DEAL::  q_point=0.789 0.211 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.211 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.789 0.789 0.211
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131
+DEAL::    div= 0.00 0.00 0.622
+DEAL::  q_point=0.211 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.0352
+DEAL::    div= 0.00 0.00 0.0447
+DEAL::  q_point=0.789 0.211 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.211 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.131
+DEAL::    div= 0.00 0.00 0.167
+DEAL::  q_point=0.789 0.789 0.789
+DEAL::    value= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.491
+DEAL::    div= 0.00 0.00 0.622
diff --git a/tests/deal.II/fe_values_view_25_nonsymmetric/fe_values_view_09/cmp/generic b/tests/deal.II/fe_values_view_25_nonsymmetric/fe_values_view_09/cmp/generic
new file mode 100644 (file)
index 0000000..977609c
--- /dev/null
@@ -0,0 +1,91 @@
+
+DEAL::FE=FESystem<2>[FE_Q<2>(1)-FE_RaviartThomas<2>(1)-FE_Nedelec<2>(1)]
+DEAL::component=0
+DEAL::0.807 4.36
+DEAL::0.807 5.29
+DEAL::1.31 4.07
+DEAL::1.31 5.59
+DEAL::component=1
+DEAL::18.2 352.
+DEAL::-18.2 337.
+DEAL::333. 180.
+DEAL::-295. 296.
+DEAL::component=2
+DEAL::43.2 -31.7
+DEAL::43.2 118.
+DEAL::125. -61.1
+DEAL::125. -40.8
+DEAL::component=3
+DEAL::2.58e-16 5.84
+DEAL::2.58e-16 5.84
+DEAL::3.58e-16 9.52
+DEAL::3.58e-16 9.52
+DEAL::component=4
+DEAL::5.84 -6.75
+DEAL::5.84 6.75
+DEAL::9.52 -11.0
+DEAL::9.52 11.0
+DEAL::FE=FESystem<3>[FE_Q<3>(1)-FE_RaviartThomas<3>(1)-FE_Nedelec<3>(1)]
+DEAL::component=0
+DEAL::2.37 4.73 9.46
+DEAL::2.37 4.73 9.46
+DEAL::2.37 4.73 9.46
+DEAL::2.37 4.73 9.46
+DEAL::2.37 4.73 9.46
+DEAL::2.37 4.73 9.46
+DEAL::2.37 4.73 9.46
+DEAL::2.37 4.73 9.46
+DEAL::component=1
+DEAL::8.18e-13 -275. -138.
+DEAL::-5.71e-13 -275. -138.
+DEAL::-367. -275. 4.71e+03
+DEAL::367. -275. 4.92e+03
+DEAL::-184. 4.57e+03 -138.
+DEAL::184. 4.79e+03 -138.
+DEAL::4.80e+03 4.57e+03 4.71e+03
+DEAL::-4.38e+03 4.79e+03 4.92e+03
+DEAL::component=2
+DEAL::-275. 1.30e-12 -138.
+DEAL::-275. -459. 4.80e+03
+DEAL::-275. -1.40e-12 -138.
+DEAL::-275. 459. 5.02e+03
+DEAL::4.67e+03 -91.8 -138.
+DEAL::4.67e+03 3.52e+03 4.80e+03
+DEAL::4.88e+03 91.8 -138.
+DEAL::4.88e+03 -3.09e+03 5.02e+03
+DEAL::component=3
+DEAL::-275. -138. 1.36e-12
+DEAL::-275. 4.90e+03 -367.
+DEAL::4.76e+03 -138. -184.
+DEAL::4.76e+03 4.90e+03 2.23e+03
+DEAL::-275. -138. -3.64e-12
+DEAL::-275. 5.11e+03 367.
+DEAL::4.97e+03 -138. 184.
+DEAL::4.97e+03 5.11e+03 -1.81e+03
+DEAL::component=4
+DEAL::-7.70e-16 5.60 22.4
+DEAL::-7.70e-16 5.60 22.4
+DEAL::-6.69e-16 5.60 22.4
+DEAL::-6.69e-16 5.60 22.4
+DEAL::1.58e-16 5.60 22.4
+DEAL::1.58e-16 5.60 22.4
+DEAL::9.52e-16 5.60 22.4
+DEAL::9.52e-16 5.60 22.4
+DEAL::component=5
+DEAL::5.60 5.78e-16 22.4
+DEAL::5.60 -1.39e-15 22.4
+DEAL::5.60 5.78e-16 22.4
+DEAL::5.60 -1.39e-15 22.4
+DEAL::5.60 -6.60e-16 22.4
+DEAL::5.60 3.53e-16 22.4
+DEAL::5.60 -6.60e-16 22.4
+DEAL::5.60 3.53e-16 22.4
+DEAL::component=6
+DEAL::5.60 11.2 -1.38e-15
+DEAL::5.60 11.2 1.96e-16
+DEAL::5.60 11.2 1.23e-17
+DEAL::5.60 11.2 1.50e-15
+DEAL::5.60 11.2 -1.38e-15
+DEAL::5.60 11.2 1.96e-16
+DEAL::5.60 11.2 1.23e-17
+DEAL::5.60 11.2 1.50e-15
index 09b0f2cc2e295c7f9994d3255455197a71200872..e3196d4155e9fc33466d105f841d52af5505cca3 100644 (file)
@@ -136,9 +136,12 @@ void do_project (const Triangulation<dim> &triangulation,
          
       if (q<=p-order_difference)
        if (error.l2_norm() > 1e-10*projection.l2_norm())
+         {
        deallog << "Projection failed with relative error "
                << error.l2_norm() / projection.l2_norm()
                << std::endl;
+       Assert (false, ExcInternalError());
+         }
     }
 }
 
index e2ea51574e86ab89c114a3fdf8c81d4a757894f4..607cb56221a41d777c7484031de357099ccd0820 100644 (file)
 
 // check that VectorTools::project works for Q elements correctly
 
+#include "../tests.h"
+#include <deal.II/base/function.h>
+#include <deal.II/base/logstream.h>
+#include <deal.II/base/quadrature_lib.h>
+#include <deal.II/lac/vector.h>
+
+#include <deal.II/grid/tria.h>
+#include <deal.II/dofs/dof_handler.h>
+#include <deal.II/lac/constraint_matrix.h>
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/grid_refinement.h>
+#include <deal.II/grid/tria_accessor.h>
+#include <deal.II/grid/tria_iterator.h>
+#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/dofs/dof_accessor.h>
+#include <deal.II/dofs/dof_tools.h>
+#include <deal.II/numerics/vector_tools.h>
+#include <deal.II/fe/fe_abf.h>
+#include <deal.II/fe/fe_dgp.h>
+#include <deal.II/fe/fe_dgp_monomial.h>
+#include <deal.II/fe/fe_dgp_nonparametric.h>
+#include <deal.II/fe/fe_dgq.h>
+#include <deal.II/fe/fe_nedelec.h>
+#include <deal.II/fe/fe_q.h>
+#include <deal.II/fe/fe_q_hierarchical.h>
+#include <deal.II/fe/fe_raviart_thomas.h>
+#include <deal.II/fe/fe_system.h>
+
+#include <fstream>
+#include <vector>
+
 char logname[] = "project_q_01/output";
 
 
-#include "project_common.h"
+template <int dim>
+class F :  public Function<dim>
+{
+  public:
+    F (const unsigned int q,
+       const unsigned int n_components)
+                    :
+                    Function<dim>(n_components),
+                    q(q)
+      {}
+
+    virtual double value (const Point<dim> &p,
+                          const unsigned int component) const
+      {
+        Assert ((component == 0) && (this->n_components == 1),
+                ExcInternalError());
+        double val = 0;
+        for (unsigned int d=0; d<dim; ++d)
+          for (unsigned int i=0; i<=q; ++i)
+            val += (d+1)*(i+1)*std::pow (p[d], 1.*i);
+        return val;
+      }
+
+
+    virtual void vector_value (const Point<dim> &p,
+                               Vector<double>   &v) const
+      {
+        for (unsigned int c=0; c<v.size(); ++c)
+          {
+            v(c) = 0;
+            for (unsigned int d=0; d<dim; ++d)
+              for (unsigned int i=0; i<=q; ++i)
+                v(c) += (d+1)*(i+1)*std::pow (p[d], 1.*i)+c;
+          }
+      }
+
+  private:
+    const unsigned int q;
+};
+
+
+template <int dim>
+void do_project (const Triangulation<dim> &triangulation,
+                 const FiniteElement<dim> &fe,
+                 const unsigned int        p,
+                 const unsigned int        order_difference)
+{
+  DoFHandler<dim>        dof_handler(triangulation);
+  dof_handler.distribute_dofs (fe);
+
+  deallog << "n_dofs=" << dof_handler.n_dofs() << std::endl;
+
+  ConstraintMatrix constraints;
+  DoFTools::make_hanging_node_constraints (dof_handler,
+                                           constraints);
+  constraints.close ();
+
+  Vector<double> projection (dof_handler.n_dofs());
+  Vector<float>  error (triangulation.n_active_cells());
+  for (unsigned int q=0; q<=p+2-order_difference; ++q)
+    {
+                                       // project the function
+      VectorTools::project (dof_handler,
+                            constraints,
+                            QGauss<dim>(p+2),
+                            F<dim> (q, fe.n_components()),
+                            projection);
+                                       // just to make sure it doesn't get
+                                       // forgotten: handle hanging node
+                                       // constraints
+      constraints.distribute (projection);
+
+                                       // then compute the interpolation error
+      VectorTools::integrate_difference (dof_handler,
+                                         projection,
+                                         F<dim> (q, fe.n_components()),
+                                         error,
+                                         QGauss<dim>(std::max(p,q)+1),
+                                         VectorTools::L2_norm);
+      deallog << fe.get_name() << ", P_" << q
+              << ", rel. error=" << error.l2_norm() / projection.l2_norm()
+              << std::endl;
+
+      if (q<=p-order_difference)
+        if (error.l2_norm() > 1e-10*projection.l2_norm())
+          {
+        deallog << "Projection failed with relative error "
+                << error.l2_norm() / projection.l2_norm()
+                << std::endl;
+        Assert (false, ExcInternalError());
+          }
+    }
+}
+
+
+
+// check the given element of polynomial order p. the last parameter, if
+// given, denotes a gap in convergence order; for example, the Nedelec element
+// of polynomial degree p has normal components of degree p-1 and therefore
+// can only represent polynomials of degree p-1 exactly. the gap is then 1.
+template <int dim>
+void test_no_hanging_nodes (const FiniteElement<dim> &fe,
+                            const unsigned int        p,
+                            const unsigned int        order_difference = 0)
+{
+  Triangulation<dim>     triangulation;
+  GridGenerator::hyper_cube (triangulation);
+  triangulation.refine_global (3);
+
+  do_project (triangulation, fe, p, order_difference);
+}
+
 
 
 template <int dim>
@@ -26,3 +168,18 @@ void test ()
   for (unsigned int p=1; p<6-dim; ++p)
     test_no_hanging_nodes (FE_Q<dim>(p), p);
 }
+
+
+int main ()
+{
+  std::ofstream logfile(logname);
+  deallog << std::setprecision (3);
+
+  deallog.attach(logfile);
+  deallog.depth_console(0);
+  deallog.threshold_double(1.e-10);
+
+  test<1>();
+  test<2>();
+  test<3>();
+}

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.