ManifoldEdgeSetTopologyContainer.h

Go to the documentation of this file.

/******************************************************************************
*       SOFA, Simulation Open-Framework Architecture, version 1.0 beta 3      *
*                (c) 2006-2008 MGH, INRIA, USTL, UJF, CNRS                    *
*                                                                             *
* This library is free software; you can redistribute it and/or modify it     *
* under the terms of the GNU Lesser General Public License as published by    *
* the Free Software Foundation; either version 2.1 of the License, or (at     *
* your option) any later version.                                             *
*                                                                             *
* This library is distributed in the hope that it will be useful, but WITHOUT *
* ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or       *
* FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License *
* for more details.                                                           *
*                                                                             *
* You should have received a copy of the GNU Lesser General Public License    *
* along with this library; if not, write to the Free Software Foundation,     *
* Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301 USA.          *
*******************************************************************************
*                               SOFA :: Modules                               *
*                                                                             *
* Authors: The SOFA Team and external contributors (see Authors.txt)          *
*                                                                             *
* Contact information: contact@sofa-framework.org                             *
******************************************************************************/
#ifndef SOFA_COMPONENT_TOPOLOGY_MANIFOLDEDGESETTOPOLOGYCONTAINER_H
#define SOFA_COMPONENT_TOPOLOGY_MANIFOLDEDGESETTOPOLOGYCONTAINER_H

#include <sofa/component/topology/EdgeSetTopologyContainer.h>

namespace sofa
{

namespace component
{

namespace topology
{   
    class ManifoldEdgeSetTopologyModifier;
    
    using core::componentmodel::topology::BaseMeshTopology;

    typedef BaseMeshTopology::PointID           PointID;
    typedef BaseMeshTopology::EdgeID            EdgeID;
    typedef BaseMeshTopology::Edge              Edge;
    typedef BaseMeshTopology::SeqEdges          SeqEdges;
    typedef BaseMeshTopology::VertexEdges       VertexEdges;

    class ManifoldEdgeSetTopologyContainer : public EdgeSetTopologyContainer 
    {
        friend class ManifoldEdgeSetTopologyModifier;
        
    public:
        ManifoldEdgeSetTopologyContainer();

        ManifoldEdgeSetTopologyContainer(const sofa::helper::vector< Edge > &edges);

        virtual ~ManifoldEdgeSetTopologyContainer() {}

        virtual void init();

        virtual void clear();
        
        virtual bool checkTopology() const;

        virtual int getNumberConnectedComponents(sofa::helper::vector<unsigned int>& components); 

        inline friend std::ostream& operator<< (std::ostream& out, const ManifoldEdgeSetTopologyContainer& t)
        {
            out << t.m_edge.size();
            for (unsigned int i=0; i<t.m_edge.size(); i++)
                out << " " << t.m_edge[i][0] << " " << t.m_edge[i][1] ;

            return out;
        }

        inline friend std::istream& operator>>(std::istream& in, ManifoldEdgeSetTopologyContainer& t)
        {
            unsigned int s;
            in >> s;
            for (unsigned int i=0; i<s; i++)
            {
                Edge T; in >> T;
                t.m_edge.push_back(T);
            }
            return in;
        }
    protected:
        virtual void createEdgeVertexShellArray();

    private:
        // Describe each connected component, which can be seen as an oriented line
        class ConnectedComponent
        {
        public:

            ConnectedComponent(unsigned int FirstVertexIndex=-1, unsigned int LastVertexIndex=-1, unsigned int size=0,unsigned int ccIndex=0) 
            : FirstVertexIndex(FirstVertexIndex), LastVertexIndex(LastVertexIndex), size(size), ccIndex(ccIndex)
            {}

            virtual ~ConnectedComponent() {}

        public:
            unsigned int FirstVertexIndex; // index of the first vertex on the line
            unsigned int LastVertexIndex; // index of the last vertex on the line

            unsigned int size; // number of the vertices on the line

            unsigned int ccIndex; // index of the connected component stored in the m_ConnectedComponentArray           
        };  

        void resetConnectedComponent()
        {
            m_ComponentVertexArray.clear();
            m_ConnectedComponentArray.clear();
        }

        bool isvoid_ConnectedComponent()
        {           
            return m_ConnectedComponentArray.size()==0;
        }

        void computeConnectedComponent();

        virtual int getNumberOfConnectedComponents()
        {
            computeConnectedComponent();
            return m_ConnectedComponentArray.size();
        }

        virtual int getFirstVertexIndex(unsigned int i)
        {
            computeConnectedComponent();
            assert(i<m_ConnectedComponentArray.size());
            return m_ConnectedComponentArray[i].FirstVertexIndex;
        }

        virtual int getLastVertexIndex(unsigned int i)
        {
            computeConnectedComponent();
            assert(i<m_ConnectedComponentArray.size());
            return m_ConnectedComponentArray[i].LastVertexIndex;
        }

        virtual int getComponentSize(unsigned int i)
        {
            computeConnectedComponent();
            assert(i<m_ConnectedComponentArray.size());
            return m_ConnectedComponentArray[i].size;
        }

        virtual int getComponentIndex(unsigned int i)
        {
            computeConnectedComponent();
            assert(i<m_ConnectedComponentArray.size());
            return m_ConnectedComponentArray[i].ccIndex;
        }

        virtual bool isComponentClosed(unsigned int i)
        {
            computeConnectedComponent();
            assert(i<m_ConnectedComponentArray.size());
            return (m_ConnectedComponentArray[i].FirstVertexIndex == m_ConnectedComponentArray[i].LastVertexIndex);
        }

        int getNextVertex(const unsigned int i) 
        {
            assert(getEdgeVertexShell(i).size()>0);     
            if((getEdgeVertexShell(i)).size()==1 && (getEdge((getEdgeVertexShell(i))[0]))[1]==i){
                return -1;
            }else{
                return (getEdge((getEdgeVertexShell(i))[0]))[1];
            }
        }

        int getPreviousVertex(const unsigned int i) 
        {               
            assert(getEdgeVertexShell(i).size()>0);     
            if((getEdgeVertexShell(i)).size()==1 && (getEdge((getEdgeVertexShell(i))[0]))[0]==i){
                return -1;
            }else{
                return (getEdge((getEdgeVertexShell(i))[0]))[0];
            }
        }

        int getNextEdge(const unsigned int i) 
        {
            if((getEdgeVertexShell(getEdge(i)[1])).size()==1){
                return -1;
            }else{
                return (getEdgeVertexShell(getEdge(i)[1]))[1];
            }
        }

        int getPreviousEdge(const unsigned int i)
        {
            if((getEdgeVertexShell(getEdge(i)[0])).size()==1){
                return -1;
            }else{
                return (getEdgeVertexShell(getEdge(i)[0]))[0];
            }
        }

        const sofa::helper::vector< unsigned int > &getComponentVertexArray() const 
        {
            return m_ComponentVertexArray;
        }

        const sofa::helper::vector< ConnectedComponent > &getConnectedComponentArray() const 
        {
            return m_ConnectedComponentArray;
        }

    private:
        sofa::helper::vector< unsigned int > m_ComponentVertexArray;

        sofa::helper::vector< ConnectedComponent > m_ConnectedComponentArray; 
    };

} // namespace topology

} // namespace component

} // namespace sofa

#endif