diff --git a/analysis/Minkowski.cpp b/analysis/Minkowski.cpp index 9652db35..e98cfdc7 100644 --- a/analysis/Minkowski.cpp +++ b/analysis/Minkowski.cpp @@ -88,7 +88,7 @@ void Minkowski::ComputeScalar(const DoubleArray& Field, const double isovalue) // printf(" (%i,%i,%i) RP(%i,%i)={%f,%f,%f} {%f,%f,%f} a=%f l=%f \n",i,j,k,e3,object.halfedge.twin(e3),P3.x,P3.y,P3.z,P1.x,P1.y,P1.z,a3,s3); //} // Euler characteristic (half edge rule: one face - 0.5*(three edges)) - Xi -= 1.5; + Xi -= 0.5; } // Euler characteristic -- each vertex shared by four cubes Xi += 0.25*double(object.VertexCount); diff --git a/tests/TestTorus.cpp b/tests/TestTorus.cpp index d3cffa05..0cad888a 100644 --- a/tests/TestTorus.cpp +++ b/tests/TestTorus.cpp @@ -169,10 +169,10 @@ int main(int argc, char **argv) // Averages->Reduce(); Object->MeasureObject(); - double Vi = Object.V(); - double Ai = Object.A(); - double Hi = Object.H(); - double Xi = Object.X(); + double Vi = Object->V(); + double Ai = Object->A(); + double Hi = Object->H(); + double Xi = Object->X(); Vi=sumReduce( Dm->Comm, Vi); Ai=sumReduce( Dm->Comm, Ai); Hi=sumReduce( Dm->Comm, Hi);