From: bangerth Date: Sat, 10 Aug 2013 14:00:56 +0000 (+0000) Subject: Make sure the common case actually works... X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=349586343f7a9f341b569699f0b7c9320f208828;p=dealii-svn.git Make sure the common case actually works... git-svn-id: https://svn.dealii.org/trunk@30278 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/source/fe/fe_q_base.cc b/deal.II/source/fe/fe_q_base.cc index 36491b4237..7e54963284 100644 --- a/deal.II/source/fe/fe_q_base.cc +++ b/deal.II/source/fe/fe_q_base.cc @@ -1053,8 +1053,13 @@ face_to_cell_index (const unsigned int face_index, // of pictures of how all the faces can look like with the various // flips and rotations. // - // that said, the Q2 case is easy enough to implement - Assert (this->dofs_per_line <= 1, ExcNotImplemented()); + // that said, the Q2 case is easy enough to implement, as is the case + // where everything is in standard orientation + Assert ((this->dofs_per_line <= 1) || + ((face_orientation == true) && + (face_flip == false) && + (face_rotation == false)), + ExcNotImplemented()); adjusted_dof_index_on_line = dof_index_on_line; break; } @@ -1079,7 +1084,11 @@ face_to_cell_index (const unsigned int face_index, // the same is true here as above for the 3d case -- someone will // just have to draw a bunch of pictures. in the meantime, // we can implement the Q2 case in which it is simple - Assert (this->dofs_per_quad <= 1, ExcNotImplemented()); + Assert ((this->dofs_per_quad <= 1) || + ((face_orientation == true) && + (face_flip == false) && + (face_rotation == false)), + ExcNotImplemented()); return (this->first_quad_index + face * this->dofs_per_quad + index);