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

bugfix: sample dirichlet values on higher levels

[[Imported from SVN: r9853]]
parent 4184403b
No related branches found
No related tags found
No related merge requests found
......@@ -12,6 +12,7 @@
#include <dune/istl/io.hh>
#include <dune/ag-common/prolongboundarypatch.hh>
#include <dune/ag-common/sampleonbitfield.hh>
#include <dune/ag-common/readbitfield.hh>
#include <dune/ag-common/quadraticipopt.hh>
#include <dune/ag-common/iterativesolver.hh>
......@@ -91,6 +92,8 @@ int main (int argc, char *argv[]) try
}
sampleOnBitField(grid, dirichletValues, dirichletNodes);
int maxlevel = grid.maxLevel();
// Determine Dirichlet dofs
......@@ -107,6 +110,7 @@ int main (int argc, char *argv[]) try
#if defined IPOPT
QuadraticIPOptSolver<OperatorType,VectorType> solver;
solver.verbosity_ = NumProc::QUIET;
solver.dirichletNodes_ = &dirichletNodes[maxlevel];
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment