diff --git a/analysis/Minkowski.cpp b/analysis/Minkowski.cpp index 352a7464..d2ff0fec 100644 --- a/analysis/Minkowski.cpp +++ b/analysis/Minkowski.cpp @@ -206,7 +206,7 @@ void Minkowski::ComputeScalar(const DoubleArray Field, const double isovalue) a1 = object.EdgeAngle(e1); a2 = object.EdgeAngle(e2); a3 = object.EdgeAngle(e3); - Ji += 0.08333333333333*(a1*s1+a2*s2+a3*s3); + Ji += 0.16666666666666*(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)); diff --git a/analysis/decl.cpp b/analysis/decl.cpp index d6ad631f..ec4200c4 100644 --- a/analysis/decl.cpp +++ b/analysis/decl.cpp @@ -315,7 +315,6 @@ Point DECL::TriNormal(int edge) { Point P,Q,R; Point U,V,W; - double ux,uy,uz,vx,vy,vz; double nx,ny,nz,len; // at cube faces define outward normal to cube if (edge == -1){