(68b53177) · Commits · 郑智淋 / lammps

src/GRANULAR/fix_nve_gran.cpp

+2 −0

Original line number	Diff line number	Diff line
		@@ -83,6 +83,7 @@ void FixNVEGran::initial_integrate()
		x[i][0] += dtv * v[i][0];
		x[i][1] += dtv * v[i][1];
		x[i][2] += dtv * v[i][2];

		dtfm = dtfphi / (radius[i]radius[i]rmass[i]);
		omega[i][0] += dtfm * torque[i][0];
		omega[i][1] += dtfm * torque[i][1];
		@@ -115,6 +116,7 @@ void FixNVEGran::final_integrate()
		v[i][0] += dtfm * f[i][0];
		v[i][1] += dtfm * f[i][1];
		v[i][2] += dtfm * f[i][2];

		dtfm = dtfphi / (radius[i]radius[i]rmass[i]);
		omega[i][0] += dtfm * torque[i][0];
		omega[i][1] += dtfm * torque[i][1];