From 0725f408c4e1e67c8eb69834b549e360ba792a35 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Thu, 28 Nov 2013 08:33:33 +0000
Subject: [PATCH] Do not use geometric information when setting the Dirichlet
values
[[Imported from SVN: r9546]]
---
cosserat-continuum.cc | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/cosserat-continuum.cc b/cosserat-continuum.cc
index c07af92a..ad93f6f2 100644
--- a/cosserat-continuum.cc
+++ b/cosserat-continuum.cc
@@ -308,7 +308,7 @@ int main (int argc, char *argv[]) try
for (vIt=grid->leafbegin<dim>(); vIt!=vEndIt; ++vIt) {
int idx = grid->leafIndexSet().index(*vIt);
- if (dirichletNodes[idx][0] and vIt->geometry().corner(0)[0] > upper[0]-1e-3) {
+ if (dirichletNodes[idx][0]) {
// Only the positions have Dirichlet values
dirichletValues(vIt->geometry().corner(0), x[idx].r,
--
GitLab