From: hartmann <hartmann@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Fri, 15 Oct 2004 07:07:37 +0000 (+0000)
Subject: GridGenerator::laplace_transformation function: don't request all boundary points... 
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=08410bad08e7fb601e0cf224943fdd40ae02e1a9;p=dealii-svn.git

GridGenerator::laplace_transformation function: don't request all boundary points being included in new_points.


git-svn-id: https://svn.dealii.org/trunk@9712 0785d39b-7218-0410-832d-ea1e28bc413d
---

diff --git a/deal.II/deal.II/source/grid/grid_generator.cc b/deal.II/deal.II/source/grid/grid_generator.cc
index 046a404262..a11ad36fe3 100644
--- a/deal.II/deal.II/source/grid/grid_generator.cc
+++ b/deal.II/deal.II/source/grid/grid_generator.cc
@@ -1311,11 +1311,11 @@ void GridGenerator::laplace_transformation (Triangulation<dim> &tria,
 		{
 		  const unsigned int vertex_index=face->vertex_index(vertex_no);
 		  map_iter=new_points.find(vertex_index);
-		  Assert(map_iter!=map_end, ExcInternalError());
-		  
-		  for (unsigned int i=0; i<dim; ++i)
-		    m[i].insert(std::pair<unsigned int,double> (
-				  face->vertex_dof_index(vertex_no, 0), map_iter->second(i)));
+
+		  if (map_iter!=map_end)
+		    for (unsigned int i=0; i<dim; ++i)
+		      m[i].insert(std::pair<unsigned int,double> (
+				    face->vertex_dof_index(vertex_no, 0), map_iter->second(i)));
 		}
 	  }
     }