Hi all,
I am very sorry but there were some typos in the last e-mail with subject
'pressure calculation'. The routine that I wrote for pressure calculation should
be as follows:
Hi all,
I am very sorry but there were some typos in the last e-mail with subject
'pressure calculation'. The routine that I wrote for pressure calculation should
be as follows:
I suggest you look at what LAMMPS is doing and compare
to your code. See compute() in pair_lj_cut.cpp and compute_pressure.cpp
where the virial is computed, then converted to pressure.
Steve