must always initialize jtag

2022-01-29 17:09:40 -05:00
parent c9de8fca1d
commit 7b3adb3f1a
1 changed files with 1 additions and 1 deletions
--- a/src/compute_pair_local.cpp
+++ b/src/compute_pair_local.cpp
@ -226,13 +226,13 @@ int ComputePairLocal::compute_pairs(int flag)
      factor_lj = special_lj[sbmask(j)];
      factor_coul = special_coul[sbmask(j)];
      j &= NEIGHMASK;
+      jtag = tag[j];

      if (!(mask[j] & groupbit)) continue;

      // itag = jtag is possible for long cutoffs that include images of self

      if (newton_pair == 0 && j >= nlocal) {
-        jtag = tag[j];
        if (itag > jtag) {
          if ((itag+jtag) % 2 == 0) continue;
        } else if (itag < jtag) {