superellipsoid.cc
a6b7299ddc2blammm-devel

superellipsoid.cc
View Options

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

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

			// This function returns a random floating point number between 0 and 1
			double rnd() {return double(rand())/RAND_MAX;}

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

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

			// Cut the cell by 5000 random planes that are scaled to create a
			// superellipsoid
			for(int i=0;i<5000;i++) {
			x=2*rnd()-1;
			y=2*rnd()-1;
			z=2*rnd()-1;
			rsq=xxxx+yyyy+zzz*z;
			if(rsq>0.01&&rsq<1) {
			r=1/sqrt(sqrt(rsq));
			x=r;y=r;z*=r;
			v.plane(xxx,yyy,zzz,xxxx+yyyy+zzz*z);
			}
			}

			// Output the Voronoi cell to a file, in the gnuplot format
			v.draw_gnuplot(0,0,0,"superellipsoid.gnu");

			// Output the Voronoi cell to a file in POV-Ray formats
			v.draw_pov(0,0,0,"superellipsoid_v.pov");
			v.draw_pov_mesh(0,0,0,"superellipsoid_m.pov");
			}