avoid division by zero

2020-11-01 16:09:30 -05:00
parent 4e147632be
commit 0b51bba75c
1 changed files with 1 additions and 1 deletions
--- a/src/GRANULAR/pair_granular.cpp
+++ b/src/GRANULAR/pair_granular.cpp
@ -1402,7 +1402,7 @@ double PairGranular::single(int i, int j, int itype, int jtype,
  radi = radius[i];
  radj = radius[j];
  radsum = radi + radj;
-  Reff = radi*radj/radsum;
+  Reff = (radsum > 0.0) ? radi*radj/radsum : 0.0;

  bool touchflag;
  E = normal_coeffs[itype][jtype][0];