FDM can fix high-order SIDE_HIERARCIC for now

Working too slowly is better than not working.

Author: roystgnr

src/fe/fe_hierarchic_shape_3D.C

@@ -1230,6 +1230,10 @@ Real FE<3,SIDE_HIERARCHIC>::shape_deriv(const Elem * elem,
     {
     case HEX27:
       {
+        // I need to debug the p>2 case here...
+        if (totalorder > 2)
+          return fe_fdm_deriv(elem, order, i, j, p, add_p_level, FE<3,SIDE_HIERARCHIC>::shape);
+
         const unsigned int dofs_per_side = (totalorder+1u)*(totalorder+1u);
         libmesh_assert_less(i, 6*dofs_per_side);
 
@@ -1486,6 +1490,11 @@ Real FE<3,SIDE_HIERARCHIC>::shape_second_deriv(const Elem * elem,
     {
     case HEX27:
       {
+        // I need to debug the p>2 case here...
+        if (totalorder > 2)
+          return fe_fdm_second_deriv(elem, order, i, j, p, add_p_level,
+                                     FE<3,SIDE_HIERARCHIC>::shape_deriv);
+
         const unsigned int dofs_per_side = (totalorder+1u)*(totalorder+1u);
         libmesh_assert_less(i, 6*dofs_per_side);
 

src/fe/fe_interface.C

@@ -2676,8 +2676,7 @@ unsigned int FEInterface::max_order(const FEType & fe_t,
         case HEX20:
           return 0;
         case HEX27:
-          return 2; // p=3+ is still buggy in 3D
-          // return unlimited;
+          return unlimited;
         case PRISM6:
         case PRISM15:
         case PRISM18: