From: Peter Munch <peterrmuench@gmail.com>
Date: Thu, 15 Oct 2020 04:45:14 +0000 (+0200)
Subject: Add MatrixFree::get_active_fe_index and ::get_active_quadrature_index
X-Git-Tag: v9.3.0-rc1~1007^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=e41ce6b30470c082352a4952f0e825042c80da73;p=dealii.git

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 a05f3db31b..0e86aa2209 100644
--- 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,