From 25af2e45b2744eefc5f49809fa034b0cec4e7a69 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?=E5=A4=A7=E7=AA=AA=20=E8=B2=B4=E6=B4=8B?=
 <takahiro@ookubonoMac-mini.local>
Date: Thu, 17 Mar 2022 18:25:54 +0900
Subject: [PATCH] add omp critical, corrected openmp

---
 dump2sq.c | 3 +++
 1 file changed, 3 insertions(+)

diff --git a/dump2sq.c b/dump2sq.c
index a27c30b..f9de1f4 100644
--- a/dump2sq.c
+++ b/dump2sq.c
@@ -248,8 +248,11 @@ void CalcCn_pp(DATA *data, HEAD *head, TYPE *type, DATASET *cn,
       if (cn->xmin < d && d < cn->xmax){
 	p = (int)floor((d - cn->xmin)/cn->dx);
 	k = type->pairid[data[i].type_new][data[j].type_new];
+	#pragma omp critical
+	{
 	cn_[p][k] = cn_[p][k] + 1.0;
 	cn_[p][cn->m-1] = cn_[p][cn->m-1] + 1.0;
+	}
       }
     }
   }