r32662 overwrote the changes in r32640, reapplied them

git-svn-id: https://svn.dealii.org/trunk@32672 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 6d5869bb34b091d09cb45ee24fdfed63ae2cb9d1..cef2e584ce61ae297f7d5b754619f7eee3a748e6 100644 (file)
--- a/deal.II/include/deal.II/numerics/vector_tools.templates.h
+++ b/deal.II/include/deal.II/numerics/vector_tools.templates.h
@@ -3153,11 +3153,7 @@ namespace VectorTools
                           + values[q_point] (first_vector_component + 1)
                           * tangentials[q_point] (1))
                        * (fe_values[vec].value (fe.face_to_cell_index (i, face), q_point)
-                          * tangentials[q_point])
-                       / std::sqrt (jacobians[q_point][0][face_coordinate_direction[face]]
-                                    * jacobians[q_point][0][face_coordinate_direction[face]]
-                                    + jacobians[q_point][1][face_coordinate_direction[face]]
-                                    * jacobians[q_point][1][face_coordinate_direction[face]]);
+                          * tangentials[q_point]);
 
                     if (q_point == 0)
                       dofs_processed[i] = true;
@@ -4666,8 +4662,8 @@ namespace VectorTools
         // contributions. Here we do not want to average the normal
         // vectors. Since we have DIM contributions, let's assume (and
         // verify) that they are in fact all linearly independent; in that
-        // case, all vector components are constrained and we need to solve
-        // a linear equation
+        // case, all vector components are constrained and we need to set all 
+        // of them to the corresponding boundary values
         case dim:
         {
           // assert that indeed only a single cell has contributed