diff --git a/cosserat-continuum.cc b/cosserat-continuum.cc index 70b1a94e930b67edc05eabebf6f2c57dca86fb46..943be144bfc104a377124bb22f9a0ad288b53dbf 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;