Add MatrixFree::get_active_fe_index and ::get_active_quadrature_index

diff --git a/include/deal.II/matrix_free/fe_evaluation.h b/include/deal.II/matrix_free/fe_evaluation.h
index a05f3db31b121f917d0f64a6a86642cfcd44dba9..0e86aa2209acee3a8eb1bfe0f9751e1b1ccdd22c 100644 (file)
--- a/include/deal.II/matrix_free/fe_evaluation.h
+++ b/include/deal.II/matrix_free/fe_evaluation.h
@@ -245,6 +245,20 @@ public:
   unsigned int
   get_current_cell_index() const;
 
+  /**
+   * Return the active fe index for this class for efficient indexing in the hp
+   * case.
+   */
+  unsigned int
+  get_active_fe_index() const;
+
+  /**
+   * Return the active quadrature index for this class for efficient indexing in
+   * the hp case.
+   */
+  unsigned int
+  get_active_quadrature_index() const;
+
 protected:
   /**
    * Constructor. Made protected to prevent users from directly using this
@@ -3837,6 +3851,26 @@ FEEvaluationBaseData<dim, Number, is_face, VectorizedArrayType>::
 }
 
 
+
+template <int dim, typename Number, bool is_face, typename VectorizedArrayType>
+inline unsigned int
+FEEvaluationBaseData<dim, Number, is_face, VectorizedArrayType>::
+  get_active_fe_index() const
+{
+  return active_fe_index;
+}
+
+
+
+template <int dim, typename Number, bool is_face, typename VectorizedArrayType>
+inline unsigned int
+FEEvaluationBaseData<dim, Number, is_face, VectorizedArrayType>::
+  get_active_quadrature_index() const
+{
+  return active_quad_index;
+}
+
+
 /*----------------------- FEEvaluationBase ----------------------------------*/
 
 template <int dim,