From 8d19d5b15a53e47482d67be6c3930d9cb98f064e Mon Sep 17 00:00:00 2001
From: Matt Jackson <mattrjackson@gmail.com>
Date: Wed, 20 Mar 2024 12:22:11 -0400
Subject: [PATCH] Avoid excess allocations in buildUpdateMap

---
 SparseGrids/tsgGridLocalPolynomial.cpp | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/SparseGrids/tsgGridLocalPolynomial.cpp b/SparseGrids/tsgGridLocalPolynomial.cpp
index aaf8ce9a1..8e4afb36e 100644
--- 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;