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

Continuing to refactor Euler characteristic

Browse Source
This commit is contained in:
James E McClure
2015-08-30 13:51:00 -04:00
parent df0a3f5635
commit 97084f4db7
2 changed files with 2 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
common/TwoPhase.cpp
View File

@@ -782,7 +782,7 @@ void TwoPhase::ComponentAverages()
// Compute the Euler characteristic
n_nw_pts=n_nw_tris=0;
geomavg_MarchingCubes(SDn,fluid_isovalue,i,j,k,nw _pts,n_nw_pts,nw_tris,n_nw_tris);
geomavg_MarchingCubes(SDn,fluid_isovalue,i,j,k,nw_pts,n_nw_pts,nw_tris,n_nw_tris);
if (n_nw_pts > 0 ){
int nvert = n_nw_pts;
int nside = 2*nvert-3;

2
common/pmmc.h
View File

@@ -4357,7 +4357,7 @@ inline double geomavg_EulerCharacteristic(DTMutableList<Point> &Points, IntArray
// Exclude edges and vertices shared with between multiple cubes
double EulerChar;
int nvert=npts;
int nside=2*vert-3;
int nside=2*nvert-3;
int nface=nvert-2;
//if (ntris != nface){
// nface = ntris;