]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Add vector and tensor properties
authorRene Gassmoeller <rene.gassmoeller@mailbox.org>
Thu, 28 May 2020 00:28:09 +0000 (17:28 -0700)
committerRene Gassmoeller <rene.gassmoeller@mailbox.org>
Thu, 28 May 2020 00:34:51 +0000 (17:34 -0700)
doc/news/changes/minor/20200527Gassmoeller [new file with mode: 0644]
include/deal.II/particles/data_out.h
source/base/data_out_base.cc
source/particles/data_out.cc
tests/particles/data_out_03.cc [new file with mode: 0644]
tests/particles/data_out_03.output [new file with mode: 0644]
tests/particles/data_out_04.cc [new file with mode: 0644]
tests/particles/data_out_04.output [new file with mode: 0644]

diff --git a/doc/news/changes/minor/20200527Gassmoeller b/doc/news/changes/minor/20200527Gassmoeller
new file mode 100644 (file)
index 0000000..8464738
--- /dev/null
@@ -0,0 +1,5 @@
+New: The class Particles::DataOut can now output particle properties as
+scalars, vectors, or tensors, depending on the arguments handed over to the
+Particles::DataOut::build_patches() function.
+<br>
+(Rene Gassmoeller, 2020/05/27)
index cce7a0718e979747f7ca3ca690bf6f75b561d775..cfbc789f56dd4e47078ccb775294caf97db74287 100644 (file)
@@ -64,15 +64,22 @@ namespace Particles
      * A dim=0 patch is built for each particle. The position of the particle is
      * used to build the node position and the ID of the particle is added as a
      * single data element.
+     * @param [in] data_component_names An optional vector of strings that
+     * describe the properties of each particle. Particle properties will only
+     * be written if this vector
+     * is provided.
+     * @param [in] data_component_interpretations An optional vector that
+     * controls if the particle properties are interpreted as scalars, vectors,
+     * or tensors. Has to be of the same length as @p data_component_names.
      *
      * @author Bruno Blais, Luca Heltai 2019
      */
     void
     build_patches(const Particles::ParticleHandler<dim, spacedim> &particles,
-                  const std::vector<std::string> &                 names = {},
+                  const std::vector<std::string> &data_component_names = {},
                   const std::vector<
                     DataComponentInterpretation::DataComponentInterpretation>
-                    &data_component_interpretation = {});
+                    &data_component_interpretations = {});
 
   protected:
     /**
@@ -89,18 +96,40 @@ namespace Particles
     virtual std::vector<std::string>
     get_dataset_names() const override;
 
+
+    /**
+     * Overload of the respective DataOutInterface::get_nonscalar_data_ranges()
+     * function. See there for a more extensive documentation.
+     * This function is a reimplementation of the function
+     * DataOut_DoFData::get_nonscalar_data_ranges().
+     */
+    virtual std::vector<
+      std::tuple<unsigned int,
+                 unsigned int,
+                 std::string,
+                 DataComponentInterpretation::DataComponentInterpretation>>
+    get_nonscalar_data_ranges() const override;
+
   private:
     /**
-     * This is a list of patches that is created each time build_patches() is
+     * This is a vector of patches that is created each time build_patches() is
      * called. These patches are used in the output routines of the base
      * classes.
      */
     std::vector<DataOutBase::Patch<0, spacedim>> patches;
 
     /**
-     * A list of field names for all data components stored in patches.
+     * A vector of field names for all data components stored in patches.
      */
     std::vector<std::string> dataset_names;
+
+    /**
+     * A vector that for each of the data components of the
+     * current data set indicates whether they are scalar fields, parts of a
+     * vector-field, or any of the other supported kinds of data.
+     */
+    std::vector<DataComponentInterpretation::DataComponentInterpretation>
+      data_component_interpretations;
   };
 
 } // namespace Particles
index e35adfab5687335f27e45b5f755a22e198e6e19d..a2642d1ef36ef782377098bfa16bcac01c41f9a0 100644 (file)
@@ -5421,7 +5421,7 @@ namespace DataOutBase
                     for (unsigned int c = 0; c < size; ++c)
                       {
                         const auto ind =
-                          Tensor<2, dim>::unrolled_to_component_indices(c);
+                          Tensor<2, spacedim>::unrolled_to_component_indices(c);
                         vtk_data[ind[0]][ind[1]] =
                           data_vectors(first_component + c, n);
                       }
index 3ff8ab4baef3004b0148ee0b34a8038da30c1f45..99d20548edda8aa75e2d99ede56dfd79e5ef70d9 100644 (file)
@@ -16,6 +16,9 @@
 #include <deal.II/particles/data_out.h>
 #include <deal.II/particles/particle_handler.h>
 
+// We use some exceptions declared in this header
+#include <deal.II/numerics/data_out_dof_data.h>
+
 DEAL_II_NAMESPACE_OPEN
 
 namespace Particles
@@ -26,10 +29,10 @@ namespace Particles
     const Particles::ParticleHandler<dim, spacedim> &particles,
     const std::vector<std::string> &                 data_component_names,
     const std::vector<DataComponentInterpretation::DataComponentInterpretation>
-      &data_component_interpretation)
+      &data_component_interpretations_)
   {
-    AssertThrow(
-      data_component_names.size() == data_component_interpretation.size(),
+    Assert(
+      data_component_names.size() == data_component_interpretations_.size(),
       ExcMessage(
         "When calling Particles::DataOut::build_patches with data component "
         "names and interpretations you need to provide as many data component "
@@ -38,7 +41,7 @@ namespace Particles
 
     if (data_component_names.size() > 0)
       {
-        AssertThrow(
+        Assert(
           data_component_names.size() ==
             particles.begin()->get_properties().size(),
           ExcMessage(
@@ -53,6 +56,14 @@ namespace Particles
                          data_component_names.begin(),
                          data_component_names.end());
 
+    data_component_interpretations.clear();
+    data_component_interpretations.emplace_back(
+      DataComponentInterpretation::component_is_scalar);
+    data_component_interpretations.insert(
+      data_component_interpretations.end(),
+      data_component_interpretations_.begin(),
+      data_component_interpretations_.end());
+
     const unsigned int n_property_components = data_component_names.size();
     const unsigned int n_data_components     = dataset_names.size();
 
@@ -100,6 +111,140 @@ namespace Particles
     return dataset_names;
   }
 
+
+
+  template <int dim, int spacedim>
+  std::vector<
+    std::tuple<unsigned int,
+               unsigned int,
+               std::string,
+               DataComponentInterpretation::DataComponentInterpretation>>
+  DataOut<dim, spacedim>::get_nonscalar_data_ranges() const
+  {
+    std::vector<
+      std::tuple<unsigned int,
+                 unsigned int,
+                 std::string,
+                 DataComponentInterpretation::DataComponentInterpretation>>
+      ranges;
+
+    // Make sure the data structures were set up correctly. Since they
+    // can only be filled by build_patches() above, they should have
+    // been checked already.
+    Assert(dataset_names.size() == data_component_interpretations.size(),
+           ExcInternalError());
+
+    // collect the ranges of particle data
+    const unsigned int n_output_components =
+      data_component_interpretations.size();
+    unsigned int output_component = 0;
+    for (unsigned int i = 0; i < n_output_components; /* i is updated below */)
+      // see what kind of data we have here. note that for the purpose of the
+      // current function all we care about is vector data
+      switch (data_component_interpretations[i])
+        {
+          case DataComponentInterpretation::component_is_scalar:
+            {
+              // Just move component forward by one
+              ++i;
+              ++output_component;
+
+              break;
+            }
+          case DataComponentInterpretation::component_is_part_of_vector:
+            {
+              // ensure that there is a continuous number of next space_dim
+              // components that all deal with vectors
+              Assert(
+                i + spacedim <= n_output_components,
+                Exceptions::DataOutImplementation::ExcInvalidVectorDeclaration(
+                  i, dataset_names[i]));
+              for (unsigned int dd = 1; dd < spacedim; ++dd)
+                Assert(
+                  data_component_interpretations[i + dd] ==
+                    DataComponentInterpretation::component_is_part_of_vector,
+                  Exceptions::DataOutImplementation::
+                    ExcInvalidVectorDeclaration(i, dataset_names[i]));
+
+              // all seems right, so figure out whether there is a common
+              // name to these components. if not, leave the name empty and
+              // let the output format writer decide what to do here
+              std::string name = dataset_names[i];
+              for (unsigned int dd = 1; dd < spacedim; ++dd)
+                if (name != dataset_names[i + dd])
+                  {
+                    name = "";
+                    break;
+                  }
+
+              // Finally add a corresponding range.
+              //
+              // This sort of logic is also explained in some detail in
+              //   DataOut::build_one_patch().
+              ranges.emplace_back(std::forward_as_tuple(
+                output_component,
+                output_component + spacedim - 1,
+                name,
+                DataComponentInterpretation::component_is_part_of_vector));
+
+              // increase the 'component' counter by the appropriate amount,
+              // same for 'i', since we have already dealt with all these
+              // components
+              output_component += spacedim;
+              i += spacedim;
+
+              break;
+            }
+
+          case DataComponentInterpretation::component_is_part_of_tensor:
+            {
+              const unsigned int size = spacedim * spacedim;
+              // ensure that there is a continuous number of next
+              // spacedim*spacedim components that all deal with tensors
+              Assert(
+                i + size <= n_output_components,
+                Exceptions::DataOutImplementation::ExcInvalidTensorDeclaration(
+                  i, dataset_names[i]));
+              for (unsigned int dd = 1; dd < size; ++dd)
+                Assert(
+                  data_component_interpretations[i + dd] ==
+                    DataComponentInterpretation::component_is_part_of_tensor,
+                  Exceptions::DataOutImplementation::
+                    ExcInvalidTensorDeclaration(i, dataset_names[i]));
+
+              // all seems right, so figure out whether there is a common
+              // name to these components. if not, leave the name empty and
+              // let the output format writer decide what to do here
+              std::string name = dataset_names[i];
+              for (unsigned int dd = 1; dd < size; ++dd)
+                if (name != dataset_names[i + dd])
+                  {
+                    name = "";
+                    break;
+                  }
+
+              // Finally add a corresponding range.
+              ranges.emplace_back(std::forward_as_tuple(
+                output_component,
+                output_component + size - 1,
+                name,
+                DataComponentInterpretation::component_is_part_of_tensor));
+
+              // increase the 'component' counter by the appropriate amount,
+              // same for 'i', since we have already dealt with all these
+              // components
+              output_component += size;
+              i += size;
+              break;
+            }
+
+          default:
+            Assert(false, ExcNotImplemented());
+        }
+
+    return ranges;
+  }
+
 } // namespace Particles
 
 #include "data_out.inst"
diff --git a/tests/particles/data_out_03.cc b/tests/particles/data_out_03.cc
new file mode 100644 (file)
index 0000000..4644f34
--- /dev/null
@@ -0,0 +1,148 @@
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2017 - 2019 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.md at
+// the top level directory of deal.II.
+//
+// ---------------------------------------------------------------------
+
+
+
+// Create a particle handler then generate an output.
+// Tests that vector particle properties are output correctly.
+
+#include <deal.II/base/data_out_base.h>
+
+#include <deal.II/distributed/tria.h>
+
+#include <deal.II/fe/mapping_q.h>
+
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/grid_tools.h>
+
+#include <deal.II/particles/data_out.h>
+#include <deal.II/particles/particle_handler.h>
+
+#include "../tests.h"
+
+template <int dim, int spacedim>
+void
+create_regular_particle_distribution(
+  Particles::ParticleHandler<dim, spacedim> &particle_handler,
+  const Triangulation<dim, spacedim> &       tr,
+  const unsigned int                         particles_per_direction = 3)
+{
+  for (unsigned int i = 0; i < particles_per_direction; ++i)
+    for (unsigned int j = 0; j < particles_per_direction; ++j)
+      {
+        Point<spacedim> position;
+        Point<dim>      reference_position;
+        unsigned int    id = i * particles_per_direction + j;
+
+        position[0] = static_cast<double>(i) /
+                      static_cast<double>(particles_per_direction - 1);
+        position[1] = static_cast<double>(j) /
+                      static_cast<double>(particles_per_direction - 1);
+
+        if (dim > 2)
+          for (unsigned int k = 0; k < particles_per_direction; ++k)
+            {
+              position[2] = static_cast<double>(j) /
+                            static_cast<double>(particles_per_direction - 1);
+              id = i * particles_per_direction * particles_per_direction +
+                   j * particles_per_direction + k;
+              Particles::Particle<dim, spacedim> particle(position,
+                                                          reference_position,
+                                                          id);
+
+              typename Triangulation<dim, spacedim>::active_cell_iterator cell =
+                GridTools::find_active_cell_around_point(
+                  tr, particle.get_location());
+
+              Particles::ParticleIterator<dim, spacedim> pit =
+                particle_handler.insert_particle(particle, cell);
+
+              for (unsigned int i = 0; i < spacedim; ++i)
+                pit->get_properties()[i] = pit->get_location()[i];
+            }
+        else
+          {
+            Particles::Particle<dim, spacedim> particle(position,
+                                                        reference_position,
+                                                        id);
+
+            typename Triangulation<dim, spacedim>::active_cell_iterator cell =
+              GridTools::find_active_cell_around_point(tr,
+                                                       particle.get_location());
+
+            Particles::ParticleIterator<dim, spacedim> pit =
+              particle_handler.insert_particle(particle, cell);
+
+            for (unsigned int i = 0; i < spacedim; ++i)
+              pit->get_properties()[i] = pit->get_location()[i];
+          }
+      }
+}
+
+template <int dim, int spacedim>
+void
+test()
+{
+  {
+    Triangulation<dim, spacedim> tr;
+
+    GridGenerator::hyper_cube(tr);
+    tr.refine_global(2);
+    MappingQ<dim, spacedim> mapping(1);
+
+    // Create a particle handler using two manually created particles
+    Particles::ParticleHandler<dim, spacedim> particle_handler(tr,
+                                                               mapping,
+                                                               spacedim);
+
+    create_regular_particle_distribution(particle_handler, tr, 2);
+
+    std::vector<std::string> data_names;
+    std::vector<DataComponentInterpretation::DataComponentInterpretation>
+      data_interpretations;
+
+    for (unsigned int i = 0; i < spacedim; ++i)
+      {
+        data_names.emplace_back("property_coord_" + std::to_string(i));
+        data_interpretations.emplace_back(
+          DataComponentInterpretation::component_is_part_of_vector);
+      }
+
+    Particles::DataOut<dim, spacedim> particle_output;
+    particle_output.build_patches(particle_handler,
+                                  data_names,
+                                  data_interpretations);
+    particle_output.write_vtk(deallog.get_file_stream());
+  }
+
+  deallog << "OK" << std::endl;
+}
+
+
+
+int
+main(int argc, char *argv[])
+{
+  initlog();
+  deallog.push("2d/2d");
+  test<2, 2>();
+  deallog.pop();
+  deallog.push("2d/3d");
+  test<2, 3>();
+  deallog.pop();
+  deallog.push("3d/3d");
+  test<3, 3>();
+  deallog.pop();
+}
diff --git a/tests/particles/data_out_03.output b/tests/particles/data_out_03.output
new file mode 100644 (file)
index 0000000..25107a6
--- /dev/null
@@ -0,0 +1,100 @@
+
+# vtk DataFile Version 3.0
+#This file was generated 
+ASCII
+DATASET UNSTRUCTURED_GRID
+
+POINTS 4 double
+0.00000 0.00000 0
+1.00000 0.00000 0
+0.00000 1.00000 0
+1.00000 1.00000 0
+
+CELLS 4 8
+1      0
+1      1
+1      2
+1      3
+
+CELL_TYPES 4
+ 1 1 1 1
+POINT_DATA 4
+VECTORS property_coord_0__property_coord_1 double
+0.00000 0.00000 0
+1.00000 0.00000 0
+0.00000 1.00000 0
+1.00000 1.00000 0
+SCALARS id double 1
+LOOKUP_TABLE default
+0.00000 2.00000 1.00000 3.00000 
+DEAL:2d/2d::OK
+# vtk DataFile Version 3.0
+#This file was generated 
+ASCII
+DATASET UNSTRUCTURED_GRID
+
+POINTS 4 double
+0.00000 0.00000 0.00000
+1.00000 0.00000 0.00000
+0.00000 1.00000 0.00000
+1.00000 1.00000 0.00000
+
+CELLS 4 8
+1      0
+1      1
+1      2
+1      3
+
+CELL_TYPES 4
+ 1 1 1 1
+POINT_DATA 4
+VECTORS property_coord_0__property_coord_1__property_coord_2 double
+0.00000 0.00000 0.00000
+1.00000 0.00000 0.00000
+0.00000 1.00000 0.00000
+1.00000 1.00000 0.00000
+SCALARS id double 1
+LOOKUP_TABLE default
+0.00000 2.00000 1.00000 3.00000 
+DEAL:2d/3d::OK
+# vtk DataFile Version 3.0
+#This file was generated 
+ASCII
+DATASET UNSTRUCTURED_GRID
+
+POINTS 8 double
+0.00000 0.00000 0.00000
+0.00000 0.00000 0.00000
+1.00000 0.00000 0.00000
+1.00000 0.00000 0.00000
+0.00000 1.00000 1.00000
+0.00000 1.00000 1.00000
+1.00000 1.00000 1.00000
+1.00000 1.00000 1.00000
+
+CELLS 8 16
+1      0
+1      1
+1      2
+1      3
+1      4
+1      5
+1      6
+1      7
+
+CELL_TYPES 8
+ 1 1 1 1 1 1 1 1
+POINT_DATA 8
+VECTORS property_coord_0__property_coord_1__property_coord_2 double
+0.00000 0.00000 0.00000
+0.00000 0.00000 0.00000
+1.00000 0.00000 0.00000
+1.00000 0.00000 0.00000
+0.00000 1.00000 1.00000
+0.00000 1.00000 1.00000
+1.00000 1.00000 1.00000
+1.00000 1.00000 1.00000
+SCALARS id double 1
+LOOKUP_TABLE default
+0.00000 1.00000 4.00000 5.00000 2.00000 3.00000 6.00000 7.00000 
+DEAL:3d/3d::OK
diff --git a/tests/particles/data_out_04.cc b/tests/particles/data_out_04.cc
new file mode 100644 (file)
index 0000000..42c3043
--- /dev/null
@@ -0,0 +1,160 @@
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2017 - 2019 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.md at
+// the top level directory of deal.II.
+//
+// ---------------------------------------------------------------------
+
+
+
+// Create a particle handler then generate an output.
+// Tests that tensor particle properties are output correctly.
+
+#include <deal.II/base/data_out_base.h>
+
+#include <deal.II/distributed/tria.h>
+
+#include <deal.II/fe/mapping_q.h>
+
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/grid_tools.h>
+
+#include <deal.II/particles/data_out.h>
+#include <deal.II/particles/particle_handler.h>
+
+#include "../tests.h"
+
+template <int dim, int spacedim>
+void
+create_regular_particle_distribution(
+  Particles::ParticleHandler<dim, spacedim> &particle_handler,
+  const Triangulation<dim, spacedim> &       tr,
+  const unsigned int                         particles_per_direction = 3)
+{
+  for (unsigned int i = 0; i < particles_per_direction; ++i)
+    for (unsigned int j = 0; j < particles_per_direction; ++j)
+      {
+        Point<spacedim> position;
+        Point<dim>      reference_position;
+        unsigned int    id = i * particles_per_direction + j;
+
+        position[0] = static_cast<double>(i) /
+                      static_cast<double>(particles_per_direction - 1);
+        position[1] = static_cast<double>(j) /
+                      static_cast<double>(particles_per_direction - 1);
+
+        if (dim > 2)
+          for (unsigned int k = 0; k < particles_per_direction; ++k)
+            {
+              position[2] = static_cast<double>(j) /
+                            static_cast<double>(particles_per_direction - 1);
+              id = i * particles_per_direction * particles_per_direction +
+                   j * particles_per_direction + k;
+              Particles::Particle<dim, spacedim> particle(position,
+                                                          reference_position,
+                                                          id);
+
+              typename Triangulation<dim, spacedim>::active_cell_iterator cell =
+                GridTools::find_active_cell_around_point(
+                  tr, particle.get_location());
+
+              Particles::ParticleIterator<dim, spacedim> pit =
+                particle_handler.insert_particle(particle, cell);
+
+              for (unsigned int i = 0; i < spacedim; ++i)
+                pit->get_properties()[i] = pit->get_location()[i];
+            }
+        else
+          {
+            Particles::Particle<dim, spacedim> particle(position,
+                                                        reference_position,
+                                                        id);
+
+            typename Triangulation<dim, spacedim>::active_cell_iterator cell =
+              GridTools::find_active_cell_around_point(tr,
+                                                       particle.get_location());
+
+            Particles::ParticleIterator<dim, spacedim> pit =
+              particle_handler.insert_particle(particle, cell);
+
+            for (unsigned int i = 0; i < spacedim * spacedim; ++i)
+              pit->get_properties()[i] = pit->get_location()[i % spacedim];
+          }
+      }
+}
+
+template <int dim, int spacedim>
+void
+test()
+{
+  {
+    Triangulation<dim, spacedim> tr;
+
+    GridGenerator::hyper_cube(tr);
+    tr.refine_global(2);
+    MappingQ<dim, spacedim> mapping(1);
+
+    // Create a particle handler using two manually created particles
+    Particles::ParticleHandler<dim, spacedim> particle_handler(
+      tr, mapping, spacedim * spacedim);
+
+    create_regular_particle_distribution(particle_handler, tr, 2);
+
+    std::vector<std::string> data_names;
+    std::vector<DataComponentInterpretation::DataComponentInterpretation>
+      data_interpretations;
+
+    for (unsigned int i = 0; i < spacedim * spacedim; ++i)
+      {
+        data_names.emplace_back("property_coord");
+        data_interpretations.emplace_back(
+          DataComponentInterpretation::component_is_part_of_tensor);
+      }
+
+    Particles::DataOut<dim, spacedim> particle_output;
+    particle_output.build_patches(particle_handler,
+                                  data_names,
+                                  data_interpretations);
+
+    // Write gnuplot to check the component output
+    particle_output.write_gnuplot(deallog.get_file_stream());
+
+    // Write to vtu to check that vtu can output tensors
+    std::stringstream test_stream;
+    particle_output.write_vtu(test_stream);
+
+    // Write pvtu to check the component interpretation
+    // We can not write vtk like in data_out_03, because vtk
+    // only supports vectors, but not tensors (at the moment of writing this
+    // test)
+    particle_output.write_pvtu_record(deallog.get_file_stream(),
+                                      {"solution-01"});
+  }
+
+  deallog << "OK" << std::endl;
+}
+
+
+
+int
+main(int argc, char *argv[])
+{
+  initlog();
+  deallog.push("2d/2d");
+  test<2, 2>();
+  deallog.pop();
+  deallog.push("2d/3d");
+  test<2, 3>();
+  deallog.pop();
+  deallog.push("3d/3d");
+  test<3, 3>();
+  deallog.pop();
+}
diff --git a/tests/particles/data_out_04.output b/tests/particles/data_out_04.output
new file mode 100644 (file)
index 0000000..8a69758
--- /dev/null
@@ -0,0 +1,105 @@
+
+# This file was generated by the deal.II library.
+
+
+#
+# For a description of the GNUPLOT format see the GNUPLOT manual.
+#
+# <x> <y> <id> <property_coord> <property_coord> <property_coord> <property_coord> 
+0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 
+
+1.00000 0.00000 2.00000 1.00000 0.00000 1.00000 0.00000 
+
+0.00000 1.00000 1.00000 0.00000 1.00000 0.00000 1.00000 
+
+1.00000 1.00000 3.00000 1.00000 1.00000 1.00000 1.00000 
+
+<?xml version="1.0"?>
+<!--
+#This file was generated 
+-->
+<VTKFile type="PUnstructuredGrid" version="0.1" byte_order="LittleEndian">
+  <PUnstructuredGrid GhostLevel="0">
+    <PPointData Scalars="scalars">
+    <PDataArray type="Float32" Name="property_coord" NumberOfComponents="9" format="ascii"/>
+    <PDataArray type="Float32" Name="id" format="ascii"/>
+    </PPointData>
+    <PPoints>
+      <PDataArray type="Float32" NumberOfComponents="3"/>
+    </PPoints>
+    <Piece Source="solution-01"/>
+  </PUnstructuredGrid>
+</VTKFile>
+DEAL:2d/2d::OK
+# This file was generated by the deal.II library.
+
+
+#
+# For a description of the GNUPLOT format see the GNUPLOT manual.
+#
+# <x> <y> <z> <id> <property_coord> <property_coord> <property_coord> <property_coord> <property_coord> <property_coord> <property_coord> <property_coord> <property_coord> 
+0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 
+
+1.00000 0.00000 0.00000 2.00000 1.00000 0.00000 0.00000 1.00000 0.00000 0.00000 1.00000 0.00000 0.00000 
+
+0.00000 1.00000 0.00000 1.00000 0.00000 1.00000 0.00000 0.00000 1.00000 0.00000 0.00000 1.00000 0.00000 
+
+1.00000 1.00000 0.00000 3.00000 1.00000 1.00000 0.00000 1.00000 1.00000 0.00000 1.00000 1.00000 0.00000 
+
+<?xml version="1.0"?>
+<!--
+#This file was generated 
+-->
+<VTKFile type="PUnstructuredGrid" version="0.1" byte_order="LittleEndian">
+  <PUnstructuredGrid GhostLevel="0">
+    <PPointData Scalars="scalars">
+    <PDataArray type="Float32" Name="property_coord" NumberOfComponents="9" format="ascii"/>
+    <PDataArray type="Float32" Name="id" format="ascii"/>
+    </PPointData>
+    <PPoints>
+      <PDataArray type="Float32" NumberOfComponents="3"/>
+    </PPoints>
+    <Piece Source="solution-01"/>
+  </PUnstructuredGrid>
+</VTKFile>
+DEAL:2d/3d::OK
+# This file was generated by the deal.II library.
+
+
+#
+# For a description of the GNUPLOT format see the GNUPLOT manual.
+#
+# <x> <y> <z> <id> <property_coord> <property_coord> <property_coord> <property_coord> <property_coord> <property_coord> <property_coord> <property_coord> <property_coord> 
+0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 
+
+0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 
+
+1.00000 0.00000 0.00000 4.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 
+
+1.00000 0.00000 0.00000 5.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 
+
+0.00000 1.00000 1.00000 2.00000 0.00000 1.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 
+
+0.00000 1.00000 1.00000 3.00000 0.00000 1.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 
+
+1.00000 1.00000 1.00000 6.00000 1.00000 1.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 
+
+1.00000 1.00000 1.00000 7.00000 1.00000 1.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 
+
+<?xml version="1.0"?>
+<!--
+#This file was generated 
+-->
+<VTKFile type="PUnstructuredGrid" version="0.1" byte_order="LittleEndian">
+  <PUnstructuredGrid GhostLevel="0">
+    <PPointData Scalars="scalars">
+    <PDataArray type="Float32" Name="property_coord" NumberOfComponents="9" format="ascii"/>
+    <PDataArray type="Float32" Name="id" format="ascii"/>
+    </PPointData>
+    <PPoints>
+      <PDataArray type="Float32" NumberOfComponents="3"/>
+    </PPoints>
+    <Piece Source="solution-01"/>
+  </PUnstructuredGrid>
+</VTKFile>
+DEAL:3d/3d::OK

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.