shape_3rd_derivative* checks for update_3rd_derivatives

The functions shape_3rd_derivative and shape_3rd_derivative_component
check that update_flags contains update_hessians, but they should check
that update_flags contain update_3rd_derivatives. Fix this.

diff --git a/include/deal.II/fe/fe_values.h b/include/deal.II/fe/fe_values.h
index 2b0e292ced3c4c3a0617f9669367da14f91c1b2b..30155cc84828516b9cb4da97c9c75b998ca2b474 100644 (file)
--- a/include/deal.II/fe/fe_values.h
+++ b/include/deal.II/fe/fe_values.h
@@ -5157,7 +5157,7 @@ FEValuesBase<dim, spacedim>::shape_3rd_derivative(const unsigned int i,
                                                   const unsigned int j) const
 {
   Assert(i < fe->dofs_per_cell, ExcIndexRange(i, 0, fe->dofs_per_cell));
-  Assert(this->update_flags & update_hessians,
+  Assert(this->update_flags & update_3rd_derivatives,
          ExcAccessToUninitializedField("update_3rd_derivatives"));
   Assert(fe->is_primitive(i), ExcShapeFunctionNotPrimitive(i));
   Assert(present_cell.get() != nullptr,
@@ -5194,7 +5194,7 @@ FEValuesBase<dim, spacedim>::shape_3rd_derivative_component(
   const unsigned int component) const
 {
   Assert(i < fe->dofs_per_cell, ExcIndexRange(i, 0, fe->dofs_per_cell));
-  Assert(this->update_flags & update_hessians,
+  Assert(this->update_flags & update_3rd_derivatives,
          ExcAccessToUninitializedField("update_3rd_derivatives"));
   Assert(component < fe->n_components(),
          ExcIndexRange(component, 0, fe->n_components()));