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

Trying to fix Euler characteristic (fail unit tests/TestTorus)

Browse Source
This commit is contained in:
James E McClure 2015-11-22 19:33:22 -05:00
parent 5950712b1a
commit b13a55b37c
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
common/pmmc.h
View File

@ -4437,9 +4437,10 @@ inline double geomavg_EulerCharacteristic(DTMutableList<Point> &Points, IntArray
} }
// Each vertex is shared by four cubes // Each vertex is shared by four cubes
nvert = 0.25*npts; nvert = 0.25*double(npts);
// Subtract shared sides to avoid double counting // Subtract shared sides to avoid double counting
nside = 3.0*ntris - ShareSideInternal - 0.5*CountSideExternal; //nside = 3.0*ntris - ShareSideInternal - 0.5*CountSideExternal;
nside = 2.0*double(npts)- 3.0 - 0.5*double(npts-1);
// Number of faces is number of triangles // Number of faces is number of triangles
nface = 1.0*ntris; nface = 1.0*ntris;