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

fix Torus

Browse Source
This commit is contained in:
JamesEMcclure 2020-06-24 14:01:10 -04:00
parent 64b12587c7
commit 1a5b46156d
2 changed files with 5 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
analysis/Minkowski.cpp
View File

@ -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); // 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)) // 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 // Euler characteristic -- each vertex shared by four cubes
Xi += 0.25*double(object.VertexCount); Xi += 0.25*double(object.VertexCount);

8
tests/TestTorus.cpp
View File

@ -169,10 +169,10 @@ int main(int argc, char **argv)
// Averages->Reduce(); // Averages->Reduce();
Object->MeasureObject(); Object->MeasureObject();
double Vi = Object.V(); double Vi = Object->V();
double Ai = Object.A(); double Ai = Object->A();
double Hi = Object.H(); double Hi = Object->H();
double Xi = Object.X(); double Xi = Object->X();
Vi=sumReduce( Dm->Comm, Vi); Vi=sumReduce( Dm->Comm, Vi);
Ai=sumReduce( Dm->Comm, Ai); Ai=sumReduce( Dm->Comm, Ai);
Hi=sumReduce( Dm->Comm, Hi); Hi=sumReduce( Dm->Comm, Hi);