From 0edd0a4b07233271d5b9fde08f7199ce9095e61f Mon Sep 17 00:00:00 2001 From: kanschat Date: Wed, 24 Oct 2007 19:35:25 +0000 Subject: [PATCH] reenter test for dof on THIS face git-svn-id: https://svn.dealii.org/trunk@15375 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/numerics/matrices.cc | 34 ++++++++++++++------- 1 file changed, 23 insertions(+), 11 deletions(-) diff --git a/deal.II/deal.II/source/numerics/matrices.cc b/deal.II/deal.II/source/numerics/matrices.cc index 327b871742..02a144c014 100644 --- a/deal.II/deal.II/source/numerics/matrices.cc +++ b/deal.II/deal.II/source/numerics/matrices.cc @@ -1037,6 +1037,10 @@ create_boundary_mass_matrix_1 (const Mapping &mapping, std::vector dofs (dofs_per_cell); std::vector dofs_on_face_vector (dofs_per_face); + // for each dof on the cell, have a + // flag whether it is on the face + std::vector dof_is_on_face(dofs_per_cell); + typename DoFHandler::active_cell_iterator cell = range.first; for (; cell!=range.second; ++cell) for (unsigned int face=0; face::faces_per_cell; ++face) @@ -1085,7 +1089,7 @@ create_boundary_mass_matrix_1 (const Mapping &mapping, for (unsigned int point=0; point &mapping, // searches. cell->get_dof_indices (dofs); cell->face(face)->get_dof_indices (dofs_on_face_vector); + for (unsigned int i=0; i