From 509f13cb7d997c5f80116e3417d570a68b7af56c Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Fri, 22 Aug 2014 10:14:56 +0000
Subject: [PATCH] Minor simplification
[[Imported from SVN: r9854]]
---
dune/gfe/riemanniantrsolver.cc | 8 +++-----
1 file changed, 3 insertions(+), 5 deletions(-)
diff --git a/dune/gfe/riemanniantrsolver.cc b/dune/gfe/riemanniantrsolver.cc
index 130704e6..532b0361 100644
--- a/dune/gfe/riemanniantrsolver.cc
+++ b/dune/gfe/riemanniantrsolver.cc
@@ -250,9 +250,7 @@ void RiemannianTrustRegionSolver<GridType,TargetSpace>::solve()
MaxNormTrustRegion<blocksize> trustRegion(guIndex_->nGlobalEntity(), initialTrustRegionRadius_);
- std::vector<std::vector<BoxConstraint<field_type,blocksize> > > trustRegionObstacles((mgStep)
- ? mgStep->numLevels()
- : 0);
+ std::vector<BoxConstraint<field_type,blocksize> > trustRegionObstacles;
// /////////////////////////////////////////////////////
// Set up the log file, if requested
@@ -331,8 +329,8 @@ void RiemannianTrustRegionSolver<GridType,TargetSpace>::solve()
{
mgStep->setProblem(stiffnessMatrix, corr_global, rhs_global);
- trustRegionObstacles.back() = trustRegion.obstacles();
- mgStep->obstacles_ = &trustRegionObstacles.back();
+ trustRegionObstacles = trustRegion.obstacles();
+ mgStep->obstacles_ = &trustRegionObstacles;
innerSolver_->preprocess();
--
GitLab