Added tests for fix for wrong degree calculation in project_curl_conforming_l2 for...

diff --git a/tests/numerics/project_bv_curl_conf_03.cc b/tests/numerics/project_bv_curl_conf_03.cc
new file mode 100755 (executable)
index 0000000..9f7b496
--- /dev/null
+++ b/tests/numerics/project_bv_curl_conf_03.cc
@@ -0,0 +1,70 @@
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2014 - 2019 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.
+//
+// ---------------------------------------------------------------------
+// The test checks that project_boundary_values_curl_conforming
+// works for FESystems with FE_Nedelec of different order.
+
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/fe/fe_system.h>
+#include <deal.II/fe/fe_q.h>
+#include <deal.II/fe/fe_nedelec.h>
+#include <deal.II/fe/mapping_q.h>
+#include <deal.II/lac/constraint_matrix.h>
+#include <deal.II/dofs/dof_handler.h>
+#include <deal.II/numerics/vector_tools.h>
+
+#include "../tests.h"
+
+using namespace dealii;
+
+int main(int argc, char **argv) {
+  initlog();
+  const unsigned int dim=3;
+
+  Triangulation<dim> triangulation;
+  GridGenerator::hyper_cube(triangulation, 0, 1, true);
+  MappingQ<dim> mapping(1);
+  ConstraintMatrix constraints;
+
+  //Testing cases with only Nedelec Elements in FE_System
+  {
+    deallog << "Equal order Nedelec degree 0 \t";
+    //Equal-order elements
+    FESystem<dim> fe_system(FE_Nedelec<dim>(0), 1, FE_Nedelec<dim>(0), 1);
+    DoFHandler<dim> dof_handler(triangulation);
+    dof_handler.distribute_dofs(fe_system);
+    
+    constraints.clear();
+    VectorTools::project_boundary_values_curl_conforming_l2(dof_handler, dim,
+                                                           Functions::ZeroFunction<dim>(dim + dim), 0, constraints, mapping);
+    constraints.close();
+    deallog << "OK" << std::endl;
+  }
+
+  {
+    deallog << "Mixed order Nedelec (degree 1,degree 0) \t";
+    //Elements of mixed order
+    FESystem<dim> fe_system(FE_Nedelec<dim>(1), 1, FE_Nedelec<dim>(0), 1);
+    DoFHandler<dim> dof_handler(triangulation);
+    dof_handler.distribute_dofs(fe_system);
+    
+    constraints.clear();
+    VectorTools::project_boundary_values_curl_conforming_l2(dof_handler, dim,
+                                                           Functions::ZeroFunction<dim>(dim + dim), 0, constraints, mapping);
+    constraints.close();
+    deallog << "OK" << std::endl;
+  }
+  
+  return 0;
+}

diff --git a/tests/numerics/project_bv_curl_conf_03.output b/tests/numerics/project_bv_curl_conf_03.output
new file mode 100644 (file)
index 0000000..f544d24
--- /dev/null
+++ b/tests/numerics/project_bv_curl_conf_03.output
@@ -0,0 +1,3 @@
+
+DEAL::Equal order Nedelec degree 0     OK
+DEAL::Mixed order Nedelec (degree 1,degree 0)  OK