singularSet.hh

Go to the documentation of this file.

 
#ifndef singularSet_hh
#define singularSet_hh
 
#include "basics/vectorsMatrices.hh"
#include "toolbox/scannerConnectors.hh"
#include "geometry/topology.hh"
#include "geometry/mesh.hh"
 
namespace hp3D {
// forward declarations
  class Hexahedron;
  class TrivialWeight;
 
  using concepts::Real;
 
  // ********************************************************** SingularEdge **
 
  class SingularEdge {
    friend std::ostream& operator<<(std::ostream& os, const SingularEdge& s);
  public:
    SingularEdge(const concepts::Connector1* edge, const concepts::Real3d vtx0,
     const concepts::Real3d vtx1) :
      edge_(edge), vtx0_(vtx0), vtx1_(vtx1) {
      diffVector_ = diffVectorNormed_ = (vtx1_-vtx0_);
      diffVectorNormed_ *= 1.0 / diffVectorNormed_.l2_2();
      diffVectorTimesVtx0Normed_ = diffVectorNormed_ * vtx0_;
    }
    const concepts::Connector1* edge() const { return edge_; }
    Real distance(const concepts::Real3d point, const Hexahedron* elm) const;
    Real distanceV0(const concepts::Real3d point, const Hexahedron* elm) const;
    Real distanceV1(const concepts::Real3d point, const Hexahedron* elm) const;
  private:
    const concepts::Connector1* edge_;
    const concepts::Real3d vtx0_, vtx1_;
    concepts::Real3d diffVector_;
    concepts::Real3d diffVectorNormed_;
    concepts::Real diffVectorTimesVtx0Normed_;
  };
 
  // ******************************************************** SingularVertex **
 
  class SingularVertex {
    friend std::ostream& operator<<(std::ostream& os, const SingularVertex& s);
  public:
    SingularVertex(const concepts::Connector0* vertex,
       const concepts::Real3d vtx0) :
      vertex_(vertex), vtx0_(vtx0) {}
    const concepts::Connector0* vertex() const { return vertex_; }
    Real distance(concepts::Real3d point) const;
    Real distance(concepts::Real3d point, const Hexahedron* elm) const;
  private:
    const concepts::Connector0* vertex_;
    const concepts::Real3d vtx0_;
  };
 
  // *********************************************************** SingularSet **
 
  class SingularSet {
    friend std::ostream& operator<<(std::ostream& os, const SingularSet& s);
  public:
    SingularSet() : vertices_(0), edges_(0) {}
    ~SingularSet();
    void add(const concepts::Attribute& attrib, concepts::Mesh3& msh);
    concepts::Scan<SingularVertex>* verticesScan() const {
      return new concepts::PListScan<SingularVertex>(*vertices_);
    }
    concepts::Scan<SingularEdge>* edgesScan() const {
      return new concepts::PListScan<SingularEdge>(*edges_);
    }
  private:
    concepts::Joiner<SingularVertex*, 1>* vertices_;
    concepts::Joiner<SingularEdge*, 1>* edges_;
  };
 
} // namespace hp3D
 
#endif // singularSet_hh