From cc4724ddb28feaedd8e3e1e578ef11bcf0fc2569 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Wed, 17 Oct 2012 15:28:23 +0200
Subject: [PATCH] Drop Neumann boundary
---
src/one-body-sample.cc | 21 +++++++++++++++------
1 file changed, 15 insertions(+), 6 deletions(-)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index a10414ec..055ac9f6 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -114,14 +114,23 @@ void setup_boundary(GridView const &gridView,
assert(false);
}
- if (coordinates[1] == upperRight[1]) // upper face
- ignoreNodes[id] = true;
- else if (coordinates[1] == lowerLeft[1]) { // lower face
+ // lower face
+ if (coordinates[1] == lowerLeft[1]) {
frictionalNodes[id] = true;
ignoreNodes[id][1] = true;
- } else if (coordinates[0] == lowerLeft[0] // left face
- || coordinates[0] == upperRight[0]) // right face
- neumannNodes[id] = true;
+ }
+
+ // upper face
+ else if (coordinates[1] == upperRight[1])
+ ignoreNodes[id] = true;
+
+ // right face (except for both corners)
+ else if (coordinates[0] == upperRight[0])
+ ;
+
+ // left face (except for both corners)
+ else if (coordinates[0] == lowerLeft[0])
+ ;
}
}
--
GitLab