From: Wolfgang Bangerth <bangerth@math.tamu.edu>
Date: Mon, 22 Jul 2013 12:56:52 +0000 (+0000)
Subject: Avoid filling a map with numeric_limit<types::boundary_id>::max() elements -- too... 
X-Git-Tag: v8.0.0~29
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=dee9589ea156718b645e3b4a16531453a37035c0;p=dealii.git

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
---

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);
           }