From: hartmann 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