Skip to content
Snippets Groups Projects
Commit b2970767 authored by Praetorius, Simon's avatar Praetorius, Simon
Browse files

evaluate marker-function only on the barycenter of the element

parent f5c04c79
No related branches found
No related tags found
No related merge requests found
...@@ -284,11 +284,10 @@ Flag GridFunctionMarker<Grid, PreGridFct>::markGrid(AdaptInfo& adaptInfo) ...@@ -284,11 +284,10 @@ Flag GridFunctionMarker<Grid, PreGridFct>::markGrid(AdaptInfo& adaptInfo)
int currentLevel = e.level(); int currentLevel = e.level();
auto geo = e.geometry(); auto geo = e.geometry();
auto const& ref = Dune::referenceElement(geo); auto const& ref = Dune::referenceElement(geo);
int targetLevel = 0;
int codim = ref.dimension; int codim = ref.dimension;
for (int i = 0; i < ref.size(codim); i++) // evaluate in the center of the element
targetLevel = std::max(targetLevel, int(std::round(localFct(ref.position(i, codim))))); int targetLevel = int(std::round(localFct(ref.position(0,0))));
int m = ((((targetLevel > currentLevel) && (currentLevel < this->maxRefineLevel_)) int m = ((((targetLevel > currentLevel) && (currentLevel < this->maxRefineLevel_))
|| (currentLevel < this->minRefineLevel_)) || (currentLevel < this->minRefineLevel_))
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment