space.hh

Go to the documentation of this file.

 
#ifndef dgSpace_hh
#define dgSpace_hh
 
#include "basics/typedefs.hh"
#include "element.hh"
#include "elementPair.hh"
#include "space/space.hh"
#include "space/elementPairs.hh"
#include "toolbox/scannerConnectors.hh"
#include "geometry/boundaryConditions.hh"
#include "geometry/mesh.hh"
#include "graphics/spaceTraits.hh"
 
#include <map>
#include <memory>
 
#define DEBUG_SPACE 0
 
namespace concepts {
 
  // *************************************************************** Scanner **
  template<>
  class Scan<linDG3D::FvdgElement> : public Scan<ElementWithCell<Real> > {
  public:
    linDG3D::FvdgElement& operator++(int) = 0;
  };
} // namespace concepts
 
namespace linDG3D {
  using concepts::Real;
 
// ****************************************************************** Spaces **
 
  class FvdgSpace : public concepts::SpaceOnCells<Real> {
  public:
    virtual ~FvdgSpace() 
    { concepts::Joiner<FvdgElement*, 1>::destructor(elm_); }
    typedef concepts::Scan<FvdgElement> Scanner;
    uint dim() const { return dim_; }
    uint nelm() const { return nelm_; }
    Scanner* scan() const {
      return new concepts::PListScan<FvdgElement>(*elm_);
    }
    const concepts::ElementPairList<Real>& boundaryElmPairList() const {
      return boundaryElmPairList_;
    }
    const concepts::ElementPairList<Real>& innerElmPairList() const {
      return innerElmPairList_;
    }
    
  protected:
    FvdgSpace(concepts::BoundaryConditions* bc) :
      dim_(0), nelm_(0), bc_(bc), elm_(0), innerElmPairList_(0) {};
    virtual std::ostream& info(std::ostream& os) const;
    uint dim_;
    uint nelm_;
    concepts::BoundaryConditions* bc_;
    concepts::Joiner<FvdgElement*, 1>* elm_;
    
    concepts::ElementPairList<Real> innerElmPairList_;
    concepts::ElementPairList<Real> boundaryElmPairList_;
    
    void createElementPairList_();
  };
  
  //************************************************************ FvdgSpaceP0 **
  class FvdgSpaceP0 : public FvdgSpace {
  public:
    FvdgSpaceP0(concepts::Mesh3& msh, concepts::BoundaryConditions* bc = 0);
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  };
 
  //************************************************************ FvdgSpaceP1 **
  class FvdgSpaceP1 : public FvdgSpace {
  public:
    FvdgSpaceP1(concepts::Mesh3& msh, concepts::BoundaryConditions* bc = 0);
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  };
 
} // namespace linDG3D
 
// *********************************************** Space Traits for Graphics **
 
namespace graphics {
  template<>
  struct spaceTraits<linDG3D::FvdgSpace> {
    static inline uint dim() { return 3; }
    static inline bool positionConnection() { return false; }
    static inline bool boundaryElements() { return false; }
  };
  template<>
  struct spaceTraits<linDG3D::FvdgSpaceP0> {
    static inline uint dim() { return 3; }
    static inline bool positionConnection() { return false; }
    static inline bool boundaryElements() { return false; }
  };
  template<>
  struct spaceTraits<linDG3D::FvdgSpaceP1> {
    static inline uint dim() { return 3; }
    static inline bool positionConnection() { return false; }
    static inline bool boundaryElements() { return false; }
  };
} // namespace graphics
 
#endif // dgSpace_hh