From: bangerth Date: Thu, 27 Jul 2006 03:58:56 +0000 (+0000) Subject: Fix a few thinkos in the code. Remove debug output X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=936b64abf090467553f4f9565fbed4b603808d28;p=dealii-svn.git Fix a few thinkos in the code. Remove debug output git-svn-id: https://svn.dealii.org/trunk@13441 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/source/dofs/hp_dof_handler.cc b/deal.II/deal.II/source/dofs/hp_dof_handler.cc index df5bcc7b82..a322dbfe49 100644 --- a/deal.II/deal.II/source/dofs/hp_dof_handler.cc +++ b/deal.II/deal.II/source/dofs/hp_dof_handler.cc @@ -2036,7 +2036,8 @@ namespace hp // if the finite element tells us // that they should have the same // value - std::vector new_dof_indices (used_dofs); + std::vector new_dof_indices (used_dofs, + deal_II_numbers::invalid_unsigned_int); // Step 2a: do so for vertices { @@ -2044,7 +2045,8 @@ namespace hp std::vector > Identities; Table<2,boost::shared_ptr > - vertex_dof_identities; + vertex_dof_identities (finite_elements->size(), + finite_elements->size()); // loop over all vertices and // see which one we need to @@ -2094,32 +2096,41 @@ namespace hp // constrained to // anything else, // except for to - // each other + // each other. use + // the rule that we + // will always + // constrain the + // dof with the + // higher fe + // index to the + // one with the + // lower, to avoid + // circular + // reasoning. Identities &identities = *vertex_dof_identities[first_fe_index][other_fe_index]; for (unsigned int i=0; i