From d2407ffbe6e1433a47e1c3dbd1e319b81f4c19db Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Wed, 19 Mar 2014 16:39:09 +0000
Subject: [PATCH] allow for finer grids in the boundary-handling code
[[Imported from SVN: r9649]]
---
cosserat-continuum.cc | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/cosserat-continuum.cc b/cosserat-continuum.cc
index 70b1a94e..943be144 100644
--- a/cosserat-continuum.cc
+++ b/cosserat-continuum.cc
@@ -186,7 +186,7 @@ int main (int argc, char *argv[]) try
neumannNodes[grid->leafIndexSet().index(*vIt)][0] = true;
#endif
#if 1 // Boundary conditions for the shearing/wrinkling example
- if (vIt->geometry().corner(0)[1] < 1e-3 or vIt->geometry().corner(0)[1] > upper[1]-1e-3 ) {
+ if (vIt->geometry().corner(0)[1] < 1e-4 or vIt->geometry().corner(0)[1] > upper[1]-1e-4 ) {
// Only translation dofs are Dirichlet
for (int j=0; j<3; j++)
dirichletNodes[grid->leafIndexSet().index(*vIt)][j] = true;
--
GitLab