]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Fix subdivided_parallelepiped
authorTimo Heister <timo.heister@gmail.com>
Tue, 13 Oct 2015 18:01:45 +0000 (14:01 -0400)
committerTimo Heister <timo.heister@gmail.com>
Tue, 13 Oct 2015 18:04:04 +0000 (14:04 -0400)
- missing reorder_cells() would generate unconnected cells
- faces had wrong boundary indicators (if colorizing)

source/grid/grid_generator.cc
tests/grid/grid_parallelepiped_04.cc [new file with mode: 0644]
tests/grid/grid_parallelepiped_04.output [new file with mode: 0644]

index b0ce6ad61dfed4f15d0ef4a0d5dd1c9b20c3def7..2173a50d25171f361f970d4a5f7876cd64e33a38 100644 (file)
@@ -983,11 +983,26 @@ namespace GridGenerator
       }
 
     // Create triangulation
+    // reorder the cells to ensure that they satisfy the convention for
+    // edge and face directions
+    GridReordering<dim>::reorder_cells(cells, true);
     tria.create_triangulation (points, cells, SubCellData());
 
     // Finally assign boundary indicators according to hyper_rectangle
     if (colorize)
-      colorize_hyper_rectangle (tria);
+      {
+        typename Triangulation<dim>::active_cell_iterator
+        cell = tria.begin_active(),
+        endc = tria.end();
+        for (; cell!=endc; ++cell)
+          {
+            for (unsigned int face=0; face<GeometryInfo<dim>::faces_per_cell; ++face)
+              {
+                if (cell->face(face)->at_boundary())
+                  cell->face(face)->set_boundary_id(face);
+              }
+          }
+      }
   }
 
 
diff --git a/tests/grid/grid_parallelepiped_04.cc b/tests/grid/grid_parallelepiped_04.cc
new file mode 100644 (file)
index 0000000..d672c26
--- /dev/null
@@ -0,0 +1,168 @@
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2013 - 2014 by the deal.II authors
+//
+// This file is part of the deal.II library.
+//
+// The deal.II library is free software; you can use it, redistribute
+// it, and/or modify it under the terms of the GNU Lesser General
+// Public License as published by the Free Software Foundation; either
+// version 2.1 of the License, or (at your option) any later version.
+// The full text of the license can be found in the file LICENSE at
+// the top level of the deal.II distribution.
+//
+// ---------------------------------------------------------------------
+
+// check boundary indicators of colored subdivided_parallelepiped
+
+#include "../tests.h"
+#include <deal.II/grid/tria_accessor.h>
+#include <deal.II/grid/tria_iterator.h>
+#include <deal.II/grid/tria.h>
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/grid_out.h>
+#include <deal.II/base/point.h>
+#include <deal.II/base/tensor.h>
+#include <deal.II/base/logstream.h>
+
+#include <deal.II/dofs/dof_handler.h>
+#include <deal.II/lac/constraint_matrix.h>
+#include <deal.II/dofs/dof_accessor.h>
+#include <deal.II/dofs/dof_tools.h>
+#include <deal.II/fe/fe_q.h>
+#include <deal.II/fe/fe_values.h>
+#include <deal.II/numerics/data_out.h>
+#include <deal.II/numerics/vector_tools.h>
+
+
+
+#include <cmath>
+#include <cstdlib>
+
+#include <fstream>
+#include <iostream>
+#include <iomanip>
+#include <map>
+
+template<int dim>
+void check_parallelepiped (bool colorize, bool log, const unsigned int (&subd)[dim])
+{
+  deallog << "* checking dim=" << dim
+          << " subd=";
+  for (unsigned int i=0; i<dim; ++i)
+    deallog << subd[i] << " ";
+  deallog << std::endl;
+
+  // Data structure defining dim coordinates that make up a
+  // parallelepiped.
+  Point<dim> (corners) [dim];
+
+  // build corners for this particular dim:
+  switch (dim)
+    {
+    case 1:
+      corners[0] = Point<dim> (0.5);
+      break;
+
+    case 2:
+      corners[0] = Point<dim> (0.0, 0.5);
+      corners[1] = Point<dim> (0.5, 0.0);
+      break;
+
+    case 3:
+      corners[0] = Point<dim> (0.0, 0.3, 0.5);
+      corners[1] = Point<dim> (0.4, 0.0, 0.5);
+      corners[2] = Point<dim> (0.4, 0.3, 0.0);
+      break;
+
+    default:
+      Assert (false, ExcInternalError ());
+    }
+
+  Triangulation<dim> triangulation;
+
+  GridGenerator::subdivided_parallelepiped (triangulation, subd, corners, colorize);
+
+  {
+
+    std::map<unsigned int, unsigned int> boundary_count;
+    typename Triangulation<dim>::active_cell_iterator
+    cell = triangulation.begin_active(),
+    endc = triangulation.end();
+    for (; cell!=endc; ++cell)
+      {
+        for (unsigned int face=0; face<GeometryInfo<dim>::faces_per_cell; ++face)
+          {
+            if (cell->face(face)->at_boundary())
+              {
+                boundary_count[cell->face(face)->boundary_id()]++;
+                deallog << " center: " << cell->face(face)->center()
+                        << " id: " << (int)cell->face(face)->boundary_id()
+                        << std::endl;
+              }
+
+          }
+      }
+
+    deallog << " boundary indicators: ";
+    for (std::map<unsigned int, unsigned int>::iterator it=boundary_count.begin();
+         it!=boundary_count.end();
+         ++it)
+      {
+        deallog << it->first << "(" << it->second << " times) ";
+      }
+    deallog << std::endl;
+  }
+
+  GridOut grid_out;
+
+  if (log)
+    {
+      FE_Q<dim> fe(2);
+      DoFHandler<dim> dh (triangulation);
+      dh.distribute_dofs (fe);
+      DataOut<dim> d_o;
+      d_o.attach_dof_handler (dh);
+      Vector<double> vec (dh.n_dofs());
+      ConstraintMatrix constraints;
+      for (unsigned int c=0; c<6; ++c)
+        VectorTools::interpolate_boundary_values (dh,
+                                                  c,
+                                                  ConstantFunction<dim>(c),
+                                                  constraints);
+      constraints.close();
+      constraints.distribute(vec);
+
+      d_o.add_data_vector (vec, "v");
+      d_o.build_patches (2);
+      char fname[]="0.vtk";
+      static int counter=0;
+      fname[0]+=counter;
+      ++counter;
+
+      std::ofstream ss(fname);
+
+      // d_o.write_vtk(ss);
+      d_o.write_vtk(deallog.get_file_stream());
+    }
+
+}
+
+int main ()
+{
+  initlog(true);
+
+  //check_parallelepiped<1> (false, true);
+  //check_parallelepiped<2> (false, true);
+  for (unsigned int subd=1; subd<=3; ++subd)
+    {
+      unsigned int subdivisions[3]= {subd,subd,subd};
+      check_parallelepiped<3> (true, subd==2, subdivisions);
+    }
+  for (unsigned int subd=1; subd<=3; ++subd)
+    {
+      unsigned int subdivisions[3]= {1,2,subd};
+      check_parallelepiped<3> (true, false, subdivisions);
+    }
+
+}
diff --git a/tests/grid/grid_parallelepiped_04.output b/tests/grid/grid_parallelepiped_04.output
new file mode 100644 (file)
index 0000000..6975d35
--- /dev/null
@@ -0,0 +1,440 @@
+
+DEAL::* checking dim=3 subd=1 1 1 
+DEAL:: center: 0.400000 0.150000 0.250000 id: 0
+DEAL:: center: 0.400000 0.450000 0.750000 id: 1
+DEAL:: center: 0.200000 0.300000 0.250000 id: 2
+DEAL:: center: 0.600000 0.300000 0.750000 id: 3
+DEAL:: center: 0.200000 0.150000 0.500000 id: 4
+DEAL:: center: 0.600000 0.450000 0.500000 id: 5
+DEAL:: boundary indicators: 0(1 times) 1(1 times) 2(1 times) 3(1 times) 4(1 times) 5(1 times) 
+DEAL::* checking dim=3 subd=2 2 2 
+DEAL:: center: 0.200000 0.0750000 0.125000 id: 0
+DEAL:: center: 0.100000 0.150000 0.125000 id: 2
+DEAL:: center: 0.100000 0.0750000 0.250000 id: 4
+DEAL:: center: 0.200000 0.375000 0.625000 id: 1
+DEAL:: center: 0.100000 0.300000 0.375000 id: 2
+DEAL:: center: 0.100000 0.225000 0.500000 id: 4
+DEAL:: center: 0.400000 0.0750000 0.375000 id: 0
+DEAL:: center: 0.500000 0.150000 0.625000 id: 3
+DEAL:: center: 0.300000 0.0750000 0.500000 id: 4
+DEAL:: center: 0.400000 0.375000 0.875000 id: 1
+DEAL:: center: 0.500000 0.300000 0.875000 id: 3
+DEAL:: center: 0.300000 0.225000 0.750000 id: 4
+DEAL:: center: 0.400000 0.225000 0.125000 id: 0
+DEAL:: center: 0.300000 0.300000 0.125000 id: 2
+DEAL:: center: 0.500000 0.375000 0.250000 id: 5
+DEAL:: center: 0.400000 0.525000 0.625000 id: 1
+DEAL:: center: 0.300000 0.450000 0.375000 id: 2
+DEAL:: center: 0.500000 0.525000 0.500000 id: 5
+DEAL:: center: 0.600000 0.225000 0.375000 id: 0
+DEAL:: center: 0.700000 0.300000 0.625000 id: 3
+DEAL:: center: 0.700000 0.375000 0.500000 id: 5
+DEAL:: center: 0.600000 0.525000 0.875000 id: 1
+DEAL:: center: 0.700000 0.450000 0.875000 id: 3
+DEAL:: center: 0.700000 0.525000 0.750000 id: 5
+DEAL:: boundary indicators: 0(4 times) 1(4 times) 2(4 times) 3(4 times) 4(4 times) 5(4 times) 
+# vtk DataFile Version 3.0
+#This file was generated 
+ASCII
+DATASET UNSTRUCTURED_GRID
+
+POINTS 216 double
+0.00000 0.00000 0.00000
+0.00000 0.0750000 0.125000
+0.00000 0.150000 0.250000
+0.100000 0.00000 0.125000
+0.100000 0.0750000 0.250000
+0.100000 0.150000 0.375000
+0.200000 0.00000 0.250000
+0.200000 0.0750000 0.375000
+0.200000 0.150000 0.500000
+0.100000 0.0750000 0.00000
+0.100000 0.150000 0.125000
+0.100000 0.225000 0.250000
+0.200000 0.0750000 0.125000
+0.200000 0.150000 0.250000
+0.200000 0.225000 0.375000
+0.300000 0.0750000 0.250000
+0.300000 0.150000 0.375000
+0.300000 0.225000 0.500000
+0.200000 0.150000 0.00000
+0.200000 0.225000 0.125000
+0.200000 0.300000 0.250000
+0.300000 0.150000 0.125000
+0.300000 0.225000 0.250000
+0.300000 0.300000 0.375000
+0.400000 0.150000 0.250000
+0.400000 0.225000 0.375000
+0.400000 0.300000 0.500000
+0.00000 0.150000 0.250000
+0.00000 0.225000 0.375000
+0.00000 0.300000 0.500000
+0.100000 0.150000 0.375000
+0.100000 0.225000 0.500000
+0.100000 0.300000 0.625000
+0.200000 0.150000 0.500000
+0.200000 0.225000 0.625000
+0.200000 0.300000 0.750000
+0.100000 0.225000 0.250000
+0.100000 0.300000 0.375000
+0.100000 0.375000 0.500000
+0.200000 0.225000 0.375000
+0.200000 0.300000 0.500000
+0.200000 0.375000 0.625000
+0.300000 0.225000 0.500000
+0.300000 0.300000 0.625000
+0.300000 0.375000 0.750000
+0.200000 0.300000 0.250000
+0.200000 0.375000 0.375000
+0.200000 0.450000 0.500000
+0.300000 0.300000 0.375000
+0.300000 0.375000 0.500000
+0.300000 0.450000 0.625000
+0.400000 0.300000 0.500000
+0.400000 0.375000 0.625000
+0.400000 0.450000 0.750000
+0.200000 0.00000 0.250000
+0.200000 0.0750000 0.375000
+0.200000 0.150000 0.500000
+0.300000 0.00000 0.375000
+0.300000 0.0750000 0.500000
+0.300000 0.150000 0.625000
+0.400000 0.00000 0.500000
+0.400000 0.0750000 0.625000
+0.400000 0.150000 0.750000
+0.300000 0.0750000 0.250000
+0.300000 0.150000 0.375000
+0.300000 0.225000 0.500000
+0.400000 0.0750000 0.375000
+0.400000 0.150000 0.500000
+0.400000 0.225000 0.625000
+0.500000 0.0750000 0.500000
+0.500000 0.150000 0.625000
+0.500000 0.225000 0.750000
+0.400000 0.150000 0.250000
+0.400000 0.225000 0.375000
+0.400000 0.300000 0.500000
+0.500000 0.150000 0.375000
+0.500000 0.225000 0.500000
+0.500000 0.300000 0.625000
+0.600000 0.150000 0.500000
+0.600000 0.225000 0.625000
+0.600000 0.300000 0.750000
+0.200000 0.150000 0.500000
+0.200000 0.225000 0.625000
+0.200000 0.300000 0.750000
+0.300000 0.150000 0.625000
+0.300000 0.225000 0.750000
+0.300000 0.300000 0.875000
+0.400000 0.150000 0.750000
+0.400000 0.225000 0.875000
+0.400000 0.300000 1.00000
+0.300000 0.225000 0.500000
+0.300000 0.300000 0.625000
+0.300000 0.375000 0.750000
+0.400000 0.225000 0.625000
+0.400000 0.300000 0.750000
+0.400000 0.375000 0.875000
+0.500000 0.225000 0.750000
+0.500000 0.300000 0.875000
+0.500000 0.375000 1.00000
+0.400000 0.300000 0.500000
+0.400000 0.375000 0.625000
+0.400000 0.450000 0.750000
+0.500000 0.300000 0.625000
+0.500000 0.375000 0.750000
+0.500000 0.450000 0.875000
+0.600000 0.300000 0.750000
+0.600000 0.375000 0.875000
+0.600000 0.450000 1.00000
+0.200000 0.150000 0.00000
+0.200000 0.225000 0.125000
+0.200000 0.300000 0.250000
+0.300000 0.150000 0.125000
+0.300000 0.225000 0.250000
+0.300000 0.300000 0.375000
+0.400000 0.150000 0.250000
+0.400000 0.225000 0.375000
+0.400000 0.300000 0.500000
+0.300000 0.225000 0.00000
+0.300000 0.300000 0.125000
+0.300000 0.375000 0.250000
+0.400000 0.225000 0.125000
+0.400000 0.300000 0.250000
+0.400000 0.375000 0.375000
+0.500000 0.225000 0.250000
+0.500000 0.300000 0.375000
+0.500000 0.375000 0.500000
+0.400000 0.300000 0.00000
+0.400000 0.375000 0.125000
+0.400000 0.450000 0.250000
+0.500000 0.300000 0.125000
+0.500000 0.375000 0.250000
+0.500000 0.450000 0.375000
+0.600000 0.300000 0.250000
+0.600000 0.375000 0.375000
+0.600000 0.450000 0.500000
+0.200000 0.300000 0.250000
+0.200000 0.375000 0.375000
+0.200000 0.450000 0.500000
+0.300000 0.300000 0.375000
+0.300000 0.375000 0.500000
+0.300000 0.450000 0.625000
+0.400000 0.300000 0.500000
+0.400000 0.375000 0.625000
+0.400000 0.450000 0.750000
+0.300000 0.375000 0.250000
+0.300000 0.450000 0.375000
+0.300000 0.525000 0.500000
+0.400000 0.375000 0.375000
+0.400000 0.450000 0.500000
+0.400000 0.525000 0.625000
+0.500000 0.375000 0.500000
+0.500000 0.450000 0.625000
+0.500000 0.525000 0.750000
+0.400000 0.450000 0.250000
+0.400000 0.525000 0.375000
+0.400000 0.600000 0.500000
+0.500000 0.450000 0.375000
+0.500000 0.525000 0.500000
+0.500000 0.600000 0.625000
+0.600000 0.450000 0.500000
+0.600000 0.525000 0.625000
+0.600000 0.600000 0.750000
+0.400000 0.150000 0.250000
+0.400000 0.225000 0.375000
+0.400000 0.300000 0.500000
+0.500000 0.150000 0.375000
+0.500000 0.225000 0.500000
+0.500000 0.300000 0.625000
+0.600000 0.150000 0.500000
+0.600000 0.225000 0.625000
+0.600000 0.300000 0.750000
+0.500000 0.225000 0.250000
+0.500000 0.300000 0.375000
+0.500000 0.375000 0.500000
+0.600000 0.225000 0.375000
+0.600000 0.300000 0.500000
+0.600000 0.375000 0.625000
+0.700000 0.225000 0.500000
+0.700000 0.300000 0.625000
+0.700000 0.375000 0.750000
+0.600000 0.300000 0.250000
+0.600000 0.375000 0.375000
+0.600000 0.450000 0.500000
+0.700000 0.300000 0.375000
+0.700000 0.375000 0.500000
+0.700000 0.450000 0.625000
+0.800000 0.300000 0.500000
+0.800000 0.375000 0.625000
+0.800000 0.450000 0.750000
+0.400000 0.300000 0.500000
+0.400000 0.375000 0.625000
+0.400000 0.450000 0.750000
+0.500000 0.300000 0.625000
+0.500000 0.375000 0.750000
+0.500000 0.450000 0.875000
+0.600000 0.300000 0.750000
+0.600000 0.375000 0.875000
+0.600000 0.450000 1.00000
+0.500000 0.375000 0.500000
+0.500000 0.450000 0.625000
+0.500000 0.525000 0.750000
+0.600000 0.375000 0.625000
+0.600000 0.450000 0.750000
+0.600000 0.525000 0.875000
+0.700000 0.375000 0.750000
+0.700000 0.450000 0.875000
+0.700000 0.525000 1.00000
+0.600000 0.450000 0.500000
+0.600000 0.525000 0.625000
+0.600000 0.600000 0.750000
+0.700000 0.450000 0.625000
+0.700000 0.525000 0.750000
+0.700000 0.600000 0.875000
+0.800000 0.450000 0.750000
+0.800000 0.525000 0.875000
+0.800000 0.600000 1.00000
+
+CELLS 64 576
+8      0       1       4       3       9       10      13      12
+8      1       2       5       4       10      11      14      13
+8      3       4       7       6       12      13      16      15
+8      4       5       8       7       13      14      17      16
+8      9       10      13      12      18      19      22      21
+8      10      11      14      13      19      20      23      22
+8      12      13      16      15      21      22      25      24
+8      13      14      17      16      22      23      26      25
+8      27      28      31      30      36      37      40      39
+8      28      29      32      31      37      38      41      40
+8      30      31      34      33      39      40      43      42
+8      31      32      35      34      40      41      44      43
+8      36      37      40      39      45      46      49      48
+8      37      38      41      40      46      47      50      49
+8      39      40      43      42      48      49      52      51
+8      40      41      44      43      49      50      53      52
+8      54      55      58      57      63      64      67      66
+8      55      56      59      58      64      65      68      67
+8      57      58      61      60      66      67      70      69
+8      58      59      62      61      67      68      71      70
+8      63      64      67      66      72      73      76      75
+8      64      65      68      67      73      74      77      76
+8      66      67      70      69      75      76      79      78
+8      67      68      71      70      76      77      80      79
+8      81      82      85      84      90      91      94      93
+8      82      83      86      85      91      92      95      94
+8      84      85      88      87      93      94      97      96
+8      85      86      89      88      94      95      98      97
+8      90      91      94      93      99      100     103     102
+8      91      92      95      94      100     101     104     103
+8      93      94      97      96      102     103     106     105
+8      94      95      98      97      103     104     107     106
+8      108     109     112     111     117     118     121     120
+8      109     110     113     112     118     119     122     121
+8      111     112     115     114     120     121     124     123
+8      112     113     116     115     121     122     125     124
+8      117     118     121     120     126     127     130     129
+8      118     119     122     121     127     128     131     130
+8      120     121     124     123     129     130     133     132
+8      121     122     125     124     130     131     134     133
+8      135     136     139     138     144     145     148     147
+8      136     137     140     139     145     146     149     148
+8      138     139     142     141     147     148     151     150
+8      139     140     143     142     148     149     152     151
+8      144     145     148     147     153     154     157     156
+8      145     146     149     148     154     155     158     157
+8      147     148     151     150     156     157     160     159
+8      148     149     152     151     157     158     161     160
+8      162     163     166     165     171     172     175     174
+8      163     164     167     166     172     173     176     175
+8      165     166     169     168     174     175     178     177
+8      166     167     170     169     175     176     179     178
+8      171     172     175     174     180     181     184     183
+8      172     173     176     175     181     182     185     184
+8      174     175     178     177     183     184     187     186
+8      175     176     179     178     184     185     188     187
+8      189     190     193     192     198     199     202     201
+8      190     191     194     193     199     200     203     202
+8      192     193     196     195     201     202     205     204
+8      193     194     197     196     202     203     206     205
+8      198     199     202     201     207     208     211     210
+8      199     200     203     202     208     209     212     211
+8      201     202     205     204     210     211     214     213
+8      202     203     206     205     211     212     215     214
+
+CELL_TYPES 64
+ 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12
+POINT_DATA 216
+SCALARS v double 1
+LOOKUP_TABLE default
+0.00000 2.00000 2.00000 0.00000 4.00000 4.00000 0.00000 4.00000 4.00000 0.00000 2.00000 2.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 2.00000 2.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 2.00000 2.00000 1.00000 4.00000 4.00000 1.00000 4.00000 4.00000 1.00000 2.00000 2.00000 1.00000 0.00000 0.00000 1.00000 0.00000 0.00000 1.00000 2.00000 2.00000 1.00000 0.00000 0.00000 1.00000 0.00000 0.00000 1.00000 0.00000 4.00000 4.00000 0.00000 4.00000 4.00000 0.00000 3.00000 3.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 3.00000 3.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 3.00000 3.00000 4.00000 4.00000 1.00000 4.00000 4.00000 1.00000 3.00000 3.00000 1.00000 0.00000 0.00000 1.00000 0.00000 0.00000 1.00000 3.00000 3.00000 1.00000 0.00000 0.00000 1.00000 0.00000 0.00000 1.00000 3.00000 3.00000 1.00000 0.00000 2.00000 2.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 2.00000 2.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 2.00000 2.00000 0.00000 5.00000 5.00000 0.00000 5.00000 5.00000 2.00000 2.00000 1.00000 0.00000 0.00000 1.00000 0.00000 0.00000 1.00000 2.00000 2.00000 1.00000 0.00000 0.00000 1.00000 0.00000 0.00000 1.00000 2.00000 2.00000 1.00000 5.00000 5.00000 1.00000 5.00000 5.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 3.00000 3.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 3.00000 3.00000 0.00000 5.00000 5.00000 0.00000 5.00000 5.00000 0.00000 3.00000 3.00000 0.00000 0.00000 1.00000 0.00000 0.00000 1.00000 3.00000 3.00000 1.00000 0.00000 0.00000 1.00000 0.00000 0.00000 1.00000 3.00000 3.00000 1.00000 5.00000 5.00000 1.00000 5.00000 5.00000 1.00000 3.00000 3.00000 1.00000 
+DEAL::* checking dim=3 subd=3 3 3 
+DEAL:: center: 0.133333 0.0500000 0.0833333 id: 0
+DEAL:: center: 0.0666667 0.100000 0.0833333 id: 2
+DEAL:: center: 0.0666667 0.0500000 0.166667 id: 4
+DEAL:: center: 0.0666667 0.200000 0.250000 id: 2
+DEAL:: center: 0.0666667 0.150000 0.333333 id: 4
+DEAL:: center: 0.133333 0.350000 0.583333 id: 1
+DEAL:: center: 0.0666667 0.300000 0.416667 id: 2
+DEAL:: center: 0.0666667 0.250000 0.500000 id: 4
+DEAL:: center: 0.266667 0.0500000 0.250000 id: 0
+DEAL:: center: 0.200000 0.0500000 0.333333 id: 4
+DEAL:: center: 0.200000 0.150000 0.500000 id: 4
+DEAL:: center: 0.266667 0.350000 0.750000 id: 1
+DEAL:: center: 0.200000 0.250000 0.666667 id: 4
+DEAL:: center: 0.400000 0.0500000 0.416667 id: 0
+DEAL:: center: 0.466667 0.100000 0.583333 id: 3
+DEAL:: center: 0.333333 0.0500000 0.500000 id: 4
+DEAL:: center: 0.466667 0.200000 0.750000 id: 3
+DEAL:: center: 0.333333 0.150000 0.666667 id: 4
+DEAL:: center: 0.400000 0.350000 0.916667 id: 1
+DEAL:: center: 0.466667 0.300000 0.916667 id: 3
+DEAL:: center: 0.333333 0.250000 0.833333 id: 4
+DEAL:: center: 0.266667 0.150000 0.0833333 id: 0
+DEAL:: center: 0.200000 0.200000 0.0833333 id: 2
+DEAL:: center: 0.200000 0.300000 0.250000 id: 2
+DEAL:: center: 0.266667 0.450000 0.583333 id: 1
+DEAL:: center: 0.200000 0.400000 0.416667 id: 2
+DEAL:: center: 0.400000 0.150000 0.250000 id: 0
+DEAL:: center: 0.400000 0.450000 0.750000 id: 1
+DEAL:: center: 0.533333 0.150000 0.416667 id: 0
+DEAL:: center: 0.600000 0.200000 0.583333 id: 3
+DEAL:: center: 0.600000 0.300000 0.750000 id: 3
+DEAL:: center: 0.533333 0.450000 0.916667 id: 1
+DEAL:: center: 0.600000 0.400000 0.916667 id: 3
+DEAL:: center: 0.400000 0.250000 0.0833333 id: 0
+DEAL:: center: 0.333333 0.300000 0.0833333 id: 2
+DEAL:: center: 0.466667 0.350000 0.166667 id: 5
+DEAL:: center: 0.333333 0.400000 0.250000 id: 2
+DEAL:: center: 0.466667 0.450000 0.333333 id: 5
+DEAL:: center: 0.400000 0.550000 0.583333 id: 1
+DEAL:: center: 0.333333 0.500000 0.416667 id: 2
+DEAL:: center: 0.466667 0.550000 0.500000 id: 5
+DEAL:: center: 0.533333 0.250000 0.250000 id: 0
+DEAL:: center: 0.600000 0.350000 0.333333 id: 5
+DEAL:: center: 0.600000 0.450000 0.500000 id: 5
+DEAL:: center: 0.533333 0.550000 0.750000 id: 1
+DEAL:: center: 0.600000 0.550000 0.666667 id: 5
+DEAL:: center: 0.666667 0.250000 0.416667 id: 0
+DEAL:: center: 0.733333 0.300000 0.583333 id: 3
+DEAL:: center: 0.733333 0.350000 0.500000 id: 5
+DEAL:: center: 0.733333 0.400000 0.750000 id: 3
+DEAL:: center: 0.733333 0.450000 0.666667 id: 5
+DEAL:: center: 0.666667 0.550000 0.916667 id: 1
+DEAL:: center: 0.733333 0.500000 0.916667 id: 3
+DEAL:: center: 0.733333 0.550000 0.833333 id: 5
+DEAL:: boundary indicators: 0(9 times) 1(9 times) 2(9 times) 3(9 times) 4(9 times) 5(9 times) 
+DEAL::* checking dim=3 subd=1 2 1 
+DEAL:: center: 0.300000 0.150000 0.125000 id: 0
+DEAL:: center: 0.300000 0.450000 0.625000 id: 1
+DEAL:: center: 0.200000 0.300000 0.250000 id: 2
+DEAL:: center: 0.100000 0.150000 0.375000 id: 4
+DEAL:: center: 0.500000 0.450000 0.375000 id: 5
+DEAL:: center: 0.500000 0.150000 0.375000 id: 0
+DEAL:: center: 0.500000 0.450000 0.875000 id: 1
+DEAL:: center: 0.600000 0.300000 0.750000 id: 3
+DEAL:: center: 0.300000 0.150000 0.625000 id: 4
+DEAL:: center: 0.700000 0.450000 0.625000 id: 5
+DEAL:: boundary indicators: 0(2 times) 1(2 times) 2(1 times) 3(1 times) 4(2 times) 5(2 times) 
+DEAL::* checking dim=3 subd=1 2 2 
+DEAL:: center: 0.200000 0.0750000 0.125000 id: 0
+DEAL:: center: 0.200000 0.375000 0.625000 id: 1
+DEAL:: center: 0.100000 0.225000 0.250000 id: 2
+DEAL:: center: 0.100000 0.150000 0.375000 id: 4
+DEAL:: center: 0.400000 0.0750000 0.375000 id: 0
+DEAL:: center: 0.400000 0.375000 0.875000 id: 1
+DEAL:: center: 0.500000 0.225000 0.750000 id: 3
+DEAL:: center: 0.300000 0.150000 0.625000 id: 4
+DEAL:: center: 0.400000 0.225000 0.125000 id: 0
+DEAL:: center: 0.400000 0.525000 0.625000 id: 1
+DEAL:: center: 0.300000 0.375000 0.250000 id: 2
+DEAL:: center: 0.500000 0.450000 0.375000 id: 5
+DEAL:: center: 0.600000 0.225000 0.375000 id: 0
+DEAL:: center: 0.600000 0.525000 0.875000 id: 1
+DEAL:: center: 0.700000 0.375000 0.750000 id: 3
+DEAL:: center: 0.700000 0.450000 0.625000 id: 5
+DEAL:: boundary indicators: 0(4 times) 1(4 times) 2(2 times) 3(2 times) 4(2 times) 5(2 times) 
+DEAL::* checking dim=3 subd=1 2 3 
+DEAL:: center: 0.166667 0.0500000 0.125000 id: 0
+DEAL:: center: 0.166667 0.350000 0.625000 id: 1
+DEAL:: center: 0.0666667 0.200000 0.250000 id: 2
+DEAL:: center: 0.100000 0.150000 0.375000 id: 4
+DEAL:: center: 0.366667 0.0500000 0.375000 id: 0
+DEAL:: center: 0.366667 0.350000 0.875000 id: 1
+DEAL:: center: 0.466667 0.200000 0.750000 id: 3
+DEAL:: center: 0.300000 0.150000 0.625000 id: 4
+DEAL:: center: 0.300000 0.150000 0.125000 id: 0
+DEAL:: center: 0.300000 0.450000 0.625000 id: 1
+DEAL:: center: 0.200000 0.300000 0.250000 id: 2
+DEAL:: center: 0.500000 0.150000 0.375000 id: 0
+DEAL:: center: 0.500000 0.450000 0.875000 id: 1
+DEAL:: center: 0.600000 0.300000 0.750000 id: 3
+DEAL:: center: 0.433333 0.250000 0.125000 id: 0
+DEAL:: center: 0.433333 0.550000 0.625000 id: 1
+DEAL:: center: 0.333333 0.400000 0.250000 id: 2
+DEAL:: center: 0.500000 0.450000 0.375000 id: 5
+DEAL:: center: 0.633333 0.250000 0.375000 id: 0
+DEAL:: center: 0.633333 0.550000 0.875000 id: 1
+DEAL:: center: 0.733333 0.400000 0.750000 id: 3
+DEAL:: center: 0.700000 0.450000 0.625000 id: 5
+DEAL:: boundary indicators: 0(6 times) 1(6 times) 2(3 times) 3(3 times) 4(2 times) 5(2 times) 

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.