From: bangerth Date: Fri, 4 Apr 2008 17:24:33 +0000 (+0000) Subject: Don't try to unify dofs on lines when there are none. This happens for example for... X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=88cebd5ecb547a50bc57469612b09a8ba7d46026;p=dealii-svn.git Don't try to unify dofs on lines when there are none. This happens for example for Raviart-Thomas elements in 3d. git-svn-id: https://svn.dealii.org/trunk@15984 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/source/hp/dof_handler.cc b/deal.II/deal.II/source/hp/dof_handler.cc index 081c9d397d..381747f7f6 100644 --- a/deal.II/deal.II/source/hp/dof_handler.cc +++ b/deal.II/deal.II/source/hp/dof_handler.cc @@ -2,7 +2,7 @@ // $Id$ // Version: $Name$ // -// Copyright (C) 2003, 2004, 2006, 2007 by the deal.II authors +// Copyright (C) 2003, 2004, 2006, 2007, 2008 by the deal.II authors // // This file is subject to QPL and may not be distributed // without copyright and license information. Please refer @@ -2378,9 +2378,11 @@ namespace hp const unsigned int fe_index_1 = line->nth_active_fe_index (f), fe_index_2 = line->nth_active_fe_index (g); - if ((*finite_elements)[fe_index_1].dofs_per_line - == - (*finite_elements)[fe_index_2].dofs_per_line) + if (((*finite_elements)[fe_index_1].dofs_per_line + == + (*finite_elements)[fe_index_2].dofs_per_line) + && + ((*finite_elements)[fe_index_1].dofs_per_line > 0)) { internal::hp::ensure_existence_of_dof_identities ((*finite_elements)[fe_index_1],