traces.hh

Go to the documentation of this file.

 
#ifndef HP3D_TRACES_HH
#define HP3D_TRACES_HH
 
#include "space/space.hh"
#include "geometry/normalVector.hh"
#include "hp2D/quad.hh"
#include "hp3D/element.hh"
 
namespace concepts {
  // forward declarations
  template <class F>
  class Set; 
}
 
namespace hp3D {
  
  // forward declarations
  class Hexahedron;
  class HexahedronFaceBase;
  
  using concepts::Real;
  
  // ******************************************************** hp3D::TraceSpace **
  
  class TraceSpace : public concepts::SpaceOnCells<Real> {
  public:
    typedef concepts::Scan<hp2D::Element<Real> > Scan;
    typedef concepts::ElementAndFacette<hp3D::Element<Real> > UnderlyingElement;
    
    enum traceTypes { FIRST, MEAN, JUMP }; 
    // for the moment only FIRST available
    
    TraceSpace(const concepts::SpaceOnCells<Real> &spc,
               const concepts::Set<uint> faceAttr,
               enum traceTypes type = FIRST,
               const concepts::FaceNormalVectorRule&
               normalVectorRule = concepts::FaceNormalVectorRule());
    
    virtual ~TraceSpace();
    
    virtual uint dim() const { return dim_; }
    
    virtual uint nelm() const { return nelm_; }
    
    virtual Scan* scan() const 
    {
      return new concepts::PListScan<hp2D::Element<Real> >(*elm_);
    }
    
    const concepts::Sequence<UnderlyingElement>
    uelm(const concepts::Quad &quad) const
    {
        concepts::HashMap<concepts::Sequence<UnderlyingElement> >::
            const_iterator i = uelm_.find(quad.key());
        if( i == uelm_.end() )
          return concepts::Sequence<UnderlyingElement>();
        return i->second;
    }
 
    const concepts::HashMap<concepts::Sequence<UnderlyingElement> > uelm() const
    {
        return uelm_;
    }
 
    virtual void recomputeShapeFunctions();
    
    concepts::Set<uint> getDofIds() const;
    
  protected:
    virtual std::ostream& info(std::ostream& os) const;
    
  private:
    const uint dim_;
    
    uint nelm_;
    
    concepts::Joiner<hp2D::Element<Real>* , 1>* elm_;
 
    concepts::HashMap<concepts::Sequence<UnderlyingElement> >  uelm_;
    
    concepts::HashMap<hp2D::Quad<Real>*> faces_;
    
    std::unique_ptr<HexahedronFaceBase> faceTransfer_;
    
    void setType_(enum traceTypes type,
                  const concepts::FaceNormalVectorRule &normalVectorRule);
    
    template<class F>
    bool build_(const Hexahedron *elm, F condition);
  };
  
}
 
#endif // HP3D_TRACES_HH