From: heister <heister@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Tue, 4 Mar 2014 14:10:02 +0000 (+0000)
Subject: add test
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=aa9e1a6f17565c3f8a29b01dba9649f433378471;p=dealii-svn.git

add test

git-svn-id: https://svn.dealii.org/trunk@32606 0785d39b-7218-0410-832d-ea1e28bc413d
---

diff --git a/tests/mpi/codim_01.cc b/tests/mpi/codim_01.cc
new file mode 100644
index 0000000000..3c77850ef2
--- /dev/null
+++ b/tests/mpi/codim_01.cc
@@ -0,0 +1,117 @@
+// ---------------------------------------------------------------------
+// $Id: 2d_refinement_01.cc 31349 2013-10-20 19:07:06Z maier $
+//
+// Copyright (C) 2008 - 2013 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.
+//
+// ---------------------------------------------------------------------
+
+
+
+// Test Tria<2,3> and DataOutput.
+
+#include "../tests.h"
+#include "coarse_grid_common.h"
+#include <deal.II/base/logstream.h>
+#include <deal.II/base/tensor.h>
+#include <deal.II/grid/tria.h>
+#include <deal.II/distributed/tria.h>
+#include <deal.II/grid/tria_accessor.h>
+#include <deal.II/grid/tria_iterator.h>
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/grid_out.h>
+#include <deal.II/numerics/data_out.h>
+
+#include <fstream>
+#include <deal.II/fe/fe_q.h>
+
+template <int dim, int spacedim=dim>
+void write_vtk (const parallel::distributed::Triangulation<dim,spacedim> &tria,
+                const char                                *filename)
+{
+  deallog << "Checksum: "
+          << tria.get_checksum ()
+          << std::endl;
+
+  tria.write_mesh_vtk (filename);
+
+  // copy the .pvtu and .vtu files
+  // into the logstream
+  {
+    std::ifstream in((std::string(filename) + ".pvtu").c_str());
+    while (in)
+      {
+        std::string s;
+        std::getline(in, s);
+        deallog.get_file_stream() << s << "\n";
+      }
+  }
+
+  {
+    std::ifstream in((std::string(filename) + "_0000.vtu").c_str());
+    while (in)
+      {
+        std::string s;
+        std::getline(in, s);
+        deallog.get_file_stream() << s << "\n";
+      }
+  }
+}
+
+template<int dim, int spacedim>
+void test(std::ostream & /*out*/)
+{
+  parallel::distributed::Triangulation<dim,spacedim> tr(MPI_COMM_WORLD);
+
+  GridGenerator::torus(tr, 1, 0.2);
+  tr.refine_global();
+  
+  tr.begin_active()->set_refine_flag();
+  (++(tr.begin_active()))->set_refine_flag();
+
+  tr.execute_coarsening_and_refinement ();
+
+  write_vtk(tr, "file");
+
+  FE_Q<dim,spacedim> fe(1);
+  DoFHandler<dim,spacedim> dh(tr);
+  dh.distribute_dofs(fe);
+  deallog << "dofs " << dh.n_dofs() << std::endl;
+
+  DataOut<dim,DoFHandler<dim,spacedim> > data_out;
+    data_out.attach_triangulation (tr);
+    Vector<float> subdomain (tr.n_active_cells());
+   for (unsigned int i=0; i<subdomain.size(); ++i)
+      subdomain(i) = tr.locally_owned_subdomain();
+   data_out.add_data_vector (subdomain, "subdomain");
+
+   std::string name = "f0.vtu";
+   name[1] += tr.locally_owned_subdomain();
+
+   {
+     std::ofstream file(name.c_str());
+     data_out.build_patches ();
+     data_out.write_vtu (file);
+   }
+   cat_file(name.c_str());
+}
+
+
+int main(int argc, char *argv[])
+{
+  Utilities::MPI::MPI_InitFinalize mpi_initialization(argc, argv, 1);
+
+  MPILogInitAll init;
+
+  deallog.push("2-3");
+  test<2,3>(deallog.get_file_stream());  
+  deallog.pop();
+}
diff --git a/tests/mpi/codim_01.mpirun=3.debug.output b/tests/mpi/codim_01.mpirun=3.debug.output
new file mode 100644
index 0000000000..c93a2f0615
--- /dev/null
+++ b/tests/mpi/codim_01.mpirun=3.debug.output
@@ -0,0 +1,597 @@
+
+DEAL:0:2-3::Checksum: 296355976
+<?xml version="1.0"?>
+<VTKFile type="PUnstructuredGrid" version="0.1" byte_order="LittleEndian">
+  <PUnstructuredGrid GhostLevel="0">
+    <PPoints>
+      <PDataArray type="Float32" Name="Position" NumberOfComponents="3" format="ascii"/>
+    </PPoints>
+    <PCellData Scalars="mpirank,treeid">
+      <PDataArray type="Int32" Name="mpirank" format="ascii"/>
+      <PDataArray type="Int32" Name="treeid" format="ascii"/>
+    </PCellData>
+    <PPointData>
+    </PPointData>
+    <Piece Source="file_0000.vtu"/>
+    <Piece Source="file_0001.vtu"/>
+    <Piece Source="file_0002.vtu"/>
+  </PUnstructuredGrid>
+</VTKFile>
+
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian">
+  <UnstructuredGrid>
+    <Piece NumberOfPoints="92" NumberOfCells="23">
+      <Points>
+        <DataArray type="Float32" Name="Position" NumberOfComponents="3" format="ascii">
+            7.96242177e-01   1.24999997e-03   5.00781229e-03
+            6.05945289e-01   1.24999997e-03   1.95304692e-01
+            8.43445301e-01   4.87499982e-02   5.30468766e-03
+            6.41867161e-01   4.87499982e-02   2.06882820e-01
+            5.95929682e-01   1.24999997e-03   2.05320314e-01
+            4.05632824e-01   1.24999997e-03   3.95617187e-01
+            6.31257832e-01   4.87499982e-02   2.17492193e-01
+            4.29679692e-01   4.87499982e-02   4.19070303e-01
+            8.45929682e-01   5.12499996e-02   5.32031246e-03
+            6.43757820e-01   5.12499996e-02   2.07492188e-01
+            8.93132806e-01   9.87500027e-02   5.61718736e-03
+            6.79679692e-01   9.87500027e-02   2.19070315e-01
+            6.33117199e-01   5.12499996e-02   2.18132809e-01
+            4.30945307e-01   5.12499996e-02   4.20304686e-01
+            6.68445289e-01   9.87500027e-02   2.30304688e-01
+            4.54992175e-01   9.87500027e-02   4.43757802e-01
+            3.91218752e-01   2.49999994e-03   4.11281258e-01
+            1.00312503e-02   2.49999994e-03   7.92468727e-01
+            4.37531263e-01   9.74999964e-02   4.59968746e-01
+            1.12187499e-02   9.74999964e-02   8.86281252e-01
+            8.91218722e-01   1.02499999e-01   1.12812500e-02
+            4.62531239e-01   1.02499999e-01   4.39968765e-01
+            9.85031247e-01   1.97500005e-01   1.24687497e-02
+            5.11218727e-01   1.97500005e-01   4.86281246e-01
+            4.39968765e-01   1.02499999e-01   4.62531239e-01
+            1.12812500e-02   1.02499999e-01   8.91218722e-01
+            4.86281246e-01   1.97500005e-01   5.11218727e-01
+            1.24687497e-02   1.97500005e-01   9.85031247e-01
+            9.85031247e-01  -1.97500005e-01   1.24687497e-02
+            5.11218727e-01  -1.97500005e-01   4.86281246e-01
+            8.91218722e-01  -1.02499999e-01   1.12812500e-02
+            4.62531239e-01  -1.02499999e-01   4.39968765e-01
+            4.86281246e-01  -1.97500005e-01   5.11218727e-01
+            1.24687497e-02  -1.97500005e-01   9.85031247e-01
+            4.39968765e-01  -1.02499999e-01   4.62531239e-01
+            1.12812500e-02  -1.02499999e-01   8.91218722e-01
+            8.86281252e-01  -9.74999964e-02   1.12187499e-02
+            4.59968746e-01  -9.74999964e-02   4.37531263e-01
+            7.92468727e-01  -2.49999994e-03   1.00312503e-02
+            4.11281258e-01  -2.49999994e-03   3.91218752e-01
+            4.37531263e-01  -9.74999964e-02   4.59968746e-01
+            1.12187499e-02  -9.74999964e-02   8.86281252e-01
+            3.91218752e-01  -2.49999994e-03   4.11281258e-01
+            1.00312503e-02  -2.49999994e-03   7.92468727e-01
+            9.89968777e-01   1.97500005e-01   1.25312498e-02
+            5.13781250e-01   1.97500005e-01   4.88718748e-01
+            1.08378124e+00   1.02499999e-01   1.37187503e-02
+            5.62468767e-01   1.02499999e-01   5.35031259e-01
+            4.88718748e-01   1.97500005e-01   5.13781250e-01
+            1.25312498e-02   1.97500005e-01   9.89968777e-01
+            5.35031259e-01   1.02499999e-01   5.62468767e-01
+            1.37187503e-02   1.02499999e-01   1.08378124e+00
+            1.08871877e+00   9.74999964e-02   1.37812505e-02
+            5.65031230e-01   9.74999964e-02   5.37468731e-01
+            1.18253124e+00   2.49999994e-03   1.49687501e-02
+            6.13718748e-01   2.49999994e-03   5.83781242e-01
+            5.37468731e-01   9.74999964e-02   5.65031230e-01
+            1.37812505e-02   9.74999964e-02   1.08871877e+00
+            5.83781242e-01   2.49999994e-03   6.13718748e-01
+            1.49687501e-02   2.49999994e-03   1.18253124e+00
+           -1.00312503e-02   2.49999994e-03   7.92468727e-01
+           -3.91218752e-01   2.49999994e-03   4.11281258e-01
+           -1.12187499e-02   9.74999964e-02   8.86281252e-01
+           -4.37531263e-01   9.74999964e-02   4.59968746e-01
+           -4.11281258e-01   2.49999994e-03   3.91218752e-01
+           -7.92468727e-01   2.49999994e-03   1.00312503e-02
+           -4.59968746e-01   9.74999964e-02   4.37531263e-01
+           -8.86281252e-01   9.74999964e-02   1.12187499e-02
+           -1.12812500e-02   1.02499999e-01   8.91218722e-01
+           -4.39968765e-01   1.02499999e-01   4.62531239e-01
+           -1.24687497e-02   1.97500005e-01   9.85031247e-01
+           -4.86281246e-01   1.97500005e-01   5.11218727e-01
+           -4.62531239e-01   1.02499999e-01   4.39968765e-01
+           -8.91218722e-01   1.02499999e-01   1.12812500e-02
+           -5.11218727e-01   1.97500005e-01   4.86281246e-01
+           -9.85031247e-01   1.97500005e-01   1.24687497e-02
+            1.00312503e-02   2.49999994e-03  -7.92468727e-01
+            3.91218752e-01   2.49999994e-03  -4.11281258e-01
+            1.12187499e-02   9.74999964e-02  -8.86281252e-01
+            4.37531263e-01   9.74999964e-02  -4.59968746e-01
+            4.11281258e-01   2.49999994e-03  -3.91218752e-01
+            7.92468727e-01   2.49999994e-03  -1.00312503e-02
+            4.59968746e-01   9.74999964e-02  -4.37531263e-01
+            8.86281252e-01   9.74999964e-02  -1.12187499e-02
+            1.12812500e-02   1.02499999e-01  -8.91218722e-01
+            4.39968765e-01   1.02499999e-01  -4.62531239e-01
+            1.24687497e-02   1.97500005e-01  -9.85031247e-01
+            4.86281246e-01   1.97500005e-01  -5.11218727e-01
+            4.62531239e-01   1.02499999e-01  -4.39968765e-01
+            8.91218722e-01   1.02499999e-01  -1.12812500e-02
+            5.11218727e-01   1.97500005e-01  -4.86281246e-01
+            9.85031247e-01   1.97500005e-01  -1.24687497e-02
+        </DataArray>
+      </Points>
+      <Cells>
+        <DataArray type="Int32" Name="connectivity" format="ascii">
+          0 1 2 3
+          4 5 6 7
+          8 9 10 11
+          12 13 14 15
+          16 17 18 19
+          20 21 22 23
+          24 25 26 27
+          28 29 30 31
+          32 33 34 35
+          36 37 38 39
+          40 41 42 43
+          44 45 46 47
+          48 49 50 51
+          52 53 54 55
+          56 57 58 59
+          60 61 62 63
+          64 65 66 67
+          68 69 70 71
+          72 73 74 75
+          76 77 78 79
+          80 81 82 83
+          84 85 86 87
+          88 89 90 91
+        </DataArray>
+        <DataArray type="Int32" Name="offsets" format="ascii">
+          4 8 12 16 20 24 28 32
+          36 40 44 48 52 56 60 64
+          68 72 76 80 84 88 92
+        </DataArray>
+        <DataArray type="UInt8" Name="types" format="ascii">
+          8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8
+          8 8 8
+        </DataArray>
+      </Cells>
+      <CellData Scalars="mpirank,treeid">
+        <DataArray type="Int32" Name="mpirank" format="ascii">
+          0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0
+        </DataArray>
+        <DataArray type="Int32" Name="treeid" format="ascii">
+          0 0 0 0 0 0 0 1 1 1 1 2 2 2 2 3 3 3 3 4
+          4 4 4
+        </DataArray>
+      </CellData>
+      <PointData>
+      </PointData>
+    </Piece>
+  </UnstructuredGrid>
+</VTKFile>
+
+DEAL:0:2-3::dofs 69
+<?xml version="1.0" ?> 
+<!-- 
+# vtk DataFile Version 3.0
+#This file was generated 
+-->
+<VTKFile type="UnstructuredGrid" version="0.1" compressor="vtkZLibDataCompressor" byte_order="LittleEndian">
+<UnstructuredGrid>
+<Piece NumberOfPoints="92" NumberOfCells="23" >
+  <Points>
+    <DataArray type="Float64" NumberOfComponents="3" format="binary">
+AQAAAKAIAACgCAAAJAEAAA==eNqdlTEOwjAMRbkT7LkNnInegCMw/bmiHAC1QxEwIBhAvQAqSYzi9AdMFuRP+H5xHFOtx3Vws7AqFaf61e2bcXXOx1uJ031Rv8tn1OO+ss8j6LXi6JXem3mYT7q/Fg7FA8KDlEfOC3JepJy5j+JB5Cj7yHlB6v+VJ97z9P3nfWLlqYz9Zu0HHz+zfvPxUfQYW/uh7PPpt5Tjd575ew1ZPbx+cjq28jCf6XwDvZdwbyj3h54DHfQ7ZT4hD3R+5hPyQNfZylOeAz3I/DHy5D4hD6b7PKs/yDsCOS/InKQ+bA4wHzbnrTyk/iDvFP/x5D5szrA6LNS7ib/z+i34bOT71XJcZ9H9vp1jPmxOln1apbdmHlY3H19UnRszD/Nh/+MvbR2P5g==
+    </DataArray>
+  </Points>
+
+  <Cells>
+    <DataArray type="Int32" Name="connectivity" format="binary">
+AQAAAHABAABwAQAAlwAAAA==eNoNw4c2QgEAANBXKoRESCUjq4XsmZQZiczi//+je8+5QRAEIUcMGzHqqDHHHHfCuJNOOW3CpDOmnHXOedMuuGjGnFmXzLvisquuuW7BDTfddsuiJSuWrbrjnrvW3PfQA4889tQTzzz30guvrNvw2hub3tryznsfffDJti8+2/HVN7u+2/PTD7/89tcf+w78988h7c0QWw==
+    </DataArray>
+    <DataArray type="Int32" Name="offsets" format="binary">
+AQAAAFwAAABcAAAAQAAAAA==eNoNw5ESgEAUAMAHB8FBEATBQRAEQRAEQRAcBEHQ/39HuzObIqIx29rZO1gcnZxdXN3cPTy9rN4+vn7+kBwEUQ==
+    </DataArray>
+    <DataArray type="UInt8" Name="types" format="binary">
+AQAAABcAAAAXAAAACwAAAA==eNrj5MQKAAnLANA=
+    </DataArray>
+  </Cells>
+  <PointData Scalars="scalars">
+    <DataArray type="Float64" Name="subdomain" format="binary">
+AQAAAOACAADgAgAADwAAAA==eNpjYBgFo2D4AgAC4AAB    </DataArray>
+  </PointData>
+ </Piece>
+ </UnstructuredGrid>
+</VTKFile>
+
+
+DEAL:1:2-3::Checksum: 0
+<?xml version="1.0"?>
+<VTKFile type="PUnstructuredGrid" version="0.1" byte_order="LittleEndian">
+  <PUnstructuredGrid GhostLevel="0">
+    <PPoints>
+      <PDataArray type="Float32" Name="Position" NumberOfComponents="3" format="ascii"/>
+    </PPoints>
+    <PCellData Scalars="mpirank,treeid">
+      <PDataArray type="Int32" Name="mpirank" format="ascii"/>
+      <PDataArray type="Int32" Name="treeid" format="ascii"/>
+    </PCellData>
+    <PPointData>
+
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian">
+  <UnstructuredGrid>
+    <Piece NumberOfPoints="92" NumberOfCells="23">
+      <Points>
+        <DataArray type="Float32" Name="Position" NumberOfComponents="3" format="ascii">
+            7.96242177e-01   1.24999997e-03   5.00781229e-03
+            6.05945289e-01   1.24999997e-03   1.95304692e-01
+            8.43445301e-01   4.87499982e-02   5.30468766e-03
+            6.41867161e-01   4.87499982e-02   2.06882820e-01
+            5.95929682e-01   1.24999997e-03   2.05320314e-01
+            4.05632824e-01   1.24999997e-03   3.95617187e-01
+            6.31257832e-01   4.87499982e-02   2.17492193e-01
+            4.29679692e-01   4.87499982e-02   4.19070303e-01
+            8.45929682e-01   5.12499996e-02   5.32031246e-03
+            6.43757820e-01   5.12499996e-02   2.07492188e-01
+            8.93132806e-01   9.87500027e-02   5.61718736e-03
+            6.79679692e-01   9.87500027e-02   2.19070315e-01
+            6.33117199e-01   5.12499996e-02   2.18132809e-01
+            4.30945307e-01   5.12499996e-02   4.20304686e-01
+            6.68445289e-01   9.87500027e-02   2.30304688e-01
+            4.54992175e-01   9.87500027e-02   4.43757802e-01
+            3.91218752e-01   2.49999994e-03   4.11281258e-01
+            1.00312503e-02   2.49999994e-03   7.92468727e-01
+            4.37531263e-01   9.74999964e-02   4.59968746e-01
+            1.12187499e-02   9.74999964e-02   8.86281252e-01
+            8.91218722e-01   1.02499999e-01   1.12812500e-02
+            4.62531239e-01   1.02499999e-01   4.39968765e-01
+            9.85031247e-01   1.97500005e-01   1.24687497e-02
+            5.11218727e-01   1.97500005e-01   4.86281246e-01
+            4.39968765e-01   1.02499999e-01   4.62531239e-01
+            1.12812500e-02   1.02499999e-01   8.91218722e-01
+            4.86281246e-01   1.97500005e-01   5.11218727e-01
+            1.24687497e-02   1.97500005e-01   9.85031247e-01
+            9.85031247e-01  -1.97500005e-01   1.24687497e-02
+            5.11218727e-01  -1.97500005e-01   4.86281246e-01
+            8.91218722e-01  -1.02499999e-01   1.12812500e-02
+            4.62531239e-01  -1.02499999e-01   4.39968765e-01
+            4.86281246e-01  -1.97500005e-01   5.11218727e-01
+            1.24687497e-02  -1.97500005e-01   9.85031247e-01
+            4.39968765e-01  -1.02499999e-01   4.62531239e-01
+            1.12812500e-02  -1.02499999e-01   8.91218722e-01
+            8.86281252e-01  -9.74999964e-02   1.12187499e-02
+            4.59968746e-01  -9.74999964e-02   4.37531263e-01
+            7.92468727e-01  -2.49999994e-03   1.00312503e-02
+            4.11281258e-01  -2.49999994e-03   3.91218752e-01
+            4.37531263e-01  -9.74999964e-02   4.59968746e-01
+            1.12187499e-02  -9.74999964e-02   8.86281252e-01
+            3.91218752e-01  -2.49999994e-03   4.11281258e-01
+            1.00312503e-02  -2.49999994e-03   7.92468727e-01
+            9.89968777e-01   1.97500005e-01   1.25312498e-02
+            5.13781250e-01   1.97500005e-01   4.88718748e-01
+            1.08378124e+00   1.02499999e-01   1.37187503e-02
+            5.62468767e-01   1.02499999e-01   5.35031259e-01
+            4.88718748e-01   1.97500005e-01   5.13781250e-01
+            1.25312498e-02   1.97500005e-01   9.89968777e-01
+            5.35031259e-01   1.02499999e-01   5.62468767e-01
+            1.37187503e-02   1.02499999e-01   1.08378124e+00
+            1.08871877e+00   9.74999964e-02   1.37812505e-02
+            5.65031230e-01   9.74999964e-02   5.37468731e-01
+            1.18253124e+00   2.49999994e-03   1.49687501e-02
+            6.13718748e-01   2.49999994e-03   5.83781242e-01
+            5.37468731e-01   9.74999964e-02   5.65031230e-01
+            1.37812505e-02   9.74999964e-02   1.08871877e+00
+            5.83781242e-01   2.49999994e-03   6.13718748e-01
+            1.49687501e-02   2.49999994e-03   1.18253124e+00
+           -1.00312503e-02   2.49999994e-03   7.92468727e-01
+           -3.91218752e-01   2.49999994e-03   4.11281258e-01
+           -1.12187499e-02   9.74999964e-02   8.86281252e-01
+           -4.37531263e-01   9.74999964e-02   4.59968746e-01
+           -4.11281258e-01   2.49999994e-03   3.91218752e-01
+           -7.92468727e-01   2.49999994e-03   1.00312503e-02
+           -4.59968746e-01   9.74999964e-02   4.37531263e-01
+           -8.86281252e-01   9.74999964e-02   1.12187499e-02
+           -1.12812500e-02   1.02499999e-01   8.91218722e-01
+           -4.39968765e-01   1.02499999e-01   4.62531239e-01
+           -1.24687497e-02   1.97500005e-01   9.85031247e-01
+           -4.86281246e-01   1.97500005e-01   5.11218727e-01
+           -4.62531239e-01   1.02499999e-01   4.39968765e-01
+           -8.91218722e-01   1.02499999e-01   1.12812500e-02
+           -5.11218727e-01   1.97500005e-01   4.86281246e-01
+           -9.85031247e-01   1.97500005e-01   1.24687497e-02
+            1.00312503e-02   2.49999994e-03  -7.92468727e-01
+            3.91218752e-01   2.49999994e-03  -4.11281258e-01
+            1.12187499e-02   9.74999964e-02  -8.86281252e-01
+            4.37531263e-01   9.74999964e-02  -4.59968746e-01
+            4.11281258e-01   2.49999994e-03  -3.91218752e-01
+            7.92468727e-01   2.49999994e-03  -1.00312503e-02
+            4.59968746e-01   9.74999964e-02  -4.37531263e-01
+            8.86281252e-01   9.74999964e-02  -1.12187499e-02
+            1.12812500e-02   1.02499999e-01  -8.91218722e-01
+            4.39968765e-01   1.02499999e-01  -4.62531239e-01
+            1.24687497e-02   1.97500005e-01  -9.85031247e-01
+            4.86281246e-01   1.97500005e-01  -5.11218727e-01
+            4.62531239e-01   1.02499999e-01  -4.39968765e-01
+            8.91218722e-01   1.02499999e-01  -1.12812500e-02
+            5.11218727e-01   1.97500005e-01  -4.86281246e-01
+            9.85031247e-01   1.97500005e-01  -1.24687497e-02
+        </DataArray>
+      </Points>
+      <Cells>
+        <DataArray type="Int32" Name="connectivity" format="ascii">
+          0 1 2 3
+          4 5 6 7
+          8 9 10 11
+          12 13 14 15
+          16 17 18 19
+          20 21 22 23
+          24 25 26 27
+          28 29 30 31
+          32 33 34 35
+          36 37 38 39
+          40 41 42 43
+          44 45 46 47
+          48 49 50 51
+          52 53 54 55
+          56 57 58 59
+          60 61 62 63
+          64 65 66 67
+          68 69 70 71
+          72 73 74 75
+          76 77 78 79
+          80 81 82 83
+          84 85 86 87
+          88 89 90 91
+        </DataArray>
+        <DataArray type="Int32" Name="offsets" format="ascii">
+          4 8 12 16 20 24 28 32
+          36 40 44 48 52 56 60 64
+          68 72 76 80 84 88 92
+        </DataArray>
+        <DataArray type="UInt8" Name="types" format="ascii">
+          8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8
+          8 8 8
+        </DataArray>
+      </Cells>
+      <CellData Scalars="mpirank,treeid">
+        <DataArray type="Int32" Name="mpirank" format="ascii">
+          0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0
+        </DataArray>
+        <DataArray type="Int32" Name="treeid" format="ascii">
+          0 0 0 0 0 0 0 1 1 1 1 2 2 2 2 3 3 3 3 4
+          4 4 4
+        </DataArray>
+      </CellData>
+      <PointData>
+      </PointData>
+    </Piece>
+  </UnstructuredGrid>
+</VTKFile>
+
+DEAL:1:2-3::dofs 69
+<?xml version="1.0" ?> 
+<!-- 
+# vtk DataFile Version 3.0
+#This file was generated 
+-->
+<VTKFile type="UnstructuredGrid" version="0.1" compressor="vtkZLibDataCompressor" byte_order="LittleEndian">
+<UnstructuredGrid>
+<Piece NumberOfPoints="80" NumberOfCells="20" >
+  <Points>
+    <DataArray type="Float64" NumberOfComponents="3" format="binary">
+AQAAAIAHAACABwAABQEAAA==eNqdVMENgzAM7E4M4DnbTXj5XakdAIkHqC+kPoJYoKLYkXLhmqT+RLnA+WSf3XV7bHKB+MLdS/B+u+6x2tmrvx/3OeJ+Zzzn+TZhPGkex9dmPcdrMPyuqYoJ8KlZD+NJv3c8FPBJf/P0+nzs8ZbjHDXFxwJPUKwD47E8inX+Tw/+t1i9BujfIK16GI/lyXRD/YXUX0j9wX+zpvhc4AmAR99mPJYn4qT+lXryubO51PN5r9fDeCwP8c1C8KyPCn0X4lvwzUD9gToYD9bZdbTqYX0nPmnWw3iIz5XsHyV7TMm8x7mEOS3w5HuY8fheqdw/BT15331PnPutXg/j8T2H+T+XphyK
+    </DataArray>
+  </Points>
+
+  <Cells>
+    <DataArray type="Int32" Name="connectivity" format="binary">
+AQAAAEABAABAAQAAhQAAAA==eNoNw4c2QgEAANBHpSHRUAhF00pWQyVNlPj/v+nec24QBMGOIXcNGzHqnjHj7psw6YGHpjwybdaMOY8tmPfEU4ueee6FJS8te2XFa6vWbFi36Y133nrvg4+2bPvki8+++mbXjj37Dnx36MixH346cebUuQu/XPrtj2tX/rrx3z+3NvsMWQ==
+    </DataArray>
+    <DataArray type="Int32" Name="offsets" format="binary">
+AQAAAFAAAABQAAAAOwAAAA==eNpjYWBg4ABiHiAWAGIRIJYAYhkgVgBiFSDWAGIdIDYAYhMgtgBiGyB2AGIXIPYAYh8gDgBiAGCQA0k=
+    </DataArray>
+    <DataArray type="UInt8" Name="types" format="binary">
+AQAAABQAAAAUAAAACwAAAA==eNrj5MQEAAd2ALU=
+    </DataArray>
+  </Cells>
+  <PointData Scalars="scalars">
+    <DataArray type="Float64" Name="subdomain" format="binary">
+AQAAAIACAACAAgAAFAAAAA==eNpjYACBD/YMo/QoPQA0AJRLXrE=    </DataArray>
+  </PointData>
+ </Piece>
+ </UnstructuredGrid>
+</VTKFile>
+
+
+
+DEAL:2:2-3::Checksum: 0
+
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian">
+  <UnstructuredGrid>
+    <Piece NumberOfPoints="92" NumberOfCells="23">
+      <Points>
+        <DataArray type="Float32" Name="Position" NumberOfComponents="3" format="ascii">
+            7.96242177e-01   1.24999997e-03   5.00781229e-03
+            6.05945289e-01   1.24999997e-03   1.95304692e-01
+            8.43445301e-01   4.87499982e-02   5.30468766e-03
+            6.41867161e-01   4.87499982e-02   2.06882820e-01
+            5.95929682e-01   1.24999997e-03   2.05320314e-01
+            4.05632824e-01   1.24999997e-03   3.95617187e-01
+            6.31257832e-01   4.87499982e-02   2.17492193e-01
+            4.29679692e-01   4.87499982e-02   4.19070303e-01
+            8.45929682e-01   5.12499996e-02   5.32031246e-03
+            6.43757820e-01   5.12499996e-02   2.07492188e-01
+            8.93132806e-01   9.87500027e-02   5.61718736e-03
+            6.79679692e-01   9.87500027e-02   2.19070315e-01
+            6.33117199e-01   5.12499996e-02   2.18132809e-01
+            4.30945307e-01   5.12499996e-02   4.20304686e-01
+            6.68445289e-01   9.87500027e-02   2.30304688e-01
+            4.54992175e-01   9.87500027e-02   4.43757802e-01
+            3.91218752e-01   2.49999994e-03   4.11281258e-01
+            1.00312503e-02   2.49999994e-03   7.92468727e-01
+            4.37531263e-01   9.74999964e-02   4.59968746e-01
+            1.12187499e-02   9.74999964e-02   8.86281252e-01
+            8.91218722e-01   1.02499999e-01   1.12812500e-02
+            4.62531239e-01   1.02499999e-01   4.39968765e-01
+            9.85031247e-01   1.97500005e-01   1.24687497e-02
+            5.11218727e-01   1.97500005e-01   4.86281246e-01
+            4.39968765e-01   1.02499999e-01   4.62531239e-01
+            1.12812500e-02   1.02499999e-01   8.91218722e-01
+            4.86281246e-01   1.97500005e-01   5.11218727e-01
+            1.24687497e-02   1.97500005e-01   9.85031247e-01
+            9.85031247e-01  -1.97500005e-01   1.24687497e-02
+            5.11218727e-01  -1.97500005e-01   4.86281246e-01
+            8.91218722e-01  -1.02499999e-01   1.12812500e-02
+            4.62531239e-01  -1.02499999e-01   4.39968765e-01
+            4.86281246e-01  -1.97500005e-01   5.11218727e-01
+            1.24687497e-02  -1.97500005e-01   9.85031247e-01
+            4.39968765e-01  -1.02499999e-01   4.62531239e-01
+            1.12812500e-02  -1.02499999e-01   8.91218722e-01
+            8.86281252e-01  -9.74999964e-02   1.12187499e-02
+            4.59968746e-01  -9.74999964e-02   4.37531263e-01
+            7.92468727e-01  -2.49999994e-03   1.00312503e-02
+            4.11281258e-01  -2.49999994e-03   3.91218752e-01
+            4.37531263e-01  -9.74999964e-02   4.59968746e-01
+            1.12187499e-02  -9.74999964e-02   8.86281252e-01
+            3.91218752e-01  -2.49999994e-03   4.11281258e-01
+            1.00312503e-02  -2.49999994e-03   7.92468727e-01
+            9.89968777e-01   1.97500005e-01   1.25312498e-02
+            5.13781250e-01   1.97500005e-01   4.88718748e-01
+            1.08378124e+00   1.02499999e-01   1.37187503e-02
+            5.62468767e-01   1.02499999e-01   5.35031259e-01
+            4.88718748e-01   1.97500005e-01   5.13781250e-01
+            1.25312498e-02   1.97500005e-01   9.89968777e-01
+            5.35031259e-01   1.02499999e-01   5.62468767e-01
+            1.37187503e-02   1.02499999e-01   1.08378124e+00
+            1.08871877e+00   9.74999964e-02   1.37812505e-02
+            5.65031230e-01   9.74999964e-02   5.37468731e-01
+            1.18253124e+00   2.49999994e-03   1.49687501e-02
+            6.13718748e-01   2.49999994e-03   5.83781242e-01
+            5.37468731e-01   9.74999964e-02   5.65031230e-01
+            1.37812505e-02   9.74999964e-02   1.08871877e+00
+            5.83781242e-01   2.49999994e-03   6.13718748e-01
+            1.49687501e-02   2.49999994e-03   1.18253124e+00
+           -1.00312503e-02   2.49999994e-03   7.92468727e-01
+           -3.91218752e-01   2.49999994e-03   4.11281258e-01
+           -1.12187499e-02   9.74999964e-02   8.86281252e-01
+           -4.37531263e-01   9.74999964e-02   4.59968746e-01
+           -4.11281258e-01   2.49999994e-03   3.91218752e-01
+           -7.92468727e-01   2.49999994e-03   1.00312503e-02
+           -4.59968746e-01   9.74999964e-02   4.37531263e-01
+           -8.86281252e-01   9.74999964e-02   1.12187499e-02
+           -1.12812500e-02   1.02499999e-01   8.91218722e-01
+           -4.39968765e-01   1.02499999e-01   4.62531239e-01
+           -1.24687497e-02   1.97500005e-01   9.85031247e-01
+           -4.86281246e-01   1.97500005e-01   5.11218727e-01
+           -4.62531239e-01   1.02499999e-01   4.39968765e-01
+           -8.91218722e-01   1.02499999e-01   1.12812500e-02
+           -5.11218727e-01   1.97500005e-01   4.86281246e-01
+           -9.85031247e-01   1.97500005e-01   1.24687497e-02
+            1.00312503e-02   2.49999994e-03  -7.92468727e-01
+            3.91218752e-01   2.49999994e-03  -4.11281258e-01
+            1.12187499e-02   9.74999964e-02  -8.86281252e-01
+            4.37531263e-01   9.74999964e-02  -4.59968746e-01
+            4.11281258e-01   2.49999994e-03  -3.91218752e-01
+            7.92468727e-01   2.49999994e-03  -1.00312503e-02
+            4.59968746e-01   9.74999964e-02  -4.37531263e-01
+            8.86281252e-01   9.74999964e-02  -1.12187499e-02
+            1.12812500e-02   1.02499999e-01  -8.91218722e-01
+            4.39968765e-01   1.02499999e-01  -4.62531239e-01
+            1.24687497e-02   1.97500005e-01  -9.85031247e-01
+            4.86281246e-01   1.97500005e-01  -5.11218727e-01
+            4.62531239e-01   1.02499999e-01  -4.39968765e-01
+            8.91218722e-01   1.02499999e-01  -1.12812500e-02
+            5.11218727e-01   1.97500005e-01  -4.86281246e-01
+            9.85031247e-01   1.97500005e-01  -1.24687497e-02
+        </DataArray>
+      </Points>
+      <Cells>
+        <DataArray type="Int32" Name="connectivity" format="ascii">
+          0 1 2 3
+          4 5 6 7
+          8 9 10 11
+          12 13 14 15
+          16 17 18 19
+          20 21 22 23
+          24 25 26 27
+          28 29 30 31
+          32 33 34 35
+          36 37 38 39
+          40 41 42 43
+          44 45 46 47
+          48 49 50 51
+          52 53 54 55
+          56 57 58 59
+          60 61 62 63
+          64 65 66 67
+          68 69 70 71
+          72 73 74 75
+          76 77 78 79
+          80 81 82 83
+          84 85 86 87
+          88 89 90 91
+        </DataArray>
+        <DataArray type="Int32" Name="offsets" format="ascii">
+          4 8 12 16 20 24 28 32
+          36 40 44 48 52 56 60 64
+          68 72 76 80 84 88 92
+        </DataArray>
+        <DataArray type="UInt8" Name="types" format="ascii">
+          8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8
+          8 8 8
+        </DataArray>
+      </Cells>
+      <CellData Scalars="mpirank,treeid">
+        <DataArray type="Int32" Name="mpirank" format="ascii">
+          0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0
+        </DataArray>
+        <DataArray type="Int32" Name="treeid" format="ascii">
+          0 0 0 0 0 0 0 1 1 1 1 2 2 2 2 3 3 3 3 4
+          4 4 4
+        </DataArray>
+      </CellData>
+      <PointData>
+      </PointData>
+    </Piece>
+  </UnstructuredGrid>
+</VTKFile>
+
+DEAL:2:2-3::dofs 69
+<?xml version="1.0" ?> 
+<!-- 
+# vtk DataFile Version 3.0
+#This file was generated 
+-->
+<VTKFile type="UnstructuredGrid" version="0.1" compressor="vtkZLibDataCompressor" byte_order="LittleEndian">
+<UnstructuredGrid>
+<Piece NumberOfPoints="96" NumberOfCells="24" >
+  <Points>
+    <DataArray type="Float64" NumberOfComponents="3" format="binary">
+AQAAAAAJAAAACQAAGgEAAA==eNqdlTEOwjAMRXsnDuBzwk2YPCPBAZAYWnWqxFDUC6CSOMh2fhs3S5QveP5Jftyu0+P0GwuleWStjyTry3kdV07zh9I8sNYHyMl12JTvECfXKXr9f/t+9O9vmTNnvWet93TMj+fkOkaXgXTH4cd9HW9XN+mvossaccA58Dbnf87iI+onrSdwf0+266gfxKnXmxjlEOSW0b2DnDDi1OstjDgg52E/0bxF/UTzZvwQ8GPeZ9kvgf2S9uk5xo95n4hT9kvg/Hf9HMtbu59o3pBu+73l1HMofdHl33HkOwP6v+Povr7b/6Gf7RxKv3f5D/rxHPnO1N+dz6HhgNzaHBafNv+AM5N9d4hj9gvurdVPe96ifqJ5+wICvYbI
+    </DataArray>
+  </Points>
+
+  <Cells>
+    <DataArray type="Int32" Name="connectivity" format="binary">
+AQAAAIABAACAAQAAnQAAAA==eNoNw4c2QgEAANCnSJJVRsmmYUV2SLISQigz/v8v3HvODYIg6DFsyF777Ddi1AEHjRl3yBGHHXXMpAnHnXDKSVOmzTjtjLPOO+eCiy675IpZ8+YsuOq6a2646ZZFty2564577nvogUeWPfHYUytWPfPcmpdeeOW1dW+8teG9dz7Y9MlHn2356otvtn2344effvvlj13//PUfAicR0Q==
+    </DataArray>
+    <DataArray type="Int32" Name="offsets" format="binary">
+AQAAAGAAAABgAAAAQgAAAA==eNoNw5ESgEAUAMAHQRAEQRAcBEEQBMFBcBAEQRAE/f9XtDuzVUTUNrZ29g4mRydnF1c3s7vFw9PL28fXzx+i4ASx
+    </DataArray>
+    <DataArray type="UInt8" Name="types" format="binary">
+AQAAABgAAAAYAAAACwAAAA==eNrj5MQOAAqkANk=
+    </DataArray>
+  </Cells>
+  <PointData Scalars="scalars">
+    <DataArray type="Float64" Name="subdomain" format="binary">
+AQAAAAADAAAAAwAAEwAAAA==eNpjYAADB4ZRepQegTQAvQ0YAQ==    </DataArray>
+  </PointData>
+ </Piece>
+ </UnstructuredGrid>
+</VTKFile>
+
+