geometry/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 << "CellToCellMapping()";
    }
  };

  // ********************************************************* 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 << "RelativeCells()";
    }
  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