From d4ca5635745cbd872a3b7675694e07642c20efb4 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Thu, 9 Apr 2015 12:00:14 +0200
Subject: [PATCH] Use correct constructors for the Basis and FufemBasis
---
dune/gfe/riemanniantrsolver.cc | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/dune/gfe/riemanniantrsolver.cc b/dune/gfe/riemanniantrsolver.cc
index 849e9447..835d3039 100644
--- a/dune/gfe/riemanniantrsolver.cc
+++ b/dune/gfe/riemanniantrsolver.cc
@@ -121,7 +121,7 @@ setup(const GridType& grid,
// //////////////////////////////////////////////////////////////////////////////////////
typedef DuneFunctionsBasis<Basis> FufemBasis;
- FufemBasis basis(grid.leafGridView());
+ FufemBasis basis(assembler_->basis_);
OperatorAssembler<FufemBasis,FufemBasis> operatorAssembler(basis, basis);
LaplaceAssembler<GridType, typename FufemBasis::LocalFiniteElement, typename FufemBasis::LocalFiniteElement> laplaceStiffness;
@@ -313,7 +313,7 @@ void RiemannianTrustRegionSolver<Basis,TargetSpace>::solve()
#if HAVE_MPI
MaxNormTrustRegion<blocksize> trustRegion(globalMapper_->size(), initialTrustRegionRadius_);
#else
- Basis basis(grid_->leafGridView());
+ const Basis& basis = assembler_->basis_;
MaxNormTrustRegion<blocksize> trustRegion(basis.indexSet().size(), initialTrustRegionRadius_);
#endif
trustRegion.set(initialTrustRegionRadius_, scaling_);
--
GitLab