(no subject)

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}"