From 2ad1be44599cb539ff378c381fee9d3b9fbb4bf8 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Tue, 10 Jan 2012 17:01:12 +0100
Subject: [PATCH] Do not pass end_points to setup_boundary
---
src/one-body-sample.cc | 8 +++-----
1 file changed, 3 insertions(+), 5 deletions(-)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 0cc55e51..8002b890 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -58,7 +58,6 @@ int const dim = 3;
template <class GridView>
void setup_boundary(GridView const &gridView,
- Dune::FieldVector<double, dim> const &end_points,
Dune::BitSetVector<dim> &ignoreNodes,
Dune::BitSetVector<1> &neumannNodes,
Dune::BitSetVector<1> &frictionalNodes) {
@@ -75,7 +74,7 @@ void setup_boundary(GridView const &gridView,
it != gridView.template end<dim>(); ++it) {
assert(it->geometry().corners() == 1);
Dune::FieldVector<double, dim> const coordinates = it->geometry().corner(0);
- if (coordinates[1] == end_points[0]) {
+ if (coordinates[1] == 1) {
++dirichlet_nodes;
size_t const id = myVertexMapper.map(*it);
ignoreNodes[id] = true;
@@ -84,7 +83,7 @@ void setup_boundary(GridView const &gridView,
size_t const id = myVertexMapper.map(*it);
frictionalNodes[id] = true;
ignoreNodes[id][1] = true; // Zero displacement in direction y
- } else if (coordinates[0] == 0 || coordinates[0] == end_points[0]) {
+ } else if (coordinates[0] == 0 || coordinates[0] == 1) {
++neumann_nodes;
size_t const id = myVertexMapper.map(*it);
neumannNodes[id] = true;
@@ -238,8 +237,7 @@ int main(int argc, char *argv[]) {
false);
Dune::BitSetVector<1> frictionalNodes(grid->size(grid->maxLevel(), dim),
false);
- setup_boundary(leafView, end_points, ignoreNodes, neumannNodes,
- frictionalNodes);
+ setup_boundary(leafView, ignoreNodes, neumannNodes, frictionalNodes);
typedef MyConvexProblem<OperatorType, VectorType> MyConvexProblemType;
typedef MyBlockProblem<MyConvexProblemType> MyBlockProblemType;
--
GitLab