From: wolf Date: Thu, 13 Nov 2003 15:05:58 +0000 (+0000) Subject: Remove face orientation table. Was wrong anyway, and not used anywhere. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=9495d9cb1b626d8b8ac8e25047c0ba9055a6c06b;p=dealii-svn.git Remove face orientation table. Was wrong anyway, and not used anywhere. git-svn-id: https://svn.dealii.org/trunk@8181 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/include/fe/fe_values.h b/deal.II/deal.II/include/fe/fe_values.h index 585a5b81c4..583dd43251 100644 --- a/deal.II/deal.II/include/fe/fe_values.h +++ b/deal.II/deal.II/include/fe/fe_values.h @@ -985,25 +985,6 @@ class FEFaceValuesBase : public FEValuesBase * @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 @@ -1038,21 +1019,6 @@ class FEFaceValuesBase : public FEValuesBase * formula here. */ const Quadrature 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::faces_per_cell]; }; @@ -1519,7 +1485,6 @@ FEFaceValuesBase::normal_vector (const unsigned int i) const } - template inline const Tensor<1,dim> & @@ -1535,15 +1500,6 @@ FEFaceValuesBase::boundary_form (const unsigned int i) const -template -inline -bool -FEFaceValuesBase::orientation () const -{ - return my_orientation; -} - - template inline const Quadrature & diff --git a/deal.II/deal.II/source/fe/fe_values.cc b/deal.II/deal.II/source/fe/fe_values.cc index 6ae722cd54..7ac4f7d1c0 100644 --- a/deal.II/deal.II/source/fe/fe_values.cc +++ b/deal.II/deal.II/source/fe/fe_values.cc @@ -621,36 +621,6 @@ FEValues::memory_consumption () const /*------------------------------- 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 FEFaceValuesBase::FEFaceValuesBase (const unsigned int n_q_points, const unsigned int dofs_per_cell, @@ -769,7 +739,6 @@ void FEFaceValues::reinit (const typename DoFHandler::cell_iterator &c typename FEValuesBase::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, @@ -871,7 +840,6 @@ void FESubfaceValues::reinit (const typename DoFHandler::cell_iterator "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,