From: Wolfgang Bangerth Date: Wed, 20 Oct 2004 12:52:29 +0000 (+0000) Subject: Extend a comment. X-Git-Tag: v8.0.0~14758 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=043be2e48c4721eb9909a66ad8d0a6b57ba457fb;p=dealii.git Extend a comment. git-svn-id: https://svn.dealii.org/trunk@9714 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/source/grid/tria_boundary.cc b/deal.II/deal.II/source/grid/tria_boundary.cc index a0e35d5fb2..37cd46bbc1 100644 --- a/deal.II/deal.II/source/grid/tria_boundary.cc +++ b/deal.II/deal.II/source/grid/tria_boundary.cc @@ -125,6 +125,21 @@ get_new_point_on_quad (const Triangulation<3>::quad_iterator &quad) const // situation by using existing edge // midpoints if available, or center() if // not + // + // note that this situation can not happen + // during mesh refinement, as there the + // edges are refined first and only then + // the face. thus, the check whether a line + // has children does not lead to the + // situation where the new face midpoints + // have different positions depending on + // which of the two cells is refined first. + // + // the situation where the edges aren't + // refined happens when the a higher order + // MappingQ requests the midpoint of a + // face, though, and it is for these cases + // that we need to have the check available return (quad->vertex(0) + quad->vertex(1) + quad->vertex(2) + quad->vertex(3) + (quad->line(0)->has_children() ?