jianmu
/
lammps
forked from lijiext/lammps
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@13111 f3b2605a-c512-4ea7-a41b-209d697bcdaa

This commit is contained in:
sjplimp 2015-02-13 21:01:20 +00:00
parent 195cbef33a
commit 0b86f84588
1 changed files with 1 additions and 1 deletions

2
src/fix_langevin.cpp
View File

@ -681,7 +681,7 @@ void FixLangevin::omega_thermostat()
double inertiaone;
for (int i = 0; i < nlocal; i++) {
if (mask[i] & groupbit) {
if ((mask[i] & groupbit) && (radius[i] > 0.0)) {
inertiaone = SINERTIA*radius[i]*radius[i]*rmass[i];
if (tstyle == ATOM) tsqrt = sqrt(tforce[i]);
gamma1 = -tendivthree*inertiaone / t_period / ftm2v;