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

debugging decl

Browse Source
This commit is contained in:
James E McClure 2018-09-13 21:23:30 -04:00
parent c2b5605e67
commit 33eeda579a
2 changed files with 10 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
analysis/Minkowski.cpp
View File

@ -207,6 +207,14 @@ void Minkowski::ComputeScalar(const DoubleArray Field, const double isovalue)
a2 = object.EdgeAngle(e2);
a3 = object.EdgeAngle(e3);
Ji += 0.08333333333333*(a1*s1+a2*s2+a3*s3);
if (0.08333333333333*(a1*s1+a2*s2+a3*s3) < 0.f){
double intcurv=0.08333333333333*(a1*s1+a2*s2+a3*s3);
double surfarea=sqrt(s*(s-s1)*(s-s2)*(s-s3));
printf("%f, %f: s1=%f,s2=%f,s3=%f, a1=%f,a2=%f,a3=%f\n",surfarea,intcurv,s1,s2,s3,a1,a2,a3);
printf(" P={%f,%f,%f}\n",P1.x,P1.y,P1.z);
printf(" Q={%f,%f,%f}\n",P2.x,P2.y,P2.z);
printf(" R={%f,%f,%f}\n",P3.x,P3.y,P3.z);
}
// Euler characteristic (half edge rule: one face - 0.5*(three edges))
Xi -= 0.5;
}

4
analysis/decl.cpp
View File

@ -388,11 +388,11 @@ double DECL::EdgeAngle(int edge)
angle *= 0.5; // half edge value
}
//1.570796326794897
if (fabs(angle) > 0.f){
/*if (angle < 0.f){
printf("%f, %f (Edge=%i, twin=%i): U={%f, %f, %f}, V={%f, %f, %f}\n",angle,dotprod,edge,halfedge.twin(edge),U.x,U.y,U.z,V.x,V.y,V.z);
printf(" P={%f, %f, %f}, Q={%f, %f, %f}, R={%f, %f, %f} \n",P.x,P.y,P.z,Q.x,Q.y,Q.z,R.x,R.y,R.z);
}
*/
return angle;
}