From e37abc65fea5e9883952e97e290b35c8069ce580 Mon Sep 17 00:00:00 2001
From: phankl <philipp.kloza@gmail.com>
Date: Thu, 18 Aug 2022 15:43:28 +0200
Subject: [PATCH] fixed bug in neigh_common with missing self-inclusion of
 ghost bond atoms

---
 src/MESONT/pair_mesocnt.cpp | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/src/MESONT/pair_mesocnt.cpp b/src/MESONT/pair_mesocnt.cpp
index 31ca825676..def7ee5e83 100644
--- a/src/MESONT/pair_mesocnt.cpp
+++ b/src/MESONT/pair_mesocnt.cpp
@@ -1257,6 +1257,7 @@ void PairMesoCNT::neigh_common(int i1, int i2, int &numred, int *redlist)
   else if (i1 > nlocal - 1) {
     numneigh2 = numneigh[i2];
     neighlist2 = firstneigh[i2];
+    redlist[numred++] = i2;
     for (int j = 0; j < numneigh2; j++)
       redlist[numred++] = neighlist2[j] & NEIGHMASK;
     return;
@@ -1264,6 +1265,7 @@ void PairMesoCNT::neigh_common(int i1, int i2, int &numred, int *redlist)
   else if (i2 > nlocal - 1) {
     numneigh1 = numneigh[i1];
     neighlist1 = firstneigh[i1];
+    redlist[numred++] = i1;
     for (int j = 0; j < numneigh1; j++)
       redlist[numred++] = neighlist1[j] & NEIGHMASK;
     return;