space2D.hh

Go to the documentation of this file.

 
#ifndef linSpace2D_hh
#define linSpace2D_hh
#include "toolbox/hashMap.hh"
#include "element2D.hh"
#include "basics/typedefs.hh"
#include "space/space.hh"
#include "space/element.hh"
#include "toolbox/scannerConnectors.hh"
 
namespace concepts {
 
  // forward declarations
  class Mesh2;
  class BoundaryConditions;
 
  // ************************************************************** Scanners **
 
  template<>
    class Scan<linearFEM::Triangle> : public Scan<linearFEM::Element> {
      public:
        virtual linearFEM::Triangle& operator++(int) = 0;
    };
 
  template<>
    class Scan<linearFEM::Quad> : public Scan<linearFEM::Element> {
      public:
        virtual linearFEM::Quad& operator++(int) = 0;
    };
 
} // namespace concepts
 
namespace linearFEM {
 
  // ************************************************************** Linear2d **
 
  class Linear2d : public concepts::Space<Real> {
    public:
      typedef concepts::Scan<Element> Scan;
      typedef void (*SMap)(const concepts::TColumn<Real>&,
          concepts::TColumn<Real>&);
 
      Linear2d(concepts::Mesh2& msh, uint level,
          concepts::BoundaryConditions* bc = 0);
      virtual ~Linear2d();
      virtual uint dim() const { return dim_; }
      virtual uint nelm() const { return nelm_; }
      virtual Scan* scan() const;
      void rebuild();
      void adjust(const Quad& elm, const short level);
 
      const SMap& S(uint i) const;
 
      static void S0(const concepts::TColumn<Real>& src,
          concepts::TColumn<Real>& dst); // bottom left
      static void S1(const concepts::TColumn<Real>& src,
          concepts::TColumn<Real>& dst); // bottom right
      static void S2(const concepts::TColumn<Real>& src,
          concepts::TColumn<Real>& dst); // top right
      static void S3(const concepts::TColumn<Real>& src,
          concepts::TColumn<Real>& dst); // top left
 
    protected:
      virtual std::ostream& info(std::ostream& os) const;
    private:
      uint dim_;
      uint nelm_;
      bool rebuild_;
      concepts::Mesh2& msh_;
      concepts::BoundaryConditions* bc_;
      concepts::Joiner<Element*, 1>* elm_;
      std::unordered_map<uint, short> adj_;
      std::unordered_map<uint, short> ctrl0_;
      std::unordered_map<uint, short> ctrl0L_;
      std::unordered_map<uint, short> ctrl1_;
      std::unordered_map<uint, short> ctrl1L_;
      std::unordered_map<uint, short> ctrl2_;
      std::unordered_map<uint, uint> dof_;
 
      static SMap S_[4];
 
      void rebuild0_(concepts::Connector2& cntr, const int l, int& L);
 
      void rebuild1_(concepts::Quad2d& cell, concepts::TColumn<Real>* T0 = 0);
 
      void rebuild1_(concepts::Triangle2d& cell,
          concepts::TColumn<Real>* T0 = 0);
      void deactivate_(concepts::Connector1& edg);
  };
 
} // namespace linearFEM
 
#endif // linSpace2D_hh