From 0a24d717c50086bb5384e08da732f1b4ccc9c2a9 Mon Sep 17 00:00:00 2001 From: bangerth Date: Fri, 19 Jan 2007 20:26:13 +0000 Subject: [PATCH] Fix a bug: the array size was already nonzero (the base class constructor sets it correctly, so we only have to copy data, not insert it. git-svn-id: https://svn.dealii.org/trunk@14349 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/fe/fe_system.cc | 20 ++++++++++++++------ 1 file changed, 14 insertions(+), 6 deletions(-) diff --git a/deal.II/deal.II/source/fe/fe_system.cc b/deal.II/deal.II/source/fe/fe_system.cc index b6db8a9000..ef70cb174b 100644 --- a/deal.II/deal.II/source/fe/fe_system.cc +++ b/deal.II/deal.II/source/fe/fe_system.cc @@ -2,7 +2,7 @@ // $Id$ // Version: $Name$ // -// Copyright (C) 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006 by the deal.II authors +// Copyright (C) 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007 by the deal.II authors // // This file is subject to QPL and may not be distributed // without copyright and license information. Please refer @@ -1864,18 +1864,26 @@ template <> void FESystem<3>::initialize_quad_dof_index_permutation () { + // the array into which we want to write + // should have the correct size already. + Assert (adjust_quad_dof_index_for_face_orientation_table.size()==this->dofs_per_quad, + ExcInternalError()); // to obtain the shifts for this composed // element, concatenate the shift vectors of // the base elements + unsigned int index = 0; for (unsigned int b=0; b &temp=this->base_element(b).adjust_quad_dof_index_for_face_orientation_table; + const std::vector &temp + = this->base_element(b).adjust_quad_dof_index_for_face_orientation_table; for (unsigned int c=0; cdofs_per_quad, + Assert (index == this->dofs_per_quad, ExcInternalError()); } -- 2.39.5