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

fix hypotenuse

Browse Source
This commit is contained in:
James E McClure 2018-09-13 11:15:17 -04:00
parent ba8a12ce5b
commit d7341d03c7
1 changed files with 5 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
analysis/decl.cpp
View File

@ -369,12 +369,14 @@ double DECL::EdgeAngle(int edge)
if (dotprod > 1.f) dotprod=1.f;
angle = acos(dotprod);
// triangle corners
int e2 = object.halfedge.next(edge);
int e3 = object.halfedge.next(edge);
P=vertex.coords(halfedge.v1(edge));
Q=vertex.coords(halfedge.v2(edge));
R=vertex.coords(halfedge.v2(halfedge.next(edge)));
Q=vertex.coords(halfedge.v1(e2));
R=vertex.coords(halfedge.v1(e3));
// 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
double hypotenuse = sqrt(W.x*W.x+W.y*W.y+W.z*W.z + V.x*V.x+V.y*V.y+V.z*V.z); // hypotenuse of right triangle
double length = sqrt((W.x+V.x)*(W.x+V.x) + (W.y+V.y)*(W.y+V.y) + (W.z+V.z)*(W.z+V.z));
if (length > hypotenuse){
// concave