From a8482112fa77d5326863a01a2b78b6a6630d7e13 Mon Sep 17 00:00:00 2001 From: "daniel.arndt" Date: Thu, 20 Mar 2014 09:45:16 +0000 Subject: [PATCH] r32662 overwrote the changes in r32640, reapplied them git-svn-id: https://svn.dealii.org/trunk@32672 0785d39b-7218-0410-832d-ea1e28bc413d --- .../include/deal.II/numerics/vector_tools.templates.h | 10 +++------- 1 file changed, 3 insertions(+), 7 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 6d5869bb34..cef2e584ce 100644 --- 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 -- 2.39.5