when computing only rotational temperature, we must not subtract the default n-dim extra DOFs

This commit is contained in:
Axel Kohlmeyer 2017-04-07 15:11:26 -04:00
parent 1725832b6c
commit 62b9fa22b8
2 changed files with 10 additions and 0 deletions

View File

@ -73,6 +73,11 @@ ComputeTempAsphere::ComputeTempAsphere(LAMMPS *lmp, int narg, char **arg) :
} else error->all(FLERR,"Illegal compute temp/asphere command");
}
// when computing only the rotational temperature,
// do not remove DOFs for translation as set by default
if (mode == ROTATE) extra_dof = 0;
vector = new double[6];
}

View File

@ -67,6 +67,11 @@ ComputeTempSphere::ComputeTempSphere(LAMMPS *lmp, int narg, char **arg) :
} else error->all(FLERR,"Illegal compute temp/sphere command");
}
// when computing only the rotational temperature,
// do not remove DOFs for translation as set by default
if (mode == ROTATE) extra_dof = 0;
vector = new double[6];
// error checks