Skip to content
Snippets Groups Projects

Use PARMG as the inner solver in Trust Region

Closed lisa_julia.nebel_at_tu-dresden.de requested to merge (removed):master into master
1 file
+ 1
1
Compare changes
  • Side-by-side
  • Inline
@@ -248,11 +248,11 @@ void TrustRegionSolver<BasisType,VectorType>::solve()
@@ -248,11 +248,11 @@ void TrustRegionSolver<BasisType,VectorType>::solve()
BasisType basis(grid_->levelGridView(grid_->maxLevel()));
BasisType basis(grid_->levelGridView(grid_->maxLevel()));
BasisType coarseBasis(grid_->levelGridView(0));
BasisType coarseBasis(grid_->levelGridView(0));
std::vector<BoxConstraint<typename VectorType::field_type, blocksize>> coarseTrustRegionObstacles(coarseBasis.size());
std::vector<BoxConstraint<typename VectorType::field_type, blocksize>> coarseTrustRegionObstacles(coarseBasis.size());
MaxNormTrustRegion<blocksize> trustRegion(basis.size(), initialTrustRegionRadius_);
int numLevels = grid_->maxLevel()+1;
int numLevels = grid_->maxLevel()+1;
auto& levelOp = mgSetup_->levelOps_;
auto& levelOp = mgSetup_->levelOps_;
#endif
#endif
 
MaxNormTrustRegion<blocksize> trustRegion(basis.size(), initialTrustRegionRadius_);
std::vector<BoxConstraint<field_type,blocksize> > trustRegionObstacles;
std::vector<BoxConstraint<field_type,blocksize> > trustRegionObstacles;
Loading