From: wolf Date: Tue, 11 Jun 2002 08:58:44 +0000 (+0000) Subject: Yeah, fix that bug. Constraint matrices are annoying!!! The TODO was X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=058f96eb16d5e70d12675b3dbe8b71d53bd29c57;p=dealii-svn.git Yeah, fix that bug. Constraint matrices are annoying!!! The TODO was long done, by the way. git-svn-id: https://svn.dealii.org/trunk@6054 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/source/fe/fe_system.cc b/deal.II/deal.II/source/fe/fe_system.cc index 76614cb894..9551529c03 100644 --- a/deal.II/deal.II/source/fe/fe_system.cc +++ b/deal.II/deal.II/source/fe/fe_system.cc @@ -1173,8 +1173,6 @@ FESystem::build_face_tables() template void FESystem::build_interface_constraints () { -//TODO: Review for base elements with more than one component. do it in the same way as for building the cell restriction and interpolation tables, i.e. top-down - // the layout of the constraints // matrix is described in the // FiniteElement class. you may @@ -1356,7 +1354,7 @@ void FESystem::build_interface_constraints () 4*base_element(m_index.first.first).dofs_per_vertex; Assert (tmp2 < base_element(m_index.first.first).dofs_per_line, ExcInternalError()); - m_index.second = base_element(m_index.first.first).dofs_per_vertex + + m_index.second = 5*base_element(m_index.first.first).dofs_per_vertex + base_element(m_index.first.first).dofs_per_line*sub_line + tmp2; } @@ -1385,7 +1383,7 @@ void FESystem::build_interface_constraints () 4*base_element(m_index.first.first).dofs_per_line; Assert (tmp2 < base_element(m_index.first.first).dofs_per_quad, ExcInternalError()); - m_index.second = base_element(m_index.first.first).dofs_per_vertex + + m_index.second = 5*base_element(m_index.first.first).dofs_per_vertex + 12*base_element(m_index.first.first).dofs_per_line + base_element(m_index.first.first).dofs_per_quad*sub_quad + tmp2;