space1D.hh

Go to the documentation of this file.

 
#ifndef linSpace1D_hh
#define linSpace1D_hh
 
#include "element1D.hh"
#include "basics/typedefs.hh"
#include "space/space.hh"
#include "space/element.hh"
#include "toolbox/scannerConnectors.hh"
#include "toolbox/dynArray.hh"
 
namespace concepts {
 
  // forward declarations
  class Mesh1;
  class BoundaryConditions;
 
  // ************************************************************** Scanners **
 
  template<>
  class Scan<linearFEM::Line> : public Scan<linearFEM::Element> {
  public:
    virtual linearFEM::Line& operator++(int) = 0;
  };
 
} // namespace concepts
 
namespace linearFEM {
 
  // ************************************************************** Linear1d **
 
  class Linear1d : public concepts::Space<Real> {
  public:
    typedef concepts::Scan<Line> Scan;
    typedef void (*SMap)(const concepts::TColumn<Real>&,
       concepts::TColumn<Real>&);
    Linear1d(concepts::Mesh1& msh, uint level,
       concepts::BoundaryConditions* bc = 0);
    virtual ~Linear1d();
    virtual uint dim() const { return dim_; }
    virtual uint nelm() const { return nelm_; }
    virtual Scan* scan() const { return new concepts::PListScan<Line>(*elm_); }
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    uint dim_;
    uint nelm_;
    concepts::Mesh1& msh_;
    concepts::BoundaryConditions* bc_;
    concepts::Joiner<Line*, 1>* elm_;
    void enforceBC_(concepts::Cell1& cell, uint level,
        concepts::DynArray<uint>& bcMap);
    void constructor_(concepts::Cell1& cell, uint level,
          concepts::DynArray<uint>& dofMap,
          const concepts::DynArray<uint>& bcMap);
  };
 
} // namespace linearFEM
 
#endif // linSpace1D_hh