Skip to content
Snippets Groups Projects
Commit 71f6cc1b authored by Oliver Sander's avatar Oliver Sander Committed by sander@FU-BERLIN.DE
Browse files

Bugfix for the case a dof is Dirichlet in a few directions, and Neumann in others

[[Imported from SVN: r7428]]
parent 88fab69e
No related branches found
No related tags found
No related merge requests found
......@@ -226,10 +226,11 @@ void RiemannianTrustRegionSolver<GridType,TargetSpace>::solve()
for (; cIt!=cEndIt; ++cIt) {
for (int k=0; k<blocksize; k++) {
if (ignoreNodes_->operator[](j)[k])
if (ignoreNodes_->operator[](j)[k]) {
(*cIt)[k] = 0;
if (j==cIt.index())
(*cIt)[k][k] = 1;
if (j==cIt.index())
(*cIt)[k][k] = 1;
}
}
}
......
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