]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Implement using the correct transformation also for subfaces
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 25 Aug 2006 14:48:03 +0000 (14:48 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 25 Aug 2006 14:48:03 +0000 (14:48 +0000)
git-svn-id: https://svn.dealii.org/trunk@13743 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/fe/fe_poly_tensor.cc
tests/fe/rt_14.cc
tests/fe/rt_14/cmp/generic
tests/fe/rt_15.cc

index cf97b95d1e0e4601542d64c06f64a6edbe1dfb56..3a980be316ddfcaf22ba066050a1c384325628da 100644 (file)
@@ -506,11 +506,11 @@ FE_PolyTensor<POLY,dim>::fill_fe_face_values (
 //      ExcNotImplemented());
 //TODO: Size assertions
 
-  // Create table with sign changes, due to the special structure of the RT elements.
-  // TODO: Preliminary hack to demonstrate the overall prinicple!
+// Create table with sign changes, due to the special structure of the RT elements.
+// TODO: Preliminary hack to demonstrate the overall prinicple!
 
-  // Compute eventual sign changes depending on the neighborhood
-  // between two faces.
+                                  // Compute eventual sign changes depending
+                                  // on the neighborhood between two faces.
   std::vector<double> sign_change (n_dofs, 1.0);
   get_face_sign_change (cell, sign_change);
   
@@ -524,37 +524,54 @@ FE_PolyTensor<POLY,dim>::fill_fe_face_values (
            {
              case independent:
              case independent_on_cartesian:
+             {
                for (unsigned int k=0; k<n_quad; ++k)
                  for (unsigned int d=0;d<dim;++d)
                    data.shape_values(first+d,k) = fe_data.shape_values[i][k+offset][d];
                break;
+             }
+             
              case covariant:
              case contravariant:
-               if (true)
-                 {
-                                                    // Use auxiliary vector for transformation
-                   std::vector<Tensor<1,dim> > shape_values (n_quad);
-                   if (mapping_type == covariant)
-                     mapping.transform_covariant(fe_data.shape_values[i], offset,
+             {
+                                                // Use auxiliary vector
+                                                // for transformation
+               std::vector<Tensor<1,dim> > shape_values (n_quad);
+               if (mapping_type == covariant)
+                 mapping.transform_covariant(fe_data.shape_values[i], offset,
+                                             shape_values, mapping_data);
+               else
+                 mapping.transform_contravariant(fe_data.shape_values[i], offset,
                                                  shape_values, mapping_data);
-                   else
-                     mapping.transform_contravariant(fe_data.shape_values[i], offset,
-                                                     shape_values, mapping_data);
+               
+                                                // then copy over to target:
+               for (unsigned int k=0; k<n_quad; ++k)
+                 {
+                                                    // recompute
+                                                    // determinant. note that
+                                                    // we have to take the
+                                                    // *cell* transformation
+                                                    // (cell_JxW_values)
+                                                    // because the usual
+                                                    // JxW_values field
+                                                    // contains the
+                                                    // transformation of the
+                                                    // face in
+                                                    // fill_fe_face_values
+                   const double
+                     J = (mapping_type == contravariant ?
+                          data.cell_JxW_values[k] / quadrature.weight(k)
+                          :
+                          1.0);
                    
-                                                    // then copy over to target:
-                   for (unsigned int k=0; k<n_quad; ++k)
-                     {
-                       // Recompute determinant
-//TODO: Note that the Jacobian we compute here is that of the transformation from the reference FACE to the real FACE. However, what we need is the Jacobian of the transformation of the reference to the real CELL, which differs by one power of h. Consequently, the tests fe/rt_{11,13,14,15} presently compute the wrong results
-                       double J = 1.0;
-                       if (mapping_type == contravariant)
-                         J = data.cell_JxW_values[k] / quadrature.weight(k);
-
-                       for (unsigned int d=0; d<dim; ++d)
-                         data.shape_values(first+d,k) = sign_change[i] * (shape_values[k][d] / J);
-                     }
+                   for (unsigned int d=0; d<dim; ++d)
+                     data.shape_values(first+d,k)
+                       = sign_change[i] * (shape_values[k][d] / J);
                  }
+               
                break;
+             }
+             
              default:
                Assert(false, ExcNotImplemented());
            }
@@ -564,17 +581,22 @@ FE_PolyTensor<POLY,dim>::fill_fe_face_values (
        {
          std::vector<Tensor<2,dim> > shape_grads1 (n_quad);
          std::vector<Tensor<2,dim> > shape_grads2 (n_quad);
+
          switch(mapping_type)
            {
              case independent:
              case independent_on_cartesian:
+             {
                mapping.transform_covariant(fe_data.shape_grads[i], offset,
                                            shape_grads1, mapping_data);
                for (unsigned int k=0; k<n_quad; ++k)
                  for (unsigned int d=0;d<dim;++d)
                    data.shape_gradients[first+d][k] = shape_grads1[k][d];
                break;
+             }
+             
              case covariant:
+             {
                mapping.transform_covariant(fe_data.shape_grads[i], offset,
                                            shape_grads1,
                                            mapping_data);
@@ -591,8 +613,10 @@ FE_PolyTensor<POLY,dim>::fill_fe_face_values (
                  for (unsigned int d=0;d<dim;++d)
                    data.shape_gradients[first+d][k] = shape_grads1[k][d];
                break;
-
+             }
+             
              case contravariant:
+             {
                mapping.transform_covariant(fe_data.shape_grads[i], offset,
                                            shape_grads1,
                                            mapping_data);
@@ -604,15 +628,17 @@ FE_PolyTensor<POLY,dim>::fill_fe_face_values (
                for (unsigned int k=0; k<n_quad; ++k)
                  for (unsigned int d=0;d<dim;++d)
                    {
-                     // Recompute determinant
-//TODO: Note that the Jacobian we compute here is that of the transformation from the reference FACE to the real FACE. However, what we need is the Jacobian of the transformation of the reference to the real CELL, which differs by one power of h. Consequently, the tests fe/rt_{11,13,14,15} presently compute the wrong results
-
-                     double J = data.cell_JxW_values[k] / quadrature.weight(k);
-                     data.shape_gradients[first+d][k] = sign_change[i] * 
-                       shape_grads2[k][d] / J;
+                                                      // recompute
+                                                      // determinant in the
+                                                      // same way as above
+                     const double
+                       J = data.cell_JxW_values[k] / quadrature.weight(k);
+                     data.shape_gradients[first+d][k]
+                       = sign_change[i] * shape_grads2[k][d] / J;
                    }
                break;
-               
+             }
+             
              default:
                Assert(false, ExcNotImplemented());
            }
@@ -674,29 +700,54 @@ FE_PolyTensor<POLY,dim>::fill_fe_subface_values (
            {
              case independent:
              case independent_on_cartesian:
+             {
                for (unsigned int k=0; k<n_quad; ++k)
                  for (unsigned int d=0;d<dim;++d)
                    data.shape_values(first+d,k) = fe_data.shape_values[i][k+offset][d];
                break;
+             }
+             
              case covariant:
              case contravariant:
-               if (true)
-                 {
-                                                    // Use auxiliary vector for transformation
-                   std::vector<Tensor<1,dim> > shape_values (n_quad);
-                   if (mapping_type == covariant)
-                     mapping.transform_covariant(fe_data.shape_values[i], offset,
+             {
+                                                // Use auxiliary vector for
+                                                // transformation
+               std::vector<Tensor<1,dim> > shape_values (n_quad);
+               if (mapping_type == covariant)
+                 mapping.transform_covariant(fe_data.shape_values[i], offset,
+                                             shape_values, mapping_data);
+               else
+                 mapping.transform_contravariant(fe_data.shape_values[i], offset,
                                                  shape_values, mapping_data);
-                   else
-                     mapping.transform_contravariant(fe_data.shape_values[i], offset,
-                                                     shape_values, mapping_data);
                    
-                                                    // then copy over to target:
-                   for (unsigned int k=0; k<n_quad; ++k)
-                     for (unsigned int d=0; d<dim; ++d)
-                       data.shape_values(first+d,k) = shape_values[k][d];
+                                                // then copy over to target:
+               for (unsigned int k=0; k<n_quad; ++k)
+                 {
+                                                    // recompute
+                                                    // determinant. note that
+                                                    // we have to take the
+                                                    // *cell* transformation
+                                                    // (cell_JxW_values)
+                                                    // because the usual
+                                                    // JxW_values field
+                                                    // contains the
+                                                    // transformation of the
+                                                    // face in
+                                                    // fill_fe_face_values
+                   const double
+                     J = (mapping_type == contravariant ?
+                          data.cell_JxW_values[k] / quadrature.weight(k)
+                          :
+                          1.0);
+                   
+                   for (unsigned int d=0; d<dim; ++d)
+                     data.shape_values(first+d,k)
+                       = /* sign_change[i] * */ (shape_values[k][d] / J);
                  }
+
                break;
+             }
+             
              default:
                Assert(false, ExcNotImplemented());
            }
@@ -706,17 +757,22 @@ FE_PolyTensor<POLY,dim>::fill_fe_subface_values (
        {
          std::vector<Tensor<2,dim> > shape_grads1 (n_quad);
          std::vector<Tensor<2,dim> > shape_grads2 (n_quad);
+
          switch(mapping_type)
            {
              case independent:
              case independent_on_cartesian:
+             {
                mapping.transform_covariant(fe_data.shape_grads[i], offset,
                                            shape_grads1, mapping_data);
                for (unsigned int k=0; k<n_quad; ++k)
                  for (unsigned int d=0;d<dim;++d)
                    data.shape_gradients[first+d][k] = shape_grads1[k][d];
                break;
+             }
+             
              case covariant:
+             {
                mapping.transform_covariant(fe_data.shape_grads[i], offset,
                                            shape_grads1,
                                            mapping_data);
@@ -733,6 +789,32 @@ FE_PolyTensor<POLY,dim>::fill_fe_subface_values (
                  for (unsigned int d=0;d<dim;++d)
                    data.shape_gradients[first+d][k] = shape_grads1[k][d];
                break;
+             }
+             
+             case contravariant:
+             {
+               mapping.transform_covariant(fe_data.shape_grads[i], offset,
+                                           shape_grads1,
+                                           mapping_data);
+
+               mapping.transform_contravariant(shape_grads1, 0,
+                                               shape_grads2,
+                                               mapping_data);
+
+               for (unsigned int k=0; k<n_quad; ++k)
+                 for (unsigned int d=0;d<dim;++d)
+                   {
+                                                      // recompute
+                                                      // determinant in the
+                                                      // same way as above
+                     const double
+                       J = data.cell_JxW_values[k] / quadrature.weight(k);
+                     data.shape_gradients[first+d][k]
+                       = /* sign_change[i] * */ shape_grads2[k][d] / J;
+                   }
+               break;
+             }
+
              default:
                Assert(false, ExcNotImplemented());
            }
index 3d163826dd326f59ff14b22de32a44c03da3e360..3aeff2cff1afaf16d2a6281164bbf32dba82468e 100644 (file)
@@ -13,7 +13,7 @@
 
 // like rt_11, but use FESubfaceValues
 //
-// the test presently fails because of the issue with computing the
+// the test used to fail because of the issue with computing the
 // normals using FEFaceValue, where FEFaceValue by accident uses the
 // *face* mapping, not the *cell* mapping to compute the Piola
 // transform (leading to a missing power of h in the determinant)
index c927bbf7b67830d7dc0c9fd4d4b42737b52383b1..ba571c1c8dd1cc448583da0a8f09c0ee371eb817 100644 (file)
 
 DEAL::Degree=0
 DEAL::  h=1.00
-DEAL::    Quadrature point 0: [-1.00 0.00 0.00 0.00 ][1.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -1.00 ][0.00 0.00 0.00 1.00 ]
-DEAL::    Quadrature point 1: [-1.00 0.00 0.00 0.00 ][1.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -1.00 ][0.00 0.00 0.00 1.00 ]
-DEAL::    Quadrature point 2: [-1.00 0.00 0.00 0.00 ][1.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -1.00 ][0.00 0.00 0.00 1.00 ]
-DEAL::    Quadrature point 3: [-1.00 0.00 0.00 0.00 ][1.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -1.00 ][0.00 0.00 0.00 1.00 ]
+DEAL::    Quadrature point 0: [1.00 0 ][0 0 ][0 1.00 ][0 0 ]
+DEAL::    Quadrature point 1: [1.00 0 ][0 0 ][0 0.50 ][0 0.50 ]
 DEAL::  h=0.50
-DEAL::    Quadrature point 0: [-4.00 0.00 0.00 0.00 ][4.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -4.00 ][0.00 0.00 0.00 4.00 ]
-DEAL::    Quadrature point 1: [-4.00 0.00 0.00 0.00 ][4.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -4.00 ][0.00 0.00 0.00 4.00 ]
-DEAL::    Quadrature point 2: [-4.00 0.00 0.00 0.00 ][4.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -4.00 ][0.00 0.00 0.00 4.00 ]
-DEAL::    Quadrature point 3: [-4.00 0.00 0.00 0.00 ][4.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -4.00 ][0.00 0.00 0.00 4.00 ]
+DEAL::    Quadrature point 0: [2.00 0 ][0 0 ][0 2.00 ][0 0 ]
+DEAL::    Quadrature point 1: [2.00 0 ][0 0 ][0 1.00 ][0 1.00 ]
 DEAL::  h=0.25
-DEAL::    Quadrature point 0: [-16.00 0.00 0.00 0.00 ][16.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -16.00 ][0.00 0.00 0.00 16.00 ]
-DEAL::    Quadrature point 1: [-16.00 0.00 0.00 0.00 ][16.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -16.00 ][0.00 0.00 0.00 16.00 ]
-DEAL::    Quadrature point 2: [-16.00 0.00 0.00 0.00 ][16.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -16.00 ][0.00 0.00 0.00 16.00 ]
-DEAL::    Quadrature point 3: [-16.00 0.00 0.00 0.00 ][16.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -16.00 ][0.00 0.00 0.00 16.00 ]
+DEAL::    Quadrature point 0: [4.00 0 ][0 0 ][0 4.00 ][0 0 ]
+DEAL::    Quadrature point 1: [4.00 0 ][0 0 ][0 2.00 ][0 2.00 ]
 DEAL::  h=0.12
-DEAL::    Quadrature point 0: [-64.00 0.00 0.00 0.00 ][64.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -64.00 ][0.00 0.00 0.00 64.00 ]
-DEAL::    Quadrature point 1: [-64.00 0.00 0.00 0.00 ][64.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -64.00 ][0.00 0.00 0.00 64.00 ]
-DEAL::    Quadrature point 2: [-64.00 0.00 0.00 0.00 ][64.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -64.00 ][0.00 0.00 0.00 64.00 ]
-DEAL::    Quadrature point 3: [-64.00 0.00 0.00 0.00 ][64.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -64.00 ][0.00 0.00 0.00 64.00 ]
+DEAL::    Quadrature point 0: [8.00 0 ][0 0 ][0 8.00 ][0 0 ]
+DEAL::    Quadrature point 1: [8.00 0 ][0 0 ][0 4.00 ][0 4.00 ]
 DEAL::  h=0.06
-DEAL::    Quadrature point 0: [-256.00 0.00 0.00 0.00 ][256.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -256.00 ][0.00 0.00 0.00 256.00 ]
-DEAL::    Quadrature point 1: [-256.00 0.00 0.00 0.00 ][256.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -256.00 ][0.00 0.00 0.00 256.00 ]
-DEAL::    Quadrature point 2: [-256.00 0.00 0.00 0.00 ][256.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -256.00 ][0.00 0.00 0.00 256.00 ]
-DEAL::    Quadrature point 3: [-256.00 0.00 0.00 0.00 ][256.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -256.00 ][0.00 0.00 0.00 256.00 ]
+DEAL::    Quadrature point 0: [16.00 0 ][0 0 ][0 16.00 ][0 0 ]
+DEAL::    Quadrature point 1: [16.00 0 ][0 0 ][0 8.00 ][0 8.00 ]
 DEAL::  h=0.03
-DEAL::    Quadrature point 0: [-1024.00 0.00 0.00 0.00 ][1024.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -1024.00 ][0.00 0.00 0.00 1024.00 ]
-DEAL::    Quadrature point 1: [-1024.00 0.00 0.00 0.00 ][1024.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -1024.00 ][0.00 0.00 0.00 1024.00 ]
-DEAL::    Quadrature point 2: [-1024.00 0.00 0.00 0.00 ][1024.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -1024.00 ][0.00 0.00 0.00 1024.00 ]
-DEAL::    Quadrature point 3: [-1024.00 0.00 0.00 0.00 ][1024.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -1024.00 ][0.00 0.00 0.00 1024.00 ]
+DEAL::    Quadrature point 0: [32.00 0 ][0 0 ][0 32.00 ][0 0 ]
+DEAL::    Quadrature point 1: [32.00 0 ][0 0 ][0 16.00 ][0 16.00 ]
 DEAL::  h=0.02
-DEAL::    Quadrature point 0: [-4096.00 0.00 0.00 0.00 ][4096.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -4096.00 ][0.00 0.00 0.00 4096.00 ]
-DEAL::    Quadrature point 1: [-4096.00 0.00 0.00 0.00 ][4096.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -4096.00 ][0.00 0.00 0.00 4096.00 ]
-DEAL::    Quadrature point 2: [-4096.00 0.00 0.00 0.00 ][4096.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -4096.00 ][0.00 0.00 0.00 4096.00 ]
-DEAL::    Quadrature point 3: [-4096.00 0.00 0.00 0.00 ][4096.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -4096.00 ][0.00 0.00 0.00 4096.00 ]
+DEAL::    Quadrature point 0: [64.00 0 ][0 0 ][0 64.00 ][0 0 ]
+DEAL::    Quadrature point 1: [64.00 0 ][0 0 ][0 32.00 ][0 32.00 ]
 DEAL::Degree=1
 DEAL::  h=1.00
-DEAL::    Quadrature point 0: [-4.00 0.00 0.00 0.00 ][6.93 3.46 0.00 0.00 ][-2.00 0.00 0.00 0.00 ][3.46 0.00 0.00 0.00 ][0.00 0.00 0.00 -4.00 ][0.00 0.00 3.46 6.93 ][0.00 0.00 0.00 -2.00 ][0.00 0.00 0.00 3.46 ][6.00 0.00 0.00 0.00 ][0.00 0.00 0.00 6.00 ][-10.39 0.00 0.00 0.00 ][0.00 0.00 0.00 -10.39 ]
-DEAL::    Quadrature point 1: [2.00 0.00 0.00 0.00 ][-3.46 0.00 0.00 0.00 ][4.00 0.00 0.00 0.00 ][-6.93 3.46 0.00 0.00 ][0.00 0.00 0.00 -4.00 ][0.00 0.00 3.46 -6.93 ][0.00 0.00 0.00 -2.00 ][0.00 0.00 0.00 -3.46 ][-6.00 0.00 0.00 0.00 ][0.00 0.00 0.00 6.00 ][10.39 0.00 0.00 0.00 ][0.00 0.00 0.00 10.39 ]
-DEAL::    Quadrature point 2: [-4.00 0.00 0.00 0.00 ][-6.93 3.46 0.00 0.00 ][-2.00 0.00 0.00 0.00 ][-3.46 0.00 0.00 0.00 ][0.00 0.00 0.00 2.00 ][0.00 0.00 0.00 -3.46 ][0.00 0.00 0.00 4.00 ][0.00 0.00 3.46 -6.93 ][6.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -6.00 ][10.39 0.00 0.00 0.00 ][0.00 0.00 0.00 10.39 ]
-DEAL::    Quadrature point 3: [2.00 0.00 0.00 0.00 ][3.46 0.00 0.00 0.00 ][4.00 0.00 0.00 0.00 ][6.93 3.46 0.00 0.00 ][0.00 0.00 0.00 2.00 ][0.00 0.00 0.00 3.46 ][0.00 0.00 0.00 4.00 ][0.00 0.00 3.46 6.93 ][-6.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -6.00 ][-10.39 0.00 0.00 0.00 ][0.00 0.00 0.00 -10.39 ]
+DEAL::    Quadrature point 0: [1.00 0 ][-1.73 0 ][0 0 ][0 0 ][0 1.00 ][0 -1.73 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL::    Quadrature point 1: [1.00 0 ][0 0 ][0 0 ][0 0 ][0 -0.25 ][0 0.43 ][0 -0.25 ][0 0.43 ][0 0 ][0 1.50 ][0 0 ][0 -2.60 ]
 DEAL::  h=0.50
-DEAL::    Quadrature point 0: [-16.00 0.00 0.00 0.00 ][27.71 13.86 0.00 0.00 ][-8.00 0.00 0.00 0.00 ][13.86 0.00 0.00 0.00 ][0.00 0.00 0.00 -16.00 ][0.00 0.00 13.86 27.71 ][0.00 0.00 0.00 -8.00 ][0.00 0.00 0.00 13.86 ][24.00 0.00 0.00 0.00 ][0.00 0.00 0.00 24.00 ][-41.57 0.00 0.00 0.00 ][0.00 0.00 0.00 -41.57 ]
-DEAL::    Quadrature point 1: [8.00 0.00 0.00 0.00 ][-13.86 0.00 0.00 0.00 ][16.00 0.00 0.00 0.00 ][-27.71 13.86 0.00 0.00 ][0.00 0.00 0.00 -16.00 ][0.00 0.00 13.86 -27.71 ][0.00 0.00 0.00 -8.00 ][0.00 0.00 0.00 -13.86 ][-24.00 0.00 0.00 0.00 ][0.00 0.00 0.00 24.00 ][41.57 0.00 0.00 0.00 ][0.00 0.00 0.00 41.57 ]
-DEAL::    Quadrature point 2: [-16.00 0.00 0.00 0.00 ][-27.71 13.86 0.00 0.00 ][-8.00 0.00 0.00 0.00 ][-13.86 0.00 0.00 0.00 ][0.00 0.00 0.00 8.00 ][0.00 0.00 0.00 -13.86 ][0.00 0.00 0.00 16.00 ][0.00 0.00 13.86 -27.71 ][24.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -24.00 ][41.57 0.00 0.00 0.00 ][0.00 0.00 0.00 41.57 ]
-DEAL::    Quadrature point 3: [8.00 0.00 0.00 0.00 ][13.86 0.00 0.00 0.00 ][16.00 0.00 0.00 0.00 ][27.71 13.86 0.00 0.00 ][0.00 0.00 0.00 8.00 ][0.00 0.00 0.00 13.86 ][0.00 0.00 0.00 16.00 ][0.00 0.00 13.86 27.71 ][-24.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -24.00 ][-41.57 0.00 0.00 0.00 ][0.00 0.00 0.00 -41.57 ]
+DEAL::    Quadrature point 0: [2.00 0 ][-3.46 0 ][0 0 ][0 0 ][0 2.00 ][0 -3.46 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL::    Quadrature point 1: [2.00 0 ][0 0 ][0 0 ][0 0 ][0 -0.50 ][0 0.87 ][0 -0.50 ][0 0.87 ][0 0 ][0 3.00 ][0 0 ][0 -5.20 ]
 DEAL::  h=0.25
-DEAL::    Quadrature point 0: [-64.00 0.00 0.00 0.00 ][110.85 55.43 0.00 0.00 ][-32.00 0.00 0.00 0.00 ][55.43 0.00 0.00 0.00 ][0.00 0.00 0.00 -64.00 ][0.00 0.00 55.43 110.85 ][0.00 0.00 0.00 -32.00 ][0.00 0.00 0.00 55.43 ][96.00 0.00 0.00 0.00 ][0.00 0.00 0.00 96.00 ][-166.28 0.00 0.00 0.00 ][0.00 0.00 0.00 -166.28 ]
-DEAL::    Quadrature point 1: [32.00 0.00 0.00 0.00 ][-55.43 0.00 0.00 0.00 ][64.00 0.00 0.00 0.00 ][-110.85 55.43 0.00 0.00 ][0.00 0.00 0.00 -64.00 ][0.00 0.00 55.43 -110.85 ][0.00 0.00 0.00 -32.00 ][0.00 0.00 0.00 -55.43 ][-96.00 0.00 0.00 0.00 ][0.00 0.00 0.00 96.00 ][166.28 0.00 0.00 0.00 ][0.00 0.00 0.00 166.28 ]
-DEAL::    Quadrature point 2: [-64.00 0.00 0.00 0.00 ][-110.85 55.43 0.00 0.00 ][-32.00 0.00 0.00 0.00 ][-55.43 0.00 0.00 0.00 ][0.00 0.00 0.00 32.00 ][0.00 0.00 0.00 -55.43 ][0.00 0.00 0.00 64.00 ][0.00 0.00 55.43 -110.85 ][96.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -96.00 ][166.28 0.00 0.00 0.00 ][0.00 0.00 0.00 166.28 ]
-DEAL::    Quadrature point 3: [32.00 0.00 0.00 0.00 ][55.43 0.00 0.00 0.00 ][64.00 0.00 0.00 0.00 ][110.85 55.43 0.00 0.00 ][0.00 0.00 0.00 32.00 ][0.00 0.00 0.00 55.43 ][0.00 0.00 0.00 64.00 ][0.00 0.00 55.43 110.85 ][-96.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -96.00 ][-166.28 0.00 0.00 0.00 ][0.00 0.00 0.00 -166.28 ]
+DEAL::    Quadrature point 0: [4.00 0 ][-6.93 0 ][0 0 ][0 0 ][0 4.00 ][0 -6.93 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL::    Quadrature point 1: [4.00 0 ][0 0 ][0 0 ][0 0 ][0 -1.00 ][0 1.73 ][0 -1.00 ][0 1.73 ][0 0 ][0 6.00 ][0 0 ][0 -10.39 ]
 DEAL::  h=0.12
-DEAL::    Quadrature point 0: [-256.00 0.00 0.00 0.00 ][443.41 221.70 0.00 0.00 ][-128.00 0.00 0.00 0.00 ][221.70 0.00 0.00 0.00 ][0.00 0.00 0.00 -256.00 ][0.00 0.00 221.70 443.41 ][0.00 0.00 0.00 -128.00 ][0.00 0.00 0.00 221.70 ][384.00 0.00 0.00 0.00 ][0.00 0.00 0.00 384.00 ][-665.11 0.00 0.00 0.00 ][0.00 0.00 0.00 -665.11 ]
-DEAL::    Quadrature point 1: [128.00 0.00 0.00 0.00 ][-221.70 0.00 0.00 0.00 ][256.00 0.00 0.00 0.00 ][-443.41 221.70 0.00 0.00 ][0.00 0.00 0.00 -256.00 ][0.00 0.00 221.70 -443.41 ][0.00 0.00 0.00 -128.00 ][0.00 0.00 0.00 -221.70 ][-384.00 0.00 0.00 0.00 ][0.00 0.00 0.00 384.00 ][665.11 0.00 0.00 0.00 ][0.00 0.00 0.00 665.11 ]
-DEAL::    Quadrature point 2: [-256.00 0.00 0.00 0.00 ][-443.41 221.70 0.00 0.00 ][-128.00 0.00 0.00 0.00 ][-221.70 0.00 0.00 0.00 ][0.00 0.00 0.00 128.00 ][0.00 0.00 0.00 -221.70 ][0.00 0.00 0.00 256.00 ][0.00 0.00 221.70 -443.41 ][384.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -384.00 ][665.11 0.00 0.00 0.00 ][0.00 0.00 0.00 665.11 ]
-DEAL::    Quadrature point 3: [128.00 0.00 0.00 0.00 ][221.70 0.00 0.00 0.00 ][256.00 0.00 0.00 0.00 ][443.41 221.70 0.00 0.00 ][0.00 0.00 0.00 128.00 ][0.00 0.00 0.00 221.70 ][0.00 0.00 0.00 256.00 ][0.00 0.00 221.70 443.41 ][-384.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -384.00 ][-665.11 0.00 0.00 0.00 ][0.00 0.00 0.00 -665.11 ]
+DEAL::    Quadrature point 0: [8.00 0 ][-13.86 0 ][0 0 ][0 0 ][0 8.00 ][0 -13.86 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL::    Quadrature point 1: [8.00 0 ][0 0 ][0 0 ][0 0 ][0 -2.00 ][0 3.46 ][0 -2.00 ][0 3.46 ][0 0 ][0 12.00 ][0 0 ][0 -20.78 ]
 DEAL::  h=0.06
-DEAL::    Quadrature point 0: [-1024.00 0.00 0.00 0.00 ][1773.62 886.81 0.00 0.00 ][-512.00 0.00 0.00 0.00 ][886.81 0.00 0.00 0.00 ][0.00 0.00 0.00 -1024.00 ][0.00 0.00 886.81 1773.62 ][0.00 0.00 0.00 -512.00 ][0.00 0.00 0.00 886.81 ][1536.00 0.00 0.00 0.00 ][0.00 0.00 0.00 1536.00 ][-2660.43 0.00 0.00 0.00 ][0.00 0.00 0.00 -2660.43 ]
-DEAL::    Quadrature point 1: [512.00 0.00 0.00 0.00 ][-886.81 0.00 0.00 0.00 ][1024.00 0.00 0.00 0.00 ][-1773.62 886.81 0.00 0.00 ][0.00 0.00 0.00 -1024.00 ][0.00 0.00 886.81 -1773.62 ][0.00 0.00 0.00 -512.00 ][0.00 0.00 0.00 -886.81 ][-1536.00 0.00 0.00 0.00 ][0.00 0.00 0.00 1536.00 ][2660.43 0.00 0.00 0.00 ][0.00 0.00 0.00 2660.43 ]
-DEAL::    Quadrature point 2: [-1024.00 0.00 0.00 0.00 ][-1773.62 886.81 0.00 0.00 ][-512.00 0.00 0.00 0.00 ][-886.81 0.00 0.00 0.00 ][0.00 0.00 0.00 512.00 ][0.00 0.00 0.00 -886.81 ][0.00 0.00 0.00 1024.00 ][0.00 0.00 886.81 -1773.62 ][1536.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -1536.00 ][2660.43 0.00 0.00 0.00 ][0.00 0.00 0.00 2660.43 ]
-DEAL::    Quadrature point 3: [512.00 0.00 0.00 0.00 ][886.81 0.00 0.00 0.00 ][1024.00 0.00 0.00 0.00 ][1773.62 886.81 0.00 0.00 ][0.00 0.00 0.00 512.00 ][0.00 0.00 0.00 886.81 ][0.00 0.00 0.00 1024.00 ][0.00 0.00 886.81 1773.62 ][-1536.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -1536.00 ][-2660.43 0.00 0.00 0.00 ][0.00 0.00 0.00 -2660.43 ]
+DEAL::    Quadrature point 0: [16.00 0 ][-27.71 0 ][0 0 ][0 0 ][0 16.00 ][0 -27.71 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL::    Quadrature point 1: [16.00 0 ][0 0 ][0 0 ][0 0 ][0 -4.00 ][0 6.93 ][0 -4.00 ][0 6.93 ][0 0 ][0 24.00 ][0 0 ][0 -41.57 ]
 DEAL::  h=0.03
-DEAL::    Quadrature point 0: [-4096.00 0.00 0.00 0.00 ][7094.48 3547.24 0.00 0.00 ][-2048.00 0.00 0.00 0.00 ][3547.24 0.00 0.00 0.00 ][0.00 0.00 0.00 -4096.00 ][0.00 0.00 3547.24 7094.48 ][0.00 0.00 0.00 -2048.00 ][0.00 0.00 0.00 3547.24 ][6144.00 0.00 0.00 0.00 ][0.00 0.00 0.00 6144.00 ][-10641.72 0.00 0.00 0.00 ][0.00 0.00 0.00 -10641.72 ]
-DEAL::    Quadrature point 1: [2048.00 0.00 0.00 0.00 ][-3547.24 0.00 0.00 0.00 ][4096.00 0.00 0.00 0.00 ][-7094.48 3547.24 0.00 0.00 ][0.00 0.00 0.00 -4096.00 ][0.00 0.00 3547.24 -7094.48 ][0.00 0.00 0.00 -2048.00 ][0.00 0.00 0.00 -3547.24 ][-6144.00 0.00 0.00 0.00 ][0.00 0.00 0.00 6144.00 ][10641.72 0.00 0.00 0.00 ][0.00 0.00 0.00 10641.72 ]
-DEAL::    Quadrature point 2: [-4096.00 0.00 0.00 0.00 ][-7094.48 3547.24 0.00 0.00 ][-2048.00 0.00 0.00 0.00 ][-3547.24 0.00 0.00 0.00 ][0.00 0.00 0.00 2048.00 ][0.00 0.00 0.00 -3547.24 ][0.00 0.00 0.00 4096.00 ][0.00 0.00 3547.24 -7094.48 ][6144.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -6144.00 ][10641.72 0.00 0.00 0.00 ][0.00 0.00 0.00 10641.72 ]
-DEAL::    Quadrature point 3: [2048.00 0.00 0.00 0.00 ][3547.24 0.00 0.00 0.00 ][4096.00 0.00 0.00 0.00 ][7094.48 3547.24 0.00 0.00 ][0.00 0.00 0.00 2048.00 ][0.00 0.00 0.00 3547.24 ][0.00 0.00 0.00 4096.00 ][0.00 0.00 3547.24 7094.48 ][-6144.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -6144.00 ][-10641.72 0.00 0.00 0.00 ][0.00 0.00 0.00 -10641.72 ]
+DEAL::    Quadrature point 0: [32.00 0 ][-55.43 0 ][0 0 ][0 0 ][0 32.00 ][0 -55.43 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL::    Quadrature point 1: [32.00 0 ][0 0 ][0 0 ][0 0 ][0 -8.00 ][0 13.86 ][0 -8.00 ][0 13.86 ][0 0 ][0 48.00 ][0 0 ][0 -83.14 ]
 DEAL::  h=0.02
-DEAL::    Quadrature point 0: [-16384.00 0.00 0.00 0.00 ][28377.92 14188.96 0.00 0.00 ][-8192.00 0.00 0.00 0.00 ][14188.96 0.00 0.00 0.00 ][0.00 0.00 0.00 -16384.00 ][0.00 0.00 14188.96 28377.92 ][0.00 0.00 0.00 -8192.00 ][0.00 0.00 0.00 14188.96 ][24576.00 0.00 0.00 0.00 ][0.00 0.00 0.00 24576.00 ][-42566.88 0.00 0.00 0.00 ][0.00 0.00 0.00 -42566.88 ]
-DEAL::    Quadrature point 1: [8192.00 0.00 0.00 0.00 ][-14188.96 0.00 0.00 0.00 ][16384.00 0.00 0.00 0.00 ][-28377.92 14188.96 0.00 0.00 ][0.00 0.00 0.00 -16384.00 ][0.00 0.00 14188.96 -28377.92 ][0.00 0.00 0.00 -8192.00 ][0.00 0.00 0.00 -14188.96 ][-24576.00 0.00 0.00 0.00 ][0.00 0.00 0.00 24576.00 ][42566.88 0.00 0.00 0.00 ][0.00 0.00 0.00 42566.88 ]
-DEAL::    Quadrature point 2: [-16384.00 0.00 0.00 0.00 ][-28377.92 14188.96 0.00 0.00 ][-8192.00 0.00 0.00 0.00 ][-14188.96 0.00 0.00 0.00 ][0.00 0.00 0.00 8192.00 ][0.00 0.00 0.00 -14188.96 ][0.00 0.00 0.00 16384.00 ][0.00 0.00 14188.96 -28377.92 ][24576.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -24576.00 ][42566.88 0.00 0.00 0.00 ][0.00 0.00 0.00 42566.88 ]
-DEAL::    Quadrature point 3: [8192.00 0.00 0.00 0.00 ][14188.96 0.00 0.00 0.00 ][16384.00 0.00 0.00 0.00 ][28377.92 14188.96 0.00 0.00 ][0.00 0.00 0.00 8192.00 ][0.00 0.00 0.00 14188.96 ][0.00 0.00 0.00 16384.00 ][0.00 0.00 14188.96 28377.92 ][-24576.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -24576.00 ][-42566.88 0.00 0.00 0.00 ][0.00 0.00 0.00 -42566.88 ]
+DEAL::    Quadrature point 0: [64.00 0 ][-110.85 0 ][0 0 ][0 0 ][0 64.00 ][0 -110.85 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL::    Quadrature point 1: [64.00 0 ][0 0 ][0 0 ][0 0 ][0 -16.00 ][0 27.71 ][0 -16.00 ][0 27.71 ][0 0 ][0 96.00 ][0 0 ][0 -166.28 ]
 DEAL::Degree=2
 DEAL::  h=1.00
-DEAL::    Quadrature point 0: [-9.00 0.00 0.00 0.00 ][15.59 3.46 0.00 0.00 ][-20.12 -13.42 0.00 0.00 ][3.00 0.00 0.00 0.00 ][-5.20 0.00 0.00 0.00 ][6.71 0.00 0.00 0.00 ][0.00 0.00 0.00 -9.00 ][0.00 0.00 3.46 15.59 ][0.00 0.00 -13.42 -20.12 ][0.00 0.00 0.00 3.00 ][0.00 0.00 0.00 -5.20 ][0.00 0.00 0.00 6.71 ][6.00 0.00 0.00 0.00 ][0.00 0.00 0.00 6.00 ][-17.32 0.00 0.00 0.00 ][0.00 0.00 0.00 -10.39 ][-10.39 0.00 0.00 0.00 ][0.00 0.00 0.00 13.42 ][30.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -17.32 ][13.42 0.00 0.00 0.00 ][0.00 0.00 0.00 30.00 ][-38.73 0.00 0.00 0.00 ][0.00 0.00 0.00 -38.73 ]
-DEAL::    Quadrature point 1: [-3.00 0.00 0.00 0.00 ][5.20 0.00 0.00 0.00 ][-6.71 0.00 0.00 0.00 ][9.00 0.00 0.00 0.00 ][-15.59 3.46 0.00 0.00 ][20.12 -13.42 0.00 0.00 ][0.00 0.00 0.00 -9.00 ][0.00 0.00 3.46 -15.59 ][0.00 0.00 13.42 -20.12 ][0.00 0.00 0.00 3.00 ][0.00 0.00 0.00 5.20 ][0.00 0.00 0.00 6.71 ][-6.00 0.00 0.00 0.00 ][0.00 0.00 0.00 6.00 ][-17.32 0.00 0.00 0.00 ][0.00 0.00 0.00 10.39 ][10.39 0.00 0.00 0.00 ][0.00 0.00 0.00 13.42 ][30.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -17.32 ][-13.42 0.00 0.00 0.00 ][0.00 0.00 0.00 -30.00 ][-38.73 0.00 0.00 0.00 ][0.00 0.00 0.00 -38.73 ]
-DEAL::    Quadrature point 2: [-9.00 0.00 0.00 0.00 ][-15.59 3.46 0.00 0.00 ][-20.12 13.42 0.00 0.00 ][3.00 0.00 0.00 0.00 ][5.20 0.00 0.00 0.00 ][6.71 0.00 0.00 0.00 ][0.00 0.00 0.00 -3.00 ][0.00 0.00 0.00 5.20 ][0.00 0.00 0.00 -6.71 ][0.00 0.00 0.00 9.00 ][0.00 0.00 3.46 -15.59 ][0.00 0.00 -13.42 20.12 ][6.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -6.00 ][-17.32 0.00 0.00 0.00 ][0.00 0.00 0.00 10.39 ][10.39 0.00 0.00 0.00 ][0.00 0.00 0.00 -13.42 ][-30.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -17.32 ][13.42 0.00 0.00 0.00 ][0.00 0.00 0.00 30.00 ][-38.73 0.00 0.00 0.00 ][0.00 0.00 0.00 -38.73 ]
-DEAL::    Quadrature point 3: [-3.00 0.00 0.00 0.00 ][-5.20 0.00 0.00 0.00 ][-6.71 0.00 0.00 0.00 ][9.00 0.00 0.00 0.00 ][15.59 3.46 0.00 0.00 ][20.12 13.42 0.00 0.00 ][0.00 0.00 0.00 -3.00 ][0.00 0.00 0.00 -5.20 ][0.00 0.00 0.00 -6.71 ][0.00 0.00 0.00 9.00 ][0.00 0.00 3.46 15.59 ][0.00 0.00 13.42 20.12 ][-6.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -6.00 ][-17.32 0.00 0.00 0.00 ][0.00 0.00 0.00 -10.39 ][-10.39 0.00 0.00 0.00 ][0.00 0.00 0.00 -13.42 ][-30.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -17.32 ][-13.42 0.00 0.00 0.00 ][0.00 0.00 0.00 -30.00 ][-38.73 0.00 0.00 0.00 ][0.00 0.00 0.00 -38.73 ]
+DEAL::    Quadrature point 0: [1.00 0 ][-1.73 0 ][2.24 0 ][0 0 ][0 0 ][0 0 ][0 1.00 ][0 -1.73 ][0 2.24 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL::    Quadrature point 1: [1.00 0 ][0 0 ][-1.12 0 ][0 0 ][0 0 ][0 0 ][0 -0.25 ][0 0.43 ][0 -0.56 ][0 -0.25 ][0 0.43 ][0 -0.56 ][0 0 ][0 1.50 ][0 0 ][0 -2.60 ][0 0 ][0 3.35 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
 DEAL::  h=0.50
-DEAL::    Quadrature point 0: [-36.00 0.00 0.00 0.00 ][62.35 13.86 0.00 0.00 ][-80.50 -53.67 0.00 0.00 ][12.00 0.00 0.00 0.00 ][-20.78 0.00 0.00 0.00 ][26.83 0.00 0.00 0.00 ][0.00 0.00 0.00 -36.00 ][0.00 0.00 13.86 62.35 ][0.00 0.00 -53.67 -80.50 ][0.00 0.00 0.00 12.00 ][0.00 0.00 0.00 -20.78 ][0.00 0.00 0.00 26.83 ][24.00 0.00 0.00 0.00 ][0.00 0.00 0.00 24.00 ][-69.28 0.00 0.00 0.00 ][0.00 0.00 0.00 -41.57 ][-41.57 0.00 0.00 0.00 ][0.00 0.00 0.00 53.67 ][120.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -69.28 ][53.67 0.00 0.00 0.00 ][0.00 0.00 0.00 120.00 ][-154.92 0.00 0.00 0.00 ][0.00 0.00 0.00 -154.92 ]
-DEAL::    Quadrature point 1: [-12.00 0.00 0.00 0.00 ][20.78 0.00 0.00 0.00 ][-26.83 0.00 0.00 0.00 ][36.00 0.00 0.00 0.00 ][-62.35 13.86 0.00 0.00 ][80.50 -53.67 0.00 0.00 ][0.00 0.00 0.00 -36.00 ][0.00 0.00 13.86 -62.35 ][0.00 0.00 53.67 -80.50 ][0.00 0.00 0.00 12.00 ][0.00 0.00 0.00 20.78 ][0.00 0.00 0.00 26.83 ][-24.00 0.00 0.00 0.00 ][0.00 0.00 0.00 24.00 ][-69.28 0.00 0.00 0.00 ][0.00 0.00 0.00 41.57 ][41.57 0.00 0.00 0.00 ][0.00 0.00 0.00 53.67 ][120.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -69.28 ][-53.67 0.00 0.00 0.00 ][0.00 0.00 0.00 -120.00 ][-154.92 0.00 0.00 0.00 ][0.00 0.00 0.00 -154.92 ]
-DEAL::    Quadrature point 2: [-36.00 0.00 0.00 0.00 ][-62.35 13.86 0.00 0.00 ][-80.50 53.67 0.00 0.00 ][12.00 0.00 0.00 0.00 ][20.78 0.00 0.00 0.00 ][26.83 0.00 0.00 0.00 ][0.00 0.00 0.00 -12.00 ][0.00 0.00 0.00 20.78 ][0.00 0.00 0.00 -26.83 ][0.00 0.00 0.00 36.00 ][0.00 0.00 13.86 -62.35 ][0.00 0.00 -53.67 80.50 ][24.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -24.00 ][-69.28 0.00 0.00 0.00 ][0.00 0.00 0.00 41.57 ][41.57 0.00 0.00 0.00 ][0.00 0.00 0.00 -53.67 ][-120.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -69.28 ][53.67 0.00 0.00 0.00 ][0.00 0.00 0.00 120.00 ][-154.92 0.00 0.00 0.00 ][0.00 0.00 0.00 -154.92 ]
-DEAL::    Quadrature point 3: [-12.00 0.00 0.00 0.00 ][-20.78 0.00 0.00 0.00 ][-26.83 0.00 0.00 0.00 ][36.00 0.00 0.00 0.00 ][62.35 13.86 0.00 0.00 ][80.50 53.67 0.00 0.00 ][0.00 0.00 0.00 -12.00 ][0.00 0.00 0.00 -20.78 ][0.00 0.00 0.00 -26.83 ][0.00 0.00 0.00 36.00 ][0.00 0.00 13.86 62.35 ][0.00 0.00 53.67 80.50 ][-24.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -24.00 ][-69.28 0.00 0.00 0.00 ][0.00 0.00 0.00 -41.57 ][-41.57 0.00 0.00 0.00 ][0.00 0.00 0.00 -53.67 ][-120.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -69.28 ][-53.67 0.00 0.00 0.00 ][0.00 0.00 0.00 -120.00 ][-154.92 0.00 0.00 0.00 ][0.00 0.00 0.00 -154.92 ]
+DEAL::    Quadrature point 0: [2.00 0 ][-3.46 0 ][4.47 0 ][0 0 ][0 0 ][0 0 ][0 2.00 ][0 -3.46 ][0 4.47 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL::    Quadrature point 1: [2.00 0 ][0 0 ][-2.24 0 ][0 0 ][0 0 ][0 0 ][0 -0.50 ][0 0.87 ][0 -1.12 ][0 -0.50 ][0 0.87 ][0 -1.12 ][0 0 ][0 3.00 ][0 0 ][0 -5.20 ][0 0 ][0 6.71 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
 DEAL::  h=0.25
-DEAL::    Quadrature point 0: [-144.00 0.00 0.00 0.00 ][249.42 55.43 0.00 0.00 ][-321.99 -214.66 0.00 0.00 ][48.00 0.00 0.00 0.00 ][-83.14 0.00 0.00 0.00 ][107.33 0.00 0.00 0.00 ][0.00 0.00 0.00 -144.00 ][0.00 0.00 55.43 249.42 ][0.00 0.00 -214.66 -321.99 ][0.00 0.00 0.00 48.00 ][0.00 0.00 0.00 -83.14 ][0.00 0.00 0.00 107.33 ][96.00 0.00 0.00 0.00 ][0.00 0.00 0.00 96.00 ][-277.13 0.00 0.00 0.00 ][0.00 0.00 0.00 -166.28 ][-166.28 0.00 0.00 0.00 ][0.00 0.00 0.00 214.66 ][480.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -277.13 ][214.66 0.00 0.00 0.00 ][0.00 0.00 0.00 480.00 ][-619.68 0.00 0.00 0.00 ][0.00 0.00 0.00 -619.68 ]
-DEAL::    Quadrature point 1: [-48.00 0.00 0.00 0.00 ][83.14 0.00 0.00 0.00 ][-107.33 0.00 0.00 0.00 ][144.00 0.00 0.00 0.00 ][-249.42 55.43 0.00 0.00 ][321.99 -214.66 0.00 0.00 ][0.00 0.00 0.00 -144.00 ][0.00 0.00 55.43 -249.42 ][0.00 0.00 214.66 -321.99 ][0.00 0.00 0.00 48.00 ][0.00 0.00 0.00 83.14 ][0.00 0.00 0.00 107.33 ][-96.00 0.00 0.00 0.00 ][0.00 0.00 0.00 96.00 ][-277.13 0.00 0.00 0.00 ][0.00 0.00 0.00 166.28 ][166.28 0.00 0.00 0.00 ][0.00 0.00 0.00 214.66 ][480.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -277.13 ][-214.66 0.00 0.00 0.00 ][0.00 0.00 0.00 -480.00 ][-619.68 0.00 0.00 0.00 ][0.00 0.00 0.00 -619.68 ]
-DEAL::    Quadrature point 2: [-144.00 0.00 0.00 0.00 ][-249.42 55.43 0.00 0.00 ][-321.99 214.66 0.00 0.00 ][48.00 0.00 0.00 0.00 ][83.14 0.00 0.00 0.00 ][107.33 0.00 0.00 0.00 ][0.00 0.00 0.00 -48.00 ][0.00 0.00 0.00 83.14 ][0.00 0.00 0.00 -107.33 ][0.00 0.00 0.00 144.00 ][0.00 0.00 55.43 -249.42 ][0.00 0.00 -214.66 321.99 ][96.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -96.00 ][-277.13 0.00 0.00 0.00 ][0.00 0.00 0.00 166.28 ][166.28 0.00 0.00 0.00 ][0.00 0.00 0.00 -214.66 ][-480.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -277.13 ][214.66 0.00 0.00 0.00 ][0.00 0.00 0.00 480.00 ][-619.68 0.00 0.00 0.00 ][0.00 0.00 0.00 -619.68 ]
-DEAL::    Quadrature point 3: [-48.00 0.00 0.00 0.00 ][-83.14 0.00 0.00 0.00 ][-107.33 0.00 0.00 0.00 ][144.00 0.00 0.00 0.00 ][249.42 55.43 0.00 0.00 ][321.99 214.66 0.00 0.00 ][0.00 0.00 0.00 -48.00 ][0.00 0.00 0.00 -83.14 ][0.00 0.00 0.00 -107.33 ][0.00 0.00 0.00 144.00 ][0.00 0.00 55.43 249.42 ][0.00 0.00 214.66 321.99 ][-96.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -96.00 ][-277.13 0.00 0.00 0.00 ][0.00 0.00 0.00 -166.28 ][-166.28 0.00 0.00 0.00 ][0.00 0.00 0.00 -214.66 ][-480.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -277.13 ][-214.66 0.00 0.00 0.00 ][0.00 0.00 0.00 -480.00 ][-619.68 0.00 0.00 0.00 ][0.00 0.00 0.00 -619.68 ]
+DEAL::    Quadrature point 0: [4.00 0 ][-6.93 0 ][8.94 0 ][0 0 ][0 0 ][0 0 ][0 4.00 ][0 -6.93 ][0 8.94 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL::    Quadrature point 1: [4.00 0 ][0 0 ][-4.47 0 ][0 0 ][0 0 ][0 0 ][0 -1.00 ][0 1.73 ][0 -2.24 ][0 -1.00 ][0 1.73 ][0 -2.24 ][0 0 ][0 6.00 ][0 0 ][0 -10.39 ][0 0 ][0 13.42 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
 DEAL::  h=0.12
-DEAL::    Quadrature point 0: [-576.00 0.00 0.00 0.00 ][997.66 221.70 0.00 0.00 ][-1287.98 -858.65 0.00 0.00 ][192.00 0.00 0.00 0.00 ][-332.55 0.00 0.00 0.00 ][429.33 0.00 0.00 0.00 ][0.00 0.00 0.00 -576.00 ][0.00 0.00 221.70 997.66 ][0.00 0.00 -858.65 -1287.98 ][0.00 0.00 0.00 192.00 ][0.00 0.00 0.00 -332.55 ][0.00 0.00 0.00 429.33 ][384.00 0.00 0.00 0.00 ][0.00 0.00 0.00 384.00 ][-1108.51 0.00 0.00 0.00 ][0.00 0.00 0.00 -665.11 ][-665.11 0.00 0.00 0.00 ][0.00 0.00 0.00 858.65 ][1920.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -1108.51 ][858.65 0.00 0.00 0.00 ][0.00 0.00 0.00 1920.00 ][-2478.71 0.00 0.00 0.00 ][0.00 0.00 0.00 -2478.71 ]
-DEAL::    Quadrature point 1: [-192.00 0.00 0.00 0.00 ][332.55 0.00 0.00 0.00 ][-429.33 0.00 0.00 0.00 ][576.00 0.00 0.00 0.00 ][-997.66 221.70 0.00 0.00 ][1287.98 -858.65 0.00 0.00 ][0.00 0.00 0.00 -576.00 ][0.00 0.00 221.70 -997.66 ][0.00 0.00 858.65 -1287.98 ][0.00 0.00 0.00 192.00 ][0.00 0.00 0.00 332.55 ][0.00 0.00 0.00 429.33 ][-384.00 0.00 0.00 0.00 ][0.00 0.00 0.00 384.00 ][-1108.51 0.00 0.00 0.00 ][0.00 0.00 0.00 665.11 ][665.11 0.00 0.00 0.00 ][0.00 0.00 0.00 858.65 ][1920.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -1108.51 ][-858.65 0.00 0.00 0.00 ][0.00 0.00 0.00 -1920.00 ][-2478.71 0.00 0.00 0.00 ][0.00 0.00 0.00 -2478.71 ]
-DEAL::    Quadrature point 2: [-576.00 0.00 0.00 0.00 ][-997.66 221.70 0.00 0.00 ][-1287.98 858.65 0.00 0.00 ][192.00 0.00 0.00 0.00 ][332.55 0.00 0.00 0.00 ][429.33 0.00 0.00 0.00 ][0.00 0.00 0.00 -192.00 ][0.00 0.00 0.00 332.55 ][0.00 0.00 0.00 -429.33 ][0.00 0.00 0.00 576.00 ][0.00 0.00 221.70 -997.66 ][0.00 0.00 -858.65 1287.98 ][384.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -384.00 ][-1108.51 0.00 0.00 0.00 ][0.00 0.00 0.00 665.11 ][665.11 0.00 0.00 0.00 ][0.00 0.00 0.00 -858.65 ][-1920.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -1108.51 ][858.65 0.00 0.00 0.00 ][0.00 0.00 0.00 1920.00 ][-2478.71 0.00 0.00 0.00 ][0.00 0.00 0.00 -2478.71 ]
-DEAL::    Quadrature point 3: [-192.00 0.00 0.00 0.00 ][-332.55 0.00 0.00 0.00 ][-429.33 0.00 0.00 0.00 ][576.00 0.00 0.00 0.00 ][997.66 221.70 0.00 0.00 ][1287.98 858.65 0.00 0.00 ][0.00 0.00 0.00 -192.00 ][0.00 0.00 0.00 -332.55 ][0.00 0.00 0.00 -429.33 ][0.00 0.00 0.00 576.00 ][0.00 0.00 221.70 997.66 ][0.00 0.00 858.65 1287.98 ][-384.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -384.00 ][-1108.51 0.00 0.00 0.00 ][0.00 0.00 0.00 -665.11 ][-665.11 0.00 0.00 0.00 ][0.00 0.00 0.00 -858.65 ][-1920.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -1108.51 ][-858.65 0.00 0.00 0.00 ][0.00 0.00 0.00 -1920.00 ][-2478.71 0.00 0.00 0.00 ][0.00 0.00 0.00 -2478.71 ]
+DEAL::    Quadrature point 0: [8.00 0 ][-13.86 0 ][17.89 0 ][0 0 ][0 0 ][0 0 ][0 8.00 ][0 -13.86 ][0 17.89 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL::    Quadrature point 1: [8.00 0 ][0 0 ][-8.94 0 ][0 0 ][0 0 ][0 0 ][0 -2.00 ][0 3.46 ][0 -4.47 ][0 -2.00 ][0 3.46 ][0 -4.47 ][0 0 ][0 12.00 ][0 0 ][0 -20.78 ][0 0 ][0 26.83 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
 DEAL::  h=0.06
-DEAL::    Quadrature point 0: [-2304.00 0.00 0.00 0.00 ][3990.65 886.81 0.00 0.00 ][-5151.90 -3434.60 0.00 0.00 ][768.00 0.00 0.00 0.00 ][-1330.22 0.00 0.00 0.00 ][1717.30 0.00 0.00 0.00 ][0.00 0.00 0.00 -2304.00 ][0.00 0.00 886.81 3990.65 ][0.00 0.00 -3434.60 -5151.90 ][0.00 0.00 0.00 768.00 ][0.00 0.00 0.00 -1330.22 ][0.00 0.00 0.00 1717.30 ][1536.00 0.00 0.00 0.00 ][0.00 0.00 0.00 1536.00 ][-4434.05 0.00 0.00 0.00 ][0.00 0.00 0.00 -2660.43 ][-2660.43 0.00 0.00 0.00 ][0.00 0.00 0.00 3434.60 ][7680.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -4434.05 ][3434.60 0.00 0.00 0.00 ][0.00 0.00 0.00 7680.00 ][-9914.84 0.00 0.00 0.00 ][0.00 0.00 0.00 -9914.84 ]
-DEAL::    Quadrature point 1: [-768.00 0.00 0.00 0.00 ][1330.22 0.00 0.00 0.00 ][-1717.30 0.00 0.00 0.00 ][2304.00 0.00 0.00 0.00 ][-3990.65 886.81 0.00 0.00 ][5151.90 -3434.60 0.00 0.00 ][0.00 0.00 0.00 -2304.00 ][0.00 0.00 886.81 -3990.65 ][0.00 0.00 3434.60 -5151.90 ][0.00 0.00 0.00 768.00 ][0.00 0.00 0.00 1330.22 ][0.00 0.00 0.00 1717.30 ][-1536.00 0.00 0.00 0.00 ][0.00 0.00 0.00 1536.00 ][-4434.05 0.00 0.00 0.00 ][0.00 0.00 0.00 2660.43 ][2660.43 0.00 0.00 0.00 ][0.00 0.00 0.00 3434.60 ][7680.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -4434.05 ][-3434.60 0.00 0.00 0.00 ][0.00 0.00 0.00 -7680.00 ][-9914.84 0.00 0.00 0.00 ][0.00 0.00 0.00 -9914.84 ]
-DEAL::    Quadrature point 2: [-2304.00 0.00 0.00 0.00 ][-3990.65 886.81 0.00 0.00 ][-5151.90 3434.60 0.00 0.00 ][768.00 0.00 0.00 0.00 ][1330.22 0.00 0.00 0.00 ][1717.30 0.00 0.00 0.00 ][0.00 0.00 0.00 -768.00 ][0.00 0.00 0.00 1330.22 ][0.00 0.00 0.00 -1717.30 ][0.00 0.00 0.00 2304.00 ][0.00 0.00 886.81 -3990.65 ][0.00 0.00 -3434.60 5151.90 ][1536.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -1536.00 ][-4434.05 0.00 0.00 0.00 ][0.00 0.00 0.00 2660.43 ][2660.43 0.00 0.00 0.00 ][0.00 0.00 0.00 -3434.60 ][-7680.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -4434.05 ][3434.60 0.00 0.00 0.00 ][0.00 0.00 0.00 7680.00 ][-9914.84 0.00 0.00 0.00 ][0.00 0.00 0.00 -9914.84 ]
-DEAL::    Quadrature point 3: [-768.00 0.00 0.00 0.00 ][-1330.22 0.00 0.00 0.00 ][-1717.30 0.00 0.00 0.00 ][2304.00 0.00 0.00 0.00 ][3990.65 886.81 0.00 0.00 ][5151.90 3434.60 0.00 0.00 ][0.00 0.00 0.00 -768.00 ][0.00 0.00 0.00 -1330.22 ][0.00 0.00 0.00 -1717.30 ][0.00 0.00 0.00 2304.00 ][0.00 0.00 886.81 3990.65 ][0.00 0.00 3434.60 5151.90 ][-1536.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -1536.00 ][-4434.05 0.00 0.00 0.00 ][0.00 0.00 0.00 -2660.43 ][-2660.43 0.00 0.00 0.00 ][0.00 0.00 0.00 -3434.60 ][-7680.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -4434.05 ][-3434.60 0.00 0.00 0.00 ][0.00 0.00 0.00 -7680.00 ][-9914.84 0.00 0.00 0.00 ][0.00 0.00 0.00 -9914.84 ]
+DEAL::    Quadrature point 0: [16.00 0 ][-27.71 0 ][35.78 0 ][0 0 ][0 0 ][0 0 ][0 16.00 ][0 -27.71 ][0 35.78 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL::    Quadrature point 1: [16.00 0 ][0 0 ][-17.89 0 ][0 0 ][0 0 ][0 0 ][0 -4.00 ][0 6.93 ][0 -8.94 ][0 -4.00 ][0 6.93 ][0 -8.94 ][0 0 ][0 24.00 ][0 0 ][0 -41.57 ][0 0 ][0 53.67 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
 DEAL::  h=0.03
-DEAL::    Quadrature point 0: [-9216.00 0.00 0.00 0.00 ][15962.58 3547.24 0.00 0.00 ][-20607.60 -13738.40 0.00 0.00 ][3072.00 0.00 0.00 0.00 ][-5320.86 0.00 0.00 0.00 ][6869.20 0.00 0.00 0.00 ][0.00 0.00 0.00 -9216.00 ][0.00 0.00 3547.24 15962.58 ][0.00 0.00 -13738.40 -20607.60 ][0.00 0.00 0.00 3072.00 ][0.00 0.00 0.00 -5320.86 ][0.00 0.00 0.00 6869.20 ][6144.00 0.00 0.00 0.00 ][0.00 0.00 0.00 6144.00 ][-17736.20 0.00 0.00 0.00 ][0.00 0.00 0.00 -10641.72 ][-10641.72 0.00 0.00 0.00 ][0.00 0.00 0.00 13738.40 ][30720.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -17736.20 ][13738.40 0.00 0.00 0.00 ][0.00 0.00 0.00 30720.00 ][-39659.35 0.00 0.00 0.00 ][0.00 0.00 0.00 -39659.35 ]
-DEAL::    Quadrature point 1: [-3072.00 0.00 0.00 0.00 ][5320.86 0.00 0.00 0.00 ][-6869.20 0.00 0.00 0.00 ][9216.00 0.00 0.00 0.00 ][-15962.58 3547.24 0.00 0.00 ][20607.60 -13738.40 0.00 0.00 ][0.00 0.00 0.00 -9216.00 ][0.00 0.00 3547.24 -15962.58 ][0.00 0.00 13738.40 -20607.60 ][0.00 0.00 0.00 3072.00 ][0.00 0.00 0.00 5320.86 ][0.00 0.00 0.00 6869.20 ][-6144.00 0.00 0.00 0.00 ][0.00 0.00 0.00 6144.00 ][-17736.20 0.00 0.00 0.00 ][0.00 0.00 0.00 10641.72 ][10641.72 0.00 0.00 0.00 ][0.00 0.00 0.00 13738.40 ][30720.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -17736.20 ][-13738.40 0.00 0.00 0.00 ][0.00 0.00 0.00 -30720.00 ][-39659.35 0.00 0.00 0.00 ][0.00 0.00 0.00 -39659.35 ]
-DEAL::    Quadrature point 2: [-9216.00 0.00 0.00 0.00 ][-15962.58 3547.24 0.00 0.00 ][-20607.60 13738.40 0.00 0.00 ][3072.00 0.00 0.00 0.00 ][5320.86 0.00 0.00 0.00 ][6869.20 0.00 0.00 0.00 ][0.00 0.00 0.00 -3072.00 ][0.00 0.00 0.00 5320.86 ][0.00 0.00 0.00 -6869.20 ][0.00 0.00 0.00 9216.00 ][0.00 0.00 3547.24 -15962.58 ][0.00 0.00 -13738.40 20607.60 ][6144.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -6144.00 ][-17736.20 0.00 0.00 0.00 ][0.00 0.00 0.00 10641.72 ][10641.72 0.00 0.00 0.00 ][0.00 0.00 0.00 -13738.40 ][-30720.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -17736.20 ][13738.40 0.00 0.00 0.00 ][0.00 0.00 0.00 30720.00 ][-39659.35 0.00 0.00 0.00 ][0.00 0.00 0.00 -39659.35 ]
-DEAL::    Quadrature point 3: [-3072.00 0.00 0.00 0.00 ][-5320.86 0.00 0.00 0.00 ][-6869.20 0.00 0.00 0.00 ][9216.00 0.00 0.00 0.00 ][15962.58 3547.24 0.00 0.00 ][20607.60 13738.40 0.00 0.00 ][0.00 0.00 0.00 -3072.00 ][0.00 0.00 0.00 -5320.86 ][0.00 0.00 0.00 -6869.20 ][0.00 0.00 0.00 9216.00 ][0.00 0.00 3547.24 15962.58 ][0.00 0.00 13738.40 20607.60 ][-6144.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -6144.00 ][-17736.20 0.00 0.00 0.00 ][0.00 0.00 0.00 -10641.72 ][-10641.72 0.00 0.00 0.00 ][0.00 0.00 0.00 -13738.40 ][-30720.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -17736.20 ][-13738.40 0.00 0.00 0.00 ][0.00 0.00 0.00 -30720.00 ][-39659.35 0.00 0.00 0.00 ][0.00 0.00 0.00 -39659.35 ]
+DEAL::    Quadrature point 0: [32.00 0 ][-55.43 0 ][71.55 0 ][0 0 ][0 0 ][0 0 ][0 32.00 ][0 -55.43 ][0 71.55 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL::    Quadrature point 1: [32.00 0 ][0 0 ][-35.78 0 ][0 0 ][0 0 ][0 0 ][0 -8.00 ][0 13.86 ][0 -17.89 ][0 -8.00 ][0 13.86 ][0 -17.89 ][0 0 ][0 48.00 ][0 0 ][0 -83.14 ][0 0 ][0 107.33 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
 DEAL::  h=0.02
-DEAL::    Quadrature point 0: [-36864.00 0.00 0.00 0.00 ][63850.32 14188.96 0.00 0.00 ][-82430.41 -54953.61 0.00 0.00 ][12288.00 0.00 0.00 0.00 ][-21283.44 0.00 0.00 0.00 ][27476.80 0.00 0.00 0.00 ][0.00 0.00 0.00 -36864.00 ][0.00 0.00 14188.96 63850.32 ][0.00 0.00 -54953.61 -82430.41 ][0.00 0.00 0.00 12288.00 ][0.00 0.00 0.00 -21283.44 ][0.00 0.00 0.00 27476.80 ][24576.00 0.00 0.00 0.00 ][0.00 0.00 0.00 24576.00 ][-70944.80 0.00 0.00 0.00 ][0.00 0.00 0.00 -42566.88 ][-42566.88 0.00 0.00 0.00 ][0.00 0.00 0.00 54953.61 ][122880.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -70944.80 ][54953.61 0.00 0.00 0.00 ][0.00 0.00 0.00 122880.00 ][-158637.40 0.00 0.00 0.00 ][0.00 0.00 0.00 -158637.40 ]
-DEAL::    Quadrature point 1: [-12288.00 0.00 0.00 0.00 ][21283.44 0.00 0.00 0.00 ][-27476.80 0.00 0.00 0.00 ][36864.00 0.00 0.00 0.00 ][-63850.32 14188.96 0.00 0.00 ][82430.41 -54953.61 0.00 0.00 ][0.00 0.00 0.00 -36864.00 ][0.00 0.00 14188.96 -63850.32 ][0.00 0.00 54953.61 -82430.41 ][0.00 0.00 0.00 12288.00 ][0.00 0.00 0.00 21283.44 ][0.00 0.00 0.00 27476.80 ][-24576.00 0.00 0.00 0.00 ][0.00 0.00 0.00 24576.00 ][-70944.80 0.00 0.00 0.00 ][0.00 0.00 0.00 42566.88 ][42566.88 0.00 0.00 0.00 ][0.00 0.00 0.00 54953.61 ][122880.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -70944.80 ][-54953.61 0.00 0.00 0.00 ][0.00 0.00 0.00 -122880.00 ][-158637.40 0.00 0.00 0.00 ][0.00 0.00 0.00 -158637.40 ]
-DEAL::    Quadrature point 2: [-36864.00 0.00 0.00 0.00 ][-63850.32 14188.96 0.00 0.00 ][-82430.41 54953.61 0.00 0.00 ][12288.00 0.00 0.00 0.00 ][21283.44 0.00 0.00 0.00 ][27476.80 0.00 0.00 0.00 ][0.00 0.00 0.00 -12288.00 ][0.00 0.00 0.00 21283.44 ][0.00 0.00 0.00 -27476.80 ][0.00 0.00 0.00 36864.00 ][0.00 0.00 14188.96 -63850.32 ][0.00 0.00 -54953.61 82430.41 ][24576.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -24576.00 ][-70944.80 0.00 0.00 0.00 ][0.00 0.00 0.00 42566.88 ][42566.88 0.00 0.00 0.00 ][0.00 0.00 0.00 -54953.61 ][-122880.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -70944.80 ][54953.61 0.00 0.00 0.00 ][0.00 0.00 0.00 122880.00 ][-158637.40 0.00 0.00 0.00 ][0.00 0.00 0.00 -158637.40 ]
-DEAL::    Quadrature point 3: [-12288.00 0.00 0.00 0.00 ][-21283.44 0.00 0.00 0.00 ][-27476.80 0.00 0.00 0.00 ][36864.00 0.00 0.00 0.00 ][63850.32 14188.96 0.00 0.00 ][82430.41 54953.61 0.00 0.00 ][0.00 0.00 0.00 -12288.00 ][0.00 0.00 0.00 -21283.44 ][0.00 0.00 0.00 -27476.80 ][0.00 0.00 0.00 36864.00 ][0.00 0.00 14188.96 63850.32 ][0.00 0.00 54953.61 82430.41 ][-24576.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -24576.00 ][-70944.80 0.00 0.00 0.00 ][0.00 0.00 0.00 -42566.88 ][-42566.88 0.00 0.00 0.00 ][0.00 0.00 0.00 -54953.61 ][-122880.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -70944.80 ][-54953.61 0.00 0.00 0.00 ][0.00 0.00 0.00 -122880.00 ][-158637.40 0.00 0.00 0.00 ][0.00 0.00 0.00 -158637.40 ]
+DEAL::    Quadrature point 0: [64.00 0 ][-110.85 0 ][143.11 0 ][0 0 ][0 0 ][0 0 ][0 64.00 ][0 -110.85 ][0 143.11 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL::    Quadrature point 1: [64.00 0 ][0 0 ][-71.55 0 ][0 0 ][0 0 ][0 0 ][0 -16.00 ][0 27.71 ][0 -35.78 ][0 -16.00 ][0 27.71 ][0 -35.78 ][0 0 ][0 96.00 ][0 0 ][0 -166.28 ][0 0 ][0 214.66 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
 DEAL::Degree=3
 DEAL::  h=1.00
-DEAL::    Quadrature point 0: [-16.00 0.00 0.00 0.00 ][27.71 3.46 0.00 0.00 ][-35.78 -13.42 0.00 0.00 ][42.33 31.75 0.00 0.00 ][-4.00 0.00 0.00 0.00 ][6.93 0.00 0.00 0.00 ][-8.94 0.00 0.00 0.00 ][10.58 0.00 0.00 0.00 ][0.00 0.00 0.00 -16.00 ][0.00 0.00 3.46 27.71 ][0.00 0.00 -13.42 -35.78 ][0.00 0.00 31.75 42.33 ][0.00 0.00 0.00 -4.00 ][0.00 0.00 0.00 6.93 ][0.00 0.00 0.00 -8.94 ][0.00 0.00 0.00 10.58 ][20.00 0.00 0.00 0.00 ][0.00 0.00 0.00 20.00 ][-17.32 0.00 0.00 0.00 ][0.00 0.00 0.00 -34.64 ][31.30 0.00 0.00 0.00 ][0.00 0.00 0.00 44.72 ][-34.64 0.00 0.00 0.00 ][0.00 0.00 0.00 -52.92 ][30.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -17.32 ][-54.22 0.00 0.00 0.00 ][0.00 0.00 0.00 30.00 ][44.72 0.00 0.00 0.00 ][0.00 0.00 0.00 -38.73 ][-38.73 0.00 0.00 0.00 ][0.00 0.00 0.00 45.83 ][70.00 0.00 0.00 0.00 ][0.00 0.00 0.00 31.30 ][-52.92 0.00 0.00 0.00 ][0.00 0.00 0.00 -54.22 ][45.83 0.00 0.00 0.00 ][0.00 0.00 0.00 70.00 ][-82.83 0.00 0.00 0.00 ][0.00 0.00 0.00 -82.83 ]
-DEAL::    Quadrature point 1: [4.00 0.00 0.00 0.00 ][-6.93 0.00 0.00 0.00 ][8.94 0.00 0.00 0.00 ][-10.58 0.00 0.00 0.00 ][16.00 0.00 0.00 0.00 ][-27.71 3.46 0.00 0.00 ][35.78 -13.42 0.00 0.00 ][-42.33 31.75 0.00 0.00 ][0.00 0.00 0.00 -16.00 ][0.00 0.00 3.46 -27.71 ][0.00 0.00 13.42 -35.78 ][0.00 0.00 31.75 -42.33 ][0.00 0.00 0.00 -4.00 ][0.00 0.00 0.00 -6.93 ][0.00 0.00 0.00 -8.94 ][0.00 0.00 0.00 -10.58 ][-20.00 0.00 0.00 0.00 ][0.00 0.00 0.00 20.00 ][-17.32 0.00 0.00 0.00 ][0.00 0.00 0.00 34.64 ][-31.30 0.00 0.00 0.00 ][0.00 0.00 0.00 44.72 ][34.64 0.00 0.00 0.00 ][0.00 0.00 0.00 52.92 ][30.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -17.32 ][54.22 0.00 0.00 0.00 ][0.00 0.00 0.00 -30.00 ][-44.72 0.00 0.00 0.00 ][0.00 0.00 0.00 -38.73 ][-38.73 0.00 0.00 0.00 ][0.00 0.00 0.00 -45.83 ][-70.00 0.00 0.00 0.00 ][0.00 0.00 0.00 31.30 ][52.92 0.00 0.00 0.00 ][0.00 0.00 0.00 54.22 ][45.83 0.00 0.00 0.00 ][0.00 0.00 0.00 70.00 ][82.83 0.00 0.00 0.00 ][0.00 0.00 0.00 82.83 ]
-DEAL::    Quadrature point 2: [-16.00 0.00 0.00 0.00 ][-27.71 3.46 0.00 0.00 ][-35.78 13.42 0.00 0.00 ][-42.33 31.75 0.00 0.00 ][-4.00 0.00 0.00 0.00 ][-6.93 0.00 0.00 0.00 ][-8.94 0.00 0.00 0.00 ][-10.58 0.00 0.00 0.00 ][0.00 0.00 0.00 4.00 ][0.00 0.00 0.00 -6.93 ][0.00 0.00 0.00 8.94 ][0.00 0.00 0.00 -10.58 ][0.00 0.00 0.00 16.00 ][0.00 0.00 3.46 -27.71 ][0.00 0.00 -13.42 35.78 ][0.00 0.00 31.75 -42.33 ][20.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -20.00 ][-17.32 0.00 0.00 0.00 ][0.00 0.00 0.00 34.64 ][31.30 0.00 0.00 0.00 ][0.00 0.00 0.00 -44.72 ][34.64 0.00 0.00 0.00 ][0.00 0.00 0.00 52.92 ][-30.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -17.32 ][54.22 0.00 0.00 0.00 ][0.00 0.00 0.00 30.00 ][44.72 0.00 0.00 0.00 ][0.00 0.00 0.00 -38.73 ][-38.73 0.00 0.00 0.00 ][0.00 0.00 0.00 45.83 ][70.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -31.30 ][52.92 0.00 0.00 0.00 ][0.00 0.00 0.00 54.22 ][-45.83 0.00 0.00 0.00 ][0.00 0.00 0.00 -70.00 ][82.83 0.00 0.00 0.00 ][0.00 0.00 0.00 82.83 ]
-DEAL::    Quadrature point 3: [4.00 0.00 0.00 0.00 ][6.93 0.00 0.00 0.00 ][8.94 0.00 0.00 0.00 ][10.58 0.00 0.00 0.00 ][16.00 0.00 0.00 0.00 ][27.71 3.46 0.00 0.00 ][35.78 13.42 0.00 0.00 ][42.33 31.75 0.00 0.00 ][0.00 0.00 0.00 4.00 ][0.00 0.00 0.00 6.93 ][0.00 0.00 0.00 8.94 ][0.00 0.00 0.00 10.58 ][0.00 0.00 0.00 16.00 ][0.00 0.00 3.46 27.71 ][0.00 0.00 13.42 35.78 ][0.00 0.00 31.75 42.33 ][-20.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -20.00 ][-17.32 0.00 0.00 0.00 ][0.00 0.00 0.00 -34.64 ][-31.30 0.00 0.00 0.00 ][0.00 0.00 0.00 -44.72 ][-34.64 0.00 0.00 0.00 ][0.00 0.00 0.00 -52.92 ][-30.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -17.32 ][-54.22 0.00 0.00 0.00 ][0.00 0.00 0.00 -30.00 ][-44.72 0.00 0.00 0.00 ][0.00 0.00 0.00 -38.73 ][-38.73 0.00 0.00 0.00 ][0.00 0.00 0.00 -45.83 ][-70.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -31.30 ][-52.92 0.00 0.00 0.00 ][0.00 0.00 0.00 -54.22 ][-45.83 0.00 0.00 0.00 ][0.00 0.00 0.00 -70.00 ][-82.83 0.00 0.00 0.00 ][0.00 0.00 0.00 -82.83 ]
+DEAL::    Quadrature point 0: [1.00 0 ][-1.73 0 ][2.24 0 ][-2.65 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 1.00 ][0 -1.73 ][0 2.24 ][0 -2.65 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL::    Quadrature point 1: [1.00 0 ][0 0 ][-1.12 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0.19 ][0 -0.32 ][0 0.42 ][0 -0.50 ][0 0.19 ][0 -0.32 ][0 0.42 ][0 -0.50 ][0 0 ][0 0.63 ][0 0 ][0 -1.08 ][0 0 ][0 1.40 ][0 0 ][0 -1.65 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 -1.96 ][0 0 ][0 3.39 ][0 0 ][0 -4.38 ][0 0 ][0 5.18 ]
 DEAL::  h=0.50
-DEAL::    Quadrature point 0: [-64.00 0.00 0.00 0.00 ][110.85 13.86 0.00 0.00 ][-143.11 -53.67 0.00 0.00 ][169.33 127.00 0.00 0.00 ][-16.00 0.00 0.00 0.00 ][27.71 0.00 0.00 0.00 ][-35.78 0.00 0.00 0.00 ][42.33 0.00 0.00 0.00 ][0.00 0.00 0.00 -64.00 ][0.00 0.00 13.86 110.85 ][0.00 0.00 -53.67 -143.11 ][0.00 0.00 127.00 169.33 ][0.00 0.00 0.00 -16.00 ][0.00 0.00 0.00 27.71 ][0.00 0.00 0.00 -35.78 ][0.00 0.00 0.00 42.33 ][80.00 0.00 0.00 0.00 ][0.00 0.00 0.00 80.00 ][-69.28 0.00 0.00 0.00 ][0.00 0.00 0.00 -138.56 ][125.22 0.00 0.00 0.00 ][0.00 0.00 0.00 178.89 ][-138.56 0.00 0.00 0.00 ][0.00 0.00 0.00 -211.66 ][120.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -69.28 ][-216.89 0.00 0.00 0.00 ][0.00 0.00 0.00 120.00 ][178.89 0.00 0.00 0.00 ][0.00 0.00 0.00 -154.92 ][-154.92 0.00 0.00 0.00 ][0.00 0.00 0.00 183.30 ][280.00 0.00 0.00 0.00 ][0.00 0.00 0.00 125.22 ][-211.66 0.00 0.00 0.00 ][0.00 0.00 0.00 -216.89 ][183.30 0.00 0.00 0.00 ][0.00 0.00 0.00 280.00 ][-331.30 0.00 0.00 0.00 ][0.00 0.00 0.00 -331.30 ]
-DEAL::    Quadrature point 1: [16.00 0.00 0.00 0.00 ][-27.71 0.00 0.00 0.00 ][35.78 0.00 0.00 0.00 ][-42.33 0.00 0.00 0.00 ][64.00 0.00 0.00 0.00 ][-110.85 13.86 0.00 0.00 ][143.11 -53.67 0.00 0.00 ][-169.33 127.00 0.00 0.00 ][0.00 0.00 0.00 -64.00 ][0.00 0.00 13.86 -110.85 ][0.00 0.00 53.67 -143.11 ][0.00 0.00 127.00 -169.33 ][0.00 0.00 0.00 -16.00 ][0.00 0.00 0.00 -27.71 ][0.00 0.00 0.00 -35.78 ][0.00 0.00 0.00 -42.33 ][-80.00 0.00 0.00 0.00 ][0.00 0.00 0.00 80.00 ][-69.28 0.00 0.00 0.00 ][0.00 0.00 0.00 138.56 ][-125.22 0.00 0.00 0.00 ][0.00 0.00 0.00 178.89 ][138.56 0.00 0.00 0.00 ][0.00 0.00 0.00 211.66 ][120.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -69.28 ][216.89 0.00 0.00 0.00 ][0.00 0.00 0.00 -120.00 ][-178.89 0.00 0.00 0.00 ][0.00 0.00 0.00 -154.92 ][-154.92 0.00 0.00 0.00 ][0.00 0.00 0.00 -183.30 ][-280.00 0.00 0.00 0.00 ][0.00 0.00 0.00 125.22 ][211.66 0.00 0.00 0.00 ][0.00 0.00 0.00 216.89 ][183.30 0.00 0.00 0.00 ][0.00 0.00 0.00 280.00 ][331.30 0.00 0.00 0.00 ][0.00 0.00 0.00 331.30 ]
-DEAL::    Quadrature point 2: [-64.00 0.00 0.00 0.00 ][-110.85 13.86 0.00 0.00 ][-143.11 53.67 0.00 0.00 ][-169.33 127.00 0.00 0.00 ][-16.00 0.00 0.00 0.00 ][-27.71 0.00 0.00 0.00 ][-35.78 0.00 0.00 0.00 ][-42.33 0.00 0.00 0.00 ][0.00 0.00 0.00 16.00 ][0.00 0.00 0.00 -27.71 ][0.00 0.00 0.00 35.78 ][0.00 0.00 0.00 -42.33 ][0.00 0.00 0.00 64.00 ][0.00 0.00 13.86 -110.85 ][0.00 0.00 -53.67 143.11 ][0.00 0.00 127.00 -169.33 ][80.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -80.00 ][-69.28 0.00 0.00 0.00 ][0.00 0.00 0.00 138.56 ][125.22 0.00 0.00 0.00 ][0.00 0.00 0.00 -178.89 ][138.56 0.00 0.00 0.00 ][0.00 0.00 0.00 211.66 ][-120.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -69.28 ][216.89 0.00 0.00 0.00 ][0.00 0.00 0.00 120.00 ][178.89 0.00 0.00 0.00 ][0.00 0.00 0.00 -154.92 ][-154.92 0.00 0.00 0.00 ][0.00 0.00 0.00 183.30 ][280.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -125.22 ][211.66 0.00 0.00 0.00 ][0.00 0.00 0.00 216.89 ][-183.30 0.00 0.00 0.00 ][0.00 0.00 0.00 -280.00 ][331.30 0.00 0.00 0.00 ][0.00 0.00 0.00 331.30 ]
-DEAL::    Quadrature point 3: [16.00 0.00 0.00 0.00 ][27.71 0.00 0.00 0.00 ][35.78 0.00 0.00 0.00 ][42.33 0.00 0.00 0.00 ][64.00 0.00 0.00 0.00 ][110.85 13.86 0.00 0.00 ][143.11 53.67 0.00 0.00 ][169.33 127.00 0.00 0.00 ][0.00 0.00 0.00 16.00 ][0.00 0.00 0.00 27.71 ][0.00 0.00 0.00 35.78 ][0.00 0.00 0.00 42.33 ][0.00 0.00 0.00 64.00 ][0.00 0.00 13.86 110.85 ][0.00 0.00 53.67 143.11 ][0.00 0.00 127.00 169.33 ][-80.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -80.00 ][-69.28 0.00 0.00 0.00 ][0.00 0.00 0.00 -138.56 ][-125.22 0.00 0.00 0.00 ][0.00 0.00 0.00 -178.89 ][-138.56 0.00 0.00 0.00 ][0.00 0.00 0.00 -211.66 ][-120.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -69.28 ][-216.89 0.00 0.00 0.00 ][0.00 0.00 0.00 -120.00 ][-178.89 0.00 0.00 0.00 ][0.00 0.00 0.00 -154.92 ][-154.92 0.00 0.00 0.00 ][0.00 0.00 0.00 -183.30 ][-280.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -125.22 ][-211.66 0.00 0.00 0.00 ][0.00 0.00 0.00 -216.89 ][-183.30 0.00 0.00 0.00 ][0.00 0.00 0.00 -280.00 ][-331.30 0.00 0.00 0.00 ][0.00 0.00 0.00 -331.30 ]
+DEAL::    Quadrature point 0: [2.00 0 ][-3.46 0 ][4.47 0 ][-5.29 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 2.00 ][0 -3.46 ][0 4.47 ][0 -5.29 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL::    Quadrature point 1: [2.00 0 ][0 0 ][-2.24 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0.38 ][0 -0.65 ][0 0.84 ][0 -0.99 ][0 0.38 ][0 -0.65 ][0 0.84 ][0 -0.99 ][0 0 ][0 1.25 ][0 0 ][0 -2.17 ][0 0 ][0 2.80 ][0 0 ][0 -3.31 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 -3.91 ][0 0 ][0 6.78 ][0 0 ][0 -8.75 ][0 0 ][0 10.35 ]
 DEAL::  h=0.25
-DEAL::    Quadrature point 0: [-256.00 0.00 0.00 0.00 ][443.41 55.43 0.00 0.00 ][-572.43 -214.66 0.00 0.00 ][677.31 507.98 0.00 0.00 ][-64.00 0.00 0.00 0.00 ][110.85 0.00 0.00 0.00 ][-143.11 0.00 0.00 0.00 ][169.33 0.00 0.00 0.00 ][0.00 0.00 0.00 -256.00 ][0.00 0.00 55.43 443.41 ][0.00 0.00 -214.66 -572.43 ][0.00 0.00 507.98 677.31 ][0.00 0.00 0.00 -64.00 ][0.00 0.00 0.00 110.85 ][0.00 0.00 0.00 -143.11 ][0.00 0.00 0.00 169.33 ][320.00 0.00 0.00 0.00 ][0.00 0.00 0.00 320.00 ][-277.13 0.00 0.00 0.00 ][0.00 0.00 0.00 -554.26 ][500.88 0.00 0.00 0.00 ][0.00 0.00 0.00 715.54 ][-554.26 0.00 0.00 0.00 ][0.00 0.00 0.00 -846.64 ][480.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -277.13 ][-867.55 0.00 0.00 0.00 ][0.00 0.00 0.00 480.00 ][715.54 0.00 0.00 0.00 ][0.00 0.00 0.00 -619.68 ][-619.68 0.00 0.00 0.00 ][0.00 0.00 0.00 733.21 ][1120.00 0.00 0.00 0.00 ][0.00 0.00 0.00 500.88 ][-846.64 0.00 0.00 0.00 ][0.00 0.00 0.00 -867.55 ][733.21 0.00 0.00 0.00 ][0.00 0.00 0.00 1120.00 ][-1325.20 0.00 0.00 0.00 ][0.00 0.00 0.00 -1325.20 ]
-DEAL::    Quadrature point 1: [64.00 0.00 0.00 0.00 ][-110.85 0.00 0.00 0.00 ][143.11 0.00 0.00 0.00 ][-169.33 0.00 0.00 0.00 ][256.00 0.00 0.00 0.00 ][-443.41 55.43 0.00 0.00 ][572.43 -214.66 0.00 0.00 ][-677.31 507.98 0.00 0.00 ][0.00 0.00 0.00 -256.00 ][0.00 0.00 55.43 -443.41 ][0.00 0.00 214.66 -572.43 ][0.00 0.00 507.98 -677.31 ][0.00 0.00 0.00 -64.00 ][0.00 0.00 0.00 -110.85 ][0.00 0.00 0.00 -143.11 ][0.00 0.00 0.00 -169.33 ][-320.00 0.00 0.00 0.00 ][0.00 0.00 0.00 320.00 ][-277.13 0.00 0.00 0.00 ][0.00 0.00 0.00 554.26 ][-500.88 0.00 0.00 0.00 ][0.00 0.00 0.00 715.54 ][554.26 0.00 0.00 0.00 ][0.00 0.00 0.00 846.64 ][480.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -277.13 ][867.55 0.00 0.00 0.00 ][0.00 0.00 0.00 -480.00 ][-715.54 0.00 0.00 0.00 ][0.00 0.00 0.00 -619.68 ][-619.68 0.00 0.00 0.00 ][0.00 0.00 0.00 -733.21 ][-1120.00 0.00 0.00 0.00 ][0.00 0.00 0.00 500.88 ][846.64 0.00 0.00 0.00 ][0.00 0.00 0.00 867.55 ][733.21 0.00 0.00 0.00 ][0.00 0.00 0.00 1120.00 ][1325.20 0.00 0.00 0.00 ][0.00 0.00 0.00 1325.20 ]
-DEAL::    Quadrature point 2: [-256.00 0.00 0.00 0.00 ][-443.41 55.43 0.00 0.00 ][-572.43 214.66 0.00 0.00 ][-677.31 507.98 0.00 0.00 ][-64.00 0.00 0.00 0.00 ][-110.85 0.00 0.00 0.00 ][-143.11 0.00 0.00 0.00 ][-169.33 0.00 0.00 0.00 ][0.00 0.00 0.00 64.00 ][0.00 0.00 0.00 -110.85 ][0.00 0.00 0.00 143.11 ][0.00 0.00 0.00 -169.33 ][0.00 0.00 0.00 256.00 ][0.00 0.00 55.43 -443.41 ][0.00 0.00 -214.66 572.43 ][0.00 0.00 507.98 -677.31 ][320.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -320.00 ][-277.13 0.00 0.00 0.00 ][0.00 0.00 0.00 554.26 ][500.88 0.00 0.00 0.00 ][0.00 0.00 0.00 -715.54 ][554.26 0.00 0.00 0.00 ][0.00 0.00 0.00 846.64 ][-480.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -277.13 ][867.55 0.00 0.00 0.00 ][0.00 0.00 0.00 480.00 ][715.54 0.00 0.00 0.00 ][0.00 0.00 0.00 -619.68 ][-619.68 0.00 0.00 0.00 ][0.00 0.00 0.00 733.21 ][1120.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -500.88 ][846.64 0.00 0.00 0.00 ][0.00 0.00 0.00 867.55 ][-733.21 0.00 0.00 0.00 ][0.00 0.00 0.00 -1120.00 ][1325.20 0.00 0.00 0.00 ][0.00 0.00 0.00 1325.20 ]
-DEAL::    Quadrature point 3: [64.00 0.00 0.00 0.00 ][110.85 0.00 0.00 0.00 ][143.11 0.00 0.00 0.00 ][169.33 0.00 0.00 0.00 ][256.00 0.00 0.00 0.00 ][443.41 55.43 0.00 0.00 ][572.43 214.66 0.00 0.00 ][677.31 507.98 0.00 0.00 ][0.00 0.00 0.00 64.00 ][0.00 0.00 0.00 110.85 ][0.00 0.00 0.00 143.11 ][0.00 0.00 0.00 169.33 ][0.00 0.00 0.00 256.00 ][0.00 0.00 55.43 443.41 ][0.00 0.00 214.66 572.43 ][0.00 0.00 507.98 677.31 ][-320.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -320.00 ][-277.13 0.00 0.00 0.00 ][0.00 0.00 0.00 -554.26 ][-500.88 0.00 0.00 0.00 ][0.00 0.00 0.00 -715.54 ][-554.26 0.00 0.00 0.00 ][0.00 0.00 0.00 -846.64 ][-480.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -277.13 ][-867.55 0.00 0.00 0.00 ][0.00 0.00 0.00 -480.00 ][-715.54 0.00 0.00 0.00 ][0.00 0.00 0.00 -619.68 ][-619.68 0.00 0.00 0.00 ][0.00 0.00 0.00 -733.21 ][-1120.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -500.88 ][-846.64 0.00 0.00 0.00 ][0.00 0.00 0.00 -867.55 ][-733.21 0.00 0.00 0.00 ][0.00 0.00 0.00 -1120.00 ][-1325.20 0.00 0.00 0.00 ][0.00 0.00 0.00 -1325.20 ]
+DEAL::    Quadrature point 0: [4.00 0 ][-6.93 0 ][8.94 0 ][-10.58 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 4.00 ][0 -6.93 ][0 8.94 ][0 -10.58 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL::    Quadrature point 1: [4.00 0 ][0 0 ][-4.47 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0.75 ][0 -1.30 ][0 1.68 ][0 -1.98 ][0 0.75 ][0 -1.30 ][0 1.68 ][0 -1.98 ][0 0 ][0 2.50 ][0 0 ][0 -4.33 ][0 0 ][0 5.59 ][0 0 ][0 -6.61 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 -7.83 ][0 0 ][0 13.56 ][0 0 ][0 -17.50 ][0 0 ][0 20.71 ]
 DEAL::  h=0.12
-DEAL::    Quadrature point 0: [-1024.00 0.00 0.00 0.00 ][1773.62 221.70 0.00 0.00 ][-2289.73 -858.65 0.00 0.00 ][2709.25 2031.94 0.00 0.00 ][-256.00 0.00 0.00 0.00 ][443.41 0.00 0.00 0.00 ][-572.43 0.00 0.00 0.00 ][677.31 0.00 0.00 0.00 ][0.00 0.00 0.00 -1024.00 ][0.00 0.00 221.70 1773.62 ][0.00 0.00 -858.65 -2289.73 ][0.00 0.00 2031.94 2709.25 ][0.00 0.00 0.00 -256.00 ][0.00 0.00 0.00 443.41 ][0.00 0.00 0.00 -572.43 ][0.00 0.00 0.00 677.31 ][1280.00 0.00 0.00 0.00 ][0.00 0.00 0.00 1280.00 ][-1108.51 0.00 0.00 0.00 ][0.00 0.00 0.00 -2217.03 ][2003.52 0.00 0.00 0.00 ][0.00 0.00 0.00 2862.17 ][-2217.03 0.00 0.00 0.00 ][0.00 0.00 0.00 -3386.56 ][1920.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -1108.51 ][-3470.19 0.00 0.00 0.00 ][0.00 0.00 0.00 1920.00 ][2862.17 0.00 0.00 0.00 ][0.00 0.00 0.00 -2478.71 ][-2478.71 0.00 0.00 0.00 ][0.00 0.00 0.00 2932.85 ][4480.00 0.00 0.00 0.00 ][0.00 0.00 0.00 2003.52 ][-3386.56 0.00 0.00 0.00 ][0.00 0.00 0.00 -3470.19 ][2932.85 0.00 0.00 0.00 ][0.00 0.00 0.00 4480.00 ][-5300.81 0.00 0.00 0.00 ][0.00 0.00 0.00 -5300.81 ]
-DEAL::    Quadrature point 1: [256.00 0.00 0.00 0.00 ][-443.41 0.00 0.00 0.00 ][572.43 0.00 0.00 0.00 ][-677.31 0.00 0.00 0.00 ][1024.00 0.00 0.00 0.00 ][-1773.62 221.70 0.00 0.00 ][2289.73 -858.65 0.00 0.00 ][-2709.25 2031.94 0.00 0.00 ][0.00 0.00 0.00 -1024.00 ][0.00 0.00 221.70 -1773.62 ][0.00 0.00 858.65 -2289.73 ][0.00 0.00 2031.94 -2709.25 ][0.00 0.00 0.00 -256.00 ][0.00 0.00 0.00 -443.41 ][0.00 0.00 0.00 -572.43 ][0.00 0.00 0.00 -677.31 ][-1280.00 0.00 0.00 0.00 ][0.00 0.00 0.00 1280.00 ][-1108.51 0.00 0.00 0.00 ][0.00 0.00 0.00 2217.03 ][-2003.52 0.00 0.00 0.00 ][0.00 0.00 0.00 2862.17 ][2217.03 0.00 0.00 0.00 ][0.00 0.00 0.00 3386.56 ][1920.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -1108.51 ][3470.19 0.00 0.00 0.00 ][0.00 0.00 0.00 -1920.00 ][-2862.17 0.00 0.00 0.00 ][0.00 0.00 0.00 -2478.71 ][-2478.71 0.00 0.00 0.00 ][0.00 0.00 0.00 -2932.85 ][-4480.00 0.00 0.00 0.00 ][0.00 0.00 0.00 2003.52 ][3386.56 0.00 0.00 0.00 ][0.00 0.00 0.00 3470.19 ][2932.85 0.00 0.00 0.00 ][0.00 0.00 0.00 4480.00 ][5300.81 0.00 0.00 0.00 ][0.00 0.00 0.00 5300.81 ]
-DEAL::    Quadrature point 2: [-1024.00 0.00 0.00 0.00 ][-1773.62 221.70 0.00 0.00 ][-2289.73 858.65 0.00 0.00 ][-2709.25 2031.94 0.00 0.00 ][-256.00 0.00 0.00 0.00 ][-443.41 0.00 0.00 0.00 ][-572.43 0.00 0.00 0.00 ][-677.31 0.00 0.00 0.00 ][0.00 0.00 0.00 256.00 ][0.00 0.00 0.00 -443.41 ][0.00 0.00 0.00 572.43 ][0.00 0.00 0.00 -677.31 ][0.00 0.00 0.00 1024.00 ][0.00 0.00 221.70 -1773.62 ][0.00 0.00 -858.65 2289.73 ][0.00 0.00 2031.94 -2709.25 ][1280.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -1280.00 ][-1108.51 0.00 0.00 0.00 ][0.00 0.00 0.00 2217.03 ][2003.52 0.00 0.00 0.00 ][0.00 0.00 0.00 -2862.17 ][2217.03 0.00 0.00 0.00 ][0.00 0.00 0.00 3386.56 ][-1920.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -1108.51 ][3470.19 0.00 0.00 0.00 ][0.00 0.00 0.00 1920.00 ][2862.17 0.00 0.00 0.00 ][0.00 0.00 0.00 -2478.71 ][-2478.71 0.00 0.00 0.00 ][0.00 0.00 0.00 2932.85 ][4480.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -2003.52 ][3386.56 0.00 0.00 0.00 ][0.00 0.00 0.00 3470.19 ][-2932.85 0.00 0.00 0.00 ][0.00 0.00 0.00 -4480.00 ][5300.81 0.00 0.00 0.00 ][0.00 0.00 0.00 5300.81 ]
-DEAL::    Quadrature point 3: [256.00 0.00 0.00 0.00 ][443.41 0.00 0.00 0.00 ][572.43 0.00 0.00 0.00 ][677.31 0.00 0.00 0.00 ][1024.00 0.00 0.00 0.00 ][1773.62 221.70 0.00 0.00 ][2289.73 858.65 0.00 0.00 ][2709.25 2031.94 0.00 0.00 ][0.00 0.00 0.00 256.00 ][0.00 0.00 0.00 443.41 ][0.00 0.00 0.00 572.43 ][0.00 0.00 0.00 677.31 ][0.00 0.00 0.00 1024.00 ][0.00 0.00 221.70 1773.62 ][0.00 0.00 858.65 2289.73 ][0.00 0.00 2031.94 2709.25 ][-1280.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -1280.00 ][-1108.51 0.00 0.00 0.00 ][0.00 0.00 0.00 -2217.03 ][-2003.52 0.00 0.00 0.00 ][0.00 0.00 0.00 -2862.17 ][-2217.03 0.00 0.00 0.00 ][0.00 0.00 0.00 -3386.56 ][-1920.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -1108.51 ][-3470.19 0.00 0.00 0.00 ][0.00 0.00 0.00 -1920.00 ][-2862.17 0.00 0.00 0.00 ][0.00 0.00 0.00 -2478.71 ][-2478.71 0.00 0.00 0.00 ][0.00 0.00 0.00 -2932.85 ][-4480.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -2003.52 ][-3386.56 0.00 0.00 0.00 ][0.00 0.00 0.00 -3470.19 ][-2932.85 0.00 0.00 0.00 ][0.00 0.00 0.00 -4480.00 ][-5300.81 0.00 0.00 0.00 ][0.00 0.00 0.00 -5300.81 ]
+DEAL::    Quadrature point 0: [8.00 0 ][-13.86 0 ][17.89 0 ][-21.17 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 8.00 ][0 -13.86 ][0 17.89 ][0 -21.17 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL::    Quadrature point 1: [8.00 0 ][0 0 ][-8.94 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 1.50 ][0 -2.60 ][0 3.35 ][0 -3.97 ][0 1.50 ][0 -2.60 ][0 3.35 ][0 -3.97 ][0 0 ][0 5.00 ][0 0 ][0 -8.66 ][0 0 ][0 11.18 ][0 0 ][0 -13.23 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 -15.65 ][0 0 ][0 27.11 ][0 0 ][0 -35.00 ][0 0 ][0 41.41 ]
 DEAL::  h=0.06
-DEAL::    Quadrature point 0: [-4096.00 0.00 0.00 0.00 ][7094.48 886.81 0.00 0.00 ][-9158.93 -3434.60 0.00 0.00 ][10837.00 8127.75 0.00 0.00 ][-1024.00 0.00 0.00 0.00 ][1773.62 0.00 0.00 0.00 ][-2289.73 0.00 0.00 0.00 ][2709.25 0.00 0.00 0.00 ][0.00 0.00 0.00 -4096.00 ][0.00 0.00 886.81 7094.48 ][0.00 0.00 -3434.60 -9158.93 ][0.00 0.00 8127.75 10837.00 ][0.00 0.00 0.00 -1024.00 ][0.00 0.00 0.00 1773.62 ][0.00 0.00 0.00 -2289.73 ][0.00 0.00 0.00 2709.25 ][5120.00 0.00 0.00 0.00 ][0.00 0.00 0.00 5120.00 ][-4434.05 0.00 0.00 0.00 ][0.00 0.00 0.00 -8868.10 ][8014.07 0.00 0.00 0.00 ][0.00 0.00 0.00 11448.67 ][-8868.10 0.00 0.00 0.00 ][0.00 0.00 0.00 -13546.25 ][7680.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -4434.05 ][-13880.77 0.00 0.00 0.00 ][0.00 0.00 0.00 7680.00 ][11448.67 0.00 0.00 0.00 ][0.00 0.00 0.00 -9914.84 ][-9914.84 0.00 0.00 0.00 ][0.00 0.00 0.00 11731.39 ][17920.00 0.00 0.00 0.00 ][0.00 0.00 0.00 8014.07 ][-13546.25 0.00 0.00 0.00 ][0.00 0.00 0.00 -13880.77 ][11731.39 0.00 0.00 0.00 ][0.00 0.00 0.00 17920.00 ][-21203.23 0.00 0.00 0.00 ][0.00 0.00 0.00 -21203.23 ]
-DEAL::    Quadrature point 1: [1024.00 0.00 0.00 0.00 ][-1773.62 0.00 0.00 0.00 ][2289.73 0.00 0.00 0.00 ][-2709.25 0.00 0.00 0.00 ][4096.00 0.00 0.00 0.00 ][-7094.48 886.81 0.00 0.00 ][9158.93 -3434.60 0.00 0.00 ][-10837.00 8127.75 0.00 0.00 ][0.00 0.00 0.00 -4096.00 ][0.00 0.00 886.81 -7094.48 ][0.00 0.00 3434.60 -9158.93 ][0.00 0.00 8127.75 -10837.00 ][0.00 0.00 0.00 -1024.00 ][0.00 0.00 0.00 -1773.62 ][0.00 0.00 0.00 -2289.73 ][0.00 0.00 0.00 -2709.25 ][-5120.00 0.00 0.00 0.00 ][0.00 0.00 0.00 5120.00 ][-4434.05 0.00 0.00 0.00 ][0.00 0.00 0.00 8868.10 ][-8014.07 0.00 0.00 0.00 ][0.00 0.00 0.00 11448.67 ][8868.10 0.00 0.00 0.00 ][0.00 0.00 0.00 13546.25 ][7680.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -4434.05 ][13880.77 0.00 0.00 0.00 ][0.00 0.00 0.00 -7680.00 ][-11448.67 0.00 0.00 0.00 ][0.00 0.00 0.00 -9914.84 ][-9914.84 0.00 0.00 0.00 ][0.00 0.00 0.00 -11731.39 ][-17920.00 0.00 0.00 0.00 ][0.00 0.00 0.00 8014.07 ][13546.25 0.00 0.00 0.00 ][0.00 0.00 0.00 13880.77 ][11731.39 0.00 0.00 0.00 ][0.00 0.00 0.00 17920.00 ][21203.23 0.00 0.00 0.00 ][0.00 0.00 0.00 21203.23 ]
-DEAL::    Quadrature point 2: [-4096.00 0.00 0.00 0.00 ][-7094.48 886.81 0.00 0.00 ][-9158.93 3434.60 0.00 0.00 ][-10837.00 8127.75 0.00 0.00 ][-1024.00 0.00 0.00 0.00 ][-1773.62 0.00 0.00 0.00 ][-2289.73 0.00 0.00 0.00 ][-2709.25 0.00 0.00 0.00 ][0.00 0.00 0.00 1024.00 ][0.00 0.00 0.00 -1773.62 ][0.00 0.00 0.00 2289.73 ][0.00 0.00 0.00 -2709.25 ][0.00 0.00 0.00 4096.00 ][0.00 0.00 886.81 -7094.48 ][0.00 0.00 -3434.60 9158.93 ][0.00 0.00 8127.75 -10837.00 ][5120.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -5120.00 ][-4434.05 0.00 0.00 0.00 ][0.00 0.00 0.00 8868.10 ][8014.07 0.00 0.00 0.00 ][0.00 0.00 0.00 -11448.67 ][8868.10 0.00 0.00 0.00 ][0.00 0.00 0.00 13546.25 ][-7680.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -4434.05 ][13880.77 0.00 0.00 0.00 ][0.00 0.00 0.00 7680.00 ][11448.67 0.00 0.00 0.00 ][0.00 0.00 0.00 -9914.84 ][-9914.84 0.00 0.00 0.00 ][0.00 0.00 0.00 11731.39 ][17920.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -8014.07 ][13546.25 0.00 0.00 0.00 ][0.00 0.00 0.00 13880.77 ][-11731.39 0.00 0.00 0.00 ][0.00 0.00 0.00 -17920.00 ][21203.23 0.00 0.00 0.00 ][0.00 0.00 0.00 21203.23 ]
-DEAL::    Quadrature point 3: [1024.00 0.00 0.00 0.00 ][1773.62 0.00 0.00 0.00 ][2289.73 0.00 0.00 0.00 ][2709.25 0.00 0.00 0.00 ][4096.00 0.00 0.00 0.00 ][7094.48 886.81 0.00 0.00 ][9158.93 3434.60 0.00 0.00 ][10837.00 8127.75 0.00 0.00 ][0.00 0.00 0.00 1024.00 ][0.00 0.00 0.00 1773.62 ][0.00 0.00 0.00 2289.73 ][0.00 0.00 0.00 2709.25 ][0.00 0.00 0.00 4096.00 ][0.00 0.00 886.81 7094.48 ][0.00 0.00 3434.60 9158.93 ][0.00 0.00 8127.75 10837.00 ][-5120.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -5120.00 ][-4434.05 0.00 0.00 0.00 ][0.00 0.00 0.00 -8868.10 ][-8014.07 0.00 0.00 0.00 ][0.00 0.00 0.00 -11448.67 ][-8868.10 0.00 0.00 0.00 ][0.00 0.00 0.00 -13546.25 ][-7680.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -4434.05 ][-13880.77 0.00 0.00 0.00 ][0.00 0.00 0.00 -7680.00 ][-11448.67 0.00 0.00 0.00 ][0.00 0.00 0.00 -9914.84 ][-9914.84 0.00 0.00 0.00 ][0.00 0.00 0.00 -11731.39 ][-17920.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -8014.07 ][-13546.25 0.00 0.00 0.00 ][0.00 0.00 0.00 -13880.77 ][-11731.39 0.00 0.00 0.00 ][0.00 0.00 0.00 -17920.00 ][-21203.23 0.00 0.00 0.00 ][0.00 0.00 0.00 -21203.23 ]
+DEAL::    Quadrature point 0: [16.00 0 ][-27.71 0 ][35.78 0 ][-42.33 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 16.00 ][0 -27.71 ][0 35.78 ][0 -42.33 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL::    Quadrature point 1: [16.00 0 ][0 0 ][-17.89 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 3.00 ][0 -5.20 ][0 6.71 ][0 -7.94 ][0 3.00 ][0 -5.20 ][0 6.71 ][0 -7.94 ][0 0 ][0 10.00 ][0 0 ][0 -17.32 ][0 0 ][0 22.36 ][0 0 ][0 -26.46 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 -31.30 ][0 0 ][0 54.22 ][0 0 ][0 -70.00 ][0 0 ][0 82.83 ]
 DEAL::  h=0.03
-DEAL::    Quadrature point 0: [-16384.00 0.00 0.00 0.00 ][28377.92 3547.24 0.00 0.00 ][-36635.74 -13738.40 0.00 0.00 ][43347.99 32510.99 0.00 0.00 ][-4096.00 0.00 0.00 0.00 ][7094.48 0.00 0.00 0.00 ][-9158.93 0.00 0.00 0.00 ][10837.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -16384.00 ][0.00 0.00 3547.24 28377.92 ][0.00 0.00 -13738.40 -36635.74 ][0.00 0.00 32510.99 43347.99 ][0.00 0.00 0.00 -4096.00 ][0.00 0.00 0.00 7094.48 ][0.00 0.00 0.00 -9158.93 ][0.00 0.00 0.00 10837.00 ][20480.00 0.00 0.00 0.00 ][0.00 0.00 0.00 20480.00 ][-17736.20 0.00 0.00 0.00 ][0.00 0.00 0.00 -35472.40 ][32056.27 0.00 0.00 0.00 ][0.00 0.00 0.00 45794.67 ][-35472.40 0.00 0.00 0.00 ][0.00 0.00 0.00 -54184.99 ][30720.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -17736.20 ][-55523.09 0.00 0.00 0.00 ][0.00 0.00 0.00 30720.00 ][45794.67 0.00 0.00 0.00 ][0.00 0.00 0.00 -39659.35 ][-39659.35 0.00 0.00 0.00 ][0.00 0.00 0.00 46925.58 ][71680.00 0.00 0.00 0.00 ][0.00 0.00 0.00 32056.27 ][-54184.99 0.00 0.00 0.00 ][0.00 0.00 0.00 -55523.09 ][46925.58 0.00 0.00 0.00 ][0.00 0.00 0.00 71680.00 ][-84812.92 0.00 0.00 0.00 ][0.00 0.00 0.00 -84812.92 ]
-DEAL::    Quadrature point 1: [4096.00 0.00 0.00 0.00 ][-7094.48 0.00 0.00 0.00 ][9158.93 0.00 0.00 0.00 ][-10837.00 0.00 0.00 0.00 ][16384.00 0.00 0.00 0.00 ][-28377.92 3547.24 0.00 0.00 ][36635.74 -13738.40 0.00 0.00 ][-43347.99 32510.99 0.00 0.00 ][0.00 0.00 0.00 -16384.00 ][0.00 0.00 3547.24 -28377.92 ][0.00 0.00 13738.40 -36635.74 ][0.00 0.00 32510.99 -43347.99 ][0.00 0.00 0.00 -4096.00 ][0.00 0.00 0.00 -7094.48 ][0.00 0.00 0.00 -9158.93 ][0.00 0.00 0.00 -10837.00 ][-20480.00 0.00 0.00 0.00 ][0.00 0.00 0.00 20480.00 ][-17736.20 0.00 0.00 0.00 ][0.00 0.00 0.00 35472.40 ][-32056.27 0.00 0.00 0.00 ][0.00 0.00 0.00 45794.67 ][35472.40 0.00 0.00 0.00 ][0.00 0.00 0.00 54184.99 ][30720.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -17736.20 ][55523.09 0.00 0.00 0.00 ][0.00 0.00 0.00 -30720.00 ][-45794.67 0.00 0.00 0.00 ][0.00 0.00 0.00 -39659.35 ][-39659.35 0.00 0.00 0.00 ][0.00 0.00 0.00 -46925.58 ][-71680.00 0.00 0.00 0.00 ][0.00 0.00 0.00 32056.27 ][54184.99 0.00 0.00 0.00 ][0.00 0.00 0.00 55523.09 ][46925.58 0.00 0.00 0.00 ][0.00 0.00 0.00 71680.00 ][84812.92 0.00 0.00 0.00 ][0.00 0.00 0.00 84812.92 ]
-DEAL::    Quadrature point 2: [-16384.00 0.00 0.00 0.00 ][-28377.92 3547.24 0.00 0.00 ][-36635.74 13738.40 0.00 0.00 ][-43347.99 32510.99 0.00 0.00 ][-4096.00 0.00 0.00 0.00 ][-7094.48 0.00 0.00 0.00 ][-9158.93 0.00 0.00 0.00 ][-10837.00 0.00 0.00 0.00 ][0.00 0.00 0.00 4096.00 ][0.00 0.00 0.00 -7094.48 ][0.00 0.00 0.00 9158.93 ][0.00 0.00 0.00 -10837.00 ][0.00 0.00 0.00 16384.00 ][0.00 0.00 3547.24 -28377.92 ][0.00 0.00 -13738.40 36635.74 ][0.00 0.00 32510.99 -43347.99 ][20480.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -20480.00 ][-17736.20 0.00 0.00 0.00 ][0.00 0.00 0.00 35472.40 ][32056.27 0.00 0.00 0.00 ][0.00 0.00 0.00 -45794.67 ][35472.40 0.00 0.00 0.00 ][0.00 0.00 0.00 54184.99 ][-30720.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -17736.20 ][55523.09 0.00 0.00 0.00 ][0.00 0.00 0.00 30720.00 ][45794.67 0.00 0.00 0.00 ][0.00 0.00 0.00 -39659.35 ][-39659.35 0.00 0.00 0.00 ][0.00 0.00 0.00 46925.58 ][71680.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -32056.27 ][54184.99 0.00 0.00 0.00 ][0.00 0.00 0.00 55523.09 ][-46925.58 0.00 0.00 0.00 ][0.00 0.00 0.00 -71680.00 ][84812.92 0.00 0.00 0.00 ][0.00 0.00 0.00 84812.92 ]
-DEAL::    Quadrature point 3: [4096.00 0.00 0.00 0.00 ][7094.48 0.00 0.00 0.00 ][9158.93 0.00 0.00 0.00 ][10837.00 0.00 0.00 0.00 ][16384.00 0.00 0.00 0.00 ][28377.92 3547.24 0.00 0.00 ][36635.74 13738.40 0.00 0.00 ][43347.99 32510.99 0.00 0.00 ][0.00 0.00 0.00 4096.00 ][0.00 0.00 0.00 7094.48 ][0.00 0.00 0.00 9158.93 ][0.00 0.00 0.00 10837.00 ][0.00 0.00 0.00 16384.00 ][0.00 0.00 3547.24 28377.92 ][0.00 0.00 13738.40 36635.74 ][0.00 0.00 32510.99 43347.99 ][-20480.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -20480.00 ][-17736.20 0.00 0.00 0.00 ][0.00 0.00 0.00 -35472.40 ][-32056.27 0.00 0.00 0.00 ][0.00 0.00 0.00 -45794.67 ][-35472.40 0.00 0.00 0.00 ][0.00 0.00 0.00 -54184.99 ][-30720.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -17736.20 ][-55523.09 0.00 0.00 0.00 ][0.00 0.00 0.00 -30720.00 ][-45794.67 0.00 0.00 0.00 ][0.00 0.00 0.00 -39659.35 ][-39659.35 0.00 0.00 0.00 ][0.00 0.00 0.00 -46925.58 ][-71680.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -32056.27 ][-54184.99 0.00 0.00 0.00 ][0.00 0.00 0.00 -55523.09 ][-46925.58 0.00 0.00 0.00 ][0.00 0.00 0.00 -71680.00 ][-84812.92 0.00 0.00 0.00 ][0.00 0.00 0.00 -84812.92 ]
+DEAL::    Quadrature point 0: [32.00 0 ][-55.43 0 ][71.55 0 ][-84.66 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 32.00 ][0 -55.43 ][0 71.55 ][0 -84.66 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL::    Quadrature point 1: [32.00 0 ][0 0 ][-35.78 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 6.00 ][0 -10.39 ][0 13.42 ][0 -15.87 ][0 6.00 ][0 -10.39 ][0 13.42 ][0 -15.87 ][0 0 ][0 20.00 ][0 0 ][0 -34.64 ][0 0 ][0 44.72 ][0 0 ][0 -52.92 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 -62.61 ][0 0 ][0 108.44 ][0 0 ][0 -140.00 ][0 0 ][0 165.65 ]
 DEAL::  h=0.02
-DEAL::    Quadrature point 0: [-65536.00 0.00 0.00 0.00 ][113511.68 14188.96 0.00 0.00 ][-146542.95 -54953.61 0.00 0.00 ][173391.96 130043.97 0.00 0.00 ][-16384.00 0.00 0.00 0.00 ][28377.92 0.00 0.00 0.00 ][-36635.74 0.00 0.00 0.00 ][43347.99 0.00 0.00 0.00 ][0.00 0.00 0.00 -65536.00 ][0.00 0.00 14188.96 113511.68 ][0.00 0.00 -54953.61 -146542.95 ][0.00 0.00 130043.97 173391.96 ][0.00 0.00 0.00 -16384.00 ][0.00 0.00 0.00 28377.92 ][0.00 0.00 0.00 -36635.74 ][0.00 0.00 0.00 43347.99 ][81920.00 0.00 0.00 0.00 ][0.00 0.00 0.00 81920.00 ][-70944.80 0.00 0.00 0.00 ][0.00 0.00 0.00 -141889.60 ][128225.08 0.00 0.00 0.00 ][0.00 0.00 0.00 183178.69 ][-141889.60 0.00 0.00 0.00 ][0.00 0.00 0.00 -216739.95 ][122880.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -70944.80 ][-222092.36 0.00 0.00 0.00 ][0.00 0.00 0.00 122880.00 ][183178.69 0.00 0.00 0.00 ][0.00 0.00 0.00 -158637.40 ][-158637.40 0.00 0.00 0.00 ][0.00 0.00 0.00 187702.30 ][286720.00 0.00 0.00 0.00 ][0.00 0.00 0.00 128225.08 ][-216739.95 0.00 0.00 0.00 ][0.00 0.00 0.00 -222092.36 ][187702.30 0.00 0.00 0.00 ][0.00 0.00 0.00 286720.00 ][-339251.68 0.00 0.00 0.00 ][0.00 0.00 0.00 -339251.68 ]
-DEAL::    Quadrature point 1: [16384.00 0.00 0.00 0.00 ][-28377.92 0.00 0.00 0.00 ][36635.74 0.00 0.00 0.00 ][-43347.99 0.00 0.00 0.00 ][65536.00 0.00 0.00 0.00 ][-113511.68 14188.96 0.00 0.00 ][146542.95 -54953.61 0.00 0.00 ][-173391.96 130043.97 0.00 0.00 ][0.00 0.00 0.00 -65536.00 ][0.00 0.00 14188.96 -113511.68 ][0.00 0.00 54953.61 -146542.95 ][0.00 0.00 130043.97 -173391.96 ][0.00 0.00 0.00 -16384.00 ][0.00 0.00 0.00 -28377.92 ][0.00 0.00 0.00 -36635.74 ][0.00 0.00 0.00 -43347.99 ][-81920.00 0.00 0.00 0.00 ][0.00 0.00 0.00 81920.00 ][-70944.80 0.00 0.00 0.00 ][0.00 0.00 0.00 141889.60 ][-128225.08 0.00 0.00 0.00 ][0.00 0.00 0.00 183178.69 ][141889.60 0.00 0.00 0.00 ][0.00 0.00 0.00 216739.95 ][122880.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -70944.80 ][222092.36 0.00 0.00 0.00 ][0.00 0.00 0.00 -122880.00 ][-183178.69 0.00 0.00 0.00 ][0.00 0.00 0.00 -158637.40 ][-158637.40 0.00 0.00 0.00 ][0.00 0.00 0.00 -187702.30 ][-286720.00 0.00 0.00 0.00 ][0.00 0.00 0.00 128225.08 ][216739.95 0.00 0.00 0.00 ][0.00 0.00 0.00 222092.36 ][187702.30 0.00 0.00 0.00 ][0.00 0.00 0.00 286720.00 ][339251.68 0.00 0.00 0.00 ][0.00 0.00 0.00 339251.68 ]
-DEAL::    Quadrature point 2: [-65536.00 0.00 0.00 0.00 ][-113511.68 14188.96 0.00 0.00 ][-146542.95 54953.61 0.00 0.00 ][-173391.96 130043.97 0.00 0.00 ][-16384.00 0.00 0.00 0.00 ][-28377.92 0.00 0.00 0.00 ][-36635.74 0.00 0.00 0.00 ][-43347.99 0.00 0.00 0.00 ][0.00 0.00 0.00 16384.00 ][0.00 0.00 0.00 -28377.92 ][0.00 0.00 0.00 36635.74 ][0.00 0.00 0.00 -43347.99 ][0.00 0.00 0.00 65536.00 ][0.00 0.00 14188.96 -113511.68 ][0.00 0.00 -54953.61 146542.95 ][0.00 0.00 130043.97 -173391.96 ][81920.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -81920.00 ][-70944.80 0.00 0.00 0.00 ][0.00 0.00 0.00 141889.60 ][128225.08 0.00 0.00 0.00 ][0.00 0.00 0.00 -183178.69 ][141889.60 0.00 0.00 0.00 ][0.00 0.00 0.00 216739.95 ][-122880.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -70944.80 ][222092.36 0.00 0.00 0.00 ][0.00 0.00 0.00 122880.00 ][183178.69 0.00 0.00 0.00 ][0.00 0.00 0.00 -158637.40 ][-158637.40 0.00 0.00 0.00 ][0.00 0.00 0.00 187702.30 ][286720.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -128225.08 ][216739.95 0.00 0.00 0.00 ][0.00 0.00 0.00 222092.36 ][-187702.30 0.00 0.00 0.00 ][0.00 0.00 0.00 -286720.00 ][339251.68 0.00 0.00 0.00 ][0.00 0.00 0.00 339251.68 ]
-DEAL::    Quadrature point 3: [16384.00 0.00 0.00 0.00 ][28377.92 0.00 0.00 0.00 ][36635.74 0.00 0.00 0.00 ][43347.99 0.00 0.00 0.00 ][65536.00 0.00 0.00 0.00 ][113511.68 14188.96 0.00 0.00 ][146542.95 54953.61 0.00 0.00 ][173391.96 130043.97 0.00 0.00 ][0.00 0.00 0.00 16384.00 ][0.00 0.00 0.00 28377.92 ][0.00 0.00 0.00 36635.74 ][0.00 0.00 0.00 43347.99 ][0.00 0.00 0.00 65536.00 ][0.00 0.00 14188.96 113511.68 ][0.00 0.00 54953.61 146542.95 ][0.00 0.00 130043.97 173391.96 ][-81920.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -81920.00 ][-70944.80 0.00 0.00 0.00 ][0.00 0.00 0.00 -141889.60 ][-128225.08 0.00 0.00 0.00 ][0.00 0.00 0.00 -183178.69 ][-141889.60 0.00 0.00 0.00 ][0.00 0.00 0.00 -216739.95 ][-122880.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -70944.80 ][-222092.36 0.00 0.00 0.00 ][0.00 0.00 0.00 -122880.00 ][-183178.69 0.00 0.00 0.00 ][0.00 0.00 0.00 -158637.40 ][-158637.40 0.00 0.00 0.00 ][0.00 0.00 0.00 -187702.30 ][-286720.00 0.00 0.00 0.00 ][0.00 0.00 0.00 -128225.08 ][-216739.95 0.00 0.00 0.00 ][0.00 0.00 0.00 -222092.36 ][-187702.30 0.00 0.00 0.00 ][0.00 0.00 0.00 -286720.00 ][-339251.68 0.00 0.00 0.00 ][0.00 0.00 0.00 -339251.68 ]
+DEAL::    Quadrature point 0: [64.00 0 ][-110.85 0 ][143.11 0 ][-169.33 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 64.00 ][0 -110.85 ][0 143.11 ][0 -169.33 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ]
+DEAL::    Quadrature point 1: [64.00 0 ][0 0 ][-71.55 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 12.00 ][0 -20.78 ][0 26.83 ][0 -31.75 ][0 12.00 ][0 -20.78 ][0 26.83 ][0 -31.75 ][0 0 ][0 40.00 ][0 0 ][0 -69.28 ][0 0 ][0 89.44 ][0 0 ][0 -105.83 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 -125.22 ][0 0 ][0 216.89 ][0 0 ][0 -280.00 ][0 0 ][0 331.30 ]
index 7d88cb9a186d34084ca79cc98808b27af1519a05..8a9c23358023180291e0a63c1e8cf7c45469df41 100644 (file)
@@ -11,9 +11,9 @@
 //
 //----------------------------  rt_15.cc  ---------------------------
 
-// Like rt_15, but use FESubfaceValues
+// Like rt_13, but use FESubfaceValues
 //
-// the test presently fails because of the issue with computing the
+// the test used to fail because of the issue with computing the
 // normals using FEFaceValue, where FEFaceValue by accident uses the
 // *face* mapping, not the *cell* mapping to compute the Piola
 // transform (leading to a missing power of h in the determinant)

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

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.