diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc index 464b6d0e1f540df2f2fb4762efc8fe2571d8c590..3071ba3500a719aad015bb46a4e3865d56a50503 100644 --- a/src/one-body-sample.cc +++ b/src/one-body-sample.cc @@ -198,19 +198,19 @@ int main(int argc, char *argv[]) { for (auto it = leafView.begin<dims>(); it != leafView.end<dims>(); ++it) { assert(it->geometry().corners() == 1); - Dune::FieldVector<double, dims> const coordinates = + Dune::FieldVector<double, dims> const localCoordinates = it->geometry().corner(0); size_t const id = myVertexMapper.map(*it); // Find the center of the lower face switch (dims) { case 3: - if (coordinates[2] != lowerLeft[2]) + if (localCoordinates[2] != lowerLeft[2]) break; // fallthrough case 2: - if (coordinates[1] == lowerLeft[1] && - std::abs(coordinates[0] - lowerLeft[0]) < 1e-8) + if (localCoordinates[1] == lowerLeft[1] && + std::abs(localCoordinates[0] - lowerLeft[0]) < 1e-8) specialNode = id; break; default: @@ -218,21 +218,21 @@ int main(int argc, char *argv[]) { } // lower face - if (coordinates[1] == lowerLeft[1]) { + if (localCoordinates[1] == lowerLeft[1]) { frictionalNodes[id] = true; ignoreNodes[id][1] = true; } // upper face - else if (coordinates[1] == upperRight[1]) + else if (localCoordinates[1] == upperRight[1]) ignoreNodes[id] = true; // right face (except for both corners) - else if (coordinates[0] == upperRight[0]) + else if (localCoordinates[0] == upperRight[0]) ; // left face (except for both corners) - else if (coordinates[0] == lowerLeft[0]) + else if (localCoordinates[0] == lowerLeft[0]) ; } // Make sure that specialNode was set and points to a frictional node