manolas
/
vcglib
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

some cleanup of the Taubin tensor estimator (still not working properly)

This commit is contained in:
ganovelli 2008-10-10 14:58:33 +00:00
parent d944271c2a
commit 56857ecdb7
1 changed files with 8 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
vcg/complex/trimesh/update/curvature.h
View File

@ -137,29 +137,26 @@ public:
vcg::face::JumpingPos<FaceType> pos((*vi).VFp(), central_vertex); vcg::face::JumpingPos<FaceType> pos((*vi).VFp(), central_vertex);
// firstV is the first vertex of the 1ring neighboorhood
VertexType* firstV = pos.VFlip(); VertexType* firstV = pos.VFlip();
VertexType* tempV; VertexType* tempV;
float totalDoubleAreaSize = 0.0f; float totalDoubleAreaSize = 0.0f;
if (((firstV->cP()-central_vertex->cP())^(pos.VFlip()->cP()-central_vertex->cP()))*central_vertex->cN()<=0.0f)
{
pos.Set(central_vertex->VFp(), central_vertex);
pos.FlipE();
firstV = pos.VFlip();
}
else pos.Set(central_vertex->VFp(), central_vertex);
// compute the area of each triangle around the central vertex as well as their total area // compute the area of each triangle around the central vertex as well as their total area
do do
{ {
pos.NextE(); // this bring the pos to the next triangle counterclock-wise
tempV = pos.VFlip(); pos.FlipF();
pos.FlipE();
// tempV takes the next vertex in the 1ring neighborhood
tempV = pos.VFlip();
assert(tempV!=central_vertex);
AdjVertex v; AdjVertex v;
v.isBorder = pos.IsBorder(); v.isBorder = pos.IsBorder();
v.vert = tempV; v.vert = tempV;
v.doubleArea = ((pos.F()->V(1)->cP() - pos.F()->V(0)->cP()) ^ (pos.F()->V(2)->cP()- pos.F()->V(0)->cP())).Norm();; v.doubleArea = vcg::DoubleArea(*pos.F());
totalDoubleAreaSize += v.doubleArea; totalDoubleAreaSize += v.doubleArea;
vertices.push_back(v); vertices.push_back(v);