diff --git a/src/neigh_derive.cpp b/src/neigh_derive.cpp
index 7f10d17d03..022e011880 100644
--- a/src/neigh_derive.cpp
+++ b/src/neigh_derive.cpp
@@ -138,8 +138,10 @@ void Neighbor::half_from_full_newton(NeighList *list)
       } else {
 	if (j >= nall) j %= nall;
 	if (x[j][2] < ztmp) continue;
-	if (x[j][2] == ztmp && x[j][1] < ytmp) continue;
-	if (x[j][2] == ztmp && x[j][1] == ytmp && x[j][0] < xtmp) continue;
+	if (x[j][2] == ztmp) {
+	  if (x[j][1] < ytmp) continue;
+	  if (x[j][1] == ytmp && x[j][0] < xtmp) continue;
+	}
       }
       neighptr[n++] = joriginal;
     }
diff --git a/src/neigh_respa.cpp b/src/neigh_respa.cpp
index 5d33abace7..daed5d6793 100644
--- a/src/neigh_respa.cpp
+++ b/src/neigh_respa.cpp
@@ -276,7 +276,7 @@ void Neighbor::respa_nsq_newton(NeighList *list)
 	  if ((itag+jtag) % 2 == 1) continue;
 	} else {
 	  if (x[j][2] < ztmp) continue;
-	  else if (x[j][2] == ztmp) {
+	  if (x[j][2] == ztmp) {
 	    if (x[j][1] < ytmp) continue;
 	    if (x[j][1] == ytmp && x[j][0] < xtmp) continue;
 	  }