]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Yeah, fix that bug. Constraint matrices are annoying!!! The TODO was
authorwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 11 Jun 2002 08:58:44 +0000 (08:58 +0000)
committerwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 11 Jun 2002 08:58:44 +0000 (08:58 +0000)
long done, by the way.

git-svn-id: https://svn.dealii.org/trunk@6054 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/fe/fe_system.cc

index 76614cb8941a4458d6dd93922f70409b6b2f2853..9551529c03f55b09bba4b4e1c22c518bbf80545a 100644 (file)
@@ -1173,8 +1173,6 @@ FESystem<dim>::build_face_tables()
 template <int dim>
 void FESystem<dim>::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<dim>::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<dim>::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;

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.