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

concave/convex bugfix

Browse Source
This commit is contained in:
James E McClure
2018-09-12 21:07:43 -04:00
parent 5b4f724979
commit 69e12e5f37
1 changed files with 4 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
analysis/decl.cpp
View File

@@ -360,7 +360,7 @@ Point DECL::TriNormal(int edge)
double DECL::EdgeAngle(int edge)
{
double angle;
point P,Q,R; // triangle corners
Point P,Q,R; // triangle corners
Point U,V,W; // triangle normal vectors
U = TriNormal(edge);
V = TriNormal(halfedge.twin(edge));
@@ -368,9 +368,9 @@ double DECL::EdgeAngle(int edge)
if (dotprod > 1.f) dotprod=1.f;
angle = acos(dotprod);
// triangle corners
P=halfedge.v1(edge);
Q=halfedge.v2(edge);
R=halfedge.v2(halfedge.next(edge));
P=vertex.coords(halfedge.v1(edge));
Q=vertex.coords(halfedge.v2(edge));
R=vertex.coords(halfedge.v2(halfedge.next(edge)));
// determine if angle is concave or convex based on edge normal
W = 0.5*(P+Q)-R; // vector that lies in plane of triangle
double hypotenuse = sqrt(W.x*V.x+W.y*V.y+W.z*V.z); // hypotenuse of right triangle