spaceGraph.hh

Go to the documentation of this file.

 
#ifndef space_graph_hh
#define space_graph_hh
 
#include "space/space.hh"
#include "toolbox/graph.hh"
 
namespace concepts
{
  
  // ********************************************************* SpaceGraph<F> **
  template<class F>
  class SpaceGraph
  {
  protected:
    typedef typename F::type type;
    typedef typename F::type::type t_type;
    Graph<type*> graph_;
  public:
    virtual void FillGraph(const F & spc, const int power=2, const bool verbose=false);
    virtual Graph<typename F::type*> getGraph() const;
    virtual void SplitGraph(const int nPparts);
    virtual Sequence<typename F::type*> getSubDomain(const int Part) const;
    virtual Sequence<bool> getSubDomainBool(const int Part) const;
 
    SpaceGraph(const F & spc, const int power=2, const bool verbose=false)
    {
      this->graph_.EmptyGraph();
      this->FillGraph(spc, power, verbose);
    };
 
    SpaceGraph(const std::unique_ptr<F> p_spc, const int power=2, const bool verbose=false)
    {
      this->graph_.EmptyGraph();
      this->FillGraph(*p_spc, power, verbose);
    };
 
    SpaceGraph(const F* p_spc, const int power=2, const bool verbose=false)
    {
      this->graph_.EmptyGraph();
      this->FillGraph(*p_spc, power, verbose);
    };
    
    SpaceGraph() {};
  };
 
  template<class F>
  void SpaceGraph<F>::FillGraph(const F & spc, const int power, const bool verbose)
  {
    this->graph_.EmptyGraph();
        
    std::unique_ptr<typename F::Scanner> sc(spc.scan());
    while(*sc)
      {
        type & elm = (*sc)++;
        int Weight = 1, BaseWeight = elm.T().n();
        for (int i=0; i < power; i++)
          Weight *= BaseWeight;
 
        this->graph_.addVertex(&elm, Weight);
        
        DEBUGL(verbose, "Added vertex at address " << &elm << " of weight " << Weight << " to this graph");
      }
  }
 
  template<class F>
  Graph<typename F::type*> SpaceGraph<F>::getGraph() const
  {
    return this->graph_;
  }
  
  template<class F>
  void SpaceGraph<F>::SplitGraph(const int nbParts)
  {
    this->graph_.SetSubdomains(nbParts);
  }
 
  template<class F>
  Sequence<typename F::type*> SpaceGraph<F>::getSubDomain(const int Part) const
  {
    return this->graph_.GetSubdomain(Part);
  }
 
  template<class F>
  Sequence<bool> SpaceGraph<F>::getSubDomainBool(const int Part) const
  {
    return this->graph_.GetSubdomainBool(Part);
  }
 
}
 
 
#endif // space_graph_hh