use FE info

diff --git a/source/fe/mapping_fe.cc b/source/fe/mapping_fe.cc
index b84cad10f3b98e33d16169190379f358915bceaf..7f6adef6a51317c06918b9e286435b829658e343 100644 (file)
--- a/source/fe/mapping_fe.cc
+++ b/source/fe/mapping_fe.cc
@@ -2321,9 +2321,7 @@ MappingFE<dim, spacedim>::is_compatible_with(
                     Utilities::to_string(cell_type.get_dimension()) +
                     " ) do not agree."));
 
-  // TODO: query fe->reference_cell() to compare once is exists.
-
-  return true;
+  return fe->reference_cell_type() == cell_type;
 }
 
 

diff --git a/source/fe/mapping_fe_field.cc b/source/fe/mapping_fe_field.cc
index 9d1d800389b79cf6f07fdc6a0abe26a4c1b15806..51a9a4b86d6f3a53ce3c137a87b6deabdb131b1e 100644 (file)
--- a/source/fe/mapping_fe_field.cc
+++ b/source/fe/mapping_fe_field.cc
@@ -354,11 +354,7 @@ MappingFEField<dim, spacedim, VectorType, void>::is_compatible_with(
                     Utilities::to_string(cell_type.get_dimension()) +
                     " ) do not agree."));
 
-
-  // TODO: check euler_dof_handler->get_fe()
-
-
-  return true;
+  return euler_dof_handler->get_fe().reference_cell_type() == cell_type;
 }