From 8656f52fec7a5fd821896acae3ed314bb664cdc2 Mon Sep 17 00:00:00 2001
From: VCB Ferreira <vascocbf@gmail.com>
Date: Sun, 2 Aug 2026 22:18:50 +0200
Subject: [PATCH] corrected paralilization

Next improment is to paralelize once per loop in (n) and not in (n,v)
---
 nufi/grids.h | 4 ++--
 src/main.cc  | 2 +-
 2 files changed, 3 insertions(+), 3 deletions(-)

diff --git a/nufi/grids.h b/nufi/grids.h
index fc4aa3f..7fc6b62 100644
--- a/nufi/grids.h
+++ b/nufi/grids.h
@@ -2,7 +2,6 @@
 #define GRIDS_H
 
 #include "nufi/cells.h"
-#include "nufi/parameters.h"
 
 #include <deal.II/base/exceptions.h>
 #include <deal.II/dofs/dof_handler.h>
@@ -18,6 +17,7 @@
 #include <deal.II/numerics/vector_tools.h>
 
 #include <memory>
+#include <omp.h>
 #include <unordered_map>
 #include <vector>
 
@@ -59,7 +59,7 @@ template <int dim> struct GridStructure {
       cells.push_back(locations[groups.back().front()].cell);
     }
 
-#pragma omp parallel
+#pragma omp parallel if (!omp_in_parallel())
     {
       std::vector<double> local_solution_buffer(fe->n_dofs_per_cell());
       FEPointEvaluation<dim, dim> evaluator(*mapping, *fe, update_gradients);
diff --git a/src/main.cc b/src/main.cc
index b6f9e42..11bf2a6 100644
--- a/src/main.cc
+++ b/src/main.cc
@@ -16,7 +16,7 @@ void clear_results_directory(const std::string &dir) {
 }
 
 int main() {
-  // omp_set_max_active_levels(1);
+  omp_set_max_active_levels(1);
   std::cout << "Threads: " << omp_get_max_threads() << "\n";
   try {
     clear_results_directory("results");