Refactored element type check to use string comparison as suggested in #4218.

Author: pbehne

src/systems/variational_smoother_system.C

@@ -139,8 +139,10 @@ void VariationalSmootherSystem::prepare_for_smoothing()
           target_elem_inverse_jacobian_dets[elem->type()] =
               std::vector<Real>(nq_points, 1.0);
 
+          const auto type_str = Utility::enum_to_string(elem->type());
+
           // Elems deriving from Tri
-          if (elem->type() == TRI3 || elem->type() == TRI6)
+          if (type_str.compare(0, 3, "TRI") == 0)
             {
 
               // The target element will be an equilateral triangle with area equal to
@@ -180,7 +182,7 @@ void VariationalSmootherSystem::prepare_for_smoothing()
                   }
 
                 default:
-                  libmesh_error_msg("Unsupported triangular element!");
+                  libmesh_error_msg("Unsupported triangular element: " << type_str);
                   break;
                 }
 

tests/mesh/mesh_smoother_test.C

@@ -283,7 +283,7 @@ public:
     // Get mesh dimension, determine whether type is triangular
     const auto * ref_elem = &(ReferenceElem::get(type));
     const auto dim = ref_elem->dim();
-    const bool type_is_tri = dynamic_cast<const Tri*>(ref_elem);
+    const bool type_is_tri = Utility::enum_to_string(type).compare(0, 3, "TRI") == 0;
 
     unsigned int n_elems_per_side = 5;