From 9141610436fb88eb19200732bc4b721d2819a35b Mon Sep 17 00:00:00 2001
From: Wolfgang Bangerth <bangerth@math.tamu.edu>
Date: Thu, 10 Mar 2011 00:39:14 +0000
Subject: [PATCH] Also interpolate displacement boundary.

git-svn-id: https://svn.dealii.org/trunk@23475 0785d39b-7218-0410-832d-ea1e28bc413d
---
 deal.II/examples/step-46/step-46.cc | 17 +++++++++++++----
 1 file changed, 13 insertions(+), 4 deletions(-)

diff --git a/deal.II/examples/step-46/step-46.cc b/deal.II/examples/step-46/step-46.cc
index 176fbb7e1e..3bce072653 100644
--- a/deal.II/examples/step-46/step-46.cc
+++ b/deal.II/examples/step-46/step-46.cc
@@ -220,16 +220,25 @@ void StokesProblem<dim>::setup_dofs ()
 
   {
     constraints.clear ();
-    std::vector<bool> component_mask (dim+1+dim, false);
-    for (unsigned int d=0; d<dim; ++d)
-      component_mask[d] = true;
     DoFTools::make_hanging_node_constraints (dof_handler,
 					     constraints);
+
+    std::vector<bool> velocity_mask (dim+1+dim, false);
+    for (unsigned int d=0; d<dim; ++d)
+      velocity_mask[d] = true;
     VectorTools::interpolate_boundary_values (dof_handler,
 					      1,
 					      BoundaryValues<dim>(),
 					      constraints,
-					      component_mask);
+					      velocity_mask);
+    std::vector<bool> elasticity_mask (dim+1+dim, false);
+    for (unsigned int d=dim+1; d<dim+1+dim; ++d)
+      elasticity_mask[d] = true;
+    VectorTools::interpolate_boundary_values (dof_handler,
+					      0,
+					      ZeroFunction<dim>(dim+1+dim),
+					      constraints,
+					      elasticity_mask);
   }
 
 				   // make sure velocity is zero at
-- 
2.39.5