GridGenerator::laplace_transformation function: don't request all boundary 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 046a4042627067e82e31089a589247a921b1749c..a11ad36fe3d77c1c584ee2765804fd98cd4cc606 100644 (file)
--- 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)));
                }
          }
     }