Subversion Repositories gelsvn

    FaceAttributeVector<int> ftouched;

    for(FaceIDIterator f = m.faces_begin(); f != m.faces_end(); ++f, ++i)

        ftouched[*f] = i;

    vector<Vec3f> vertices;

    vertices.resize(m.no_faces());

    vector<int> faces;

    vector<int> indices;

    // Create new vertices. Each face becomes a vertex whose position

    // is the centre of the face

    i = 0;

    for(FaceIDIterator f = m.faces_begin(); f != m.faces_end(); ++f, ++i)

        vertices[i] = centre(m, *f);

    // Create new faces. Each vertex is a new face with N=valency of vertex

    // edges.

    for(VertexIDIterator v = m.vertices_begin(); v!= m.vertices_end(); ++v){

        if(boundary(m, *v))

            continue;

        HalfEdgeWalker w = m.halfedgewalker(*v);

        // Push vertex indices for this face onto indices vector.

        // The circulator moves around the face in a clockwise fashion

        // so we just reverse the ordering.

        indices.insert(indices.end(), index_tmp.rbegin(), index_tmp.rend());

        // Insert face valency in the face vector.

        faces.push_back(w.no_steps());

    }

    // Clear the manifold before new geometry is inserted.

    m.clear();

    // And build

    m.build(    vertices.size(), 

                reinterpret_cast<float*>(&vertices[0]), 

                faces.size(), 

                &faces[0],

                &indices[0]);

    }

}