space.hh

Go to the documentation of this file.

 
#ifndef Space1D_hh
#define Space1D_hh
 
#include "basics/exceptions.hh"
#include "basics/typedefs.hh"
#include "geometry/mesh.hh"
#include "geometry/boundaryConditions.hh"
#include "space/space.hh"
#include "space/hpMethod.hh"
#include "toolbox/scannerConnectors.hh"
#include "toolbox/dynArray.hh"
#include "element.hh"
 
namespace hp1D {
  using concepts::Real;
 
  // ***************************************************************** Space **
 
  class Space :
    public concepts::SpaceOnCells<Real>, 
    public concepts::AdaptiveSpace<Real, concepts::AdaptiveAdjustP<1> > {
  public:
    typedef concepts::Scan<BaseElement<Real> > Scan;
 
    Space(concepts::Mesh1& msh, uint l, uint p,
          concepts::BoundaryConditions* bc = 0);
    virtual ~Space();
    
    inline uint dim() const;
    inline uint nelm() const;
    inline Scan* scan() const;
 
    inline uint dim();
    inline uint nelm();
    inline Scan* scan();
 
    virtual void adjust(const concepts::Element<Real>& elm,
                        const concepts::AdaptiveAdjustP<1>& a);
 
    void rebuild();
 
    void recomputeShapefunctions();
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    concepts::Mesh1& msh_;
    concepts::BoundaryConditions* bc_;
 
    bool rebuild_;
 
    uint dim_;
 
    uint nelm_;
 
    concepts::Joiner<BaseElement<Real>*, 1>* elm_;
 
    concepts::DynArray<concepts::AdaptiveControl<> > ctrl0_;
    concepts::DynArray<concepts::AdaptiveControlP<1> > ctrl1_;
 
    concepts::DynArray<concepts::AdaptiveAdjustP<1> > adj_;
 
    void meshAndPoly_(concepts::EdgeNd& cell, int l, int& L, int& P);
    void buildElements_(concepts::EdgeNd& cell,
      concepts::TColumn<Real>* T0 = 0);
    void enforceBC_(concepts::EdgeNd& cell);
  };
 
  uint Space::dim() const {
    if (rebuild_) throw conceptsException(concepts::SpaceNotBuilt());
    return dim_; 
  }
 
  uint Space::nelm() const {
    if (rebuild_) throw conceptsException(concepts::SpaceNotBuilt());
    return nelm_; 
  }
 
  Space::Scan* Space::scan() const {
    if (rebuild_) throw conceptsException(concepts::SpaceNotBuilt());
    return new concepts::PListScan<BaseElement<Real> >(*elm_);
  }
 
  uint Space::dim() {
    if (rebuild_) rebuild();
    return dim_; 
  }
 
  uint Space::nelm() {
    if (rebuild_) rebuild();
    return nelm_; 
  }
 
  Space::Scan* Space::scan() {
    if (rebuild_) rebuild();
    return new concepts::PListScan<BaseElement<Real> >(*elm_);
  }
 
} // namespace hp1D
 
#endif // Space1D_hh