From dc7b63a6553cb776b2a90ff0ec59a83285510038 Mon Sep 17 00:00:00 2001
From: kronbichler <kronbichler@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Sat, 24 Nov 2012 14:58:27 +0000
Subject: [PATCH] Simplify fixed_power.

git-svn-id: https://svn.dealii.org/trunk@27678 0785d39b-7218-0410-832d-ea1e28bc413d
---
 deal.II/source/fe/fe_q.cc | 18 +++---------------
 1 file changed, 3 insertions(+), 15 deletions(-)

diff --git a/deal.II/source/fe/fe_q.cc b/deal.II/source/fe/fe_q.cc
index 658c503146..93b8de41dc 100644
--- a/deal.II/source/fe/fe_q.cc
+++ b/deal.II/source/fe/fe_q.cc
@@ -1368,13 +1368,7 @@ template <int dim, int spacedim>
 void FE_Q<dim,spacedim>::initialize_unit_face_support_points ()
 {
   const unsigned int codim = dim-1;
-
-  // number of points: (degree+1)^codim
-  unsigned int n = this->degree+1;
-  for (unsigned int i=1; i<codim; ++i)
-    n *= this->degree+1;
-
-  this->unit_face_support_points.resize(n);
+  this->unit_face_support_points.resize(Utilities::fixed_power<codim>(this->degree+1));
 
   const std::vector<unsigned int> &face_index_map_inverse=
     FE_Q_Helper::invert_numbering(face_index_map);
@@ -1403,15 +1397,9 @@ template <int dim, int spacedim>
 void FE_Q<dim,spacedim>::initialize_unit_face_support_points (const Quadrature<1> &points)
 {
   const unsigned int codim = dim-1;
+  this->unit_face_support_points.resize(Utilities::fixed_power<codim>(this->degree+1));
 
-  // number of points: (degree+1)^codim
-  unsigned int n = this->degree+1;
-  for (unsigned int i=1; i<codim; ++i)
-    n *= this->degree+1;
-
-  this->unit_face_support_points.resize(n);
-
-  const std::vector< Point<1> > edge = points.get_points();
+  const std::vector< Point<1> > &edge = points.get_points();
 
   const std::vector<unsigned int> &face_index_map_inverse=
     FE_Q_Helper::invert_numbering(face_index_map);
-- 
2.39.5