From b6c21607af51c400911aea00f12fcb432a43cd35 Mon Sep 17 00:00:00 2001 From: Guido Kanschat Date: Tue, 27 Mar 2001 18:05:32 +0000 Subject: [PATCH] Comments added git-svn-id: https://svn.dealii.org/trunk@4313 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/fe/mapping_q1.cc | 18 +++++++++++++----- 1 file changed, 13 insertions(+), 5 deletions(-) diff --git a/deal.II/deal.II/source/fe/mapping_q1.cc b/deal.II/deal.II/source/fe/mapping_q1.cc index d877681a7d..c4c3118356 100644 --- a/deal.II/deal.II/source/fe/mapping_q1.cc +++ b/deal.II/deal.II/source/fe/mapping_q1.cc @@ -287,7 +287,18 @@ MappingQ1::update_each (const UpdateFlags in) const | update_normal_vectors)); // cerr << "Mapping-each " << hex << in << ' ' << out; - + + // The following is a little incorrect: + // If not applied on a face, + // update_boundary_forms does not + // make sense. On the other hand, + // it is necessary on a + // face. Currently, + // update_boundary_forms is simply + // ignored for the interior of a + // cell. +//TODO: Consider giving this function information on whether we are on +// a face. if (out & (update_JxW_values |update_normal_vectors)) out |= update_boundary_forms; @@ -353,10 +364,7 @@ MappingQ1::compute_face_data (UpdateFlags update_flags, const unsigned int n_original_q_points, InternalData& data) const { - if (update_flags & update_JxW_values) - update_flags |= update_boundary_forms; - - compute_data (update_flags, q, n_original_q_points, data); + compute_data (update_flags, q, n_original_q_points, data); #if (deal_II_dimension>1) if (data.update_flags & update_boundary_forms) -- 2.39.5