branches/dynamic/examples/extra/box_cut.cc

   1 // Single Voronoi cell example code
   2 //
   3 // Author   : Chris H. Rycroft (LBL / UC Berkeley)
   4 // Email    : chr@alum.mit.edu
   5 // Date     : June 1st 2008
   6
   7 #include "cell.cc"
   8
   9 // Parameters controlling the center of the test box
  10 const double cx=1.5,cy=1.5,cz=1.5;
  11
  12 int main() {
  13         double x,y,z;
  14         voronoicell v;
  15
  16         // Initialize the Voronoi cell to be a cube of side length 16, centered
  17         // on the origin
  18         v.init(-8,8,-8,8,-8,8);
  19
  20         // Cut by a grid of points in a box of width one, centered on
  21         // (cx,cy,cz)
  22         for(x=cx-0.5;x<cx+0.55;x+=0.1) {
  23                 for(y=cy-0.5;y<cy+0.55;y+=0.1) {
  24                         for(z=cz-0.5;z<cz+0.55;z+=0.1) {
  25                                 v.plane(x,y,z);
  26                         }
  27                 }
  28         }
  29
  30         // Output the Voronoi cell in gnuplot format
  31         v.draw_gnuplot("box_cut.gnu",0,0,0);
  32
  33         // Now make a small file that contains the test box
  34         v.init(cx-0.5,cx+0.5,cy-0.5,cy+0.5,cz-0.5,cz+0.5);
  35         v.draw_gnuplot("box_cut.points",0,0,0);
  36 }