OilfieldToolsNet/LBPM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Refactored Euler characteristic into pmmc.h

Browse Source
This commit is contained in:
James E McClure 2015-08-30 12:25:30 -04:00
parent 63204dfef6
commit 705abd8538
1 changed files with 8 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
common/pmmc.h
View File

@ -4042,9 +4042,13 @@ inline double geomavg_EulerCharacteristic(DTMutableList<Point> &Points, IntArray
// Compute the Euler characteristic for triangles in a cube
// Exclude edges and vertices shared with between multiple cubes
double EulerChar;
int nvert=n_nw_pts + n_ns_pts - n_nws_pts;
int nside=2*nvert-3;
int nface=nvert-2;
int nvert=npts;
int nside=2*npts-3;
int nface=npts-2;
//if (ntris != nface){
// nface = ntris;
// nside =
//}
//...........................................................
// Check that this point is not on a previously computed face
// Note direction that the marching cubes algorithm marches
@ -4083,7 +4087,7 @@ inline double geomavg_EulerCharacteristic(DTMutableList<Point> &Points, IntArray
if (!newside) nside-=1;
}
EulerChar = double(nvert - nside + nface);
return EulerChar;
}