Skip to content
Snippets Groups Projects
Commit d2407ffb authored by Oliver Sander's avatar Oliver Sander Committed by sander
Browse files

allow for finer grids in the boundary-handling code

[[Imported from SVN: r9649]]
parent 155dcaa3
No related branches found
No related tags found
No related merge requests found
...@@ -186,7 +186,7 @@ int main (int argc, char *argv[]) try ...@@ -186,7 +186,7 @@ int main (int argc, char *argv[]) try
neumannNodes[grid->leafIndexSet().index(*vIt)][0] = true; neumannNodes[grid->leafIndexSet().index(*vIt)][0] = true;
#endif #endif
#if 1 // Boundary conditions for the shearing/wrinkling example #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 // Only translation dofs are Dirichlet
for (int j=0; j<3; j++) for (int j=0; j<3; j++)
dirichletNodes[grid->leafIndexSet().index(*vIt)][j] = true; dirichletNodes[grid->leafIndexSet().index(*vIt)][j] = true;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment