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

added term to internal force in order to deny area=0 for each face.

This commit is contained in:
Nico Pietroni 2005-04-21 12:24:41 +00:00
parent dbb88409e3
commit 51040fd364
1 changed files with 15 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
apps/test/segmentation3d/segmentator.h
View File

@ -105,6 +105,7 @@ public:
public: public:
bool intersected; bool intersected;
float kdihedral; float kdihedral;
ScalarType AreaRep;
MyFace() MyFace()
{ {
@ -116,10 +117,8 @@ public:
{ {
__super::Init(k,mass); __super::Init(k,mass);
kdihedral=k_dihedral; kdihedral=k_dihedral;
//if (!intersected) AreaRep=((V(1)->RPos() - V(0)->RPos()) ^ (V(2)->RPos() - V(0)->RPos())).Norm();
SetS(); SetS();
/*if (QualityFace()<0.1f)
ClearS();*/
} }
bool IsActive() bool IsActive()
@ -204,6 +203,19 @@ public:
} }
} }
///area changing constrain penalize area goes to zero
CoordType m0=(V(0)->P()+V(1)->P())/2.f;
CoordType m1=(V(1)->P()+V(2)->P())/2.f;
CoordType m2=(V(2)->P()+V(0)->P())/2.f;
m0=(V(2)->P()-m0).Normalize();//directions
m1=(V(0)->P()-m1).Normalize();
m2=(V(1)->P()-m2).Normalize();
ScalarType FArea=AreaRep/Area()*_k;
V(0)->IntForce()+=m0*FArea/3.f;
V(1)->IntForce()+=m1*FArea/3.f;
V(2)->IntForce()+=m2*FArea/3.f;
///end area constrain
return(__super::Update()); return(__super::Update());
///new ///new
//double stretch; //double stretch;