TPZBoostGraph.h

Go to the documentation of this file.

#ifndef TPZBOOSTGRAPH_H
#define TPZBOOSTGRAPH_H
#include "TPZRenumbering.h"

#ifdef USING_BOOST

#include "pzvec.h"
#include <boost/graph/adjacency_list.hpp>
#include <boost/graph/sloan_ordering.hpp>
#include <boost/graph/cuthill_mckee_ordering.hpp>
#include <boost/graph/properties.hpp>
#include <boost/graph/bandwidth.hpp>
#include <boost/graph/profile.hpp>
#include <boost/graph/wavefront.hpp>
#include <boost/graph/properties.hpp>

typedef boost::adjacency_list<boost::setS,boost::vecS,boost::undirectedS,
boost::property<boost::vertex_color_t,boost::default_color_type,
boost::property<boost::vertex_degree_t,int,
boost::property<boost::vertex_priority_t,double > > > > Graph;

typedef boost::graph_traits<Graph>::vertex_descriptor Vertex;
typedef boost::graph_traits<Graph>::vertices_size_type size_type;

typedef std::pair<std::size_t, std::size_t> Pair;


class TPZBoostGraph : public TPZRenumbering {
    
public:
                
                enum GraphType { KMC, KMCExpensive, Sloan };
                
                
                
                TPZBoostGraph(GraphType tp) : TPZRenumbering(), fGType(tp)
                {
                }
    
                TPZBoostGraph(int64_t NElements, int64_t NNodes);

                virtual ~TPZBoostGraph();

                void CompressedResequence(TPZVec<int64_t> &perm, TPZVec<int64_t> &inverseperm);
                
                void ResequenceOld(TPZVec<int> &perm, TPZVec<int> &inverseperm);
                void Resequence(TPZVec<int64_t> &perm, TPZVec<int64_t> &inverseperm);
                void setGType(GraphType M) { fGType = M; }
                virtual void ClearDataStructures();
private:
                TPZVec<Pair> m_Edges;
                Graph m_Graph;
                //Creating two iterators over the vertices
                //graph_traits<Graph>::vertex_iterator ui, ui_end;
                
                boost::property_map<Graph,boost::vertex_degree_t>::type m_Degrees;
                //    = get(vertex_degree, G);
                //   for (boost::tie(ui, ui_end) = vertices(G); ui != ui_end; ++ui)
                //     deg[*ui] = degree(*ui, G);
                
                boost::property_map<Graph, boost::vertex_index_t>::type m_Index_map;
                // = get(vertex_index, G);  
                TPZVec<int64_t> m_Connects;
                
    GraphType fGType;
};



#endif // USING_BOOST

#endif //TPZBOOSTGRAPH_H