diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc index d9148a1b54a79dcc86598ea890fa9022cdca8a00..feb5c59f4b9b84ac40322d49be941fdac9c28b47 100644 --- a/src/one-body-sample.cc +++ b/src/one-body-sample.cc @@ -95,13 +95,13 @@ void setup_boundary(GridView const &gridView, assert(it->geometry().corners() == 1); Dune::FieldVector<double, dim> const coordinates = it->geometry().corner(0); size_t const id = myVertexMapper.map(*it); - if (coordinates[1] == upperRight[1]) + if (coordinates[1] == upperRight[1]) // upper face ignoreNodes[id] = true; - else if (coordinates[1] == lowerLeft[1]) { + else if (coordinates[1] == lowerLeft[1]) { // lower face frictionalNodes[id] = true; - ignoreNodes[id][1] = true; // Zero displacement in direction y - } else if (coordinates[0] == lowerLeft[0] || - coordinates[0] == upperRight[0]) + ignoreNodes[id][1] = true; + } else if (coordinates[0] == lowerLeft[0] // left face + || coordinates[0] == upperRight[0]) // right face neumannNodes[id] = true; } }