const std::vector<FullMatrix<double> > &dofs_subcell);
/**
- * Initialize the @p unit_support_points field of the FiniteElement class.
+ * Initialize the @p generalized_support_points field of the FiniteElement class.
* Called from the constructor.
*/
- void initialize_unit_support_points ();
+ void initialize_generalized_support_points ();
/**
- * Initialize the @p unit_face_support_points field of the FiniteElement
+ * Initialize the @p generalized_face_support_points field of the FiniteElement
* class. Called from the constructor.
*/
- void initialize_unit_face_support_points ();
+ void initialize_generalized_face_support_points ();
/**
* Mapping from lexicographic to shape function numbering on first face.
/* -------------- declaration of explicit specializations ------------- */
template <>
-void FE_Q_Hierarchical<1>::initialize_unit_face_support_points ();
+void FE_Q_Hierarchical<1>::initialize_generalized_face_support_points ();
template <>
bool
// finally fill in support points
// on cell and face
- initialize_unit_support_points ();
- initialize_unit_face_support_points ();
+ initialize_generalized_support_points ();
+ initialize_generalized_face_support_points ();
}
template <int dim>
-void FE_Q_Hierarchical<dim>::initialize_unit_support_points ()
+void FE_Q_Hierarchical<dim>::initialize_generalized_support_points ()
{
// number of points: (degree+1)^dim
unsigned int n = this->degree+1;
for (unsigned int i=1; i<dim; ++i)
n *= this->degree+1;
- this->unit_support_points.resize(n);
+ this->generalized_support_points.resize(n);
const std::vector<unsigned int> &index_map_inverse=
this->poly_space.get_numbering_inverse();
else
p(2) = .5;
}
- this->unit_support_points[index_map_inverse[k++]] = p;
+ this->generalized_support_points[index_map_inverse[k++]] = p;
};
}
template <>
-void FE_Q_Hierarchical<1>::initialize_unit_face_support_points ()
+void FE_Q_Hierarchical<1>::initialize_generalized_face_support_points ()
{
// no faces in 1d, so nothing to do
}
template <int dim>
-void FE_Q_Hierarchical<dim>::initialize_unit_face_support_points ()
+void FE_Q_Hierarchical<dim>::initialize_generalized_face_support_points ()
{
const unsigned int codim = dim-1;
for (unsigned int i=1; i<codim; ++i)
n *= this->degree+1;
- this->unit_face_support_points.resize(n);
+ this->generalized_face_support_points.resize(n);
Point<codim> p;
else
p(2) = .5;
}
- this->unit_face_support_points[face_renumber[k++]] = p;
+ this->generalized_face_support_points[face_renumber[k++]] = p;
};
}