Merge pull request #757 from mattrjackson/master

Avoid excess allocations in buildUpdateMap
ORNL · Mar 20, 2024 · ef62603 · ef62603
2 parents 73cfbbb + 8d19d5b
commit ef62603
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/SparseGrids/tsgGridLocalPolynomial.cpp b/SparseGrids/tsgGridLocalPolynomial.cpp
@@ -1266,14 +1266,14 @@ Data2D<int> GridLocalPolynomial::buildUpdateMap(double tolerance, TypeRefinement
         int max_1D_parents = rule->getMaxNumParents();
 
         HierarchyManipulations::SplitDirections split(points);
-
-        #pragma omp parallel for
+        std::vector<int> global_to_pnts(num_points);
+        #pragma omp parallel for firstprivate(global_to_pnts)
         for(int j=0; j<split.getNumJobs(); j++){ // split.getNumJobs() gives the number of 1D interpolants to construct
             int d = split.getJobDirection(j);
             int nump = split.getJobNumPoints(j);
             const int *pnts = split.getJobPoints(j);
 
-            std::vector<int> global_to_pnts(num_points);
+
             std::vector<int> levels(nump);
 
             int max_level = 0;