//-----------------------------------------------------------------------------
void SimpleMesh::Initialize(){
// Calculate and store all differentials and area
// First update all face normals and triangle areas
for(unsigned int i = 0; i < mFaces.size(); i++){
mFaces.at(i).normal = FaceNormal(i);
}
// reset the clock
timespec tS;
tS.tv_sec = 0;
tS.tv_nsec = 0;
clock_settime(CLOCK_PROCESS_CPUTIME_ID, &tS);
// Then update all vertex normals and curvature
for(unsigned int i = 0; i < mVerts.size(); i++){
// Vertex normals are just weighted averages
mVerts.at(i).normal = VertexNormal(i);
}
clock_gettime(CLOCK_PROCESS_CPUTIME_ID, &tS);
std::cout << "Time for vertnorm calc is: \n" << "sec\n" << tS.tv_nsec / 1000000000.0f << "sec\n";
std::cout << tS.tv_nsec / 1000000.0f << " milliseconds" << std::endl;
// Then update vertex curvature
for(unsigned int i = 0; i < mVerts.size(); i++){
mVerts.at(i).curvature = VertexCurvature(i);
// std::cerr << mVerts.at(i).curvature << "\n";
}
// Finally update face curvature
for(unsigned int i = 0; i < mFaces.size(); i++){
mFaces.at(i).curvature = FaceCurvature(i);
}
}
void HalfEdgeMesh::Update() {
// Calculate and store all differentials and area
// First update all face normals and triangle areas
for(unsigned int i = 0; i < GetNumFaces(); i++){
f(i).normal = FaceNormal(i);
}
// Then update all vertex normals and curvature
for(unsigned int i = 0; i < GetNumVerts(); i++){
// Vertex normals are just weighted averages
mVerts.at(i).normal = VertexNormal(i);
}
// Then update vertex curvature
for(unsigned int i = 0; i < GetNumVerts(); i++){
mVerts.at(i).curvature = VertexCurvature(i);
// std::cerr << mVerts.at(i).curvature << "\n";
}
// Finally update face curvature
for(unsigned int i = 0; i < GetNumFaces(); i++){
f(i).curvature = FaceCurvature(i);
}
std::cerr << "Area: " << Area() << ".\n";
std::cerr << "Volume: " << Volume() << ".\n";
// Update vertex and face colors
if (mVisualizationMode == CurvatureVertex) {
std::vector<Vertex>::iterator iter = mVerts.begin();
std::vector<Vertex>::iterator iend = mVerts.end();
float minCurvature = (std::numeric_limits<float>::max)();
float maxCurvature = -(std::numeric_limits<float>::max)();
while (iter != iend) {
if (minCurvature > (*iter).curvature) minCurvature = (*iter).curvature;
if (maxCurvature < (*iter).curvature) maxCurvature = (*iter).curvature;
iter++;
}
std::cerr << "Mapping color based on vertex curvature with range [" << minCurvature << "," << maxCurvature << "]" << std::endl;
iter = mVerts.begin();
while (iter != iend) {
(*iter).color = mColorMap->Map((*iter).curvature, minCurvature, maxCurvature);
iter++;
}
}
else if (mVisualizationMode == CurvatureFace) {
std::vector<Face>::iterator iter = mFaces.begin();
std::vector<Face>::iterator iend = mFaces.end();
float minCurvature = (std::numeric_limits<float>::max)();
float maxCurvature = -(std::numeric_limits<float>::max)();
while (iter != iend) {
if (minCurvature > (*iter).curvature) minCurvature = (*iter).curvature;
if (maxCurvature < (*iter).curvature) maxCurvature = (*iter).curvature;
iter++;
}
std::cerr << "Mapping color based on face curvature with range [" << minCurvature << "," << maxCurvature << "]" << std::endl;
iter = mFaces.begin();
while (iter != iend) {
(*iter).color = mColorMap->Map((*iter).curvature, minCurvature, maxCurvature);
iter++;
}
}
}