box_cut.cc
View Options

			// Box cutting example code
			//
			// Author : Chris H. Rycroft (LBL / UC Berkeley)
			// Email : chr@alum.mit.edu
			// Date : August 30th 2011

			#include "voro++.hh"
			using namespace voro;

			// Parameters controlling the center of the test box
			const double cx=1.5,cy=1.5,cz=1.5;

			int main() {
			double x,y,z;
			voronoicell v;

			// Initialize the Voronoi cell to be a cube of side length 16, centered
			// on the origin
			v.init(-8,8,-8,8,-8,8);

			// Cut by a grid of points in a box of width one, centered on
			// (cx,cy,cz)
			for(x=cx-0.5;x<cx+0.55;x+=0.1) for(y=cy-0.5;y<cy+0.55;y+=0.1)
			for(z=cz-0.5;z<cz+0.55;z+=0.1) v.plane(x,y,z);

			// Output the Voronoi cell in gnuplot format
			v.draw_gnuplot(0,0,0,"box_cut.gnu");

			// Now make a small file that contains the test box
			v.init(cx-0.5,cx+0.5,cy-0.5,cy+0.5,cz-0.5,cz+0.5);
			v.draw_gnuplot(0,0,0,"box_cut.points");
			}

lammps/lib/voro++/examples/extra/box_cut.cc706c4ccdaefdlammm-master