Skip to content
Snippets Groups Projects
Commit 9d063827 authored by sjplimp's avatar sjplimp
Browse files

git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@11019 f3b2605a-c512-4ea7-a41b-209d697bcdaa
parent 54f2a2ae
No related branches found
No related tags found
No related merge requests found
......@@ -294,6 +294,7 @@ void ComputeVoronoi::processCell(voronoicell_neighbor &c, int i)
// count only faces above area threshold
c.face_areas(narea);
have_narea = true;
voro[i][1] = 0.0;
for (j=0; j<narea.size(); ++j)
if (narea[j] > fthresh) voro[i][1] += 1.0;
} else {
......@@ -306,6 +307,7 @@ void ComputeVoronoi::processCell(voronoicell_neighbor &c, int i)
voro[i][2] = c.surface_area();
} else if (surface == VOROSURF_GROUP) {
if (!have_narea) c.face_areas(narea);
voro[i][2] = 0.0;
// loop over all faces (neighbors) and check if they are in the surface group
for (j=0; j<voro[i][1]; ++j)
if (mask[neigh[j]] & sgroupbit) voro[i][2] += narea[j];
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment