Hello,
Stan’s new per-atom kspace contribution isn’t included with the ‘virial’ option of stress/atom. Looks like it just needs a flag added here:
else if (strcmp(arg[iarg],“virial”) == 0) {
pairflag = 1;
bondflag = angleflag = dihedralflag = improperflag = 1;
fixflag = 1;
} else error->all(FLERR,“Illegal compute stress/atom command”);
happy stresses,
Tim