From 3cbd40273e36bc3636dc80f9f46f8b9d86dff47c Mon Sep 17 00:00:00 2001
From: Marco Feder <marco.feder@sissa.it>
Date: Wed, 11 Dec 2024 09:08:53 +0000
Subject: [PATCH] Remove assert in dof identities for  DGP

---
 source/fe/fe_dgp.cc | 15 ++++++---------
 1 file changed, 6 insertions(+), 9 deletions(-)

diff --git a/source/fe/fe_dgp.cc b/source/fe/fe_dgp.cc
index acd6ec9442..23c1d6bf16 100644
--- a/source/fe/fe_dgp.cc
+++ b/source/fe/fe_dgp.cc
@@ -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>>();
 }
 
 
-- 
2.39.5