]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Convince myself that the FE_DGQArbitraryNodes does what it suggests.
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 17 Oct 2011 05:35:26 +0000 (05:35 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 17 Oct 2011 05:35:26 +0000 (05:35 +0000)
git-svn-id: https://svn.dealii.org/trunk@24611 0785d39b-7218-0410-832d-ea1e28bc413d

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

diff --git a/tests/fe/mass_matrix_arbitrary_nodes.cc b/tests/fe/mass_matrix_arbitrary_nodes.cc
new file mode 100644 (file)
index 0000000..ee50b48
--- /dev/null
@@ -0,0 +1,89 @@
+//----------------------------  mass_matrix_arbitrary_nodes.cc  ---------------------------
+//    mass_matrix_arbitrary_nodes.cc,v 1.2 2004/01/23 16:34:24 wolf Exp
+//    Version:
+//
+//    Copyright (C) 2011 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  mass_matrix_arbitrary_nodes.cc  ---------------------------
+
+// verify that we get a diagonal matrix when using FE_DGQArbitraryNodes with
+// the same quadrature formula for integration as for construction of the FE.
+
+#include "../tests.h"
+#include <iostream>
+#include <fstream>
+
+#include <deal.II/base/logstream.h>
+#include <deal.II/base/quadrature_lib.h>
+#include <deal.II/grid/tria_iterator.h>
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/dofs/dof_tools.h>
+#include <deal.II/lac/sparse_matrix.h>
+#include <deal.II/fe/fe_dgq.h>
+#include <deal.II/numerics/matrices.h>
+
+
+
+
+template <int dim>
+void check ()
+{
+  deallog << dim << 'D' << std::endl;
+  for (unsigned int q_points=1; q_points<7-dim; ++q_points)
+    {
+      deallog << "q_points=" << q_points << std::endl;
+
+      Triangulation<dim> tr;
+      if (dim == 1)
+       GridGenerator::hyper_cube (tr);
+      else
+       GridGenerator::hyper_ball (tr);
+      tr.refine_global (3-dim);
+
+      FE_DGQArbitraryNodes<dim> fe ((QGauss<1>(q_points)));
+      DoFHandler<dim> dh (tr);
+      dh.distribute_dofs (fe);
+
+      SparsityPattern sp (dh.n_dofs(),
+                         dh.n_dofs(),
+                         fe.dofs_per_cell);
+      DoFTools::make_sparsity_pattern (dh, sp);
+      sp.compress();
+
+      SparseMatrix<double> mass_matrix (sp);
+      MatrixTools::create_mass_matrix (dh,
+                                      QGauss<dim>(q_points),
+                                      mass_matrix);
+
+                                      // verify that the matrix is diagonal
+      for (unsigned int i=0; i<dh.n_dofs(); ++i)
+       {
+         deallog << mass_matrix.el(i,i) << std::endl;
+         for (unsigned int j=0; j<dh.n_dofs(); ++j)
+           if (i != j)
+             Assert (std::fabs(mass_matrix.el(i,j) / mass_matrix.el(i,i)) < 1e-14,
+                     ExcInternalError());
+       }
+    }
+}
+
+
+
+
+int main ()
+{
+  std::ofstream logfile("mass_matrix_arbitrary_nodes/output");
+  deallog.attach(logfile);
+  deallog.depth_console(0);
+  deallog.threshold_double(1.e-10);
+
+  check<1> ();
+  check<2> ();
+  check<3> ();
+}
+
diff --git a/tests/fe/mass_matrix_arbitrary_nodes/cmp/generic b/tests/fe/mass_matrix_arbitrary_nodes/cmp/generic
new file mode 100644 (file)
index 0000000..f1b2d9f
--- /dev/null
@@ -0,0 +1,928 @@
+
+DEAL::1D
+DEAL::q_points=1
+DEAL::0.250000
+DEAL::0.250000
+DEAL::0.250000
+DEAL::0.250000
+DEAL::q_points=2
+DEAL::0.125000
+DEAL::0.125000
+DEAL::0.125000
+DEAL::0.125000
+DEAL::0.125000
+DEAL::0.125000
+DEAL::0.125000
+DEAL::0.125000
+DEAL::q_points=3
+DEAL::0.0694444
+DEAL::0.111111
+DEAL::0.0694444
+DEAL::0.0694444
+DEAL::0.111111
+DEAL::0.0694444
+DEAL::0.0694444
+DEAL::0.111111
+DEAL::0.0694444
+DEAL::0.0694444
+DEAL::0.111111
+DEAL::0.0694444
+DEAL::q_points=4
+DEAL::0.0434819
+DEAL::0.0815181
+DEAL::0.0815181
+DEAL::0.0434819
+DEAL::0.0434819
+DEAL::0.0815181
+DEAL::0.0815181
+DEAL::0.0434819
+DEAL::0.0434819
+DEAL::0.0815181
+DEAL::0.0815181
+DEAL::0.0434819
+DEAL::0.0434819
+DEAL::0.0815181
+DEAL::0.0815181
+DEAL::0.0434819
+DEAL::q_points=5
+DEAL::0.0296159
+DEAL::0.0598286
+DEAL::0.0711111
+DEAL::0.0598286
+DEAL::0.0296159
+DEAL::0.0296159
+DEAL::0.0598286
+DEAL::0.0711111
+DEAL::0.0598286
+DEAL::0.0296159
+DEAL::0.0296159
+DEAL::0.0598286
+DEAL::0.0711111
+DEAL::0.0598286
+DEAL::0.0296159
+DEAL::0.0296159
+DEAL::0.0598286
+DEAL::0.0711111
+DEAL::0.0598286
+DEAL::0.0296159
+DEAL::2D
+DEAL::q_points=1
+DEAL::0.125000
+DEAL::0.125000
+DEAL::0.0821068
+DEAL::0.0821068
+DEAL::0.125000
+DEAL::0.0821068
+DEAL::0.125000
+DEAL::0.0821068
+DEAL::0.0857864
+DEAL::0.0857864
+DEAL::0.0857864
+DEAL::0.0857864
+DEAL::0.125000
+DEAL::0.125000
+DEAL::0.0821068
+DEAL::0.0821068
+DEAL::0.125000
+DEAL::0.0821068
+DEAL::0.125000
+DEAL::0.0821068
+DEAL::q_points=2
+DEAL::0.0343456
+DEAL::0.0343456
+DEAL::0.0281544
+DEAL::0.0281544
+DEAL::0.0343456
+DEAL::0.0343456
+DEAL::0.0281544
+DEAL::0.0281544
+DEAL::0.0236222
+DEAL::0.0236222
+DEAL::0.0174311
+DEAL::0.0174311
+DEAL::0.0236222
+DEAL::0.0236222
+DEAL::0.0174311
+DEAL::0.0174311
+DEAL::0.0343456
+DEAL::0.0281544
+DEAL::0.0343456
+DEAL::0.0281544
+DEAL::0.0236222
+DEAL::0.0174311
+DEAL::0.0236222
+DEAL::0.0174311
+DEAL::0.0343456
+DEAL::0.0281544
+DEAL::0.0343456
+DEAL::0.0281544
+DEAL::0.0236222
+DEAL::0.0174311
+DEAL::0.0236222
+DEAL::0.0174311
+DEAL::0.0214466
+DEAL::0.0214466
+DEAL::0.0214466
+DEAL::0.0214466
+DEAL::0.0214466
+DEAL::0.0214466
+DEAL::0.0214466
+DEAL::0.0214466
+DEAL::0.0214466
+DEAL::0.0214466
+DEAL::0.0214466
+DEAL::0.0214466
+DEAL::0.0214466
+DEAL::0.0214466
+DEAL::0.0214466
+DEAL::0.0214466
+DEAL::0.0343456
+DEAL::0.0343456
+DEAL::0.0281544
+DEAL::0.0281544
+DEAL::0.0343456
+DEAL::0.0343456
+DEAL::0.0281544
+DEAL::0.0281544
+DEAL::0.0236222
+DEAL::0.0236222
+DEAL::0.0174311
+DEAL::0.0174311
+DEAL::0.0236222
+DEAL::0.0236222
+DEAL::0.0174311
+DEAL::0.0174311
+DEAL::0.0343456
+DEAL::0.0281544
+DEAL::0.0343456
+DEAL::0.0281544
+DEAL::0.0236222
+DEAL::0.0174311
+DEAL::0.0236222
+DEAL::0.0174311
+DEAL::0.0343456
+DEAL::0.0281544
+DEAL::0.0343456
+DEAL::0.0281544
+DEAL::0.0236222
+DEAL::0.0174311
+DEAL::0.0236222
+DEAL::0.0174311
+DEAL::q_points=3
+DEAL::0.0109269
+DEAL::0.0174830
+DEAL::0.0109269
+DEAL::0.0154321
+DEAL::0.0246914
+DEAL::0.0154321
+DEAL::0.00836324
+DEAL::0.0133812
+DEAL::0.00836324
+DEAL::0.0109269
+DEAL::0.0174830
+DEAL::0.0109269
+DEAL::0.0154321
+DEAL::0.0246914
+DEAL::0.0154321
+DEAL::0.00836324
+DEAL::0.0133812
+DEAL::0.00836324
+DEAL::0.00761723
+DEAL::0.0121876
+DEAL::0.00761723
+DEAL::0.0101366
+DEAL::0.0162186
+DEAL::0.0101366
+DEAL::0.00505357
+DEAL::0.00808572
+DEAL::0.00505357
+DEAL::0.00761723
+DEAL::0.0121876
+DEAL::0.00761723
+DEAL::0.0101366
+DEAL::0.0162186
+DEAL::0.0101366
+DEAL::0.00505357
+DEAL::0.00808572
+DEAL::0.00505357
+DEAL::0.0109269
+DEAL::0.0154321
+DEAL::0.00836324
+DEAL::0.0174830
+DEAL::0.0246914
+DEAL::0.0133812
+DEAL::0.0109269
+DEAL::0.0154321
+DEAL::0.00836324
+DEAL::0.00761723
+DEAL::0.0101366
+DEAL::0.00505357
+DEAL::0.0121876
+DEAL::0.0162186
+DEAL::0.00808572
+DEAL::0.00761723
+DEAL::0.0101366
+DEAL::0.00505357
+DEAL::0.0109269
+DEAL::0.0154321
+DEAL::0.00836324
+DEAL::0.0174830
+DEAL::0.0246914
+DEAL::0.0133812
+DEAL::0.0109269
+DEAL::0.0154321
+DEAL::0.00836324
+DEAL::0.00761723
+DEAL::0.0101366
+DEAL::0.00505357
+DEAL::0.0121876
+DEAL::0.0162186
+DEAL::0.00808572
+DEAL::0.00761723
+DEAL::0.0101366
+DEAL::0.00505357
+DEAL::0.00661932
+DEAL::0.0105909
+DEAL::0.00661932
+DEAL::0.0105909
+DEAL::0.0169455
+DEAL::0.0105909
+DEAL::0.00661932
+DEAL::0.0105909
+DEAL::0.00661932
+DEAL::0.00661932
+DEAL::0.0105909
+DEAL::0.00661932
+DEAL::0.0105909
+DEAL::0.0169455
+DEAL::0.0105909
+DEAL::0.00661932
+DEAL::0.0105909
+DEAL::0.00661932
+DEAL::0.00661932
+DEAL::0.0105909
+DEAL::0.00661932
+DEAL::0.0105909
+DEAL::0.0169455
+DEAL::0.0105909
+DEAL::0.00661932
+DEAL::0.0105909
+DEAL::0.00661932
+DEAL::0.00661932
+DEAL::0.0105909
+DEAL::0.00661932
+DEAL::0.0105909
+DEAL::0.0169455
+DEAL::0.0105909
+DEAL::0.00661932
+DEAL::0.0105909
+DEAL::0.00661932
+DEAL::0.0109269
+DEAL::0.0174830
+DEAL::0.0109269
+DEAL::0.0154321
+DEAL::0.0246914
+DEAL::0.0154321
+DEAL::0.00836324
+DEAL::0.0133812
+DEAL::0.00836324
+DEAL::0.0109269
+DEAL::0.0174830
+DEAL::0.0109269
+DEAL::0.0154321
+DEAL::0.0246914
+DEAL::0.0154321
+DEAL::0.00836324
+DEAL::0.0133812
+DEAL::0.00836324
+DEAL::0.00761723
+DEAL::0.0121876
+DEAL::0.00761723
+DEAL::0.0101366
+DEAL::0.0162186
+DEAL::0.0101366
+DEAL::0.00505357
+DEAL::0.00808572
+DEAL::0.00505357
+DEAL::0.00761723
+DEAL::0.0121876
+DEAL::0.00761723
+DEAL::0.0101366
+DEAL::0.0162186
+DEAL::0.0101366
+DEAL::0.00505357
+DEAL::0.00808572
+DEAL::0.00505357
+DEAL::0.0109269
+DEAL::0.0154321
+DEAL::0.00836324
+DEAL::0.0174830
+DEAL::0.0246914
+DEAL::0.0133812
+DEAL::0.0109269
+DEAL::0.0154321
+DEAL::0.00836324
+DEAL::0.00761723
+DEAL::0.0101366
+DEAL::0.00505357
+DEAL::0.0121876
+DEAL::0.0162186
+DEAL::0.00808572
+DEAL::0.00761723
+DEAL::0.0101366
+DEAL::0.00505357
+DEAL::0.0109269
+DEAL::0.0154321
+DEAL::0.00836324
+DEAL::0.0174830
+DEAL::0.0246914
+DEAL::0.0133812
+DEAL::0.0109269
+DEAL::0.0154321
+DEAL::0.00836324
+DEAL::0.00761723
+DEAL::0.0101366
+DEAL::0.00505357
+DEAL::0.0121876
+DEAL::0.0162186
+DEAL::0.00808572
+DEAL::0.00761723
+DEAL::0.0101366
+DEAL::0.00505357
+DEAL::q_points=4
+DEAL::0.00434003
+DEAL::0.00813652
+DEAL::0.00813652
+DEAL::0.00434003
+DEAL::0.00750264
+DEAL::0.0140657
+DEAL::0.0140657
+DEAL::0.00750264
+DEAL::0.00667560
+DEAL::0.0125152
+DEAL::0.0125152
+DEAL::0.00667560
+DEAL::0.00322266
+DEAL::0.00604172
+DEAL::0.00604172
+DEAL::0.00322266
+DEAL::0.00434003
+DEAL::0.00813652
+DEAL::0.00813652
+DEAL::0.00434003
+DEAL::0.00750264
+DEAL::0.0140657
+DEAL::0.0140657
+DEAL::0.00750264
+DEAL::0.00667560
+DEAL::0.0125152
+DEAL::0.0125152
+DEAL::0.00667560
+DEAL::0.00322266
+DEAL::0.00604172
+DEAL::0.00604172
+DEAL::0.00322266
+DEAL::0.00304248
+DEAL::0.00570392
+DEAL::0.00570392
+DEAL::0.00304248
+DEAL::0.00507004
+DEAL::0.00950512
+DEAL::0.00950512
+DEAL::0.00507004
+DEAL::0.00424300
+DEAL::0.00795462
+DEAL::0.00795462
+DEAL::0.00424300
+DEAL::0.00192511
+DEAL::0.00360912
+DEAL::0.00360912
+DEAL::0.00192511
+DEAL::0.00304248
+DEAL::0.00570392
+DEAL::0.00570392
+DEAL::0.00304248
+DEAL::0.00507004
+DEAL::0.00950512
+DEAL::0.00950512
+DEAL::0.00507004
+DEAL::0.00424300
+DEAL::0.00795462
+DEAL::0.00795462
+DEAL::0.00424300
+DEAL::0.00192511
+DEAL::0.00360912
+DEAL::0.00360912
+DEAL::0.00192511
+DEAL::0.00434003
+DEAL::0.00750264
+DEAL::0.00667560
+DEAL::0.00322266
+DEAL::0.00813652
+DEAL::0.0140657
+DEAL::0.0125152
+DEAL::0.00604172
+DEAL::0.00813652
+DEAL::0.0140657
+DEAL::0.0125152
+DEAL::0.00604172
+DEAL::0.00434003
+DEAL::0.00750264
+DEAL::0.00667560
+DEAL::0.00322266
+DEAL::0.00304248
+DEAL::0.00507004
+DEAL::0.00424300
+DEAL::0.00192511
+DEAL::0.00570392
+DEAL::0.00950512
+DEAL::0.00795462
+DEAL::0.00360912
+DEAL::0.00570392
+DEAL::0.00950512
+DEAL::0.00795462
+DEAL::0.00360912
+DEAL::0.00304248
+DEAL::0.00507004
+DEAL::0.00424300
+DEAL::0.00192511
+DEAL::0.00434003
+DEAL::0.00750264
+DEAL::0.00667560
+DEAL::0.00322266
+DEAL::0.00813652
+DEAL::0.0140657
+DEAL::0.0125152
+DEAL::0.00604172
+DEAL::0.00813652
+DEAL::0.0140657
+DEAL::0.0125152
+DEAL::0.00604172
+DEAL::0.00434003
+DEAL::0.00750264
+DEAL::0.00667560
+DEAL::0.00322266
+DEAL::0.00304248
+DEAL::0.00507004
+DEAL::0.00424300
+DEAL::0.00192511
+DEAL::0.00570392
+DEAL::0.00950512
+DEAL::0.00795462
+DEAL::0.00360912
+DEAL::0.00570392
+DEAL::0.00950512
+DEAL::0.00795462
+DEAL::0.00360912
+DEAL::0.00304248
+DEAL::0.00507004
+DEAL::0.00424300
+DEAL::0.00192511
+DEAL::0.00259510
+DEAL::0.00486520
+DEAL::0.00486520
+DEAL::0.00259510
+DEAL::0.00486520
+DEAL::0.00912110
+DEAL::0.00912110
+DEAL::0.00486520
+DEAL::0.00486520
+DEAL::0.00912110
+DEAL::0.00912110
+DEAL::0.00486520
+DEAL::0.00259510
+DEAL::0.00486520
+DEAL::0.00486520
+DEAL::0.00259510
+DEAL::0.00259510
+DEAL::0.00486520
+DEAL::0.00486520
+DEAL::0.00259510
+DEAL::0.00486520
+DEAL::0.00912110
+DEAL::0.00912110
+DEAL::0.00486520
+DEAL::0.00486520
+DEAL::0.00912110
+DEAL::0.00912110
+DEAL::0.00486520
+DEAL::0.00259510
+DEAL::0.00486520
+DEAL::0.00486520
+DEAL::0.00259510
+DEAL::0.00259510
+DEAL::0.00486520
+DEAL::0.00486520
+DEAL::0.00259510
+DEAL::0.00486520
+DEAL::0.00912110
+DEAL::0.00912110
+DEAL::0.00486520
+DEAL::0.00486520
+DEAL::0.00912110
+DEAL::0.00912110
+DEAL::0.00486520
+DEAL::0.00259510
+DEAL::0.00486520
+DEAL::0.00486520
+DEAL::0.00259510
+DEAL::0.00259510
+DEAL::0.00486520
+DEAL::0.00486520
+DEAL::0.00259510
+DEAL::0.00486520
+DEAL::0.00912110
+DEAL::0.00912110
+DEAL::0.00486520
+DEAL::0.00486520
+DEAL::0.00912110
+DEAL::0.00912110
+DEAL::0.00486520
+DEAL::0.00259510
+DEAL::0.00486520
+DEAL::0.00486520
+DEAL::0.00259510
+DEAL::0.00434003
+DEAL::0.00813652
+DEAL::0.00813652
+DEAL::0.00434003
+DEAL::0.00750264
+DEAL::0.0140657
+DEAL::0.0140657
+DEAL::0.00750264
+DEAL::0.00667560
+DEAL::0.0125152
+DEAL::0.0125152
+DEAL::0.00667560
+DEAL::0.00322266
+DEAL::0.00604172
+DEAL::0.00604172
+DEAL::0.00322266
+DEAL::0.00434003
+DEAL::0.00813652
+DEAL::0.00813652
+DEAL::0.00434003
+DEAL::0.00750264
+DEAL::0.0140657
+DEAL::0.0140657
+DEAL::0.00750264
+DEAL::0.00667560
+DEAL::0.0125152
+DEAL::0.0125152
+DEAL::0.00667560
+DEAL::0.00322266
+DEAL::0.00604172
+DEAL::0.00604172
+DEAL::0.00322266
+DEAL::0.00304248
+DEAL::0.00570392
+DEAL::0.00570392
+DEAL::0.00304248
+DEAL::0.00507004
+DEAL::0.00950512
+DEAL::0.00950512
+DEAL::0.00507004
+DEAL::0.00424300
+DEAL::0.00795462
+DEAL::0.00795462
+DEAL::0.00424300
+DEAL::0.00192511
+DEAL::0.00360912
+DEAL::0.00360912
+DEAL::0.00192511
+DEAL::0.00304248
+DEAL::0.00570392
+DEAL::0.00570392
+DEAL::0.00304248
+DEAL::0.00507004
+DEAL::0.00950512
+DEAL::0.00950512
+DEAL::0.00507004
+DEAL::0.00424300
+DEAL::0.00795462
+DEAL::0.00795462
+DEAL::0.00424300
+DEAL::0.00192511
+DEAL::0.00360912
+DEAL::0.00360912
+DEAL::0.00192511
+DEAL::0.00434003
+DEAL::0.00750264
+DEAL::0.00667560
+DEAL::0.00322266
+DEAL::0.00813652
+DEAL::0.0140657
+DEAL::0.0125152
+DEAL::0.00604172
+DEAL::0.00813652
+DEAL::0.0140657
+DEAL::0.0125152
+DEAL::0.00604172
+DEAL::0.00434003
+DEAL::0.00750264
+DEAL::0.00667560
+DEAL::0.00322266
+DEAL::0.00304248
+DEAL::0.00507004
+DEAL::0.00424300
+DEAL::0.00192511
+DEAL::0.00570392
+DEAL::0.00950512
+DEAL::0.00795462
+DEAL::0.00360912
+DEAL::0.00570392
+DEAL::0.00950512
+DEAL::0.00795462
+DEAL::0.00360912
+DEAL::0.00304248
+DEAL::0.00507004
+DEAL::0.00424300
+DEAL::0.00192511
+DEAL::0.00434003
+DEAL::0.00750264
+DEAL::0.00667560
+DEAL::0.00322266
+DEAL::0.00813652
+DEAL::0.0140657
+DEAL::0.0125152
+DEAL::0.00604172
+DEAL::0.00813652
+DEAL::0.0140657
+DEAL::0.0125152
+DEAL::0.00604172
+DEAL::0.00434003
+DEAL::0.00750264
+DEAL::0.00667560
+DEAL::0.00322266
+DEAL::0.00304248
+DEAL::0.00507004
+DEAL::0.00424300
+DEAL::0.00192511
+DEAL::0.00570392
+DEAL::0.00950512
+DEAL::0.00795462
+DEAL::0.00360912
+DEAL::0.00570392
+DEAL::0.00950512
+DEAL::0.00795462
+DEAL::0.00360912
+DEAL::0.00304248
+DEAL::0.00507004
+DEAL::0.00424300
+DEAL::0.00192511
+DEAL::3D
+DEAL::q_points=1
+DEAL::0.0754991
+DEAL::0.227671
+DEAL::0.227671
+DEAL::0.227671
+DEAL::0.227671
+DEAL::0.227671
+DEAL::0.227671
+DEAL::q_points=2
+DEAL::0.00943739
+DEAL::0.00943739
+DEAL::0.00943739
+DEAL::0.00943739
+DEAL::0.00943739
+DEAL::0.00943739
+DEAL::0.00943739
+DEAL::0.00943739
+DEAL::0.0457532
+DEAL::0.0457532
+DEAL::0.0457532
+DEAL::0.0457532
+DEAL::0.0152511
+DEAL::0.0152511
+DEAL::0.0152511
+DEAL::0.0152511
+DEAL::0.0457532
+DEAL::0.0457532
+DEAL::0.0152511
+DEAL::0.0152511
+DEAL::0.0457532
+DEAL::0.0457532
+DEAL::0.0152511
+DEAL::0.0152511
+DEAL::0.0457532
+DEAL::0.0457532
+DEAL::0.0152511
+DEAL::0.0152511
+DEAL::0.0457532
+DEAL::0.0457532
+DEAL::0.0152511
+DEAL::0.0152511
+DEAL::0.0457532
+DEAL::0.0152511
+DEAL::0.0457532
+DEAL::0.0152511
+DEAL::0.0457532
+DEAL::0.0152511
+DEAL::0.0457532
+DEAL::0.0152511
+DEAL::0.0457532
+DEAL::0.0457532
+DEAL::0.0152511
+DEAL::0.0152511
+DEAL::0.0457532
+DEAL::0.0457532
+DEAL::0.0152511
+DEAL::0.0152511
+DEAL::0.0457532
+DEAL::0.0152511
+DEAL::0.0457532
+DEAL::0.0152511
+DEAL::0.0457532
+DEAL::0.0152511
+DEAL::0.0457532
+DEAL::0.0152511
+DEAL::q_points=3
+DEAL::0.00161821
+DEAL::0.00258913
+DEAL::0.00161821
+DEAL::0.00258913
+DEAL::0.00414261
+DEAL::0.00258913
+DEAL::0.00161821
+DEAL::0.00258913
+DEAL::0.00161821
+DEAL::0.00258913
+DEAL::0.00414261
+DEAL::0.00258913
+DEAL::0.00414261
+DEAL::0.00662818
+DEAL::0.00414261
+DEAL::0.00258913
+DEAL::0.00414261
+DEAL::0.00258913
+DEAL::0.00161821
+DEAL::0.00258913
+DEAL::0.00161821
+DEAL::0.00258913
+DEAL::0.00414261
+DEAL::0.00258913
+DEAL::0.00161821
+DEAL::0.00258913
+DEAL::0.00161821
+DEAL::0.00901889
+DEAL::0.0144302
+DEAL::0.00901889
+DEAL::0.0144302
+DEAL::0.0230884
+DEAL::0.0144302
+DEAL::0.00901889
+DEAL::0.0144302
+DEAL::0.00901889
+DEAL::0.00780764
+DEAL::0.0124922
+DEAL::0.00780764
+DEAL::0.0124922
+DEAL::0.0199876
+DEAL::0.0124922
+DEAL::0.00780764
+DEAL::0.0124922
+DEAL::0.00780764
+DEAL::0.00200193
+DEAL::0.00320309
+DEAL::0.00200193
+DEAL::0.00320309
+DEAL::0.00512495
+DEAL::0.00320309
+DEAL::0.00200193
+DEAL::0.00320309
+DEAL::0.00200193
+DEAL::0.00901889
+DEAL::0.0144302
+DEAL::0.00901889
+DEAL::0.00780764
+DEAL::0.0124922
+DEAL::0.00780764
+DEAL::0.00200193
+DEAL::0.00320309
+DEAL::0.00200193
+DEAL::0.0144302
+DEAL::0.0230884
+DEAL::0.0144302
+DEAL::0.0124922
+DEAL::0.0199876
+DEAL::0.0124922
+DEAL::0.00320309
+DEAL::0.00512495
+DEAL::0.00320309
+DEAL::0.00901889
+DEAL::0.0144302
+DEAL::0.00901889
+DEAL::0.00780764
+DEAL::0.0124922
+DEAL::0.00780764
+DEAL::0.00200193
+DEAL::0.00320309
+DEAL::0.00200193
+DEAL::0.00901889
+DEAL::0.0144302
+DEAL::0.00901889
+DEAL::0.00780764
+DEAL::0.0124922
+DEAL::0.00780764
+DEAL::0.00200193
+DEAL::0.00320309
+DEAL::0.00200193
+DEAL::0.0144302
+DEAL::0.0230884
+DEAL::0.0144302
+DEAL::0.0124922
+DEAL::0.0199876
+DEAL::0.0124922
+DEAL::0.00320309
+DEAL::0.00512495
+DEAL::0.00320309
+DEAL::0.00901889
+DEAL::0.0144302
+DEAL::0.00901889
+DEAL::0.00780764
+DEAL::0.0124922
+DEAL::0.00780764
+DEAL::0.00200193
+DEAL::0.00320309
+DEAL::0.00200193
+DEAL::0.00901889
+DEAL::0.00780764
+DEAL::0.00200193
+DEAL::0.0144302
+DEAL::0.0124922
+DEAL::0.00320309
+DEAL::0.00901889
+DEAL::0.00780764
+DEAL::0.00200193
+DEAL::0.0144302
+DEAL::0.0124922
+DEAL::0.00320309
+DEAL::0.0230884
+DEAL::0.0199876
+DEAL::0.00512495
+DEAL::0.0144302
+DEAL::0.0124922
+DEAL::0.00320309
+DEAL::0.00901889
+DEAL::0.00780764
+DEAL::0.00200193
+DEAL::0.0144302
+DEAL::0.0124922
+DEAL::0.00320309
+DEAL::0.00901889
+DEAL::0.00780764
+DEAL::0.00200193
+DEAL::0.00901889
+DEAL::0.0144302
+DEAL::0.00901889
+DEAL::0.00780764
+DEAL::0.0124922
+DEAL::0.00780764
+DEAL::0.00200193
+DEAL::0.00320309
+DEAL::0.00200193
+DEAL::0.0144302
+DEAL::0.0230884
+DEAL::0.0144302
+DEAL::0.0124922
+DEAL::0.0199876
+DEAL::0.0124922
+DEAL::0.00320309
+DEAL::0.00512495
+DEAL::0.00320309
+DEAL::0.00901889
+DEAL::0.0144302
+DEAL::0.00901889
+DEAL::0.00780764
+DEAL::0.0124922
+DEAL::0.00780764
+DEAL::0.00200193
+DEAL::0.00320309
+DEAL::0.00200193
+DEAL::0.00901889
+DEAL::0.00780764
+DEAL::0.00200193
+DEAL::0.0144302
+DEAL::0.0124922
+DEAL::0.00320309
+DEAL::0.00901889
+DEAL::0.00780764
+DEAL::0.00200193
+DEAL::0.0144302
+DEAL::0.0124922
+DEAL::0.00320309
+DEAL::0.0230884
+DEAL::0.0199876
+DEAL::0.00512495
+DEAL::0.0144302
+DEAL::0.0124922
+DEAL::0.00320309
+DEAL::0.00901889
+DEAL::0.00780764
+DEAL::0.00200193
+DEAL::0.0144302
+DEAL::0.0124922
+DEAL::0.00320309
+DEAL::0.00901889
+DEAL::0.00780764
+DEAL::0.00200193

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.