Dear LAMMPS users,
I am working with 2D granular system with the periodic boundaries. The system has disk shape particles in two different sizes. I calculated the voronoi volume using voro++ package in lammps. But it seems like there is an issue with the output. The simulation box volume is not equal to the sum of the voronoi cell volume.
#The command I used is,
compute 3 all voronoi/atom
#I added individual voronoi volumes together with following commands,
compute 4 all reduce sum c_3[1]
#Finally I get the output,
dump 1 all custom {irunf} dump.voro.{b3} id type x y z radius c_3[1] c_3[2]
#output of the total cell volume,
variable t1 equal c_4
print "${t1}"