Skip to content
Snippets Groups Projects
Commit 46691921 authored by Elias Pipping's avatar Elias Pipping Committed by Elias Pipping
Browse files

coordinates -> localCoordinates

parent c970acaa
No related branches found
No related tags found
No related merge requests found
......@@ -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
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment