Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F71937039
common.cc
No One
Temporary
Actions
Download File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Subscribers
None
File Metadata
Details
File Info
Storage
Attached
Created
Sat, Jul 13, 20:45
Size
2 KB
Mime Type
text/x-c
Expires
Mon, Jul 15, 20:45 (1 d, 23 h)
Engine
blob
Format
Raw Data
Handle
19008898
Attached To
rLAMMPS lammps
common.cc
View Options
// Voro++, a 3D cell-based Voronoi library
//
// Author : Chris H. Rycroft (LBL / UC Berkeley)
// Email : chr@alum.mit.edu
// Date : August 30th 2011
/** \file common.cc
* \brief Implementations of the small helper functions. */
#include "common.hh"
namespace
voro
{
/** \brief Prints a vector of integers.
*
* Prints a vector of integers.
* \param[in] v the vector to print.
* \param[in] fp the file stream to print to. */
void
voro_print_vector
(
std
::
vector
<
int
>
&
v
,
FILE
*
fp
)
{
int
k
=
0
,
s
=
v
.
size
();
while
(
k
+
4
<
s
)
{
fprintf
(
fp
,
"%d %d %d %d "
,
v
[
k
],
v
[
k
+
1
],
v
[
k
+
2
],
v
[
k
+
3
]);
k
+=
4
;
}
if
(
k
+
3
<=
s
)
{
if
(
k
+
4
==
s
)
fprintf
(
fp
,
"%d %d %d %d"
,
v
[
k
],
v
[
k
+
1
],
v
[
k
+
2
],
v
[
k
+
3
]);
else
fprintf
(
fp
,
"%d %d %d"
,
v
[
k
],
v
[
k
+
1
],
v
[
k
+
2
]);
}
else
{
if
(
k
+
2
==
s
)
fprintf
(
fp
,
"%d %d"
,
v
[
k
],
v
[
k
+
1
]);
else
fprintf
(
fp
,
"%d"
,
v
[
k
]);
}
}
/** \brief Prints a vector of doubles.
*
* Prints a vector of doubles.
* \param[in] v the vector to print.
* \param[in] fp the file stream to print to. */
void
voro_print_vector
(
std
::
vector
<
double
>
&
v
,
FILE
*
fp
)
{
int
k
=
0
,
s
=
v
.
size
();
while
(
k
+
4
<
s
)
{
fprintf
(
fp
,
"%g %g %g %g "
,
v
[
k
],
v
[
k
+
1
],
v
[
k
+
2
],
v
[
k
+
3
]);
k
+=
4
;
}
if
(
k
+
3
<=
s
)
{
if
(
k
+
4
==
s
)
fprintf
(
fp
,
"%g %g %g %g"
,
v
[
k
],
v
[
k
+
1
],
v
[
k
+
2
],
v
[
k
+
3
]);
else
fprintf
(
fp
,
"%g %g %g"
,
v
[
k
],
v
[
k
+
1
],
v
[
k
+
2
]);
}
else
{
if
(
k
+
2
==
s
)
fprintf
(
fp
,
"%g %g"
,
v
[
k
],
v
[
k
+
1
]);
else
fprintf
(
fp
,
"%g"
,
v
[
k
]);
}
}
/** \brief Prints a vector a face vertex information.
*
* Prints a vector of face vertex information. A value is read, which
* corresponds to the number of vertices in the next face. The routine reads
* this number of values and prints them as a bracked list. This is repeated
* until the end of the vector is reached.
* \param[in] v the vector to interpret and print.
* \param[in] fp the file stream to print to. */
void
voro_print_face_vertices
(
std
::
vector
<
int
>
&
v
,
FILE
*
fp
)
{
int
j
,
k
=
0
,
l
;
if
(
v
.
size
()
>
0
)
{
l
=
v
[
k
++
];
if
(
l
<=
1
)
{
if
(
l
==
1
)
fprintf
(
fp
,
"(%d)"
,
v
[
k
++
]);
else
fputs
(
"()"
,
fp
);
}
else
{
j
=
k
+
l
;
fprintf
(
fp
,
"(%d"
,
v
[
k
++
]);
while
(
k
<
j
)
fprintf
(
fp
,
",%d"
,
v
[
k
++
]);
fputs
(
")"
,
fp
);
}
while
((
unsigned
int
)
k
<
v
.
size
())
{
l
=
v
[
k
++
];
if
(
l
<=
1
)
{
if
(
l
==
1
)
fprintf
(
fp
,
" (%d)"
,
v
[
k
++
]);
else
fputs
(
" ()"
,
fp
);
}
else
{
j
=
k
+
l
;
fprintf
(
fp
,
" (%d"
,
v
[
k
++
]);
while
(
k
<
j
)
fprintf
(
fp
,
",%d"
,
v
[
k
++
]);
fputs
(
")"
,
fp
);
}
}
}
}
}
Event Timeline
Log In to Comment