use times inverse_np instead

This commit is contained in:
Yifan Li
2023-10-21 00:28:33 -04:00
parent 3d480a01a5
commit 80f72e7a1d

View File

@ -1347,7 +1347,7 @@ void FixPIMDLangevin::compute_tote()
void FixPIMDLangevin::compute_t_prim() void FixPIMDLangevin::compute_t_prim()
{ {
t_prim = 1.5 * atom->natoms * np * force->boltz * temp - total_spring_energy / np; t_prim = 1.5 * atom->natoms * np * force->boltz * temp - total_spring_energy * inverse_np;
} }
/* ---------------------------------------------------------------------- */ /* ---------------------------------------------------------------------- */