Added Scalar::get_function_third_derivatives_from_local_dof_values

diff --git a/include/deal.II/fe/fe_values.h b/include/deal.II/fe/fe_values.h
index e17caaf5cc4f5c1b05fd54d240ac3380504a025e..e4f5b920e778a540eef710b96c0120dbb1a0fb22 100644 (file)
--- a/include/deal.II/fe/fe_values.h
+++ b/include/deal.II/fe/fe_values.h
@@ -463,6 +463,14 @@ namespace FEValuesViews
                                          std::vector<typename ProductType<third_derivative_type,
                                          typename InputVector::value_type>::type> &third_derivatives) const;
 
+    /**
+     * @copydoc FEValuesViews::Scalar::get_function_values_from_local_dof_values()
+     */
+    template <class InputVector>
+    void get_function_third_derivatives_from_local_dof_values (const InputVector &dof_values,
+                                                               std::vector<typename OutputType<typename InputVector::value_type>::third_derivative_type> &third_derivatives) const;
+
+
   private:
     /**
      * A pointer to the FEValuesBase object we operate on.

diff --git a/source/fe/fe_values.cc b/source/fe/fe_values.cc
index 3f67cd9e3dd406dad8594a093143501e00f300c9..c26294d7ba4103c71d86d248679b6c81473cc3a0 100644 (file)
--- a/source/fe/fe_values.cc
+++ b/source/fe/fe_values.cc
@@ -1465,6 +1465,26 @@ namespace FEValuesViews
 
 
 
+  template <int dim, int spacedim>
+  template <class InputVector>
+  void
+  Scalar<dim,spacedim>::
+  get_function_third_derivatives_from_local_dof_values(const InputVector &dof_values,
+                                                       std::vector<typename OutputType<typename InputVector::value_type>::third_derivative_type> &third_derivatives) const
+  {
+    Assert (fe_values->update_flags & update_3rd_derivatives,
+            (typename FEValuesBase<dim,spacedim>::ExcAccessToUninitializedField("update_3rd_derivatives")));
+    Assert (fe_values->present_cell.get() != nullptr,
+            ExcMessage ("FEValues object is not reinit'ed to any cell"));
+    AssertDimension (dof_values.size(), fe_values->dofs_per_cell);
+
+    internal::do_function_derivatives<3,dim,spacedim>
+    (make_array_view(dof_values.begin(), dof_values.end()),
+     fe_values->finite_element_output.shape_3rd_derivatives, shape_function_data, third_derivatives);
+  }
+
+
+
   template <int dim, int spacedim>
   template <class InputVector>
   void

diff --git a/source/fe/fe_values.inst.in b/source/fe/fe_values.inst.in
index 067db65566d3655aa9f380115653783a40452ffc..69c7931298425e8db0093fcdd18ce528d0961175 100644 (file)
--- a/source/fe/fe_values.inst.in
+++ b/source/fe/fe_values.inst.in
@@ -164,6 +164,11 @@ for (VEC : GENERAL_CONTAINER_TYPES;
     void FEValuesViews::Scalar<deal_II_dimension, deal_II_space_dimension>
     ::get_function_laplacians_from_local_dof_values<VEC<Number>>
             (const VEC<Number>&, std::vector<typename OutputType<Number>::value_type>&) const;
+
+    template
+    void FEValuesViews::Scalar<deal_II_dimension, deal_II_space_dimension>
+    ::get_function_third_derivatives_from_local_dof_values<VEC<Number>>
+            (const VEC<Number>&, std::vector<typename OutputType<Number>::third_derivative_type>&) const;
 #endif
 }