gridInfo.hh

Go to the documentation of this file.

 
#ifndef hp2DgridInfo_hh
#define hp2DgridInfo_hh
 
#include "toolbox/hashMap.hh"
#include "toolbox/sequence.hh"
#include "space.hh"
#include "geometry/connector.hh"
#include "hp2D/element.hh"
 
namespace hp2D {
 
  // *************************************************************** GridInfo **
 
  template<class F = Real>
  class GridInfo: public concepts::OutputOperator {
  public:
    typedef concepts::ElementAndFacette<hp2D::Element<F> > UnderlyingElement;
    typedef concepts::HashMap<concepts::Sequence<UnderlyingElement> > EdgeToElms;
    typedef concepts::HashMap<concepts::Sequence<uint> > CoarseToSuccessor;
 
    enum Type {
      regular = 0, cIrregular = 1, sIrregular = 2
    };
 
    GridInfo(const concepts::SpaceOnCells<F>& spc,
       const concepts::Set<uint>& edgeAttr);
 
    const EdgeToElms& get(Type type) const;
 
    const CoarseToSuccessor& get() const;
 
    const concepts::Set<uint>& attribute() const;
 
  protected:
    virtual std::ostream& info(std::ostream& os) const;
 
  private:
 
    // copy to attribute to identify active edges
    concepts::Set<uint> edgeAttr_;
 
    void build_(const concepts::SpaceOnCells<F>& spc,
    const concepts::Set<uint> edgeAttr);
 
    EdgeToElms map_reg;
 
    EdgeToElms map_cIrr;
 
    EdgeToElms map_sIrr;
 
    CoarseToSuccessor map_CtS;
 
    //test for edge-attribute is in the Set of requested edgeAttributes
    inline bool activeEdge_(const concepts::Set<uint>& edgeAttr,
          const concepts::Connector1& edge) const {
      return edgeAttr.exist(edge.attrib());
    }
 
    void get_sIrr_(const concepts::Connector1* edge, uint cE,
       const std::set<uint>& edge_1, std::set<uint>& edge_s,
       concepts::HashMap<concepts::Sequence<uint> >& cChld) const;
 
  };
 
} // namespace hp2D
 
#endif // hp2DgridInfo_hh