merge from mainline

git-svn-id: https://svn.dealii.org/branches/branch_unify_linear_algebra@29971 0785d39b-7218-0410-832d-ea1e28bc413d

diff --cc deal.II/cmake/configure/configure_1_lapack.cmake
index 6445f6b1f1e6a5771a0663ae26af2bb3a5ba9889,6445f6b1f1e6a5771a0663ae26af2bb3a5ba9889..5dbd35c580cd32f108fb6a440f98c19d62653f30
--- 1/deal.II/cmake/configure/configure_1_lapack.cmake
--- 2/deal.II/cmake/configure/configure_1_lapack.cmake
+++ b/deal.II/cmake/configure/configure_1_lapack.cmake
@@@ -52,15 -52,15 +52,6 @@@ SET(DEAL_II_LAPACK_FUNCTION
  
  MACRO(CHECK_FOR_LAPACK_FUNCTIONS)
    SET(CMAKE_REQUIRED_LIBRARIES ${LAPACK_LIBRARIES})
--  #
--  # For some static lapack versions it is necessary to link with -lm.
--  # So link with it as well if -lm is already present in our link
--  # interface:
--  #
--  IF(NOT m_lib MATCHES "-NOTFOUND")
--    LIST(APPEND CMAKE_REQUIRED_LIBRARIES ${m_lib})
--  ENDIF()
--
    ADD_FLAGS(CMAKE_REQUIRED_FLAGS "${LAPACK_LINKER_FLAGS}")
  
    FOREACH(_func ${DEAL_II_LAPACK_FUNCTIONS})
@@@ -80,6 -80,6 +71,14 @@@ MACRO(RESET_LAPACK_FUNCTIONS_CACHE
    ENDFOREACH()
  ENDMACRO()
  
++#
++# Resolve a cache invalidation problem by searching for this flag
++# unconditionally. It is used in FEATURE_LAPACK_CONFIGURE_EXTERNAL
++# depending on cached variables.
++#
++ENABLE_IF_LINKS(_dummy "-lgfortran")
++ENABLE_IF_LINKS(_dummy "-lquadmath")
++
  
  MACRO(FEATURE_LAPACK_CONFIGURE_EXTERNAL)
    #
@@@ -100,10 -100,10 +99,25 @@@
      LIST(APPEND LAPACK_LIBRARIES ${BLAS_LIBRARIES})
    ENDIF()
  
++
++  #
++  # Well, in case of static archives we have to manually pick up the
++  # complete link interface. *sigh*
++  #
++  # Do this unconditionally for the most common case:
++  # TODO: Non-GNU setups...
++  #
++  IF(NOT m_lib MATCHES "-NOTFOUND")
++    LIST(APPEND LAPACK_LIBRARIES ${m_lib})
++  ENDIF()
++
++  ENABLE_IF_LINKS(LAPACK_LINKER_FLAGS "-lgfortran")
++  ENABLE_IF_LINKS(LAPACK_LINKER_FLAGS "-lquadmath")
++
++
    ADD_FLAGS(CMAKE_SHARED_LINKER_FLAGS "${LAPACK_LINKER_FLAGS}")
    LIST(APPEND DEAL_II_EXTERNAL_LIBRARIES ${LAPACK_LIBRARIES})
  
--
    CHECK_FOR_LAPACK_FUNCTIONS()
  ENDMACRO()
  

diff --cc deal.II/cmake/configure/configure_boost.cmake
index 8e28bf83e58d6b96f179836e8b7be8c970137e92,8e28bf83e58d6b96f179836e8b7be8c970137e92..23bd875f552663676414653eed4b83d26c9823ee
--- 1/deal.II/cmake/configure/configure_boost.cmake
--- 2/deal.II/cmake/configure/configure_boost.cmake
+++ b/deal.II/cmake/configure/configure_boost.cmake
@@@ -35,13 -35,13 +35,30 @@@ SET(DEAL_II_WITH_BOOST ON # Always true
  
  
  MACRO(FEATURE_BOOST_FIND_EXTERNAL var)
--  #
--  # We require at least version 1.44
--  #
++
    IF(DEAL_II_WITH_THREADS)
--    FIND_PACKAGE(Boost 1.44 COMPONENTS serialization system thread)
++    SET(_boost_components serialization system thread)
    ELSE()
--    FIND_PACKAGE(Boost 1.44 COMPONENTS serialization system)
++    SET(_boost_components serialization system)
++  ENDIF()
++
++  #
++  # Prefer static libs if BUILD_SHARED_LIBS=OFF:
++  #
++  IF(NOT BUILD_SHARED_LIBS)
++    SET(Boost_USE_STATIC_LIBS TRUE)
++  ENDIF()
++
++  FIND_PACKAGE(Boost 1.44 COMPONENTS ${_boost_components})
++
++  #
++  # Fall back to dynamic libraries if no static libraries could be found:
++  #
++  IF( Boost_USE_STATIC_LIBS AND
++      (NOT Boost_SERIALIZATION_FOUND OR NOT Boost_SYSTEM_FOUND)
++    )
++    SET(Boost_USE_STATIC_LIBS FALSE)
++    FIND_PACKAGE(Boost 1.44 COMPONENTS ${_boost_components})
    ENDIF()
  
    IF( Boost_SERIALIZATION_FOUND AND

diff --cc deal.II/cmake/modules/FindPETSC.cmake
index b3fddfac03cb2c6faed2e68512f6c49d9f14dccd,b3fddfac03cb2c6faed2e68512f6c49d9f14dccd..7e5f4480543542b4acbfab9ec121db9934bfa659
--- 1/deal.II/cmake/modules/FindPETSC.cmake
--- 2/deal.II/cmake/modules/FindPETSC.cmake
+++ b/deal.II/cmake/modules/FindPETSC.cmake
@@@ -52,8 -52,8 +52,8 @@@ FIND_LIBRARY(PETSC_LIBRAR
  FIND_FILE(PETSC_PETSCVARIABLES
    NAMES petscvariables
    HINTS
--    ${PETSC_DIR}
      ${PETSC_DIR}/${PETSC_ARCH}
++    ${PETSC_DIR}
    PATH_SUFFIXES conf
    )
  
@@@ -75,7 -75,7 +75,7 @@@ IF(NOT PETSC_PETSCVARIABLES MATCHES "-N
        IF(NOT _token MATCHES "(petsc|stdc\\+\\+|gcc_s)")
          FIND_LIBRARY(PETSC_LIBRARY_${_token}
            NAMES ${_token}
--          HINTS ${_hintes}
++          HINTS ${_hints}
            )
          IF(NOT PETSC_LIBRARY_${_token} MATCHES "-NOTFOUND")
            LIST(APPEND _petsc_libraries ${PETSC_LIBRARY_${_token}})

diff --cc deal.II/cmake/setup_deal_ii.cmake
index 5b82eedd223bbf758b4a95081d6fba4d974ca775,5b82eedd223bbf758b4a95081d6fba4d974ca775..d517590f5de03d749aee7d089ac6809e6a8ad0d2
--- 1/deal.II/cmake/setup_deal_ii.cmake
--- 2/deal.II/cmake/setup_deal_ii.cmake
+++ b/deal.II/cmake/setup_deal_ii.cmake
@@@ -147,7 -147,7 +147,7 @@@ ENDIF(
  
  ###########################################################################
  #                                                                         #
--#  Cleanup and Includes that have to happen after the call to PROJECT():  #
++#    Cleanup and setup that has to happen after the call to PROJECT():    #
  #                                                                         #
  ###########################################################################
  
@@@ -164,6 -164,6 +164,14 @@@ FOREACH(_build ${DEAL_II_BUILD_TYPES}
      )
  ENDFOREACH()
  
++#
++# Invert the search order for libraries when BUILD_SHARED_LIBS=OFF to
++# prefer static archives instead of shared libraries:
++#
++IF(NOT BUILD_SHARED_LIBS)
++  # TODO: Does this work on a Windows or CYGWIN target?
++  LIST(REVERSE CMAKE_FIND_LIBRARY_SUFFIXES)
++ENDIF()
  
  #
  # Cross compilation stuff:

diff --cc deal.II/examples/step-5/step-5.cc
index 18dfd8036e532368881717e6900296da3685ee32,18dfd8036e532368881717e6900296da3685ee32..6e4a5c34ee62a512554b146bc06ff4faa9113bc1
--- 1/deal.II/examples/step-5/step-5.cc
--- 2/deal.II/examples/step-5/step-5.cc
+++ b/deal.II/examples/step-5/step-5.cc
@@@ -533,7 -533,7 +533,7 @@@ void Step5<dim>::output_results (const 
  
  // The second to last thing in this program is the definition of the
  // <code>run()</code> function. In contrast to the previous programs, we will
--// compute on a sequence of meshes that after each iteration is globall
++// compute on a sequence of meshes that after each iteration is globally
  // refined. The function therefore consists of a loop over 6 cycles. In each
  // cycle, we first print the cycle number, and then have to decide what to do
  // with the mesh. If this is not the first cycle, we simply refine the

diff --cc deal.II/examples/step-7/step-7.cc
index 85393fd73352df14af99825a457c1eea5ee41745,85393fd73352df14af99825a457c1eea5ee41745..3e3eaed078a57f2a0e2cf86eb7defea10c997dff
--- 1/deal.II/examples/step-7/step-7.cc
--- 2/deal.II/examples/step-7/step-7.cc
+++ b/deal.II/examples/step-7/step-7.cc
@@@ -233,7 -233,7 +233,7 @@@ namespace Step
    // share the same connotation that points have and are only objects in a
    // more abstract space than the one spanned by the coordinate
    // directions. (In fact, gradients live in `reciprocal' space, since the
--  // dimension of their components is not that of a length, but one over
++  // dimension of their components is not that of a length, but of one over
    // length).
    template <int dim>
    Tensor<1,dim> Solution<dim>::gradient (const Point<dim>   &p,

diff --cc deal.II/examples/step-9/step-9.cc
index b5589551e5764520e174c348736369414a804de5,b5589551e5764520e174c348736369414a804de5..503e65e786d26c43f211bf8895ea7ed596989d06
--- 1/deal.II/examples/step-9/step-9.cc
--- 2/deal.II/examples/step-9/step-9.cc
+++ b/deal.II/examples/step-9/step-9.cc
@@@ -142,7 -142,7 +142,7 @@@ namespace Step
    // @sect3{Equation data declaration}
  
    // Next we declare a class that describes the advection field. This, of
--  // course, is a vector field with as many compents as there are space
++  // course, is a vector field with as many components as there are space
    // dimensions. One could now use a class derived from the
    // <code>Function</code> base class, as we have done for boundary values and
    // coefficients in previous examples, but there is another possibility in
@@@ -391,7 -391,7 +391,7 @@@
    // rather than the common vectors of doubles, as the additional accuracy is
    // not necessary for estimated values.
    //
--  // In addition to these two functions, the class declares to exceptions
++  // In addition to these two functions, the class declares two exceptions
    // which are raised when a cell has no neighbors in each of the space
    // directions (in which case the matrix described in the introduction would
    // be singular and can't be inverted), while the other one is used in the
@@@ -510,11 -510,11 +510,11 @@@
      // configured for multi-threading, then the number of CPUs is set to one.)
      // However, sometimes there might be reasons to use another value. For
      // example, you might want to use less processors than there are in your
--    // system in order not to use too many computational ressources. On the
++    // system in order not to use too many computational resources. On the
      // other hand, if there are several jobs running on a computer and you
      // want to get a higher percentage of CPU time, it might be worth to start
      // more threads than there are CPUs, as most operating systems assign
--    // roughly the same CPU ressources to all threads presently running. For
++    // roughly the same CPU resources to all threads presently running. For
      // this reason, the <code>MultithreadInfo</code> class contains a
      // read-write variable <code>n_default_threads</code> which is set to
      // <code>n_cpus</code> by default, but can be set to another value. This
@@@ -576,7 -576,7 +576,7 @@@
      // differ (<code>begin_active</code> returns an
      // <code>active_iterator</code>, while <code>end</code> returns a
      // <code>raw_iterator</code>), and in this case the C++ language requires
--    // us to specify the template type explicitely. For brevity, we first
++    // us to specify the template type explicitly. For brevity, we first
      // typedef this data type to an alias.
      typedef typename DoFHandler<dim>::active_cell_iterator active_cell_iterator;
      std::vector<std::pair<active_cell_iterator,active_cell_iterator> >
@@@ -611,8 -611,8 +611,8 @@@
      // When all the threads are running, the only thing we have to do is wait
      // for them to finish. This is necessary of course, as we can't proceed
      // with our tasks before the matrix and right hand side are
--    // assemblesd. Waiting for all the threads to finish can be done using the
--    // <code>joint_all</code> function in the <code>ThreadGroup</code>
++    // assembled. Waiting for all the threads to finish can be done using the
++    // <code>join_all</code> function in the <code>ThreadGroup</code>
      // container, which just calls <code>join</code> on each of the thread
      // objects it stores.
      //
@@@ -621,7 -621,7 +621,7 @@@
      threads.join_all ();
  
  
--    // After the matrix has been assembled in parallel, we stil have to
++    // After the matrix has been assembled in parallel, we still have to
      // eliminate hanging node constraints. This is something that can't be
      // done on each of the threads separately, so we have to do it now.
      hanging_node_constraints.condense (system_matrix);
@@@ -816,7 -816,7 +816,7 @@@
          //
          // On the other hand, we would now like to write the local
          // contributions to the global system of equations into the global
--        // objects. This needs some kind of synchronisation, as if we would
++        // objects. This needs some kind of synchronization, as if we would
          // not take care of the fact that multiple threads write into the
          // matrix at the same time, we might be surprised that one threads
          // reads data from the matrix that another thread is presently

diff --cc deal.II/include/deal.II/base/tensor.h
index bb5a41193fa1b4407f1e30645f7edc1e67cdfe96,bb1e0ee70449bff0f930524bf460fbef9694b0f1..052d81aaac03614e080c140b49e48b36487f9bbf
--- 1/deal.II/include/deal.II/base/tensor.h
--- 2/deal.II/include/deal.II/base/tensor.h
+++ b/deal.II/include/deal.II/base/tensor.h
@@@ -405,7 -415,29 +415,33 @@@ Tensor<rank_,dim,Number>::operator[] (c
    return subtensor[i];
  }
  
 -    const unsigned int outer_ind = indices[rank_-1];
+ template <int rank_, int dim, typename Number>
+ inline
+ Number 
+ Tensor<rank_,dim,Number>::operator[] (const TableIndices<rank_> & indices) const
+ {
 -       indices1[i] = indices[i];
 -      return (subtensor[outer_ind])[indices1];        
++    const unsigned int inner_ind = indices[0];
++    Assert (inner_ind<dim, ExcIndexRange(inner_ind, 0, dim));
++    
+       TableIndices<rank_-1> indices1;
+       for (unsigned int i = 0; i < rank_-1;i++)
++       indices1[i] = indices[i+1];
++      return (subtensor[inner_ind])[indices1];        
+ }
  
 -    const unsigned int outer_ind = indices[rank_-1];
+ template <int rank_, int dim, typename Number>
+ inline
+ Number & 
+ Tensor<rank_,dim,Number>::operator[] (const TableIndices<rank_> & indices)
+ {
 -       indices1[i] = indices[i];
 -      return (subtensor[outer_ind])[indices1];        
++    const unsigned int inner_ind = indices[0];
++    Assert (inner_ind<dim, ExcIndexRange(inner_ind, 0, dim));
++    
+       TableIndices<rank_-1> indices1;
+       for (unsigned int i = 0; i < rank_-1;i++)
++       indices1[i] = indices[i+1];
++      return (subtensor[inner_ind])[indices1];                
+ }
  
  template <int rank_, int dim, typename Number>
  inline

diff --cc deal.II/include/deal.II/fe/fe_values.h
index 4da0529a08abf0bf02ab692b82ccc8089fa2ceca,b9da03e4df24d23827509b67dd9fd7d0d3001b31..e3edf8f049e5e25d564ac85014bdb0f8c77025e1
--- 1/deal.II/include/deal.II/fe/fe_values.h
--- 2/deal.II/include/deal.II/fe/fe_values.h
+++ b/deal.II/include/deal.II/fe/fe_values.h
@@@ -4619,6 -4619,114 +4619,114 @@@ namespace FEValuesView
        }
    }
  
 -        return_value[indices] = fe_values.shape_values(snc,q_point);//last index first [jj][ii]
+   template <int dim, int spacedim>
+   inline
+   typename Tensor<2, dim, spacedim>::value_type
+   Tensor<2, dim, spacedim>::value (const unsigned int shape_function,
+                                    const unsigned int q_point) const
+   {
+     typedef FEValuesBase<dim,spacedim> FVB;
+     Assert (shape_function < fe_values.fe->dofs_per_cell,
+             ExcIndexRange (shape_function, 0, fe_values.fe->dofs_per_cell));
+     Assert (fe_values.update_flags & update_values,
+             typename FVB::ExcAccessToUninitializedField());
+ 
+     // similar to the vector case where we
+     // have more then one index and we need
+     // to convert between unrolled and
+     // component indexing for tensors
+     const int snc
+       = shape_function_data[shape_function].single_nonzero_component;
+ 
+     if (snc == -2)
+       {
+         // shape function is zero for the
+         // selected components
+         return value_type();
+ 
+       }
+     else if (snc != -1)
+       {
+         value_type return_value;
+         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);
 -              = fe_values.shape_values(shape_function_data[shape_function].row_index[d],q_point);//last index first [jj][ii]
++        return_value[indices] = fe_values.shape_values(snc,q_point);
+         return return_value;
+       }
+     else
+       {
+         value_type return_value;
+         for (unsigned int d = 0; d < dim*dim; ++d)
+           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);
+           }
+         return return_value;
+       }
+   }
+ 
+ 
+   template <int dim, int spacedim>
+   inline
+   typename Tensor<2, dim, spacedim>::divergence_type
+   Tensor<2, dim, spacedim>::divergence(const unsigned int shape_function,
+                                        const unsigned int q_point) const
+   {
+     typedef FEValuesBase<dim,spacedim> FVB;
+     Assert (shape_function < fe_values.fe->dofs_per_cell,
+             ExcIndexRange (shape_function, 0, fe_values.fe->dofs_per_cell));
+     Assert (fe_values.update_flags & update_gradients,
+             typename FVB::ExcAccessToUninitializedField());
+ 
+     const int snc = shape_function_data[shape_function].single_nonzero_component;
+ 
+     if (snc == -2)
+       {
+         // shape function is zero for the
+         // selected components
+         return divergence_type();
+       }
+     else if (snc != -1)
+       {
+         // we have a single non-zero component
+         // when the tensor is
+         // represented in unrolled form.
+         //
+         // the divergence of a second-order tensor
+         // is a first order tensor.
+         //
+         // assume the second-order tensor is
+         // A with components A_{ij}.  
+         // divergence as:
+         // b_j := \dfrac{\partial phi_{ij}}{\partial x_i}.
+         //
+         // Now, we know the nonzero component
+         // in unrolled form: it is indicated
+         // by 'snc'. we can figure out which
+         // tensor components belong to this:
+         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);
+         const unsigned int ii = indices[0];
+         const unsigned int jj = indices[1];
+ 
+         const dealii::Tensor<1, spacedim> phi_grad = fe_values.shape_gradients[snc][q_point];
+ 
+         divergence_type return_value;
+         return_value[jj] = phi_grad[ii];
+ 
+         return return_value;
+ 
+       }
+     else
+       {
+         Assert (false, ExcNotImplemented());
+         divergence_type return_value;
+         return return_value;
+       }
+   }
  }
  
  

diff --cc deal.II/include/deal.II/grid/grid_generator.h
Simple merge

diff --cc deal.II/source/fe/fe_values.cc
index 7957a09afcc0154827ad8643fe0e745c12651c4a,23b1b8c22617d003c4ac80f4ddde090c9a5088e7..1a54152551fc137abf52a7c8d2ec193e4c9a9d16
--- 1/deal.II/source/fe/fe_values.cc
--- 2/deal.II/source/fe/fe_values.cc
+++ b/deal.II/source/fe/fe_values.cc
@@@ -1173,9 -1173,7 +1173,7 @@@ namespace FEValuesView
  
                    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][i][j] += value * *shape_value_ptr++;
 -                    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;
@@@ -1189,7 -1185,7 +1185,7 @@@
                        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][i][j] += value * *shape_value_ptr++;
 -                        values[q_point][indices] += value * *shape_value_ptr++;//last index first [j][i]
++                        values[q_point][indices] += value * *shape_value_ptr++;
                      }
              }
          }

diff --cc deal.II/source/lac/sparse_direct.cc
index c8c91b628c525bd18b736e4f32640752d3fb9dcf,c8c91b628c525bd18b736e4f32640752d3fb9dcf..d02a9c1ce2369ae983f91c9a8a4a8e8e553a9c8b
--- 1/deal.II/source/lac/sparse_direct.cc
--- 2/deal.II/source/lac/sparse_direct.cc
+++ b/deal.II/source/lac/sparse_direct.cc
@@@ -379,7 -379,7 +379,7 @@@ void SparseDirectUMFPACK::factorize (co
  
  
  void
--SparseDirectUMFPACK::solve (Vector<double> &) const
++SparseDirectUMFPACK::solve (Vector<double> &, bool) const
  {
    AssertThrow(false, ExcMessage("To call this function you need UMFPACK, but configured deal.II without passing the necessary switch to 'cmake'. Please consult the installation instructions in doc/readme.html."));
  }
@@@ -387,7 -387,7 +387,7 @@@
  
  
  void
--SparseDirectUMFPACK::solve (BlockVector<double> &) const
++SparseDirectUMFPACK::solve (BlockVector<double> &, bool) const
  {
    AssertThrow(false, ExcMessage("To call this function you need UMFPACK, but configured deal.II without passing the necessary switch to 'cmake'. Please consult the installation instructions in doc/readme.html."));
  }
@@@ -396,7 -396,7 +396,8 @@@
  template <class Matrix>
  void
  SparseDirectUMFPACK::solve (const Matrix &,
--                            Vector<double> &)
++                            Vector<double> &,
++                            bool)
  {
    AssertThrow(false, ExcMessage("To call this function you need UMFPACK, but configured deal.II without passing the necessary switch to 'cmake'. Please consult the installation instructions in doc/readme.html."));
  }
@@@ -406,7 -406,7 +407,8 @@@
  template <class Matrix>
  void
  SparseDirectUMFPACK::solve (const Matrix &,
--                            BlockVector<double> &)
++                            BlockVector<double> &,
++                            bool)
  {
    AssertThrow(false, ExcMessage("To call this function you need UMFPACK, but configured deal.II without passing the necessary switch to 'cmake'. Please consult the installation instructions in doc/readme.html."));
  }

diff --cc tests/base/tensor_25.cc
index 0000000000000000000000000000000000000000,0000000000000000000000000000000000000000..82dea93c514ed673aa8fe49a368fd4882a3296cc
new file mode 100644 (file)
--- /dev/null
--- /dev/null
+++ b/tests/base/tensor_25.cc
@@@ -1,0 -1,0 +1,64 @@@
++//----------------------------  tensor_25.cc  ---------------------------
++//    $Id$
++//    Version: $Name$
++//
++//    Copyright (C) 2006, 2008, 2009, 2013 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.
++//
++//----------------------------  tensor_25.cc  ---------------------------
++
++// check Tensor<2,dim>::operator[] (TableIndices)
++
++#include "../tests.h"
++#include <deal.II/base/tensor.h>
++#include <deal.II/base/logstream.h>
++#include <deal.II/lac/vector.h>
++#include <fstream>
++#include <iomanip>
++
++
++template <int dim>
++void check ()
++{
++  typedef Tensor<3,dim> T;
++  T t;
++  for (unsigned int i=0; i<T::n_independent_components; ++i)
++    t[T::unrolled_to_component_indices (i)] = (i+1)*(i+2);
++
++  for (unsigned int i=0; i<dim; ++i)
++    for (unsigned int j=0; j<dim; ++j)
++      for (unsigned int k=0; k<dim; ++k)
++      {
++        Assert (TableIndices<3>(i,j,k)
++                ==
++                T::unrolled_to_component_indices (i*dim*dim+j*dim+k),
++                ExcInternalError());
++        Assert (T::component_to_unrolled_index(TableIndices<3>(i,j,k))
++                ==
++                i*dim*dim+j*dim+k,
++                ExcInternalError());
++        Assert (t[TableIndices<3>(i,j,k)] == t[i][j][k],
++                ExcInternalError());
++      }
++  deallog << "OK" << std::endl;
++}
++
++
++int main ()
++{
++  std::ofstream logfile("tensor_25/output");
++  deallog << std::setprecision(3);
++  deallog.attach(logfile);
++  deallog.depth_console(0);
++  deallog.threshold_double(1.e-10);
++
++  check<1> ();
++  check<2> ();
++  check<3> ();
++
++  deallog << "OK" << std::endl;
++}

diff --cc tests/base/tensor_25/cmp/generic
index 0000000000000000000000000000000000000000,0000000000000000000000000000000000000000..5f42bb230f11ca6753755403100d254e31356269
new file mode 100644 (file)
--- /dev/null
--- /dev/null
+++ b/tests/base/tensor_25/cmp/generic
@@@ -1,0 -1,0 +1,5 @@@
++
++DEAL::OK
++DEAL::OK
++DEAL::OK
++DEAL::OK

diff --cc tests/base/work_stream_03.cc
index eb6dc88d5c8019ca1e50c889ac2745fbb2de7b61,cea42e56fd6000f3e6b4384f49c2e7b8f1ac1d96..eec7d029922d000ef744cf866cd15e266d53bdd6
--- 1/tests/base/work_stream_03.cc
--- 2/tests/base/work_stream_03.cc
+++ b/tests/base/work_stream_03.cc
@@@ -11,7 -11,7 +11,6 @@@
  //
  //----------------------------  work_stream_03.cc  ---------------------------
  
--
  // Moritz originally implemented thread local scratch objects for
  // WorkStream in r24748 but it led to failures in the testsuite. what
  // exactly went on was a mystery and this test is a first step in
@@@ -19,147 -19,151 +18,143 @@@
  // the failing tests (deal.II/project_q_01) multiple times and
  // verifying that it indeed works
  
- //#include "../tests.h"
- #include <deal.II/base/function.h>
+ #include "../tests.h"
  #include <deal.II/base/logstream.h>
  #include <deal.II/base/quadrature_lib.h>
+ #include <deal.II/base/work_stream.h>
  #include <deal.II/lac/vector.h>
++#include <deal.II/base/parallel.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 <deal.II/fe/fe_nothing.h>
+ #include <deal.II/fe/fe_values.h>
  
  #include <fstream>
  #include <vector>
  
- using namespace dealii;
- 
- 
  char logname[] = "work_stream_03/output";
  
--
--
--template <int dim>
- class F :  public Function<dim>
 -double value (const Point<dim> &p)
 -{
 -  double val = 0;
 -  for (unsigned int d=0; d<dim; ++d)
 -    for (unsigned int i=0; i<=1; ++i)
 -      val += std::pow (p[d], 1.*i);
 -  return val;
 -}
++template<int dim>
++  double
++  value(const Point<dim> &p)
++  {
++    double val = 0;
++    for (unsigned int d = 0; d < dim; ++d)
++      for (unsigned int i = 0; i <= 1; ++i)
++        val += std::pow(p[d], 1. * i);
++    return val;
++  }
+ 
+ namespace
  {
-   public:
-     F (const unsigned int q)
-                   :
-                   q(q)
-       {}
 -  template <int dim>
 -  struct Scratch
 -  {
 -    Scratch (const FiniteElement<dim> &fe,
 -             const Quadrature<dim> &quadrature)
 -      :
 -      fe_collection (fe),
 -      quadrature_collection (quadrature),
 -      x_fe_values (fe_collection,
 -                   quadrature_collection,
 -                   update_q_points),
 -      rhs_values(quadrature_collection.size())
 -    {}
--
-     virtual double value (const Point<dim> &p,
-                         const unsigned int component) const
 -    Scratch (const Scratch &data)
 -      :
 -      fe_collection (data.fe_collection),
 -      quadrature_collection (data.quadrature_collection),
 -      x_fe_values (fe_collection,
 -                   quadrature_collection,
 -                   update_q_points),
 -      rhs_values (data.rhs_values)
 -    {}
 -
 -    const FiniteElement<dim>      &fe_collection;
 -    const Quadrature<dim>                &quadrature_collection;
 -
 -    FEValues<dim> x_fe_values;
 -
 -    std::vector<double>                  rhs_values;
 -  };
++  template<int dim>
++    struct Scratch
++    {
++      Scratch(const FiniteElement<dim> &fe, const Quadrature<dim> &quadrature) :
++          fe_collection(fe), quadrature_collection(quadrature), x_fe_values(
++              fe_collection, quadrature_collection, update_q_points), rhs_values(
++              quadrature_collection.size())
++      {
++      }
++
++      Scratch(const Scratch &data) :
++          fe_collection(data.fe_collection), quadrature_collection(
++              data.quadrature_collection), x_fe_values(fe_collection,
++              quadrature_collection, update_q_points), rhs_values(
++              data.rhs_values)
 +      {
-       Assert (component == 0, 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;
 +      }
  
-   private:
-     const unsigned int q;
- };
++      const FiniteElement<dim> &fe_collection;
++      const Quadrature<dim> &quadrature_collection;
 +
++      FEValues<dim> x_fe_values;
 +
- template <int dim>
- void do_project (const Triangulation<dim> &triangulation,
-                const unsigned int        p)
++      std::vector<double> rhs_values;
++    };
+ 
+   struct CopyData
+   {
 -    Vector<double>    cell_rhs;
++    std::vector<double> cell_rhs;
+   };
+ }
+ 
 -template<int dim, typename CellIterator>
++void
++zero_subrange(const unsigned int begin, const unsigned int end,
++    std::vector<double> &dst)
 +{
-   std::cout << "Start: " << __PRETTY_FUNCTION__ << ' ' << p << std::endl;
-   FE_Q<dim> fe(p);
-   DoFHandler<dim>        dof_handler(triangulation);
-   dof_handler.distribute_dofs (fe);
- 
-   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; ++q)
-     {
-                                      // project the function
-       VectorTools::project (dof_handler,
-                           constraints,
-                           QGauss<dim>(p+2),
-                           F<dim> (q),
-                           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),
-                                        error,
-                                        QGauss<dim>(std::max(p,q)+1),
-                                        VectorTools::L2_norm);
-       Assert (error.l2_norm() <= 1e-12*projection.l2_norm(),
-             ExcInternalError());
-     }
-   std::cout << "Done: " << __PRETTY_FUNCTION__ << ' ' << p << std::endl;
++  for (unsigned int i = begin; i < end; ++i)
++    dst[i] = 0;
 +}
 +
++template<int dim>
+   void
 -  mass_assembler(const CellIterator &cell, Scratch<dim> &data,
 -      CopyData &copy_data)
++  mass_assembler(const typename Triangulation<dim>::active_cell_iterator &cell,
++      Scratch<dim> &data, CopyData &copy_data)
+   {
+     data.x_fe_values.reinit(cell);
  
 -    copy_data.cell_rhs = 0;
 -    if (cell->center().distance (data.x_fe_values.quadrature_point(0)) >= 1e-6 *cell->diameter())
 -      std::cout << '.' << std::flush;
 -    else
 -      std::cout << '*' << std::flush;
++    const Point<dim> q = data.x_fe_values.quadrature_point(0);
 +
+     // this appears to be the key: the following line overwrites some of the memory
+     // in which we store the quadrature point location. if the line is moved down,
+     // the comparison in the if() always succeeds...
- template <int dim>
- void test ()
++    parallel::apply_to_subranges(0U, copy_data.cell_rhs.size(),
++        std_cxx1x::bind(&zero_subrange, std_cxx1x::_1, std_cxx1x::_2,
++            std_cxx1x::ref(copy_data.cell_rhs)), 1);
 +
 -
 -void copy_local_to_global (const CopyData &data,
 -                           double *sum)
++    std::cout << (q != data.x_fe_values.quadrature_point(0) ? '.' : '*')
++        << std::flush;
+ 
+     copy_data.cell_rhs[0] = value(data.x_fe_values.quadrature_point(0));
+   }
+ 
++void
++copy_local_to_global(const CopyData &data, double *sum)
  {
-   Triangulation<dim>     triangulation;
+   *sum += data.cell_rhs[0];
+ }
+ 
 -
 -void do_project ()
++void
++do_project()
+ {
+   static const int dim = 3;
+ 
 -  Triangulation<dim>     triangulation;
++  Triangulation < dim > triangulation;
    GridGenerator::hyper_cube (triangulation);
--  triangulation.refine_global (2);
 -
 -  FE_Nothing<3> fe;
 -  DoFHandler<dim>        dof_handler(triangulation);
 -  dof_handler.distribute_dofs (fe);
++  triangulation.refine_global(2);
+ 
 -  QMidpoint<dim> q;
++  FE_Nothing < dim > fe;
++  QMidpoint < dim > q;
  
-   Threads::TaskGroup<> g;
-   for (unsigned int p=1; p<4; ++p)
-     g += Threads::new_task (&do_project<dim>, triangulation, p);
-   g.join_all ();
 -  for (unsigned int i=0; i<12; ++i)
++  for (unsigned int i = 0; i < 12; ++i)
+     {
+       std::vector<double> tmp;
+ 
+       double sum = 0;
 -      Scratch<dim> assembler_data (dof_handler.get_fe(), q);
++      Scratch<dim> assembler_data(fe, q);
+       CopyData copy_data;
 -      copy_data.cell_rhs.reinit (8);
 -      dealii::WorkStream::run (dof_handler.begin_active(),
 -                       dof_handler.end(),
 -                       &mass_assembler<dim, typename DoFHandler<dim>::active_cell_iterator>,
 -                       std_cxx1x::bind(&
 -                                       copy_local_to_global,
 -                                       std_cxx1x::_1, &sum),
 -                       assembler_data,
 -                       copy_data,
 -                       8,
 -                       1);
 -      printf ("\nCheck: %5.3f\n", sum);
++      copy_data.cell_rhs.resize(8);
++      WorkStream::run(triangulation.begin_active(), triangulation.end(),
++          &mass_assembler<dim>,
++          std_cxx1x::bind(&copy_local_to_global, std_cxx1x::_1, &sum),
++          assembler_data, copy_data, 8, 1);
++      printf("\nCheck: %5.3f\n", sum);
++      deallog << sum << std::endl;
+     }
  }
  
  
--int main ()
++int main()
  {
    std::ofstream logfile(logname);
--  deallog << std::setprecision (3);
++  deallog << std::setprecision(3);
  
    deallog.attach(logfile);
    deallog.depth_console(0);
    deallog.threshold_double(1.e-10);
  
- 
-   Threads::TaskGroup<> g;
-   for (unsigned int i=0; i<2; ++i)
-     {
-       g += Threads::new_task (&test<1>);
-       g += Threads::new_task (&test<2>);
-       g += Threads::new_task (&test<3>);
-     }
-   g.join_all ();
- 
-   deallog << "OK" << std::endl;
 -  do_project ();
++  do_project();
  }
- 

diff --cc tests/base/work_stream_03/cmp/generic
index 0fd8fc12f0b442283edd8868867114c242b04d11,0fd8fc12f0b442283edd8868867114c242b04d11..d1b460b43fedaa6b40cc98c5a87f30b0114883ea
--- 1/tests/base/work_stream_03/cmp/generic
--- 2/tests/base/work_stream_03/cmp/generic
+++ b/tests/base/work_stream_03/cmp/generic
@@@ -1,2 -1,2 +1,13 @@@
  
--DEAL::OK
++DEAL::288.
++DEAL::288.
++DEAL::288.
++DEAL::288.
++DEAL::288.
++DEAL::288.
++DEAL::288.
++DEAL::288.
++DEAL::288.
++DEAL::288.
++DEAL::288.
++DEAL::288.

diff --cc tests/deal.II/fe_values_view_21_nonsymmetric/cmp/generic
index 0000000000000000000000000000000000000000,80180cbeef65da584a2f61c078b1bc38912f2942..84c5609545e34aac2e5639a0f5c0e14949115d56
mode 000000,100644..100644
--- /dev/null
--- 2/tests/deal.II/fe_values_view_21_nonsymmetric/cmp/generic
+++ b/tests/deal.II/fe_values_view_21_nonsymmetric/cmp/generic
@@@ -1,0 -1,9 +1,9 @@@
 -DEAL::6.50 6.00 5.00
+ 
+ DEAL::Number of stress components: 9  Number of gamma components: 1
+ DEAL::Number of active cells: 1
+ DEAL::Total number of cells: 1
+ 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::3.50 3.50 3.50

diff --cc tests/deal.II/fe_values_view_23_nonsymmetric.cc
index 0000000000000000000000000000000000000000,0000000000000000000000000000000000000000..9adfc38074c8203b259da3d5e5c4a83c57641b0d
new file mode 100644 (file)
--- /dev/null
--- /dev/null
+++ b/tests/deal.II/fe_values_view_23_nonsymmetric.cc
@@@ -1,0 -1,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 --cc tests/deal.II/fe_values_view_23_nonsymmetric/cmp/generic
index 0000000000000000000000000000000000000000,0000000000000000000000000000000000000000..870eda4ef0250a4e3c4e8f13d92ab2e733299428
new file mode 100644 (file)
--- /dev/null
--- /dev/null
+++ b/tests/deal.II/fe_values_view_23_nonsymmetric/cmp/generic
@@@ -1,0 -1,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 --cc tests/deal.II/fe_values_view_23_nonsymmetric/fe_values_view_09/cmp/generic
index 0000000000000000000000000000000000000000,0000000000000000000000000000000000000000..977609c8f410667b0f6d6e330a9f1299c7eade55
new file mode 100644 (file)
--- /dev/null
--- /dev/null
+++ b/tests/deal.II/fe_values_view_23_nonsymmetric/fe_values_view_09/cmp/generic
@@@ -1,0 -1,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 --cc tests/deal.II/fe_values_view_24_nonsymmetric.cc
index 0000000000000000000000000000000000000000,0000000000000000000000000000000000000000..7d1335f1a97f7ed80fb14b0f90a2493345c53352
new file mode 100644 (file)
--- /dev/null
--- /dev/null
+++ b/tests/deal.II/fe_values_view_24_nonsymmetric.cc
@@@ -1,0 -1,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 --cc tests/deal.II/fe_values_view_24_nonsymmetric/cmp/generic
index 0000000000000000000000000000000000000000,0000000000000000000000000000000000000000..c48eecd5e9b78bb8b36a75343abadde7d98fe46d
new file mode 100644 (file)
--- /dev/null
--- /dev/null
+++ b/tests/deal.II/fe_values_view_24_nonsymmetric/cmp/generic
@@@ -1,0 -1,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 --cc tests/deal.II/fe_values_view_24_nonsymmetric/fe_values_view_09/cmp/generic
index 0000000000000000000000000000000000000000,0000000000000000000000000000000000000000..977609c8f410667b0f6d6e330a9f1299c7eade55
new file mode 100644 (file)
--- /dev/null
--- /dev/null
+++ b/tests/deal.II/fe_values_view_24_nonsymmetric/fe_values_view_09/cmp/generic
@@@ -1,0 -1,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 --cc tests/deal.II/fe_values_view_25_nonsymmetric.cc
index 0000000000000000000000000000000000000000,0000000000000000000000000000000000000000..04fb1e2f40b471d75e14c9a14508f80de5441618
new file mode 100644 (file)
--- /dev/null
--- /dev/null
+++ b/tests/deal.II/fe_values_view_25_nonsymmetric.cc
@@@ -1,0 -1,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 --cc tests/deal.II/fe_values_view_25_nonsymmetric/cmp/generic
index 0000000000000000000000000000000000000000,0000000000000000000000000000000000000000..440f9572cff801c77fd7471bc35e775be43d15a9
new file mode 100644 (file)
--- /dev/null
--- /dev/null
+++ b/tests/deal.II/fe_values_view_25_nonsymmetric/cmp/generic
@@@ -1,0 -1,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 --cc tests/deal.II/fe_values_view_25_nonsymmetric/fe_values_view_09/cmp/generic
index 0000000000000000000000000000000000000000,0000000000000000000000000000000000000000..977609c8f410667b0f6d6e330a9f1299c7eade55
new file mode 100644 (file)
--- /dev/null
--- /dev/null
+++ b/tests/deal.II/fe_values_view_25_nonsymmetric/fe_values_view_09/cmp/generic
@@@ -1,0 -1,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