* In 1d, this function does
* nothing.
*/
- static void build_face_renumbering (const FiniteElementData<dim-1> &fe_data,
- const unsigned int degree,
+ static void build_face_renumbering (const unsigned int degree,
std::vector<unsigned int> &numbering);
/**
poly = new TensorProductPolynomials<dim> (v);
- // do some internal book-keeping
+ // do some internal book-keeping on
+ // cells and faces. if in 1d, the
+ // face function is empty
build_renumbering (*this, degree, renumber);
-//TODO:[WB] externalize this to a proper template function
-#if deal_II_dimension > 1
- build_face_renumbering (FiniteElementData<dim-1>(FE_Q<dim-1>::get_dpo_vector(degree),1),
- degree,
- face_renumber);
-#endif
+ build_face_renumbering (degree, face_renumber);
// copy constraint matrices if they
// are defined. otherwise set them
template <int dim>
void
-FE_Q<dim>::build_face_renumbering (const FiniteElementData<dim-1> &fe_data,
- const unsigned int degree,
+FE_Q<dim>::build_face_renumbering (const unsigned int degree,
std::vector<unsigned int> &numbering)
{
+ FiniteElementData<dim-1> fe_data(FE_Q<dim-1>::get_dpo_vector(degree),1);
FE_Q<dim-1>::build_renumbering (fe_data, degree, numbering);
}
template <>
void
-FE_Q<1>::build_face_renumbering (const FiniteElementData<0> &,
- const unsigned int,
+FE_Q<1>::build_face_renumbering (const unsigned int,
std::vector<unsigned int>&)
{}