From 5a51c74b9d59adcea676490f12163897a4f2b622 Mon Sep 17 00:00:00 2001 From: Plimpton Date: Thu, 21 Jan 2021 15:41:06 -0700 Subject: [PATCH] vtally change to fix smd --- src/USER-MISC/fix_smd.cpp | 11 +++++------ 1 file changed, 5 insertions(+), 6 deletions(-) diff --git a/src/USER-MISC/fix_smd.cpp b/src/USER-MISC/fix_smd.cpp index 828ceda0ed..b4645f78cf 100644 --- a/src/USER-MISC/fix_smd.cpp +++ b/src/USER-MISC/fix_smd.cpp @@ -58,7 +58,7 @@ FixSMD::FixSMD(LAMMPS *lmp, int narg, char **arg) : extvector = 1; respa_level_support = 1; ilevel_respa = 0; - virial_flag = 1; + virial_global_flag = virial_atom_flag = 1; int argoffs=3; if (strcmp(arg[argoffs],"cvel") == 0) { @@ -182,10 +182,9 @@ void FixSMD::setup(int vflag) void FixSMD::post_force(int vflag) { - // energy and virial setup + // virial setup - if (vflag) v_setup(vflag); - else evflag = 0; + v_init(vflag); if (styleflag & SMD_TETHER) smd_tether(); else smd_couple(); @@ -276,7 +275,7 @@ void FixSMD::smd_tether() v[3] = -fx*massfrac*unwrap[1]; v[4] = -fx*massfrac*unwrap[2]; v[5] = -fy*massfrac*unwrap[2]; - v_tally(i, v); + v_tally(i,v); } } } else { @@ -297,7 +296,7 @@ void FixSMD::smd_tether() v[3] = -fx*massfrac*unwrap[1]; v[4] = -fx*massfrac*unwrap[2]; v[5] = -fy*massfrac*unwrap[2]; - v_tally(i, v); + v_tally(i,v); } } }