From dc7b63a6553cb776b2a90ff0ec59a83285510038 Mon Sep 17 00:00:00 2001 From: kronbichler 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 void FE_Q::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; idegree+1; - - this->unit_face_support_points.resize(n); + this->unit_face_support_points.resize(Utilities::fixed_power(this->degree+1)); const std::vector &face_index_map_inverse= FE_Q_Helper::invert_numbering(face_index_map); @@ -1403,15 +1397,9 @@ template void FE_Q::initialize_unit_face_support_points (const Quadrature<1> &points) { const unsigned int codim = dim-1; + this->unit_face_support_points.resize(Utilities::fixed_power(this->degree+1)); - // number of points: (degree+1)^codim - unsigned int n = this->degree+1; - for (unsigned int i=1; idegree+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 &face_index_map_inverse= FE_Q_Helper::invert_numbering(face_index_map); -- 2.39.5