I have a cubic simulation box and I use fix deform xy to tilt the
simulation box. I then unfix the fix deform and run a simple nvt with
langevin thermostat. I did thermo_style custom to get Pxx,PYY,Pzz,Pxy,Pxz
and Pyx. How is the shear tensor computed for triclinic boxes? Is it equal
Furthermore,how does lammps handle periodic boundary conditions in a
triclinic box because I want to check bond lengths from dump files.