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

Area -> doublearea/2 

added some typename
This commit is contained in:
Federico Ponchio 2006-02-27 18:02:57 +00:00
parent 524fc2574c
commit ec70371a78
1 changed files with 9 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -23,6 +23,9 @@
/**************************************************************************** /****************************************************************************
History History
$Log: not supported by cvs2svn $ $Log: not supported by cvs2svn $
Revision 1.2 2005/10/25 09:17:41 spinelli
correct IsBorder
Revision 1.1 2005/02/22 16:40:29 ganovelli Revision 1.1 2005/02/22 16:40:29 ganovelli
created. This version writes the gaussian curvature on the Q() member of created. This version writes the gaussian curvature on the Q() member of
the vertex the vertex
@ -67,8 +70,8 @@ Discrete Differential-Geometry Operators for Triangulated 2-Manifolds Mark Meyer
static void Gaussian( MeshType & m){ static void Gaussian( MeshType & m){
assert(m.HasPerVertexQuality()); assert(m.HasPerVertexQuality());
MeshType::VertexIterator vi; // iteratore vertice typename MeshType::VertexIterator vi; // iteratore vertice
MeshType::FaceIterator fi; // iteratore facce typename MeshType::FaceIterator fi; // iteratore facce
double *area; // areamix vector double *area; // areamix vector
int i; // index int i; // index
double area0, area1, area2; double area0, area1, area2;
@ -109,9 +112,9 @@ static void Gaussian( MeshType & m){
} }
else // triangolo ottuso else // triangolo ottuso
{ {
(*fi).V(0)->Q() += (*fi).Area() / 3; (*fi).V(0)->Q() += (*fi).DoubleArea() / 6;
(*fi).V(1)->Q() += (*fi).Area() / 3; (*fi).V(1)->Q() += (*fi).DoubleArea() / 6;
(*fi).V(2)->Q() += (*fi).Area() / 3; (*fi).V(2)->Q() += (*fi).DoubleArea() / 6;
} }
} }
@ -129,7 +132,7 @@ static void Gaussian( MeshType & m){
{ {
if(vcg::face::IsBorder((*fi), i)) if(vcg::face::IsBorder((*fi), i))
{ {
MeshType::CoordType e1,e2; typename MeshType::CoordType e1,e2;
vcg::face::Pos<FaceType> hp(&*fi,i,(*fi).V(i)); vcg::face::Pos<FaceType> hp(&*fi,i,(*fi).V(i));
//MeshType::hedgepos_type hp(&*fi,i,(*fi).V(i)); //MeshType::hedgepos_type hp(&*fi,i,(*fi).V(i));
vcg::face::Pos<FaceType> hp1=hp; vcg::face::Pos<FaceType> hp1=hp;