From: guido Date: Tue, 27 Mar 2001 18:05:32 +0000 (+0000) Subject: Comments added X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=9dd025bae6bf561f182474a8426fa2c17157f542;p=dealii-svn.git Comments added git-svn-id: https://svn.dealii.org/trunk@4313 0785d39b-7218-0410-832d-ea1e28bc413d --- 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)