/*! Subdivides the face at faceIndex given 2 not subdividable neighbors
*/
std::vector< std::vector<Vector3<float> > >
AdaptiveLoopSubdivisionMesh::Subdivide2(unsigned int faceIndex){
// We know we have otwo false faces
// 1. Start by orienting the triangle so that we always have the same case
// We find the start edge as the edge who shares face with the _subdividable_ neighbor
EdgeIterator eit = GetEdgeIterator( f(faceIndex).edge );
unsigned int num = 0;
while( !Subdividable(eit.Pair().GetEdgeFaceIndex()) ){
eit.Pair().Next();
assert(num++ < 3);
}
// go back to inner edge
eit.Pair();
// 2. Now find the vertices
std::vector< std::vector<Vector3<float> > > faces;
// corner vertices, labeled from current edge's origin vertex
Vector3<float> v1 = VertexRule(eit.GetEdgeVertexIndex());
Vector3<float> v2 = VertexRule(eit.Next().GetEdgeVertexIndex());
Vector3<float> v3 = VertexRule(eit.Next().GetEdgeVertexIndex());
// edge vertices, labeled from start edge
Vector3<float> e1v = EdgeRule( eit.Next().GetEdgeIndex() );
// 3. Create the 2 faces and push them on the vector
std::vector<Vector3<float> > face;
face.push_back(v1); face.push_back(e1v); face.push_back(v3);
faces.push_back(face);
face.clear();
face.push_back(e1v); face.push_back(v2); face.push_back(v3);
faces.push_back(face);
face.clear();
// 4. Return
return faces;
}
/*! Subdivides the mesh one step, depending on subdividability
*/
void AdaptiveLoopSubdivisionMesh::Subdivide()
{
// Create new mesh and copy all the attributes
HalfEdgeMesh subDivMesh;
subDivMesh.SetTransform(GetTransform());
subDivMesh.SetName(GetName());
subDivMesh.SetColorMap(GetColorMap());
subDivMesh.SetWireframe(GetWireframe());
subDivMesh.SetShowNormals(GetShowNormals());
subDivMesh.SetOpacity(GetOpacity());
if (IsHovering()) subDivMesh.Hover();
if (IsSelected()) subDivMesh.Select();
subDivMesh.mMinCMap = mMinCMap;
subDivMesh.mMaxCMap = mMaxCMap;
subDivMesh.mAutoMinMax = mAutoMinMax;
// loop over each face and create new ones
for(unsigned int i=0; i<GetNumFaces(); i++){
// find neighbor faces
unsigned int f1, f2, f3;
EdgeIterator eit = GetEdgeIterator( f(i).edge );
f1 = eit.Pair().GetEdgeFaceIndex(); eit.Pair();
f2 = eit.Next().Pair().GetEdgeFaceIndex(); eit.Pair();
f3 = eit.Next().Pair().GetEdgeFaceIndex();
unsigned int numNotSubdividable = !Subdividable(f1) + !Subdividable(f2) + !Subdividable(f3);
// Do not subdivide if "self" is not subdividable
if(!Subdividable(i)){
numNotSubdividable = 3;
}
std::vector< std::vector <Vector3<float> > > faces;
switch(numNotSubdividable){
case 0:
// normal subdivision (from LoopSubdivisionMesh)
faces = LoopSubdivisionMesh::Subdivide(i);
break;
case 1:
// special case 1
faces = Subdivide1(i);
break;
case 2:
// special case 2
faces = Subdivide2(i);
break;
case 3:
// trivial case, no subdivision, same as if subdividable(fi) == false
faces = Subdivide3(i);
break;
}
// add the faces (if any) to subDivMesh
for(unsigned int j=0; j<faces.size(); j++){
subDivMesh.AddFace(faces.at(j));
}
}
// Assign the new mesh
*this = AdaptiveLoopSubdivisionMesh(subDivMesh, ++mNumSubDivs);
Update();
}
