Unverified Commit 37a046cf authored Aug 14, 2019 by Jibril B. Coulibaly Committed by GitHub Aug 14, 2019

Update pair_granular.cpp

Modified PairGranular::single function to return the total normal force into argument fforce.

This was done for pair styles gran/* but not for the granular pari_style, resulting in the variable fforce being uninitialized.

parent a12b959f

src/GRANULAR/pair_granular.cpp

+5 −1

Original line number	Diff line number	Diff line
		@@ -1631,6 +1631,10 @@ double PairGranular::single(int i, int j, int itype, int jtype,
		}
		}

		// set force and return no energy

		fforce = Fntot*rinv;

		// set single_extra quantities

		svector[0] = fs1;

Admin message