Add a test for pyramid output with gnuplot.

diff --git a/tests/simplex/data_out_write_gnuplot_02.cc b/tests/simplex/data_out_write_gnuplot_02.cc
new file mode 100644 (file)
index 0000000..47f5969
--- /dev/null
+++ b/tests/simplex/data_out_write_gnuplot_02.cc
@@ -0,0 +1,111 @@
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2020 - 2021 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.
+//
+// ---------------------------------------------------------------------
+
+
+
+// Test DataOut::write_gnuplot() for a pyramid mesh
+
+#include <deal.II/base/quadrature_lib.h>
+
+#include <deal.II/dofs/dof_handler.h>
+
+#include <deal.II/fe/fe_pyramid_p.h>
+#include <deal.II/fe/fe_simplex_p.h>
+#include <deal.II/fe/fe_simplex_p_bubbles.h>
+#include <deal.II/fe/fe_system.h>
+#include <deal.II/fe/fe_wedge_p.h>
+#include <deal.II/fe/mapping_fe.h>
+
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/grid_in.h>
+#include <deal.II/grid/grid_out.h>
+#include <deal.II/grid/tria.h>
+
+#include <deal.II/numerics/data_out.h>
+#include <deal.II/numerics/vector_tools.h>
+
+#include "../tests.h"
+
+using namespace dealii;
+
+template <int dim>
+class RightHandSideFunction : public Function<dim>
+{
+public:
+  RightHandSideFunction(const unsigned int n_components)
+    : Function<dim>(n_components)
+  {}
+
+  virtual double
+  value(const Point<dim> &p, const unsigned int component = 0) const
+  {
+    return p[component % dim] * p[component % dim];
+  }
+};
+
+
+
+template <int dim, int spacedim = dim>
+void
+test(const FiniteElement<dim, spacedim> &fe)
+{
+  Triangulation<dim, spacedim> tria;
+  GridGenerator::reference_cell(tria, ReferenceCells::Pyramid);
+
+  DoFHandler<dim> dof_handler(tria);
+
+  dof_handler.distribute_dofs(fe);
+
+  Vector<double> solution(dof_handler.n_dofs());
+
+  MappingFE<dim> mapping(FE_PyramidP<dim>(1));
+
+  AffineConstraints<double> dummy;
+  dummy.close();
+
+  VectorTools::project(mapping,
+                       dof_handler,
+                       dummy,
+                       QGaussPyramid<dim>(fe.tensor_degree() + 1),
+                       RightHandSideFunction<dim>(1),
+                       solution);
+
+
+  DataOut<dim> data_out;
+
+  data_out.attach_dof_handler(dof_handler);
+  data_out.add_data_vector(solution, "solution");
+
+
+  data_out.build_patches(mapping, 1);
+
+#if false
+  static unsigned int counter = 0;
+      std::ofstream output("test." + std::to_string(dim) + "." +
+                           std::to_string(counter++) + ".gnuplot");
+      data_out.write_gnuplot(output);
+#else
+  data_out.write_gnuplot(deallog.get_file_stream());
+#endif
+}
+
+int
+main()
+{
+  initlog();
+
+  const int dim = 3;
+  test<dim>(FE_PyramidP<dim>(1));
+}

diff --git a/tests/simplex/data_out_write_gnuplot_02.output b/tests/simplex/data_out_write_gnuplot_02.output
new file mode 100644 (file)
index 0000000..a627d9d
--- /dev/null
+++ b/tests/simplex/data_out_write_gnuplot_02.output
@@ -0,0 +1,20 @@
+
+# This file was generated by the deal.II library.
+
+
+#
+# For a description of the GNUPLOT format see the GNUPLOT manual.
+#
+# <x> <y> <z> <solution> 
+-1.00000 -1.00000 0.00000 0.311111 
+1.00000 -1.00000 0.00000 0.311111 
+1.00000 1.00000 0.00000 0.311111 
+-1.00000 1.00000 0.00000 0.311111 
+-1.00000 -1.00000 0.00000 0.311111 
+0.00000 0.00000 1.00000 -0.133333 
+1.00000 -1.00000 0.00000 0.311111 
+
+1.00000 1.00000 0.00000 0.311111 
+0.00000 0.00000 1.00000 -0.133333 
+-1.00000 1.00000 0.00000 0.311111 
+