linDG1D/space.hh

Go to the documentation of this file.

// Space for linear DG FEM in 1D

#ifndef linDG1Dspace_hh
#define linDG1Dspace_hh

#include <iostream>

#include "basics/typedefs.hh"
#include "space/space.hh"
#include "space/elementPairs.hh"
#include "linearFEM/element1D.hh"
#include "linearFEM/space1D.hh"
#include "toolbox/scannerConnectors.hh"

namespace concepts {

  // forward declarations
  class Mesh1;

} // namespace concepts

namespace linDG1D {

  using concepts::Real;

  // ************************************************************** Linear1d **

  class Linear1d : public concepts::Space<Real> {
  public:
    typedef concepts::Scan<linearFEM::Line> Scan;
    typedef void (*SMap)(const concepts::TColumn<Real>&,
       concepts::TColumn<Real>&);
    Linear1d(concepts::Mesh1& msh, uint level);
    virtual ~Linear1d();
    virtual uint dim() const { return dim_; }
    virtual uint nelm() const { return nelm_; }
    virtual Scan* scan() const
    { return new concepts::PListScan<linearFEM::Line>(*elm_); }
    const concepts::ElementPairList<Real>& elmPairs() const
    { return elmPair_; }
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    uint dim_;
    uint nelm_;
    concepts::Mesh1& msh_;
    concepts::Joiner<linearFEM::Line*, 1>* elm_;
    concepts::ElementPairList<Real> elmPair_;
    void constructor_(concepts::Cell1& cell, uint level);

    struct TwoElements {
      TwoElements() { elms[0] = 0; elms[1] = 0; }
      void addElm(const linearFEM::Line* elm, const uint idx) {
  conceptsAssert(elms[1] == 0, concepts::Assertion());
  if (elms[0] == 0) {
    elms[0] = elm;
    idxs[0] = idx;
  }
  else {
    elms[1] = elm;
    idxs[1] = idx;
  }
      }
      const linearFEM::Line* elms[2];
      uint idxs[2];
    };
  };

} // namespace linDG1D

#endif // linDG1Dspace_hh