From: bangerth Date: Fri, 19 Jan 2007 20:28:50 +0000 (+0000) Subject: Add a comment X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=999dc3d3959068e0822fd1cc9989f4ee3e6e2725;p=dealii-svn.git Add a comment git-svn-id: https://svn.dealii.org/trunk@14350 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/include/dofs/dof_accessor.templates.h b/deal.II/deal.II/include/dofs/dof_accessor.templates.h index 124e2f9065..612dd8e571 100644 --- a/deal.II/deal.II/include/dofs/dof_accessor.templates.h +++ b/deal.II/deal.II/include/dofs/dof_accessor.templates.h @@ -2,7 +2,7 @@ // $Id$ // Version: $Name$ // -// Copyright (C) 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006 by the deal.II authors +// Copyright (C) 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007 by the deal.II authors // // This file is subject to QPL and may not be distributed // without copyright and license information. Please refer @@ -588,6 +588,18 @@ DoFObjectAccessor<3,DH>::get_dof_indices (std::vector &dof_indices for (unsigned int line=0; line<12; ++line) for (unsigned int d=0; dline(line)->dof_index(d,fe_index); + + // now copy dof numbers from the face. for + // faces with the wrong orientation, we + // have already made sure that we're ok by + // picking the correct lines and vertices + // (this happens automatically in the + // line() and vertex() functions. however, + // if the face is in wrong orientation, we + // look at it in flipped orientation and we + // will have to adjust the shape function + // indices that we see to correspond to the + // correct (cell-local) ordering. for (unsigned int quad=0; quad<6; ++quad) if (this->face_orientation(quad)) for (unsigned int d=0; d >::update_cell_dof_indices_cache () const for (unsigned int line=0; line<12; ++line) for (unsigned int d=0; dline(line)->dof_index(d); + + // now copy dof numbers from the face. for + // faces with the wrong orientation, we + // have already made sure that we're ok by + // picking the correct lines and vertices + // (this happens automatically in the + // line() and vertex() functions. however, + // if the face is in wrong orientation, we + // look at it in flipped orientation and we + // will have to adjust the shape function + // indices that we see to correspond to the + // correct (cell-local) ordering. for (unsigned int quad=0; quad<6; ++quad) if (this->face_orientation(quad)) for (unsigned int d=0; d