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

factor of two

Browse Source
This commit is contained in:
James E McClure 2018-09-15 09:47:32 -04:00
parent 1abc7b2e25
commit a3c79d126f
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
analysis/decl.cpp
View File

@ -390,7 +390,7 @@ double DECL::EdgeAngle(int edge)
dotprod=-dotprod; dotprod=-dotprod;
V.x = -V.x; V.y = -V.y; V.z = -V.z; V.x = -V.x; V.y = -V.y; V.z = -V.z;
} }
if (dotprod > 1.f) dotprod=1.f; if (dotprod > 1.f) dotprod=1.f;
angle = acos(dotprod); angle = acos(dotprod);
/* project onto plane of cube face also works /* project onto plane of cube face also works
@ -415,7 +415,7 @@ double DECL::EdgeAngle(int edge)
//length = sqrt(nx*nx+ny*ny+nz*nz); //length = sqrt(nx*nx+ny*ny+nz*nz);
Point w=0.5*(P+Q)-R; Point w=0.5*(P+Q)-R;
if (W.x*w.x + W.y*w.y + W.z*w.z < 0.f){ if (W.x*w.x + W.y*w.y + W.z*w.z < 0.f){
printf("flip edge normal \n"); printf("flip edge normal \n");
W.x = -W.x; W.x = -W.x;
W.y = -W.y; W.y = -W.y;
W.z = -W.z; W.z = -W.z;
@ -424,7 +424,7 @@ double DECL::EdgeAngle(int edge)
// concave // concave
angle = -angle; angle = -angle;
} }
printf("angle=%f,dot=%f (Edge=%i, twin=%i): P={%f, %f, %f}, Q={%f, %f, %f} U={%f, %f, %f}, V={%f, %f, %f}\n",angle,dotprod,edge,halfedge.twin(edge),P.x,P.y,P.z,Q.x,Q.y,Q.z,U.x,U.y,U.z,V.x,V.y,V.z); printf("angle=%f,dot=%f (Edge=%i, twin=%i): P={%f, %f, %f}, Q={%f, %f, %f} U={%f, %f, %f}, V={%f, %f, %f}\n",angle,dotprod,edge,halfedge.twin(edge),P.x,P.y,P.z,Q.x,Q.y,Q.z,U.x,U.y,U.z,V.x,V.y,V.z);
return angle; return angle;
} }