total volume

Hi,

Once, I tried from voro++ to create voronoi cell during my simulation in a small sequential DEM code. I can give you the different lines implemented in this code to compute the apparent volume = sum volume each voronoï cell. However you need to go on the website to understand how this powerfull code works and to check if theses following lines are correct.

Good luck !!!
JF

Ps: you can also try with CGAL, but you have to be an expert on C++.

#include "../src_voro++/voro++.cc"

// Set up constants for the container geometry
const fpoint x_min=0,x_max=1; // VORONOI
const fpoint y_min=0,y_max=1; // VORONOI
const fpoint z_min=0; // VORONOI

fpoint z_max=1; // VORONOI
fpoint *Vol_cell;// VORONOI

// Set up the number of blocks that the container is divided
const int n_x=10,n_y=10,n_z=10; // VORONOI

container_poly conp(x_min,x_max,y_min,y_max,z_min,z_max,n_x,n_y,n_z,false,false,false,8); // VORONOI

loop on the particles: {
conp.put(i,_x1(i),_x2(i),_x3(i),_R(i)); // VORONOI
}

conp.store_cell_volumes(Vol_cell);// VORONOI