]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Remove assert in dof identities for DGP 17923/head
authorMarco Feder <marco.feder@sissa.it>
Wed, 11 Dec 2024 09:08:53 +0000 (09:08 +0000)
committerMarco Feder <marco.feder@sissa.it>
Wed, 11 Dec 2024 09:08:53 +0000 (09:08 +0000)
source/fe/fe_dgp.cc

index acd6ec9442fdb7d8e0e3a66114411a14b04ed206..23c1d6bf1645a9af5b3601f30957037847eadf0e 100644 (file)
@@ -166,16 +166,13 @@ FE_DGP<dim, spacedim>::hp_constraints_are_implemented() const
 template <int dim, int spacedim>
 std::vector<std::pair<unsigned int, unsigned int>>
 FE_DGP<dim, spacedim>::hp_vertex_dof_identities(
-  const FiniteElement<dim, spacedim> &fe_other) const
+  const FiniteElement<dim, spacedim> & /*fe_other*/) const
 {
-  // there are no such constraints for DGP elements at all
-  if (dynamic_cast<const FE_DGP<dim, spacedim> *>(&fe_other) != nullptr)
-    return std::vector<std::pair<unsigned int, unsigned int>>();
-  else
-    {
-      DEAL_II_NOT_IMPLEMENTED();
-      return std::vector<std::pair<unsigned int, unsigned int>>();
-    }
+  // this element is discontinuous, so by definition there can
+  // be no identities between its dofs and those of any neighbor
+  // (of whichever type the neighbor may be -- after all, we have
+  // no face dofs on this side to begin with)
+  return std::vector<std::pair<unsigned int, unsigned int>>();
 }
 
 

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.