From: kanschat Date: Fri, 19 Oct 2007 06:03:43 +0000 (+0000) Subject: rely more on precollected dofs X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=a3d888c25b3de325e0251b85b9775d123d09d510;p=dealii-svn.git rely more on precollected dofs git-svn-id: https://svn.dealii.org/trunk@15351 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/source/numerics/matrices.cc b/deal.II/deal.II/source/numerics/matrices.cc index da9de4a1b5..327b871742 100644 --- a/deal.II/deal.II/source/numerics/matrices.cc +++ b/deal.II/deal.II/source/numerics/matrices.cc @@ -1036,11 +1036,7 @@ 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) @@ -1198,57 +1194,18 @@ create_boundary_mass_matrix_1 (const Mapping &mapping, cell->get_dof_indices (dofs); cell->face(face)->get_dof_indices (dofs_on_face_vector); - // check for each of the - // dofs on this cell - // whether it is on the - // face - for (unsigned int i=0; i max_diag_entry) - max_diag_entry = std::fabs(cell_matrix(i,i)); -#endif - // lock the matrix Threads::ThreadMutex::ScopedLock lock (mutex); for (unsigned int i=0; i