singularSet.hh

Go to the documentation of this file.

 
#ifndef singularSet2D_hh
#define singularSet2D_hh
 
#include "basics/vectorsMatrices.hh"
#include "toolbox/scannerConnectors.hh"
#include "geometry/topology.hh"
#include "geometry/mesh.hh"
#include "hp2D/quad.hh"
 
namespace hp2D {
 
  using concepts::Real;
 
  // ******************************************************** SingularVertex **
 
  class SingularVertex {
    friend std::ostream& operator<<(std::ostream& os, const SingularVertex& s);
  public:
    SingularVertex(const concepts::Connector0* vertex,
       const concepts::Real2d vtx0) :
      vertex_(vertex), vtx0_(vtx0) {}
    const concepts::Connector0* vertex() const { return vertex_; }
    Real distance(concepts::Real2d point) const;
    Real distance(concepts::Real2d point, const Quad<Real>* elm) const {
      return distance(elm->chi(point[0], point[1]));
    }
  private:
    const concepts::Connector0* vertex_;
    const concepts::Real2d vtx0_;
  };
 
  // *********************************************************** SingularSet **
 
  class SingularSet {
    friend std::ostream& operator<<(std::ostream& os, const SingularSet& s);
  public:
    SingularSet() : vertices_(0) {}
    ~SingularSet();
    void add(const concepts::Attribute& attrib, concepts::Mesh2& msh);
    concepts::Scan<SingularVertex>* verticesScan() const {
      return new concepts::PListScan<SingularVertex>(*vertices_);
    }
  private:
    concepts::Joiner<SingularVertex*, 1>* vertices_;
  };
 
} // namespace hp2D
 
#endif // singularSet2D_hh