From ee34a9048e1c1b3ca560fb9ff0492373f22dd202 Mon Sep 17 00:00:00 2001
From: "Vasco C. B. Ferreira" <vascocbf@gmail.com>
Date: Fri, 31 Jul 2026 15:15:24 +0200
Subject: [PATCH] coarsing below GLOBAL_REFINE bug fixed, test made, working

---
 nufi/poisson_problem.h | 5 +++++
 1 file changed, 5 insertions(+)

diff --git a/nufi/poisson_problem.h b/nufi/poisson_problem.h
index 0d15e69..b7cfd18 100644
--- a/nufi/poisson_problem.h
+++ b/nufi/poisson_problem.h
@@ -427,6 +427,11 @@ template <int dim> void PoissonProblem<dim>::coarse_and_refine_grid(size_t it) {
       Parameters::REFINEMENT_BOTTOM_FRACTION,
       std::numeric_limits<unsigned int>::max(), VectorTools::L2_norm);
 
+  // Avoid coarsing below GLOBAL_REFINEMENT level for CellLocator
+  for (const auto &cell : triangulation.active_cell_iterators())
+    if (cell->level() <= static_cast<int>(Parameters::GLOBAL_REFINEMENT))
+      cell->clear_coarsen_flag();
+
   triangulation.execute_coarsening_and_refinement();
 
   std::cout << "Refinement Finished..." << "\n";