correct sign of virial tally in fix smd. issue reported by Dinesh SD on lammps-users (036647c1) · Commits · 郑智淋 / lammps

src/USER-MISC/fix_smd.cpp

+12 −12

Original line number	Diff line number	Diff line
		@@ -270,12 +270,12 @@ void FixSMD::smd_tether()
		ftotal[2] -= fz*massfrac;
		if (evflag) {
		domain->unmap(x[i],image[i],unwrap);
		v[0] = fxmassfracunwrap[0];
		v[1] = fymassfracunwrap[1];
		v[2] = fzmassfracunwrap[2];
		v[3] = fxmassfracunwrap[1];
		v[4] = fxmassfracunwrap[2];
		v[5] = fymassfracunwrap[2];
		v[0] = -fxmassfracunwrap[0];
		v[1] = -fymassfracunwrap[1];
		v[2] = -fzmassfracunwrap[2];
		v[3] = -fxmassfracunwrap[1];
		v[4] = -fxmassfracunwrap[2];
		v[5] = -fymassfracunwrap[2];
		v_tally(i, v);
		}
		}
		@@ -291,12 +291,12 @@ void FixSMD::smd_tether()
		ftotal[2] -= fz*massfrac;
		if (evflag) {
		domain->unmap(x[i],image[i],unwrap);
		v[0] = fxmassfracunwrap[0];
		v[1] = fymassfracunwrap[1];
		v[2] = fzmassfracunwrap[2];
		v[3] = fxmassfracunwrap[1];
		v[4] = fxmassfracunwrap[2];
		v[5] = fymassfracunwrap[2];
		v[0] = -fxmassfracunwrap[0];
		v[1] = -fymassfracunwrap[1];
		v[2] = -fzmassfracunwrap[2];
		v[3] = -fxmassfracunwrap[1];
		v[4] = -fxmassfracunwrap[2];
		v[5] = -fymassfracunwrap[2];
		v_tally(i, v);
		}
		}

Original line number	Diff line number	Diff line
		@@ -270,12 +270,12 @@ void FixSMD::smd_tether()
		ftotal[2] -= fz*massfrac;
		if (evflag) {
		domain->unmap(x[i],image[i],unwrap);
		v[0] = fxmassfracunwrap[0];
		v[1] = fymassfracunwrap[1];
		v[2] = fzmassfracunwrap[2];
		v[3] = fxmassfracunwrap[1];
		v[4] = fxmassfracunwrap[2];
		v[5] = fymassfracunwrap[2];
		v[0] = -fxmassfracunwrap[0];
		v[1] = -fymassfracunwrap[1];
		v[2] = -fzmassfracunwrap[2];
		v[3] = -fxmassfracunwrap[1];
		v[4] = -fxmassfracunwrap[2];
		v[5] = -fymassfracunwrap[2];
		v_tally(i, v);
		}
		}
		@@ -291,12 +291,12 @@ void FixSMD::smd_tether()
		ftotal[2] -= fz*massfrac;
		if (evflag) {
		domain->unmap(x[i],image[i],unwrap);
		v[0] = fxmassfracunwrap[0];
		v[1] = fymassfracunwrap[1];
		v[2] = fzmassfracunwrap[2];
		v[3] = fxmassfracunwrap[1];
		v[4] = fxmassfracunwrap[2];
		v[5] = fymassfracunwrap[2];
		v[0] = -fxmassfracunwrap[0];
		v[1] = -fymassfracunwrap[1];
		v[2] = -fzmassfracunwrap[2];
		v[3] = -fxmassfracunwrap[1];
		v[4] = -fxmassfracunwrap[2];
		v[5] = -fymassfracunwrap[2];
		v_tally(i, v);
		}
		}