]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Check same with FEFaceValues. Results are presently wrong (no scaling, when it should...
authorkayser-herold <kayser-herold@0785d39b-7218-0410-832d-ea1e28bc413d>
Sat, 29 Jul 2006 21:59:02 +0000 (21:59 +0000)
committerkayser-herold <kayser-herold@0785d39b-7218-0410-832d-ea1e28bc413d>
Sat, 29 Jul 2006 21:59:02 +0000 (21:59 +0000)
git-svn-id: https://svn.dealii.org/trunk@13507 0785d39b-7218-0410-832d-ea1e28bc413d

tests/fe/rt_11.cc [new file with mode: 0644]
tests/fe/rt_11/cmp/generic [new file with mode: 0644]

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

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.