From: Marc Fehling <mafehling.git@gmail.com>
Date: Thu, 31 Dec 2020 04:49:02 +0000 (-0700)
Subject: Added tests for FENothing constraints.
X-Git-Tag: v9.3.0-rc1~539^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F11448%2Fhead;p=dealii.git

Added tests for FENothing constraints.
---

diff --git a/tests/hp/fe_nothing_dominates.h b/tests/hp/fe_nothing_dominates.h
new file mode 100644
index 0000000000..1836c0bac3
--- /dev/null
+++ b/tests/hp/fe_nothing_dominates.h
@@ -0,0 +1,92 @@
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 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.
+//
+// ---------------------------------------------------------------------
+
+
+
+// Check for continuity requirements of two neighboring finite elements
+// and project a given function subject to these constraints.
+
+
+#include <deal.II/base/function.h>
+
+#include <deal.II/dofs/dof_handler.h>
+#include <deal.II/dofs/dof_tools.h>
+
+#include <deal.II/grid/tria.h>
+
+#include <deal.II/hp/fe_collection.h>
+#include <deal.II/hp/q_collection.h>
+
+#include <deal.II/lac/affine_constraints.h>
+
+#include <deal.II/numerics/vector_tools_project.h>
+
+#include "../tests.h"
+
+#include "../test_grids.h"
+
+
+template <int dim>
+void
+project(const hp::FECollection<dim> &fe_collection,
+        const hp::QCollection<dim> & q_collection,
+        const Function<dim> &        function)
+{
+#ifdef DEBUG
+  Assert(fe_collection.size() == 2, ExcInternalError());
+  Assert(q_collection.size() == 2, ExcInternalError());
+  for (unsigned int f = 0; f < fe_collection.size(); ++f)
+    Assert(fe_collection[f].n_components() == function.n_components,
+           ExcInternalError());
+#endif
+
+  // setup
+  // +---+---+
+  // | 0 | 1 |
+  // +---+---+
+  Triangulation<dim> tria;
+  TestGrids::hyper_line(tria, 2);
+
+  DoFHandler<dim> dofh(tria);
+  (++(dofh.begin_active()))->set_active_fe_index(1);
+  dofh.distribute_dofs(fe_collection);
+
+  // make constraints
+  AffineConstraints<double> constraints;
+  DoFTools::make_hanging_node_constraints(dofh, constraints);
+  constraints.close();
+  deallog << "constraints:" << std::endl;
+  constraints.print(deallog.get_file_stream());
+
+  // project function with constraints
+  Vector<double> solution(dofh.n_dofs());
+  VectorTools::project(dofh, constraints, q_collection, function, solution);
+
+  // verify output
+  deallog << "dof values:" << std::endl;
+  Vector<double> cell_values;
+  for (const auto &cell : dofh.active_cell_iterators())
+    {
+      cell_values.reinit(cell->get_fe().n_dofs_per_cell());
+      cell->get_dof_values(solution, cell_values);
+
+      deallog << " cell " << cell->active_cell_index() << ":";
+      for (const auto &value : cell_values)
+        deallog << " " << value;
+      deallog << std::endl;
+    }
+
+  deallog << "OK" << std::endl;
+}
diff --git a/tests/hp/fe_nothing_dominates_01.cc b/tests/hp/fe_nothing_dominates_01.cc
new file mode 100644
index 0000000000..8ddb354370
--- /dev/null
+++ b/tests/hp/fe_nothing_dominates_01.cc
@@ -0,0 +1,84 @@
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 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.
+//
+// ---------------------------------------------------------------------
+
+
+
+// Check for continuity requirements for neighboring FE_Nothing and FE_Q
+// elements.
+//
+// Note: output corresponds to current results, which are not correct.
+
+
+#include <deal.II/base/function.h>
+#include <deal.II/base/quadrature.h>
+#include <deal.II/base/quadrature_lib.h>
+
+#include <deal.II/fe/fe_nothing.h>
+#include <deal.II/fe/fe_q.h>
+
+#include <deal.II/hp/fe_collection.h>
+#include <deal.II/hp/q_collection.h>
+
+#include "../tests.h"
+
+#include "fe_nothing_dominates.h"
+
+
+template <int dim>
+void
+test(const bool fe_nothing_dominates)
+{
+  Functions::ConstantFunction<dim> function(1.);
+
+  {
+    deallog << "(FE_Nothing, FE_Q)" << std::endl;
+    hp::FECollection<dim> fe_collection(FE_Nothing<dim>(/*n_components=*/1,
+                                                        fe_nothing_dominates),
+                                        FE_Q<dim>(1));
+    hp::QCollection<dim>  q_collection(Quadrature<dim>(1), QGauss<dim>(2));
+    project(fe_collection, q_collection, function);
+  }
+
+  {
+    deallog << "(FE_Q, FE_Nothing)" << std::endl;
+    hp::FECollection<dim> fe_collection(FE_Q<dim>(1),
+                                        FE_Nothing<dim>(/*n_components=*/1,
+                                                        fe_nothing_dominates));
+    hp::QCollection<dim>  q_collection(QGauss<dim>(2), Quadrature<dim>(1));
+    project(fe_collection, q_collection, function);
+  }
+}
+
+
+int
+main()
+{
+  initlog();
+
+  deallog << std::boolalpha;
+  for (const bool fe_nothing_dominates : {false, true})
+    {
+      deallog << "FE_Nothing dominates: " << fe_nothing_dominates << std::endl;
+      deallog.push("1d");
+      test<1>(fe_nothing_dominates);
+      deallog.pop();
+      deallog.push("2d");
+      test<2>(fe_nothing_dominates);
+      deallog.pop();
+      deallog.push("3d");
+      test<3>(fe_nothing_dominates);
+      deallog.pop();
+    }
+}
diff --git a/tests/hp/fe_nothing_dominates_01.output b/tests/hp/fe_nothing_dominates_01.output
new file mode 100644
index 0000000000..8428e72878
--- /dev/null
+++ b/tests/hp/fe_nothing_dominates_01.output
@@ -0,0 +1,75 @@
+
+DEAL::FE_Nothing dominates: false
+DEAL:1d::(FE_Nothing, FE_Q)
+DEAL:1d::constraints:
+DEAL:1d::dof values:
+DEAL:1d:: cell 0:
+DEAL:1d:: cell 1: 1 1
+DEAL:1d::OK
+DEAL:1d::(FE_Q, FE_Nothing)
+DEAL:1d::constraints:
+DEAL:1d::dof values:
+DEAL:1d:: cell 0: 1 1
+DEAL:1d:: cell 1:
+DEAL:1d::OK
+DEAL:2d::(FE_Nothing, FE_Q)
+DEAL:2d::constraints:
+DEAL:2d::dof values:
+DEAL:2d:: cell 0:
+DEAL:2d:: cell 1: 1 1 1 1
+DEAL:2d::OK
+DEAL:2d::(FE_Q, FE_Nothing)
+DEAL:2d::constraints:
+DEAL:2d::dof values:
+DEAL:2d:: cell 0: 1 1 1 1
+DEAL:2d:: cell 1:
+DEAL:2d::OK
+DEAL:3d::(FE_Nothing, FE_Q)
+DEAL:3d::constraints:
+DEAL:3d::dof values:
+DEAL:3d:: cell 0:
+DEAL:3d:: cell 1: 1 1 1 1 1 1 1 1
+DEAL:3d::OK
+DEAL:3d::(FE_Q, FE_Nothing)
+DEAL:3d::constraints:
+DEAL:3d::dof values:
+DEAL:3d:: cell 0: 1 1 1 1 1 1 1 1
+DEAL:3d:: cell 1:
+DEAL:3d::OK
+DEAL::FE_Nothing dominates: true
+DEAL:1d::(FE_Nothing, FE_Q)
+DEAL:1d::constraints:
+DEAL:1d::dof values:
+DEAL:1d:: cell 0:
+DEAL:1d:: cell 1: 1 1
+DEAL:1d::OK
+DEAL:1d::(FE_Q, FE_Nothing)
+DEAL:1d::constraints:
+DEAL:1d::dof values:
+DEAL:1d:: cell 0: 1 1
+DEAL:1d:: cell 1:
+DEAL:1d::OK
+DEAL:2d::(FE_Nothing, FE_Q)
+DEAL:2d::constraints:
+DEAL:2d::dof values:
+DEAL:2d:: cell 0:
+DEAL:2d:: cell 1: 1 1 1 1
+DEAL:2d::OK
+DEAL:2d::(FE_Q, FE_Nothing)
+DEAL:2d::constraints:
+DEAL:2d::dof values:
+DEAL:2d:: cell 0: 1 1 1 1
+DEAL:2d:: cell 1:
+DEAL:2d::OK
+DEAL:3d::(FE_Nothing, FE_Q)
+DEAL:3d::constraints:
+DEAL:3d::dof values:
+DEAL:3d:: cell 0:
+DEAL:3d:: cell 1: 1 1 1 1 1 1 1 1
+DEAL:3d::OK
+DEAL:3d::(FE_Q, FE_Nothing)
+DEAL:3d::constraints:
+DEAL:3d::dof values:
+DEAL:3d:: cell 0: 1 1 1 1 1 1 1 1
+DEAL:3d:: cell 1:
+DEAL:3d::OK
diff --git a/tests/hp/fe_nothing_dominates_02.cc b/tests/hp/fe_nothing_dominates_02.cc
new file mode 100644
index 0000000000..843b20c302
--- /dev/null
+++ b/tests/hp/fe_nothing_dominates_02.cc
@@ -0,0 +1,80 @@
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 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.
+//
+// ---------------------------------------------------------------------
+
+
+
+// Check for continuity requirements for neighboring
+// (FE_NothingxFE_Nothing) and (FE_QxFE_Q) elements.
+// The twist: only one FE_Nothing elements dominates.
+//
+// Note: output corresponds to current results, which are not correct.
+
+
+#include <deal.II/base/function.h>
+#include <deal.II/base/quadrature.h>
+#include <deal.II/base/quadrature_lib.h>
+
+#include <deal.II/fe/fe_nothing.h>
+#include <deal.II/fe/fe_q.h>
+#include <deal.II/fe/fe_system.h>
+
+#include <deal.II/hp/fe_collection.h>
+#include <deal.II/hp/q_collection.h>
+
+#include "../tests.h"
+
+#include "fe_nothing_dominates.h"
+
+
+template <int dim>
+void
+test()
+{
+  Functions::ConstantFunction<dim> function(1., 2);
+
+  {
+    deallog << "(FE_QxFE_Q, FE_Nothing(true)xFE_Nothing(false))" << std::endl;
+    hp::FECollection<dim> fe_collection(
+      FESystem<dim>(FE_Q<dim>(1), FE_Q<dim>(1)),
+      FESystem<dim>(FE_Nothing<dim>(1, true), FE_Nothing<dim>(1, false)));
+    hp::QCollection<dim> q_collection(QGauss<dim>(2), Quadrature<dim>(1));
+    project(fe_collection, q_collection, function);
+  }
+  {
+    deallog << "(FE_Nothing(true)xFE_Nothing(false), FE_QxFE_Q)" << std::endl;
+    hp::FECollection<dim> fe_collection(
+      FESystem<dim>(FE_Nothing<dim>(1, true), FE_Nothing<dim>(1, false)),
+      FESystem<dim>(FE_Q<dim>(1), FE_Q<dim>(1)));
+    hp::QCollection<dim> q_collection(Quadrature<dim>(1), QGauss<dim>(2));
+    project(fe_collection, q_collection, function);
+  }
+}
+
+
+int
+main()
+{
+  initlog();
+
+  deallog.push("1d");
+  test<1>();
+  deallog.pop();
+  deallog.push("2d");
+  test<2>();
+  deallog.pop();
+  deallog.push("3d");
+  test<3>();
+  deallog.pop();
+}
diff --git a/tests/hp/fe_nothing_dominates_02.output b/tests/hp/fe_nothing_dominates_02.output
new file mode 100644
index 0000000000..7ee5d38428
--- /dev/null
+++ b/tests/hp/fe_nothing_dominates_02.output
@@ -0,0 +1,37 @@
+
+DEAL:1d::(FE_QxFE_Q, FE_Nothing(true)xFE_Nothing(false))
+DEAL:1d::constraints:
+DEAL:1d::dof values:
+DEAL:1d:: cell 0: 1 1 1 1
+DEAL:1d:: cell 1:
+DEAL:1d::OK
+DEAL:1d::(FE_Nothing(true)xFE_Nothing(false), FE_QxFE_Q)
+DEAL:1d::constraints:
+DEAL:1d::dof values:
+DEAL:1d:: cell 0:
+DEAL:1d:: cell 1: 1 1 1 1
+DEAL:1d::OK
+DEAL:2d::(FE_QxFE_Q, FE_Nothing(true)xFE_Nothing(false))
+DEAL:2d::constraints:
+DEAL:2d::dof values:
+DEAL:2d:: cell 0: 1 1 1 1 1 1 1 1
+DEAL:2d:: cell 1:
+DEAL:2d::OK
+DEAL:2d::(FE_Nothing(true)xFE_Nothing(false), FE_QxFE_Q)
+DEAL:2d::constraints:
+DEAL:2d::dof values:
+DEAL:2d:: cell 0:
+DEAL:2d:: cell 1: 1 1 1 1 1 1 1 1
+DEAL:2d::OK
+DEAL:3d::(FE_QxFE_Q, FE_Nothing(true)xFE_Nothing(false))
+DEAL:3d::constraints:
+DEAL:3d::dof values:
+DEAL:3d:: cell 0: 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+DEAL:3d:: cell 1:
+DEAL:3d::OK
+DEAL:3d::(FE_Nothing(true)xFE_Nothing(false), FE_QxFE_Q)
+DEAL:3d::constraints:
+DEAL:3d::dof values:
+DEAL:3d:: cell 0:
+DEAL:3d:: cell 1: 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+DEAL:3d::OK
diff --git a/tests/hp/fe_nothing_dominates_03.cc b/tests/hp/fe_nothing_dominates_03.cc
new file mode 100644
index 0000000000..6e37f67671
--- /dev/null
+++ b/tests/hp/fe_nothing_dominates_03.cc
@@ -0,0 +1,79 @@
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 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.
+//
+// ---------------------------------------------------------------------
+
+
+
+// Check for continuity requirements for neighboring
+// (FE_QxFE_Nothing) and (FE_NothingxFE_Q) elements.
+// The twist: only one FE_Nothing element dominates.
+//
+// Note: output corresponds to current results, which are not correct.
+
+
+#include <deal.II/base/function.h>
+#include <deal.II/base/quadrature.h>
+#include <deal.II/base/quadrature_lib.h>
+
+#include <deal.II/fe/fe_nothing.h>
+#include <deal.II/fe/fe_q.h>
+#include <deal.II/fe/fe_system.h>
+
+#include <deal.II/hp/fe_collection.h>
+#include <deal.II/hp/q_collection.h>
+
+#include "../tests.h"
+
+#include "fe_nothing_dominates.h"
+
+
+template <int dim>
+void
+test()
+{
+  Functions::ConstantFunction<dim> function(1., 2);
+  hp::QCollection<dim>             q_collection(QGauss<dim>(2), QGauss<dim>(2));
+
+  {
+    deallog << "(FE_Nothing(true)xFE_Q, FE_QxFE_Nothing(false))" << std::endl;
+    hp::FECollection<dim> fe_collection(
+      FESystem<dim>(FE_Nothing<dim>(1, true), FE_Q<dim>(1)),
+      FESystem<dim>(FE_Q<dim>(1), FE_Nothing<dim>(1, false)));
+    project(fe_collection, q_collection, function);
+  }
+  {
+    deallog << "(FE_QxFE_Nothing(false), FE_Nothing(true)xFE_Q)" << std::endl;
+    hp::FECollection<dim> fe_collection(
+      FESystem<dim>(FE_Q<dim>(1), FE_Nothing<dim>(1, false)),
+      FESystem<dim>(FE_Nothing<dim>(1, true), FE_Q<dim>(1)));
+    project(fe_collection, q_collection, function);
+  }
+}
+
+
+int
+main()
+{
+  initlog();
+
+  deallog.push("1d");
+  test<1>();
+  deallog.pop();
+  deallog.push("2d");
+  test<2>();
+  deallog.pop();
+  deallog.push("3d");
+  test<3>();
+  deallog.pop();
+}
diff --git a/tests/hp/fe_nothing_dominates_03.output b/tests/hp/fe_nothing_dominates_03.output
new file mode 100644
index 0000000000..545c03e92c
--- /dev/null
+++ b/tests/hp/fe_nothing_dominates_03.output
@@ -0,0 +1,49 @@
+
+DEAL:1d::(FE_Nothing(true)xFE_Q, FE_QxFE_Nothing(false))
+DEAL:1d::constraints:
+DEAL:1d::dof values:
+DEAL:1d:: cell 0: 1 1
+DEAL:1d:: cell 1: 1 1
+DEAL:1d::OK
+DEAL:1d::(FE_QxFE_Nothing(false), FE_Nothing(true)xFE_Q)
+DEAL:1d::constraints:
+DEAL:1d::dof values:
+DEAL:1d:: cell 0: 1 1
+DEAL:1d:: cell 1: 1 1
+DEAL:1d::OK
+DEAL:2d::(FE_Nothing(true)xFE_Q, FE_QxFE_Nothing(false))
+DEAL:2d::constraints:
+    4 = 0
+    6 = 0
+DEAL:2d::dof values:
+DEAL:2d:: cell 0: 1 1 1 1
+DEAL:2d:: cell 1: 0 1.5 0 1.5
+DEAL:2d::OK
+DEAL:2d::(FE_QxFE_Nothing(false), FE_Nothing(true)xFE_Q)
+DEAL:2d::constraints:
+    1 = 0
+    3 = 0
+DEAL:2d::dof values:
+DEAL:2d:: cell 0: 1.5 0 1.5 0
+DEAL:2d:: cell 1: 1 1 1 1
+DEAL:2d::OK
+DEAL:3d::(FE_Nothing(true)xFE_Q, FE_QxFE_Nothing(false))
+DEAL:3d::constraints:
+    8 = 0
+    10 = 0
+    12 = 0
+    14 = 0
+DEAL:3d::dof values:
+DEAL:3d:: cell 0: 1 1 1 1 1 1 1 1
+DEAL:3d:: cell 1: 0 1.5 0 1.5 0 1.5 0 1.5
+DEAL:3d::OK
+DEAL:3d::(FE_QxFE_Nothing(false), FE_Nothing(true)xFE_Q)
+DEAL:3d::constraints:
+    1 = 0
+    3 = 0
+    5 = 0
+    7 = 0
+DEAL:3d::dof values:
+DEAL:3d:: cell 0: 1.5 0 1.5 0 1.5 0 1.5 0
+DEAL:3d:: cell 1: 1 1 1 1 1 1 1 1
+DEAL:3d::OK