]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Remove some code that didn't do anything.
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Thu, 5 Nov 2009 23:59:44 +0000 (23:59 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Thu, 5 Nov 2009 23:59:44 +0000 (23:59 +0000)
git-svn-id: https://svn.dealii.org/trunk@20051 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/dofs/dof_tools.cc

index 85d322ec01ccdd60693403b4b3c79d7f20b102a3..d7e53c31442ed0d2f92d411029f80be636a0d485 100644 (file)
@@ -5401,89 +5401,12 @@ DoFTools::make_zero_boundary_constraints (const DH<dim, spacedim> &dof,
            face_dofs.resize (fe.dofs_per_face);
            face->get_dof_indices (face_dofs, cell->active_fe_index());
 
-           if (fe_is_system)
-             {
-                                   // enter those dofs into the list
-                                   // that match the component
-                                   // signature.
-               for (unsigned int i=0; i<face_dofs.size(); ++i)
-                  {
-                    unsigned int component;
-                    if (fe.is_primitive())
-                      component = fe.face_system_to_component_index(i).first;
-                    else
-                      {
-                                    // non-primitive case. make sure
-                                    // that this particular shape
-                                    // function _is_ primitive, and
-                                    // get at it's component. use
-                                    // usual trick to transfer face
-                                    // dof index to cell dof index
-                        const unsigned int cell_i
-                          = (dim == 1 ?
-                             i
-                             :
-                             (dim == 2 ?
-                              (i<2*fe.dofs_per_vertex ? i : i+2*fe.dofs_per_vertex)
-                              :
-                              (dim == 3 ?
-                               (i<4*fe.dofs_per_vertex ?
-                                i
-                                :
-                                (i<4*fe.dofs_per_vertex+4*fe.dofs_per_line ?
-                                 i+4*fe.dofs_per_vertex
-                                 :
-                                 i+4*fe.dofs_per_vertex+8*fe.dofs_per_line))
-                               :
-                               numbers::invalid_unsigned_int)));
-                        Assert (cell_i < fe.dofs_per_cell, ExcInternalError());
-
-                                    // make sure that if this is not a
-                                    // primitive shape function, then
-                                    // all the corresponding
-                                    // components in the mask are not
-                                    // set
-                        if (!fe.is_primitive(cell_i))
-                          for (unsigned int c=0; c<n_components; ++c)
-                            if (fe.get_nonzero_components(cell_i)[c])
-                              Assert (component_mask[c] == false,
-                                      FETools::ExcFENotPrimitive());
-
-                                    // pick the first possibly of more
-                                    // than one non-zero component. if
-                                    // the shape function is
-                                    // non-primitive, then we will
-                                    // ignore the result in the
-                                    // following anyway, otherwise
-                                    // there's only one non-zero
-                                    // component which we will use
-                        component = (std::find (fe.get_nonzero_components(cell_i).begin(),
-                                                fe.get_nonzero_components(cell_i).end(),
-                                                true)
-                                     -
-                                     fe.get_nonzero_components(cell_i).begin());
-                      }
-
-                                   // cast zero boundary constraints onto
-                                   // a matrix
-                   for (unsigned int i=0; i<fe.dofs_per_vertex; ++i)
-                     if (component_mask[fe.face_system_to_component_index(i).first])
-                       zero_boundary_constraints.add_line (face_dofs[i]);
-
-                  }
-             }
-           else
-             {
-                                   // get the one component that this
-                                   // function has and cast zero
-                                   // boundary constraints onto a
-                                   // matrix
-               for (unsigned int i=0; i<fe.dofs_per_face; ++i)
-                 if (component_mask[fe.face_system_to_component_index(i).first])
-                   zero_boundary_constraints.add_line (face_dofs[i]);
-
-             }
-
+                                            // enter those dofs into the list
+                                            // that match the component
+                                            // signature.
+           for (unsigned int i=0; i<face_dofs.size(); ++i)
+             if (component_mask[fe.face_system_to_component_index(i).first])
+               zero_boundary_constraints.add_line (face_dofs[i]);
          }
       }
 }

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.