From dee9589ea156718b645e3b4a16531453a37035c0 Mon Sep 17 00:00:00 2001
From: Wolfgang Bangerth <bangerth@math.tamu.edu>
Date: Mon, 22 Jul 2013 12:56:52 +0000
Subject: [PATCH] Avoid filling a map with
 numeric_limit<types::boundary_id>::max() elements -- too many if we ever
 increase the number of bits in boundary_id.

git-svn-id: https://svn.dealii.org/trunk@30093 0785d39b-7218-0410-832d-ea1e28bc413d
---
 deal.II/include/deal.II/numerics/vector_tools.templates.h | 7 +++++--
 1 file changed, 5 insertions(+), 2 deletions(-)

diff --git a/deal.II/include/deal.II/numerics/vector_tools.templates.h b/deal.II/include/deal.II/numerics/vector_tools.templates.h
index abb0372b0b..ace4f136ab 100644
--- a/deal.II/include/deal.II/numerics/vector_tools.templates.h
+++ b/deal.II/include/deal.II/numerics/vector_tools.templates.h
@@ -535,9 +535,12 @@ namespace VectorTools
             // parts. We want the
             // function to hold on
             // all parts of the boundary
+            const std::vector<types::boundary_id>
+            used_boundary_indicators = dof.get_tria().get_boundary_indicators();
+
             typename FunctionMap<spacedim>::type boundary_functions;
-            for (types::boundary_id c=0; c<numbers::internal_face_boundary_id; ++c)
-              boundary_functions[c] = &function;
+            for (unsigned int i=0; i<used_boundary_indicators.size(); ++i)
+              boundary_functions[used_boundary_indicators[i]] = &function;
             project_boundary_values (mapping, dof, boundary_functions, q_boundary,
                                      boundary_values);
           }
-- 
2.39.5