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

Minor simplification

[[Imported from SVN: r9854]]
parent de336a2a
No related branches found
No related tags found
No related merge requests found
...@@ -250,9 +250,7 @@ void RiemannianTrustRegionSolver<GridType,TargetSpace>::solve() ...@@ -250,9 +250,7 @@ void RiemannianTrustRegionSolver<GridType,TargetSpace>::solve()
MaxNormTrustRegion<blocksize> trustRegion(guIndex_->nGlobalEntity(), initialTrustRegionRadius_); MaxNormTrustRegion<blocksize> trustRegion(guIndex_->nGlobalEntity(), initialTrustRegionRadius_);
std::vector<std::vector<BoxConstraint<field_type,blocksize> > > trustRegionObstacles((mgStep) std::vector<BoxConstraint<field_type,blocksize> > trustRegionObstacles;
? mgStep->numLevels()
: 0);
// ///////////////////////////////////////////////////// // /////////////////////////////////////////////////////
// Set up the log file, if requested // Set up the log file, if requested
...@@ -331,8 +329,8 @@ void RiemannianTrustRegionSolver<GridType,TargetSpace>::solve() ...@@ -331,8 +329,8 @@ void RiemannianTrustRegionSolver<GridType,TargetSpace>::solve()
{ {
mgStep->setProblem(stiffnessMatrix, corr_global, rhs_global); mgStep->setProblem(stiffnessMatrix, corr_global, rhs_global);
trustRegionObstacles.back() = trustRegion.obstacles(); trustRegionObstacles = trustRegion.obstacles();
mgStep->obstacles_ = &trustRegionObstacles.back(); mgStep->obstacles_ = &trustRegionObstacles;
innerSolver_->preprocess(); innerSolver_->preprocess();
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment