Commit aeaa28f3 authored by Bernd Flemisch's avatar Bernd Flemisch Committed by Katharina Heck
Browse files

[cleanup] replace fvGridGeometry_ by gridGeometry_ for private variables

parent 6a59d0f8
......@@ -73,7 +73,7 @@ public:
std::shared_ptr<const GridGeometry> gridGeometry,
std::shared_ptr<const GridVariables> gridVariables)
: problem_(problem)
, fvGridGeometry_(gridGeometry)
, gridGeometry_(gridGeometry)
, gridVariables_(gridVariables)
, minLevel_(getParamFromGroup<int>(problem->paramGroup(), "Adaptive.MinLevel"))
, maxLevel_(getParamFromGroup<int>(problem->paramGroup(), "Adaptive.MaxLevel"))
......@@ -141,11 +141,11 @@ public:
void calculate(const SolutionVector& sol)
{
//! prepare an indicator for refinement
indicatorVector_.assign(fvGridGeometry_->gridView().size(0), false);
indicatorVector_.assign(gridGeometry_->gridView().size(0), false);
for (const auto& element : elements(fvGridGeometry_->gridView()))
for (const auto& element : elements(gridGeometry_->gridView()))
{
const auto eIdx = fvGridGeometry_->elementMapper().index(element);
const auto eIdx = gridGeometry_->elementMapper().index(element);
//! refine any element being below the minimum level
if (element.level() < minLevel_)
......@@ -163,7 +163,7 @@ public:
continue;
// get the fvGeometry and elementVolVars needed for the bc and source interfaces
auto fvGeometry = localView(*fvGridGeometry_);
auto fvGeometry = localView(*gridGeometry_);
fvGeometry.bind(element);
auto elemVolVars = localView(gridVariables_->curGridVolVars());
......@@ -274,14 +274,14 @@ public:
*/
int operator() (const Element& element) const
{
if (indicatorVector_[fvGridGeometry_->elementMapper().index(element)])
if (indicatorVector_[gridGeometry_->elementMapper().index(element)])
return 1;
return 0;
}
private:
std::shared_ptr<const Problem> problem_; //!< The problem to be solved
std::shared_ptr<const GridGeometry> fvGridGeometry_; //!< The finite volume grid geometry
std::shared_ptr<const GridGeometry> gridGeometry_; //!< The finite volume grid geometry
std::shared_ptr<const GridVariables> gridVariables_; //!< The secondary variables on the grid
std::vector<bool> indicatorVector_; //!< Indicator for BCs/sources
......
......@@ -96,7 +96,7 @@ public:
* \param paramGroup The parameter group in which to look for runtime parameters first (default is "")
*/
FVProblem(std::shared_ptr<const FVGridGeometry> fvGridGeometry, const std::string& paramGroup = "")
: fvGridGeometry_(fvGridGeometry)
: gridGeometry_(fvGridGeometry)
, paramGroup_(paramGroup)
{
// set a default name for the problem
......@@ -439,7 +439,7 @@ public:
{
NumEqVector source(0);
auto scvIdx = scv.indexInElement();
auto key = std::make_pair(fvGridGeometry_->elementMapper().index(element), scvIdx);
auto key = std::make_pair(gridGeometry_->elementMapper().index(element), scvIdx);
if (pointSourceMap_.count(key))
{
// call the solDependent function. Herein the user might fill/add values to the point sources
......@@ -493,7 +493,7 @@ public:
if (!sources.empty())
{
// calculate point source locations and save them in a map
PointSourceHelper::computePointSourceMap(*fvGridGeometry_,
PointSourceHelper::computePointSourceMap(*gridGeometry_,
sources,
pointSourceMap_);
}
......@@ -586,7 +586,7 @@ public:
//! The finite volume grid geometry
const FVGridGeometry& gridGeometry() const
{ return *fvGridGeometry_; }
{ return *gridGeometry_; }
//! The parameter group in which to retrieve runtime parameters
const std::string& paramGroup() const
......@@ -607,8 +607,8 @@ private:
*/
void applyInitialSolutionImpl_(SolutionVector& sol, /*isBox=*/std::true_type) const
{
const auto numDofs = fvGridGeometry_->vertexMapper().size();
const auto numVert = fvGridGeometry_->gridView().size(dim);
const auto numDofs = gridGeometry_->vertexMapper().size();
const auto numVert = gridGeometry_->gridView().size(dim);
sol.resize(numDofs);
// if there are more dofs than vertices (enriched nodal dofs), we have to
......@@ -616,11 +616,11 @@ private:
if (numDofs != numVert)
{
std::vector<bool> dofVisited(numDofs, false);
for (const auto& element : elements(fvGridGeometry_->gridView()))
for (const auto& element : elements(gridGeometry_->gridView()))
{
for (int i = 0; i < element.subEntities(dim); ++i)
{
const auto dofIdxGlobal = fvGridGeometry_->vertexMapper().subIndex(element, i, dim);
const auto dofIdxGlobal = gridGeometry_->vertexMapper().subIndex(element, i, dim);
// forward to implementation if value at dof is not set yet
if (!dofVisited[dofIdxGlobal])
......@@ -635,9 +635,9 @@ private:
// otherwise we directly loop over the vertices
else
{
for (const auto& vertex : vertices(fvGridGeometry_->gridView()))
for (const auto& vertex : vertices(gridGeometry_->gridView()))
{
const auto dofIdxGlobal = fvGridGeometry_->vertexMapper().index(vertex);
const auto dofIdxGlobal = gridGeometry_->vertexMapper().index(vertex);
sol[dofIdxGlobal] = asImp_().initial(vertex);
}
}
......@@ -648,16 +648,16 @@ private:
*/
void applyInitialSolutionImpl_(SolutionVector& sol, /*isBox=*/std::false_type) const
{
sol.resize(fvGridGeometry_->numDofs());
for (const auto& element : elements(fvGridGeometry_->gridView()))
sol.resize(gridGeometry_->numDofs());
for (const auto& element : elements(gridGeometry_->gridView()))
{
const auto dofIdxGlobal = fvGridGeometry_->elementMapper().index(element);
const auto dofIdxGlobal = gridGeometry_->elementMapper().index(element);
sol[dofIdxGlobal] = asImp_().initial(element);
}
}
//! The finite volume grid geometry
std::shared_ptr<const FVGridGeometry> fvGridGeometry_;
std::shared_ptr<const FVGridGeometry> gridGeometry_;
//! The parameter group in which to retrieve runtime parameters
std::string paramGroup_;
......
......@@ -62,7 +62,7 @@ public:
template<class Problem>
FVGridVariables(std::shared_ptr<Problem> problem,
std::shared_ptr<const GridGeometry> fvGridGeometry)
: fvGridGeometry_(fvGridGeometry)
: gridGeometry_(fvGridGeometry)
, curGridVolVars_(*problem)
, prevGridVolVars_(*problem)
, gridFluxVarsCache_(*problem)
......@@ -73,10 +73,10 @@ public:
void init(const SolutionVector& curSol)
{
// resize and update the volVars with the initial solution
curGridVolVars_.update(*fvGridGeometry_, curSol);
curGridVolVars_.update(*gridGeometry_, curSol);
// update the flux variables caches (always force flux cache update on initialization)
gridFluxVarsCache_.update(*fvGridGeometry_, curGridVolVars_, curSol, true);
gridFluxVarsCache_.update(*gridGeometry_, curGridVolVars_, curSol, true);
// set the volvars of the previous time step in case we have an instationary problem
// note that this means some memory overhead in the case of enabled caching, however
......@@ -99,10 +99,10 @@ public:
void update(const SolutionVector& curSol, bool forceFluxCacheUpdate = false)
{
// resize and update the volVars with the initial solution
curGridVolVars_.update(*fvGridGeometry_, curSol);
curGridVolVars_.update(*gridGeometry_, curSol);
// update the flux variables caches
gridFluxVarsCache_.update(*fvGridGeometry_, curGridVolVars_, curSol, forceFluxCacheUpdate);
gridFluxVarsCache_.update(*gridGeometry_, curGridVolVars_, curSol, forceFluxCacheUpdate);
}
//! update all variables after grid adaption
......@@ -134,7 +134,7 @@ public:
curGridVolVars_ = prevGridVolVars_;
// update the flux variables caches
gridFluxVarsCache_.update(*fvGridGeometry_, curGridVolVars_, solution);
gridFluxVarsCache_.update(*gridGeometry_, curGridVolVars_, solution);
}
//! return the flux variables cache
......@@ -168,11 +168,11 @@ public:
//! return the finite volume grid geometry
const GridGeometry& gridGeometry() const
{ return *fvGridGeometry_; }
{ return *gridGeometry_; }
protected:
std::shared_ptr<const GridGeometry> fvGridGeometry_; //!< pointer to the constant grid geometry
std::shared_ptr<const GridGeometry> gridGeometry_; //!< pointer to the constant grid geometry
private:
GridVolumeVariables curGridVolVars_; //!< the current volume variables (primary and secondary variables)
......
......@@ -42,7 +42,7 @@ class GridGeometryView
public:
explicit GridGeometryView(const ActualFVGridGeometry* actualFVGridGeometry)
: fvGridGeometry_(actualFVGridGeometry) {}
: gridGeometry_(actualFVGridGeometry) {}
//! export the GridView type and the discretization method
using GridView = typename ActualFVGridGeometry::GridView;
......@@ -75,35 +75,35 @@ public:
* \brief Return the gridView this grid geometry object lives on
*/
const auto& gridView() const
{ return fvGridGeometry_->gridView(); }
{ return gridGeometry_->gridView(); }
/*!
* \brief Returns the connectivity map of which dofs have derivatives with respect
* to a given dof.
*/
const auto& connectivityMap() const // TODO return correct map
{ return fvGridGeometry_->connectivityMap(); }
{ return gridGeometry_->connectivityMap(); }
/*!
* \brief Returns the mapper for vertices to indices for possibly adaptive grids.
*/
const auto& vertexMapper() const
{ return fvGridGeometry_->vertexMapper(); }
{ return gridGeometry_->vertexMapper(); }
/*!
* \brief Returns the mapper for elements to indices for constant grids.
*/
const auto& elementMapper() const
{ return fvGridGeometry_->elementMapper(); }
{ return gridGeometry_->elementMapper(); }
/*!
* \brief Returns the actual gridGeometry we are a restriction of
*/
const ActualFVGridGeometry& actualfvGridGeometry() const
{ return *fvGridGeometry_; }
{ return *gridGeometry_; }
protected:
const ActualFVGridGeometry* fvGridGeometry_;
const ActualFVGridGeometry* gridGeometry_;
};
......@@ -129,7 +129,7 @@ public:
* \brief The total number of cell centered dofs
*/
std::size_t numDofs() const
{ return this->fvGridGeometry_->numCellCenterDofs(); }
{ return this->gridGeometry_->numCellCenterDofs(); }
};
/*!
......@@ -154,7 +154,7 @@ public:
* \brief The total number of cell centered dofs
*/
std::size_t numDofs() const
{ return this->fvGridGeometry_->numFaceDofs(); }
{ return this->gridGeometry_->numFaceDofs(); }
};
/*!
......
......@@ -95,10 +95,10 @@ public:
//! return the fv grid geometry
[[deprecated("Use gridGeometry() instead. fvGridGeometry() will be removed after 3.1!")]]
const FVGridGeometry& fvGridGeometry() const
{ return (*gridVariables_->fvGridGeometry_); }
{ return (*gridVariables_->gridGeometry_); }
//! return the fv grid geometry
const FVGridGeometry& gridGeometry() const
{ return (*gridVariables_->fvGridGeometry_); }
{ return (*gridVariables_->gridGeometry_); }
// return the actual grid variables
const ActualGridVariables& gridVariables() const
......@@ -249,7 +249,7 @@ public:
void update(const SolutionVector& curSol)
{
ParentType::update(curSol[cellCenterIdx]);
curGridFaceVariables_.update(*this->fvGridGeometry_, curSol[faceIdx]);
curGridFaceVariables_.update(*this->gridGeometry_, curSol[faceIdx]);
}
//! initialize all variables (stationary case)
......@@ -257,8 +257,8 @@ public:
void init(const SolutionVector& curSol)
{
ParentType::init(curSol[cellCenterIdx]);
curGridFaceVariables_.update(*this->fvGridGeometry_, curSol[faceIdx]);
prevGridFaceVariables_.update(*this->fvGridGeometry_, curSol[faceIdx]);
curGridFaceVariables_.update(*this->gridGeometry_, curSol[faceIdx]);
prevGridFaceVariables_.update(*this->gridGeometry_, curSol[faceIdx]);
}
//! initialize all variables (instationary case)
......
......@@ -89,7 +89,7 @@ class FVSpatialParamsOneP
public:
FVSpatialParamsOneP(std::shared_ptr<const GridGeometry> gridGeometry)
: fvGridGeometry_(gridGeometry)
: gridGeometry_(gridGeometry)
, gravity_(0.0)
{
const bool enableGravity = getParam<bool>("Problem.EnableGravity");
......@@ -328,11 +328,11 @@ public:
//! The finite volume grid geometry
[[deprecated("Use gridGeometry() instead. fvGridGeometry() will be removed after 3.1!")]]
const GridGeometry& fvGridGeometry() const
{ return *fvGridGeometry_; }
{ return *gridGeometry_; }
//! The finite volume grid geometry
const GridGeometry& gridGeometry() const
{ return *fvGridGeometry_; }
{ return *gridGeometry_; }
protected:
......@@ -343,7 +343,7 @@ protected:
{ return *static_cast<const Implementation*>(this); }
private:
std::shared_ptr<const GridGeometry> fvGridGeometry_;
std::shared_ptr<const GridGeometry> gridGeometry_;
GlobalPosition gravity_; //!< The gravity vector
Scalar forchCoeffDefault_;
};
......
......@@ -67,7 +67,7 @@ class FVSpatialParamsElastic
public:
//! The constructor
FVSpatialParamsElastic(std::shared_ptr<const GridGeometry> gridGeometry)
: fvGridGeometry_(gridGeometry)
: gridGeometry_(gridGeometry)
, gravity_(0.0)
{
const bool enableGravity = getParam<bool>("Problem.EnableGravity");
......@@ -161,11 +161,11 @@ public:
//! The finite volume grid geometry
[[deprecated("Use gridGeometry() instead. fvGridGeometry() will be removed after 3.1!")]]
const GridGeometry& fvGridGeometry() const
{ return *fvGridGeometry_; }
{ return *gridGeometry_; }
//! The finite volume grid geometry
const GridGeometry& gridGeometry() const
{ return *fvGridGeometry_; }
{ return *gridGeometry_; }
protected:
Implementation &asImp_()
......@@ -175,7 +175,7 @@ protected:
{ return *static_cast<const Implementation*>(this); }
private:
std::shared_ptr<const GridGeometry> fvGridGeometry_;
std::shared_ptr<const GridGeometry> gridGeometry_;
GlobalPosition gravity_; //!< The gravity vector
};
} // end namespace Dumuxs
......
......@@ -88,7 +88,7 @@ class FVSpatialParamsPoroElastic
public:
//! The constructor
FVSpatialParamsPoroElastic(std::shared_ptr<const GridGeometry> gridGeometry)
: fvGridGeometry_(gridGeometry)
: gridGeometry_(gridGeometry)
, gravity_(0.0)
{
const bool enableGravity = getParam<bool>("Problem.EnableGravity");
......@@ -316,11 +316,11 @@ public:
//! The finite volume grid geometry
[[deprecated("Use gridGeometry() instead. fvGridGeometry() will be removed after 3.1!")]]
const GridGeometry& fvGridGeometry() const
{ return *fvGridGeometry_; }
{ return *gridGeometry_; }
//! The finite volume grid geometry
const GridGeometry& gridGeometry() const
{ return *fvGridGeometry_; }
{ return *gridGeometry_; }
protected:
Implementation &asImp_()
......@@ -330,7 +330,7 @@ protected:
{ return *static_cast<const Implementation*>(this); }
private:
std::shared_ptr<const GridGeometry> fvGridGeometry_;
std::shared_ptr<const GridGeometry> gridGeometry_;
GlobalPosition gravity_; //!< The gravity vector
};
} // end namespace Dumux
......
......@@ -67,7 +67,7 @@ public:
*/
NewtonConvergenceWriter(const GridGeometry& gridGeometry,
const std::string& name = "newton_convergence")
: fvGridGeometry_(gridGeometry)
: gridGeometry_(gridGeometry)
, writer_(gridGeometry.gridView(), name, "", "")
{
resize();
......@@ -95,7 +95,7 @@ public:
//! Resizes the output fields. This has to be called whenever the grid changes
void resize()
{
const auto numDofs = fvGridGeometry_.numDofs();
const auto numDofs = gridGeometry_.numDofs();
// resize the output fields
for (int eqIdx = 0; eqIdx < numEq; ++eqIdx)
......@@ -135,7 +135,7 @@ private:
std::size_t id_ = 0UL;
std::size_t iteration_ = 0UL;
const GridGeometry& fvGridGeometry_;
const GridGeometry& gridGeometry_;
Dune::VTKSequenceWriter<GridView> writer_;
......
......@@ -69,7 +69,7 @@ public:
*/
StaggeredNewtonConvergenceWriter(const GridGeometry& gridGeometry,
const std::string& name = "newton_convergence")
: fvGridGeometry_(gridGeometry)
: gridGeometry_(gridGeometry)
, ccWriter_(gridGeometry.gridView(), name, "", "")
, faceWriter_(std::make_shared<PointCloudVtkWriter<Scalar, GlobalPosition>>(coordinates_))
, faceSequenceWriter_(faceWriter_, name + "-face", "","",
......@@ -89,8 +89,8 @@ public:
//! Resizes the output fields. This has to be called whenever the grid changes
void resize()
{
const auto numCellCenterDofs = fvGridGeometry_.numCellCenterDofs();
const auto numFaceDofs = fvGridGeometry_.numFaceDofs();
const auto numCellCenterDofs = gridGeometry_.numCellCenterDofs();
const auto numFaceDofs = gridGeometry_.numFaceDofs();
// resize the cell center output fields
for (int eqIdx = 0; eqIdx < numEqCellCenter; ++eqIdx)
......@@ -109,9 +109,9 @@ public:
}
coordinates_.resize(numFaceDofs);
for (auto&& facet : facets(fvGridGeometry_.gridView()))
for (auto&& facet : facets(gridGeometry_.gridView()))
{
const auto dofIdxGlobal = fvGridGeometry_.gridView().indexSet().index(facet);
const auto dofIdxGlobal = gridGeometry_.gridView().indexSet().index(facet);
coordinates_[dofIdxGlobal] = facet.geometry().center();
}
}
......@@ -163,7 +163,7 @@ private:
std::size_t id_ = 0UL;
std::size_t iteration_ = 0UL;
const GridGeometry& fvGridGeometry_;
const GridGeometry& gridGeometry_;
Dune::VTKSequenceWriter<GridView> ccWriter_;
......
......@@ -64,10 +64,10 @@ public:
* before having been calculated.
*/
TwoPGridAdaptIndicator(std::shared_ptr<const GridGeometry> gridGeometry, const std::string& paramGroup = "")
: fvGridGeometry_(gridGeometry)
: gridGeometry_(gridGeometry)
, refineBound_(std::numeric_limits<Scalar>::max())
, coarsenBound_(std::numeric_limits<Scalar>::lowest())
, maxSaturationDelta_(fvGridGeometry_->gridView().size(0), 0.0)
, maxSaturationDelta_(gridGeometry_->gridView().size(0), 0.0)
, minLevel_(getParamFromGroup<std::size_t>(paramGroup, "Adaptive.MinLevel", 0))
, maxLevel_(getParamFromGroup<std::size_t>(paramGroup, "Adaptive.MaxLevel", 0))
{}
......@@ -120,7 +120,7 @@ public:
//! Reset the indicator to a state that returns false for all elements
refineBound_ = std::numeric_limits<Scalar>::max();
coarsenBound_ = std::numeric_limits<Scalar>::lowest();
maxSaturationDelta_.assign(fvGridGeometry_->gridView().size(0), 0.0);
maxSaturationDelta_.assign(gridGeometry_->gridView().size(0), 0.0);
//! maxLevel_ must be higher than minLevel_ to allow for refinement
if (minLevel_ >= maxLevel_)
......@@ -135,15 +135,15 @@ public:
Scalar globalMin = std::numeric_limits<Scalar>::max();
//! Calculate minimum and maximum saturation
for (const auto& element : elements(fvGridGeometry_->gridView()))
for (const auto& element : elements(gridGeometry_->gridView()))
{
//! Index of the current leaf-element
const auto globalIdxI = fvGridGeometry_->elementMapper().index(element);
const auto globalIdxI = gridGeometry_->elementMapper().index(element);
//! Obtain the saturation at the center of the element
const auto geometry = element.geometry();
const auto elemSol = elementSolution(element, sol, *fvGridGeometry_);
const Scalar satI = evalSolution(element, geometry, *fvGridGeometry_, elemSol, geometry.center())[saturationIdx];
const auto elemSol = elementSolution(element, sol, *gridGeometry_);
const Scalar satI = evalSolution(element, geometry, *gridGeometry_, elemSol, geometry.center())[saturationIdx];
//! Maybe update the global minimum/maximum
using std::min;
......@@ -152,22 +152,22 @@ public:
globalMax = max(satI, globalMax);
//! Calculate maximum delta in saturation for this cell
for (const auto& intersection : intersections(fvGridGeometry_->gridView(), element))
for (const auto& intersection : intersections(gridGeometry_->gridView(), element))
{
//! Only consider internal intersections
if (intersection.neighbor())
{
//! Access neighbor
const auto outside = intersection.outside();
const auto globalIdxJ = fvGridGeometry_->elementMapper().index(outside);
const auto globalIdxJ = gridGeometry_->elementMapper().index(outside);
//! Visit intersection only once
if (element.level() > outside.level() || (element.level() == outside.level() && globalIdxI < globalIdxJ))
{
//! Obtain saturation in the neighbor
const auto outsideGeometry = outside.geometry();
const auto elemSolJ = elementSolution(outside, sol, *fvGridGeometry_);
const Scalar satJ = evalSolution(outside, outsideGeometry, *fvGridGeometry_, elemSolJ, outsideGeometry.center())[saturationIdx];
const auto elemSolJ = elementSolution(outside, sol, *gridGeometry_);
const Scalar satJ = evalSolution(outside, outsideGeometry, *gridGeometry_, elemSolJ, outsideGeometry.center())[saturationIdx];
using std::abs;
Scalar localdelta = abs(satI - satJ);
......@@ -201,7 +201,7 @@ public:
//#endif
//! check if neighbors have to be refined too
for (const auto& element : elements(fvGridGeometry_->gridView(), Dune::Partitions::interior))
for (const auto& element : elements(gridGeometry_->gridView(), Dune::Partitions::interior))
if (this->operator()(element) > 0)
checkNeighborsRefine_(element);
}
......@@ -218,12 +218,12 @@ public:
int operator() (const Element& element) const
{
if (element.hasFather()
&& maxSaturationDelta_[fvGridGeometry_->elementMapper().index(element)] < coarsenBound_)
&& maxSaturationDelta_[gridGeometry_->elementMapper().index(element)] < coarsenBound_)
{
return -1;
}
else if (element.level() < maxLevel_
&& maxSaturationDelta_[fvGridGeometry_->elementMapper().index(element)] > refineBound_)
&& maxSaturationDelta_[gridGeometry_->elementMapper().index(element)] > refineBound_)
{
return 1;
}
......@@ -245,7 +245,7 @@ private:
*/
bool checkNeighborsRefine_(const Element &element, std::size_t level = 1)
{
for(const auto& intersection : intersections(fvGridGeometry_->gridView(), element))
for(const auto& intersection : intersections(gridGeometry_->gridView(), element))
{
if(!intersection.neighbor())
continue;
......@@ -260,7 +260,7 @@ private:
if (outside.level() < maxLevel_ && outside.level() < element.level())
{
// ensure refinement for outside element
maxSaturationDelta_[fvGridGeometry_->elementMapper().index(outside)] = std::numeric_limits<Scalar>::max();
maxSaturationDelta_[gridGeometry_->elementMapper().index(outside)] = std::numeric_limits<Scalar>::max();
if(level < maxLevel_)
checkNeighborsRefine_(outside, ++level);
}
......@@ -269,7 +269,7 @@ private:
return true;
}
std::shared_ptr<const GridGeometry> fvGridGeometry_;
std::shared_ptr<const GridGeometry> gridGeometry_;
Scalar refineBound_;
Scalar coarsenBound_;
......
......@@ -119,7 +119,7 @@ public:
SolutionVector& sol)
: GridDataTransfer()
, problem_(problem)
, fvGridGeometry_(gridGeometry)
, gridGeometry_(gridGeometry)
, gridVariables_(gridVariables)
, sol_(sol)
, adaptionMap_(gridGeometry->gridView().grid(), 0)
......@@ -137,7 +137,7 @@ public:
{
adaptionMap_.resize();
const auto& grid = fvGridGeometry_->gridView().grid();
const auto& grid = gridGeometry_->gridView().grid();
for (auto level = grid.maxLevel(); level >= 0; level--)
{
for (const auto& element : elements(grid.levelGridView(level)))
......@@ -148,11 +148,11 @@ public:
// put values in the map for leaf elements