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

bugfix: honour the minLevel parameter 

[[Imported from SVN: r5682]]
parent e9ce1d2f
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 6 additions and 4 deletions
rodobstacle.cc
+ 6
− 4
View file @ 4da667b4
...@@ -97,8 +97,10 @@ int main (int argc, char *argv[]) try ...@@ -97,8 +97,10 @@ int main (int argc, char *argv[]) try
typedef OneDGrid GridType; typedef OneDGrid GridType;
GridType grid(numRodBaseElements, 0, numRodBaseElements); GridType grid(numRodBaseElements, 0, numRodBaseElements);
std::vector<std::vector<BoxConstraint<double,3> > > trustRegionObstacles(1); grid.globalRefine(minLevel);
std::vector<BitSetVector<1> > hasObstacle(1);
std::vector<std::vector<BoxConstraint<double,3> > > trustRegionObstacles(minLevel+1);
std::vector<BitSetVector<1> > hasObstacle(minLevel+1);
BitSetVector<blocksize> dirichletNodes; BitSetVector<blocksize> dirichletNodes;
// //////////////////////////////// // ////////////////////////////////
...@@ -142,7 +144,7 @@ int main (int argc, char *argv[]) try ...@@ -142,7 +144,7 @@ int main (int argc, char *argv[]) try
double trustRegionRadius = 0.1; double trustRegionRadius = 0.1;
CorrectionType rhs; CorrectionType rhs;
SolutionType x(grid.size(0,1)); SolutionType x(grid.size(1));
CorrectionType corr; CorrectionType corr;
// ////////////////////////// // //////////////////////////
...@@ -160,7 +162,7 @@ int main (int argc, char *argv[]) try ...@@ -160,7 +162,7 @@ int main (int argc, char *argv[]) try
// Refinement Loop // Refinement Loop
// ///////////////////////////////////////////////////////////////////// // /////////////////////////////////////////////////////////////////////
for (int toplevel=0; toplevel<=maxLevel; toplevel++) { for (int toplevel=minLevel; toplevel<=maxLevel; toplevel++) {
std::cout << "####################################################" << std::endl; std::cout << "####################################################" << std::endl;
std::cout << " Solving on level: " << toplevel << std::endl; std::cout << " Solving on level: " << toplevel << std::endl;
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment