diff --git a/cosserat-continuum.cc b/cosserat-continuum.cc index 64c8964617b0bda3490d4aaedaebd676e26c9cf3..e2a36d5709583af31a78e873a66673ee0804d42c 100644 --- a/cosserat-continuum.cc +++ b/cosserat-continuum.cc @@ -151,19 +151,6 @@ int main (int argc, char *argv[]) try // Read Dirichlet values // ///////////////////////////////////////// -#if 0 - BitSetVector<1> allNodes(grid.size(dim)); - allNodes.setAll(); - BoundaryPatch<GridType::LeafGridView> dirichletBoundary(grid.leafView(), allNodes); - - BitSetVector<blocksize> dirichletNodes(grid.size(dim), false); - for (int i=0; i<dirichletNodes.size(); i++) { - // Only translation dofs are Dirichlet - if (dirichletBoundary.containsVertex(i)) - for (int j=0; j<3; j++) - dirichletNodes[i][j] = true; - } -#else BitSetVector<blocksize> dirichletNodes(grid.size(dim), false); BitSetVector<1> neumannNodes(grid.size(dim), false); @@ -181,7 +168,6 @@ int main (int argc, char *argv[]) try } } -#endif ////////////////////////////////////////////////////////////////////////////// // Assemble Neumann term @@ -198,9 +184,6 @@ int main (int argc, char *argv[]) try // Initial solution // ////////////////////////// - FieldVector<double,3> yAxis(0); - yAxis[1] = 1; - vIt = grid.leafbegin<dim>(); for (; vIt!=vEndIt; ++vIt) {