* @p{n.ds}.
*/
const Tensor<1,dim> & boundary_form (const unsigned int i) const;
-
- /**
- * Get the natural orientation of
- * the face.
- *
- * If this is true, then the
- * natural normal vector of the
- * face points outwards (faces
- * 0,1 in two dimensions, faces
- * 0,2,4 in three). On the other
- * faces, the natural normal
- * vector points inwards.
- *
- * With this function, you have
- * access to a consistent and
- * unique orientation of all
- * faces.
- */
- bool orientation () const;
/**
* Return the list of outward normal
* formula here.
*/
const Quadrature<dim-1> quadrature;
-
- /**
- * Value for orientation function.
- */
- bool my_orientation;
-
- /**
- * Data table for orientation
- * values. Since the faces are
- * numbered in a creative
- * fashion, we stores these
- * values instead of computing
- * them.
- */
- static const bool orientation_table[GeometryInfo<dim>::faces_per_cell];
};
}
-
template <int dim>
inline
const Tensor<1,dim> &
-template <int dim>
-inline
-bool
-FEFaceValuesBase<dim>::orientation () const
-{
- return my_orientation;
-}
-
-
template <int dim>
inline
const Quadrature<dim-1> &
/*------------------------------- FEFaceValuesBase --------------------------*/
-#if deal_II_dimension == 1
-template <>
-const bool
-FEFaceValuesBase<1>::orientation_table[GeometryInfo<1>::faces_per_cell] =
-{
- false, true
-};
-#endif
-
-
-#if deal_II_dimension == 2
-template <>
-const bool
-FEFaceValuesBase<2>::orientation_table[GeometryInfo<2>::faces_per_cell] =
-{
- true, true, false, false
-};
-#endif
-
-
-#if deal_II_dimension == 3
-template <>
-const bool
-FEFaceValuesBase<3>::orientation_table[GeometryInfo<3>::faces_per_cell] =
-{
- true, false, true, false, true, false
-};
-#endif
-
-
template <int dim>
FEFaceValuesBase<dim>::FEFaceValuesBase (const unsigned int n_q_points,
const unsigned int dofs_per_cell,
typename FEValuesBase<dim>::ExcFEDontMatch());
this->present_cell = cell;
- this->my_orientation = this->orientation_table[face_no];
this->get_mapping().fill_fe_face_values(cell, face_no,
this->quadrature,
"instead in these cases."));
this->present_cell = cell;
- this->my_orientation = this->orientation_table[face_no];
this->get_mapping().fill_fe_subface_values(cell, face_no, subface_no,
this->quadrature,