From 043be2e48c4721eb9909a66ad8d0a6b57ba457fb Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Wed, 20 Oct 2004 12:52:29 +0000 Subject: [PATCH] Extend a comment. git-svn-id: https://svn.dealii.org/trunk@9714 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/grid/tria_boundary.cc | 15 +++++++++++++++ 1 file changed, 15 insertions(+) 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() ? -- 2.39.5