From 0b51bba75ca1e58730bc2d7aa1f46fed9b090c47 Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Sun, 1 Nov 2020 16:09:30 -0500
Subject: [PATCH] avoid division by zero

---
 src/GRANULAR/pair_granular.cpp | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/GRANULAR/pair_granular.cpp b/src/GRANULAR/pair_granular.cpp
index a11ba0bf3d..da1e7a12ad 100644
--- 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];