--- /dev/null
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2017 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.
+//
+// ---------------------------------------------------------------------
+
+
+
+// Similar to no_flux_14 for a higher order element including a high-order
+// mapping. The purpose here is to test that the result of
+// Manifold::normal_vector function inside of
+// VectorTools::compute_no_normal_flux_constraints is accurate for higher
+// order mappings.
+
+
+#include "../tests.h"
+
+#include <deal.II/grid/tria.h>
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/dofs/dof_handler.h>
+#include <deal.II/lac/constraint_matrix.h>
+#include <deal.II/fe/fe_system.h>
+#include <deal.II/fe/mapping_q.h>
+#include <deal.II/numerics/vector_tools.h>
+#include <deal.II/grid/manifold_lib.h>
+
+template <int dim>
+void
+check ()
+{
+ SphericalManifold<dim> spherical;
+ Triangulation<dim> tria;
+ GridGenerator::hyper_shell(tria, Point<dim>(), 0.5, 1., 96, true);
+ tria.set_all_manifold_ids(0);
+ tria.set_manifold(0, spherical);
+
+ ConstraintMatrix cm;
+ MappingQ<dim> mapping(4);
+
+ FESystem<dim> fe(FE_Q<dim>(2),dim);
+ DoFHandler<dim> dofh(tria);
+
+ dofh.distribute_dofs (fe);
+
+ std::set<types::boundary_id> no_normal_flux_boundaries;
+ no_normal_flux_boundaries.insert (0);
+ no_normal_flux_boundaries.insert (1);
+ VectorTools::compute_no_normal_flux_constraints (dofh, 0, no_normal_flux_boundaries, cm, mapping);
+
+ cm.print (deallog.get_file_stream ());
+}
+
+
+
+int main ()
+{
+ initlog();
+ deallog.get_file_stream().precision(8);
+
+ check<3> ();
+}
projects / dealii.git / commitdiff