hp1D/space.hh

Go to the documentation of this file.

/* class for the hp FEM space in 1D
 */

#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<Element<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<Element<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<Element<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<Element<Real> >(*elm_);
  }

} // namespace hp1D

#endif // Space1D_hh