From: hartmann <hartmann@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Fri, 15 Oct 2004 06:59:10 +0000 (+0000)
Subject: Bug fix: the StraightBoundary<3>::get_new_point_on_quad requested for children of... 
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=6fe827a0bdbfd196882da94a2f95a59ab7644c4e;p=dealii-svn.git

Bug fix: the StraightBoundary<3>::get_new_point_on_quad requested for children of lines of the quad which in general do not exist. Get around this by using the centers of the lines.


git-svn-id: https://svn.dealii.org/trunk@9709 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 48e306705c..9f49f182da 100644
--- a/deal.II/deal.II/source/grid/tria_boundary.cc
+++ b/deal.II/deal.II/source/grid/tria_boundary.cc
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$
 //
-//    Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003 by the deal.II authors
+//    Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003, 2004 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -111,10 +111,10 @@ get_new_point_on_quad (const Triangulation<3>::quad_iterator &quad) const
 {
   return (quad->vertex(0) + quad->vertex(1) +
 	  quad->vertex(2) + quad->vertex(3) +
-	  quad->line(0)->child(0)->vertex(1) +
-	  quad->line(1)->child(0)->vertex(1) +
-  	  quad->line(2)->child(0)->vertex(1) +
-  	  quad->line(3)->child(0)->vertex(1)) / 8;
+	  quad->line(0)->center() +
+	  quad->line(1)->center() +
+  	  quad->line(2)->center() +
+  	  quad->line(3)->center()) / 8;
 }
 
 #endif