* argument.
*/
UpdateFlags compute_update_flags (const UpdateFlags update_flags) const;
+
+ /**
+ * Returns reference to default
+ * @p{MappingQ1} object. Needed
+ * by constructors of derived
+ * classes that uses
+ * @p{MappingQ1} implicitely.
+ */
+ static const Mapping<dim> &get_default_mapping();
};
}
-/*------------------------ Inline functions: FEFaceValuesBase --------------------*/
+/*------------------------ Inline functions: FEValues ----------------------------*/
template <int dim>
};
-
/*------------------------ Inline functions: FEFaceValuesBase --------------------*/
#endif
-// TODO:[RH,GK] replace this by non-global object
-static const MappingQ1<deal_II_dimension> mapping_q1;
-
template <int dim>
void
FEValuesData<dim>::initialize (const unsigned int n_quadrature_points,
+template <int dim>
+const Mapping<dim> &
+FEValuesBase<dim>::get_default_mapping ()
+{
+ static const MappingQ1<dim> mapping_q1;
+ return mapping_q1;
+}
+
+
/*------------------------------- FEValues -------------------------------*/
fe.dofs_per_cell,
1,
update_default,
- mapping_q1,
+ get_default_mapping(),
fe),
quadrature (q)
{
fe.dofs_per_cell,
GeometryInfo<dim>::faces_per_cell,
update_flags,
- mapping_q1,
+ get_default_mapping(),
fe, quadrature)
{
Assert (DEAL_II_COMPAT_MAPPING, ExcCompatibility("mapping"));
GeometryInfo<dim>::faces_per_cell *
GeometryInfo<dim>::subfaces_per_face,
update_flags,
- mapping_q1,
+ get_default_mapping(),
fe, quadrature)
{
Assert (DEAL_II_COMPAT_MAPPING, ExcCompatibility("mapping"));