branches/exact/src/common.cc

   1 // Voro++, a 3D cell-based Voronoi library
   2 //
   3 // Author   : Chris H. Rycroft (LBL / UC Berkeley)
   4 // Email    : chr@alum.mit.edu
   5 // Date     : August 30th 2011
   6
   7 /** \file common.cc
   8  * \brief Implementations of the small helper functions. */
   9
  10 #include "common.hh"
  11
  12 namespace voro {
  13
  14 /** \brief Prints a vector of integers.
  15  *
  16  * Prints a vector of integers.
  17  * \param[in] v the vector to print.
  18  * \param[in] fp the file stream to print to. */
  19 void voro_print_vector(vector<int> &v,FILE *fp) {
  20         int k=0,s=v.size();
  21         while(k+4<s) {
  22                 fprintf(fp,"%d %d %d %d ",v[k],v[k+1],v[k+2],v[k+3]);
  23                 k+=4;
  24         }
  25         if(k+3<=s) {
  26                 if(k+4==s) fprintf(fp,"%d %d %d %d",v[k],v[k+1],v[k+2],v[k+3]);
  27                 else fprintf(fp,"%d %d %d",v[k],v[k+1],v[k+2]);
  28         } else {
  29                 if(k+2==s) fprintf(fp,"%d %d",v[k],v[k+1]);
  30                 else fprintf(fp,"%d",v[k]);
  31         }
  32 }
  33
  34 /** \brief Prints a vector of doubles.
  35  *
  36  * Prints a vector of doubles.
  37  * \param[in] v the vector to print.
  38  * \param[in] fp the file stream to print to. */
  39 void voro_print_vector(vector<double> &v,FILE *fp) {
  40         int k=0,s=v.size();
  41         while(k+4<s) {
  42                 fprintf(fp,"%g %g %g %g ",v[k],v[k+1],v[k+2],v[k+3]);
  43                 k+=4;
  44         }
  45         if(k+3<=s) {
  46                 if(k+4==s) fprintf(fp,"%g %g %g %g",v[k],v[k+1],v[k+2],v[k+3]);
  47                 else fprintf(fp,"%g %g %g",v[k],v[k+1],v[k+2]);
  48         } else {
  49                 if(k+2==s) fprintf(fp,"%g %g",v[k],v[k+1]);
  50                 else fprintf(fp,"%g",v[k]);
  51         }
  52 }
  53
  54 /** \brief Prints a vector a face vertex information.
  55  *
  56  * Prints a vector of face vertex information. A value is read, which
  57  * corresponds to the number of vertices in the next face. The routine reads
  58  * this number of values and prints them as a bracked list. This is repeated
  59  * until the end of the vector is reached.
  60  * \param[in] v the vector to interpret and print.
  61  * \param[in] fp the file stream to print to. */
  62 void voro_print_face_vertices(vector<int> &v,FILE *fp) {
  63         int j,k=0,l;
  64         if(v.size()>0) {
  65                 l=v[k++];
  66                 if(l<=1) {
  67                         if(l==1) fprintf(fp,"(%d)",v[k++]);
  68                         else fputs("()",fp);
  69                 } else {
  70                         j=k+l;
  71                         fprintf(fp,"(%d",v[k++]);
  72                         while(k<j) fprintf(fp,",%d",v[k++]);
  73                         fputs(")",fp);
  74                 }
  75                 while((unsigned int) k<v.size()) {
  76                         l=v[k++];
  77                         if(l<=1) {
  78                                 if(l==1) fprintf(fp," (%d)",v[k++]);
  79                                 else fputs(" ()",fp);
  80                         } else {
  81                                 j=k+l;
  82                                 fprintf(fp," (%d",v[k++]);
  83                                 while(k<j) fprintf(fp,",%d",v[k++]);
  84                                 fputs(")",fp);
  85                         }
  86                 }
  87         }
  88 }
  89
  90 }