From 999dc3d3959068e0822fd1cc9989f4ee3e6e2725 Mon Sep 17 00:00:00 2001 From: bangerth Date: Fri, 19 Jan 2007 20:28:50 +0000 Subject: [PATCH] Add a comment git-svn-id: https://svn.dealii.org/trunk@14350 0785d39b-7218-0410-832d-ea1e28bc413d --- .../deal.II/include/dofs/dof_accessor.templates.h | 14 +++++++++++++- deal.II/deal.II/source/dofs/dof_accessor.cc | 14 +++++++++++++- 2 files changed, 26 insertions(+), 2 deletions(-) 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