| 1 |
francois |
979 |
// 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(std::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(std::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(std::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 |
|
|
} |
