diff --git a/src/fe/fe_hierarchic_shape_3D.C b/src/fe/fe_hierarchic_shape_3D.C index 4b4b378e882..618829e1e51 100644 --- a/src/fe/fe_hierarchic_shape_3D.C +++ b/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); diff --git a/src/fe/fe_interface.C b/src/fe/fe_interface.C index 017073c24ca..0c9f7d0fa7d 100644 --- a/src/fe/fe_interface.C +++ b/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: diff --git a/tests/fe/fe_side_test.C b/tests/fe/fe_side_test.C index 385e6c7f732..aa95e3ae073 100644 --- a/tests/fe/fe_side_test.C +++ b/tests/fe/fe_side_test.C @@ -582,8 +582,8 @@ INSTANTIATE_FESIDETEST(FOURTH, SIDE_HIERARCHIC, TRI7); INSTANTIATE_FESIDETEST(CONSTANT, SIDE_HIERARCHIC, HEX27); INSTANTIATE_FESIDETEST(FIRST, SIDE_HIERARCHIC, HEX27); INSTANTIATE_FESIDETEST(SECOND, SIDE_HIERARCHIC, HEX27); -//INSTANTIATE_FESIDETEST(THIRD, SIDE_HIERARCHIC, HEX27); -//INSTANTIATE_FESIDETEST(FOURTH, SIDE_HIERARCHIC, HEX27); +INSTANTIATE_FESIDETEST(THIRD, SIDE_HIERARCHIC, HEX27); +INSTANTIATE_FESIDETEST(FOURTH, SIDE_HIERARCHIC, HEX27); // INSTANTIATE_FESIDETEST(CONSTANT, SIDE_HIERARCHIC, TET14); INSTANTIATE_FESIDETEST(FIRST, SIDE_HIERARCHIC, TET14);