corrected paralilization

Next improment is to paralelize once per loop in (n) and not in (n,v)

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);

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");