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

dbg decl

Browse Source
This commit is contained in:
James E McClure 2018-09-14 22:00:32 -04:00
parent a2b2108d95
commit 8dfe9478bb
1 changed files with 4 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
analysis/decl.cpp
View File

@ -385,7 +385,10 @@ double DECL::EdgeAngle(int edge)
V.x = nx/len; V.y = ny/len; V.z = nz/len; V.x = nx/len; V.y = ny/len; V.z = nz/len;
dotprod = U.x*V.x + U.y*V.y + U.z*V.z; dotprod = U.x*V.x + U.y*V.y + U.z*V.z;
if (dotprod > 1.f) dotprod=1.f; if (dotprod > 1.f) dotprod=1.f;
if (dotprod < 0.f) dotprod=-dotprod; if (dotprod < 0.f){
dotprod=-dotprod;
V.x = -V.x; V.y = -V.y; V.z = -V.z;
}
angle = acos(dotprod); angle = acos(dotprod);
/* project onto plane of cube face also works /* project onto plane of cube face also works
W = U - dotprod*V; W = U - dotprod*V;