Subversion Repositories gelsvn

    void rootCC_subdivide(Manifold& m_in, Manifold& m)

    {

        if(&m != &m_in)

            m = m_in;

        VertexAttributeVector<int> vtouched(m.allocated_vertices(), 0);

        VertexAttributeVector<Vec3d> new_pos(m.allocated_vertices(), Vec3d(0));

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

            int v = valency(m, *vid);

            double beta = (4.0 - 2.0 * cos(2.0 * M_PI / v))/(9.0*v);

            new_pos[*vid] = (1.0 - v * beta) * m.pos(*vid);

            for(Walker w = m.walker(*vid); !w.full_circle(); w = w.circulate_vertex_ccw())

            {

                new_pos[*vid] += beta * m.pos(w.vertex());

            }

        }

        vector<FaceID> faces;

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

            faces.push_back(*f);

        for(int i=0;i<faces.size(); ++i)

            vtouched[m.split_face_by_vertex(faces[i])] = 1;

        for(HalfEdgeIDIterator h = m.halfedges_begin(); h != m.halfedges_end(); ++h)

        {

            Walker w = m.walker(*h);

            if(vtouched[w.vertex()] == 0 && vtouched[w.opp().vertex()] == 0)

                m.merge_faces(w.face(), *h);

        }

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

            if(vtouched[*vid] == 0)

                m.pos(*vid) = new_pos[*vid];

    }