From fb51ffd35e5f111bbf7a679a0dffc562530c1cb1 Mon Sep 17 00:00:00 2001
From: Guilhem Lavaux <guilhem.lavaux@iap.fr>
Date: Fri, 26 Apr 2024 08:53:38 +0200
Subject: [PATCH] Fixup for sph

---
 src/sphSmooth.tcc | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/src/sphSmooth.tcc b/src/sphSmooth.tcc
index c50fb7e..6643dd3 100644
--- a/src/sphSmooth.tcc
+++ b/src/sphSmooth.tcc
@@ -40,8 +40,8 @@ namespace CosmoTool {
     uint32_t requested = newNngb;
 
     if (state != 0) {
-      state->distances = boost::shared_ptr<CoordType[]>(new CoordType[Nsph]);
-      state->ngb = boost::shared_ptr<SPHCell *[]>(new SPHCell *[Nsph]);
+      state->distances = boost::shared_ptr<CoordType[]>(new CoordType[newNngb]);
+      state->ngb = boost::shared_ptr<SPHCell *[]>(new SPHCell *[newNngb]);
     } else {
       state = &internal;
       if (requested > maxNgb) {
@@ -212,7 +212,7 @@ namespace CosmoTool {
       ComputePrecision d = state->distances[i];
       SPHCell &cell = *(state->ngb[i]);
       double kernel_value = getKernel(d / state->smoothRadius) / r3;
-#pragma omp atomic
+#pragma omp atomic update
       cell.val.weight += kernel_value;
     }
   }