add test for write_vtu_in_parallel

git-svn-id: https://svn.dealii.org/branches/branch_port_the_testsuite@31172 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/tests/mpi/parallel_vtu_01.cc b/tests/mpi/parallel_vtu_01.cc
new file mode 100644 (file)
index 0000000..d137268
--- /dev/null
+++ b/tests/mpi/parallel_vtu_01.cc
@@ -0,0 +1,96 @@
+// ---------------------------------------------------------------------
+// $Id: p4est_data_out_01.cc 30688 2013-09-13 15:16:20Z maier $
+//
+// Copyright (C) 2009 - 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.
+//
+// ---------------------------------------------------------------------
+
+
+
+// create a parallel DoFHandler and output data using the parallel vtk output
+// (uses MPI IO)
+
+#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/grid_generator.h>
+#include <deal.II/grid/grid_out.h>
+#include <deal.II/dofs/dof_handler.h>
+#include <deal.II/dofs/dof_handler.h>
+#include <deal.II/dofs/dof_tools.h>
+#include <deal.II/grid/tria_accessor.h>
+#include <deal.II/grid/tria_iterator.h>
+#include <deal.II/dofs/dof_accessor.h>
+#include <deal.II/numerics/data_out.h>
+
+#include <deal.II/fe/fe_q.h>
+#include <deal.II/lac/trilinos_vector.h>
+
+#include <fstream>
+
+
+template<int dim>
+void test()
+{
+  unsigned int myid = Utilities::MPI::this_mpi_process (MPI_COMM_WORLD);
+
+  if (Utilities::MPI::this_mpi_process (MPI_COMM_WORLD) == 0)
+    deallog << "hyper_cube" << std::endl;
+
+  parallel::distributed::Triangulation<dim> tr(MPI_COMM_WORLD);
+
+  GridGenerator::hyper_cube(tr);
+  DoFHandler<dim> dofh(tr);
+
+  static const FE_Q<dim> fe(2);
+  dofh.distribute_dofs (fe);
+
+
+  TrilinosWrappers::MPI::Vector x;
+  x.reinit(dofh.locally_owned_dofs(), MPI_COMM_WORLD);
+  x=2.0;
+
+  DataOut<dim> data_out;
+  data_out.attach_dof_handler (dofh);
+  data_out.add_data_vector (x, "x");
+  data_out.build_patches ();
+
+  data_out.write_vtu_in_parallel("output.vtu", MPI_COMM_WORLD);
+  MPI_Barrier(MPI_COMM_WORLD);
+  
+  if (myid==0)
+    {
+      std::vector<std::string> filenames;
+      filenames.push_back("output.vtu");
+      {        
+       std::ofstream master("output.pvtu");
+       data_out.write_pvtu_record (master, filenames);
+      }
+      
+      cat_file("output.vtu");
+      cat_file("output.pvtu");      
+    }
+  deallog << "OK" << std::endl;
+}
+
+
+int main(int argc, char *argv[])
+{
+  Utilities::MPI::MPI_InitFinalize mpi_initialization(argc, argv, 1);
+  MPILogInitAll log;
+  unsigned int myid = Utilities::MPI::this_mpi_process (MPI_COMM_WORLD);
+
+  test<2>();
+}

diff --git a/tests/mpi/parallel_vtu_01.with_trilinos=true.mpirun=1.output b/tests/mpi/parallel_vtu_01.with_trilinos=true.mpirun=1.output
new file mode 100644 (file)
index 0000000..81da5d4
--- /dev/null
+++ b/tests/mpi/parallel_vtu_01.with_trilinos=true.mpirun=1.output
@@ -0,0 +1,52 @@
+
+DEAL:0::hyper_cube
+<?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="4" NumberOfCells="1" >
+  <Points>
+    <DataArray type="Float64" NumberOfComponents="3" format="binary">
+AQAAAGAAAABgAAAAFAAAAA==eNpjYMAHPtgzkCQP4yPEAZbkBL0=
+    </DataArray>
+  </Points>
+
+  <Cells>
+    <DataArray type="Int32" Name="connectivity" format="binary">
+AQAAABAAAAAQAAAAEwAAAA==eNpjYGBgYARiZiBmAmIAADwABw==
+    </DataArray>
+    <DataArray type="Int32" Name="offsets" format="binary">
+AQAAAAQAAAAEAAAADAAAAA==eNpjYWBgAAAAFAAF
+    </DataArray>
+    <DataArray type="UInt8" Name="types" format="binary">
+AQAAAAEAAAABAAAACQAAAA==eNrjBAAACgAK
+    </DataArray>
+  </Cells>
+  <PointData Scalars="scalars">
+    <DataArray type="Float64" Name="x" format="binary">
+AQAAACAAAAAgAAAADwAAAA==eNpjYAADBwYcNAANIAEB    </DataArray>
+  </PointData>
+ </Piece>
+ </UnstructuredGrid>
+</VTKFile>
+
+<?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="Float64" Name="x" format="ascii"/>
+    </PPointData>
+    <PPoints>
+      <PDataArray type="Float64" NumberOfComponents="3"/>
+    </PPoints>
+    <Piece Source="output.vtu"/>
+  </PUnstructuredGrid>
+</VTKFile>
+
+DEAL:0::OK

diff --git a/tests/mpi/parallel_vtu_01.with_trilinos=true.mpirun=3.output b/tests/mpi/parallel_vtu_01.with_trilinos=true.mpirun=3.output
new file mode 100644 (file)
index 0000000..7e23298
--- /dev/null
+++ b/tests/mpi/parallel_vtu_01.with_trilinos=true.mpirun=3.output
@@ -0,0 +1,74 @@
+
+DEAL:0::hyper_cube
+<?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="0" NumberOfCells="0" >
+<Cells>
+<DataArray type="UInt8" Name="types"></DataArray>
+</Cells>
+  <PointData Scalars="scalars">
+    <DataArray type="Float64" Name="x"></DataArray>
+  </PointData>
+</Piece>
+<Piece NumberOfPoints="0" NumberOfCells="0" >
+<Cells>
+<DataArray type="UInt8" Name="types"></DataArray>
+</Cells>
+  <PointData Scalars="scalars">
+    <DataArray type="Float64" Name="x"></DataArray>
+  </PointData>
+</Piece>
+<Piece NumberOfPoints="4" NumberOfCells="1" >
+  <Points>
+    <DataArray type="Float64" NumberOfComponents="3" format="binary">
+AQAAAGAAAABgAAAAFAAAAA==eNpjYMAHPtgzkCQP4yPEAZbkBL0=
+    </DataArray>
+  </Points>
+
+  <Cells>
+    <DataArray type="Int32" Name="connectivity" format="binary">
+AQAAABAAAAAQAAAAEwAAAA==eNpjYGBgYARiZiBmAmIAADwABw==
+    </DataArray>
+    <DataArray type="Int32" Name="offsets" format="binary">
+AQAAAAQAAAAEAAAADAAAAA==eNpjYWBgAAAAFAAF
+    </DataArray>
+    <DataArray type="UInt8" Name="types" format="binary">
+AQAAAAEAAAABAAAACQAAAA==eNrjBAAACgAK
+    </DataArray>
+  </Cells>
+  <PointData Scalars="scalars">
+    <DataArray type="Float64" Name="x" format="binary">
+AQAAACAAAAAgAAAADwAAAA==eNpjYAADBwYcNAANIAEB    </DataArray>
+  </PointData>
+ </Piece>
+ </UnstructuredGrid>
+</VTKFile>
+
+<?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="Float64" Name="x" format="ascii"/>
+    </PPointData>
+    <PPoints>
+      <PDataArray type="Float64" NumberOfComponents="3"/>
+    </PPoints>
+    <Piece Source="output.vtu"/>
+  </PUnstructuredGrid>
+</VTKFile>
+
+DEAL:0::OK
+
+DEAL:1::OK
+
+
+DEAL:2::OK
+