From 21e2e4b70eade344d611366f325af8169a707e3f Mon Sep 17 00:00:00 2001 From: wolf Date: Thu, 16 Aug 2001 11:31:42 +0000 Subject: [PATCH] Add comment about a case where we had a TODO previously. git-svn-id: https://svn.dealii.org/trunk@4891 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/numerics/matrices.cc | 32 +++++++++++++++++++-- 1 file changed, 29 insertions(+), 3 deletions(-) diff --git a/deal.II/deal.II/source/numerics/matrices.cc b/deal.II/deal.II/source/numerics/matrices.cc index 8d59abe53d..f8a70f1f27 100644 --- a/deal.II/deal.II/source/numerics/matrices.cc +++ b/deal.II/deal.II/source/numerics/matrices.cc @@ -671,9 +671,35 @@ create_boundary_mass_matrix_1 (const Mapping &mapping, cell_matrix(i,j)); else { -//TODO:[?] We assume here that shape functions that have support also on the boundary also have their support point on the boundary (by having their indices in dof_is_on_face set true). This is not true, however, e.g. for discontinuous elements. - // compare here for relative - // smallness + // assume that all + // shape functions + // that are nonzero + // on the boundary + // are also listed + // in the + // @p{dof_to_boundary} + // mapping. if that + // is not the case, + // then the + // boundary mass + // matrix does not + // make that much + // sense anyway, as + // it only contains + // entries for + // parts of the + // functions living + // on the boundary + // + // these, we may + // compare here for + // relative + // smallness of all + // entries in the + // local matrix + // which are not + // taken over to + // the global one Assert (std::fabs(cell_matrix(i,j)) <= 1e-10 * max_diag_entry, ExcInternalError ()); }; -- 2.39.5