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

Do not use geometric information when setting the Dirichlet values

[[Imported from SVN: r9546]]
parent 0c4778a5
Branches
Tags
No related merge requests found
...@@ -308,7 +308,7 @@ int main (int argc, char *argv[]) try ...@@ -308,7 +308,7 @@ int main (int argc, char *argv[]) try
for (vIt=grid->leafbegin<dim>(); vIt!=vEndIt; ++vIt) { for (vIt=grid->leafbegin<dim>(); vIt!=vEndIt; ++vIt) {
int idx = grid->leafIndexSet().index(*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 // Only the positions have Dirichlet values
dirichletValues(vIt->geometry().corner(0), x[idx].r, dirichletValues(vIt->geometry().corner(0), x[idx].r,
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment