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

Debuggin Euler characteristic

Browse Source
This commit is contained in:
James E McClure 2015-08-26 08:28:38 -04:00
parent 661b666406
commit f562dca4e1
1 changed files with 5 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
common/TwoPhase.cpp
View File

@ -588,13 +588,13 @@ void TwoPhase::ComponentAverages()
for (int p=0; p<n_nw_pts; p++){ for (int p=0; p<n_nw_pts; p++){
Point PT = nw_pts(p); Point PT = nw_pts(p);
// Check that this point is not on a previously computed face // Check that this point is not on a previously computed face
if (PT.x - double(i) > 1e-12 && PT.y - double(j) > 1e-12 && PT.z - double(k) > 1e-12) if (PT.x - double(i) < 1e-12 && PT.y - double(j) < 1e-12 && PT.z - double(k) < 1e-12)
nvert-=1; nvert-=1;
} }
for (int p=0; p<n_ns_pts; p++){ for (int p=0; p<n_ns_pts; p++){
Point PT = ns_pts(p); Point PT = ns_pts(p);
// Check that this point is not on a previously computed face // Check that this point is not on a previously computed face
if (PT.x - double(i) > 1e-7 && PT.y - double(j) > 1e-7 && PT.z - double(k) > 1e-7) if (PT.x - double(i) < 1e-12 && PT.y - double(j) < 1e-12 && PT.z - double(k) < 1e-12)
nvert-=1; nvert-=1;
} }
for (int p=0; p<n_nw_tris; p++){ for (int p=0; p<n_nw_tris; p++){
@ -653,6 +653,9 @@ void TwoPhase::ComponentAverages()
if (B.z - double(k) < 1e-12 && C.z-double(k) < 1e-12) newside=false; if (B.z - double(k) < 1e-12 && C.z-double(k) < 1e-12) newside=false;
if (!newside) nside-=1; if (!newside) nside-=1;
} }
if (nside < 0) printf("Negative nside \n!");
if (nface < 0) printf("Negative nface \n!");
// Add the counts for the vertices, sides and faces from the local cube // Add the counts for the vertices, sides and faces from the local cube
ComponentAverages_NWP(NVERT,LabelNWP) += nvert; ComponentAverages_NWP(NVERT,LabelNWP) += nvert;
ComponentAverages_NWP(NSIDE,LabelNWP) += nside; ComponentAverages_NWP(NSIDE,LabelNWP) += nside;