meshRelations.hh

Go to the documentation of this file.

 
#ifndef meshRelations_hh
#define meshRelations_hh
 
#include "toolbox/sequence.hh"
#include "topology.hh"
 
namespace concepts {
 
  // *********************************************************** PointInCell **
 
  template<uint dim = 2>
  class PointInCell : public OutputOperator {  
  public:
    const Connector* cntr_;
    Point<Real,dim>  p_;
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  };
 
  template<>
  class PointInCell<1> : public OutputOperator {  
  public:
    const Connector* cntr_;
    Real             p_;
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  };
 
  // ***************************************************** CellToCellMapping **
  
  template<uint sdim = 2, uint tdim = 2>
  class CellToCellMapping : public OutputOperator {
  public:
    virtual PointInCell<tdim> operator()(const PointInCell<sdim> P) const = 0;
  protected:
    virtual std::ostream& info(std::ostream& os) const {
      return os << concepts::typeOf(*this)<< "()";
    }
  };
 
  // ********************************************************* RelativeCells **
 
  class RelativeCells : public OutputOperator {
  public:
    RelativeCells(std::vector<const Connector*> cells);
 
    bool child (const PointInCell<2> P, PointInCell<2>& T) const;
    bool parent(const PointInCell<2> P, PointInCell<2>& T) const;
    PointInCell<2> finestCell  (const PointInCell<2> P) const;
    PointInCell<2> coarsestCell(const PointInCell<2> P) const;
 
    void rebuildCells();
  protected:
    virtual std::ostream& info(std::ostream& os) const {
      return os << concepts::typeOf(*this)<< "()";
    }
  private:
    mutable concepts::Sequence<const Connector*> cells_;
    
    // mapping to the parent connector
    mutable std::map<const Connector*, const Connector*> parents_;
    // which child number has the connector inside the parent connector
    mutable std::map<const Connector*, uint> childNo_;
 
    bool edgeChild_(const Edge* cntr, const Real2d xi, PointInCell<2>& P) const;
 
    bool quadChild_(const Quad* cntr, const Real2d xi, PointInCell<2>& P) const;
 
    bool edgeParent_(const Edge* cntr, const Real2d xi, uint i,
                     Real2d& Txi) const;
 
    bool quadParent_(const Quad* cntr, const Real2d xi, uint i,
                     Real2d& Txi) const;
  };
 
 
} // namespace concepts
 
 
#endif // meshRelations_hh