linearFEM/space2D.hh

Go to the documentation of this file.

// space for linear FEM in 2D

#ifndef linSpace2D_hh
#define linSpace2D_hh

#if __ICC == 800
#  include <hash_map>
#  define __gnu_cxx std
#else
#  if __GNUC__ >= 3
#    include <ext/hash_map>
#    if __GNUC_MINOR__ < 1 && __GNUC__ == 3
#      define __gnu_cxx std
#    endif
#  else
#    include <hash_map>
#    define __gnu_cxx std
#  endif
#endif

#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_;
      __gnu_cxx::hash_map<uint, short> adj_;
      __gnu_cxx::hash_map<uint, short> ctrl0_;
      __gnu_cxx::hash_map<uint, short> ctrl0L_;
      __gnu_cxx::hash_map<uint, short> ctrl1_;
      __gnu_cxx::hash_map<uint, short> ctrl1L_;
      __gnu_cxx::hash_map<uint, short> ctrl2_;
      __gnu_cxx::hash_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