linearFEM/element.hh

Go to the documentation of this file.

// elements for linear FEM

#ifndef elementLin_hh
#define elementLin_hh

#include <memory>
#include "basics/typedefs.hh"
#include "space/element.hh"
#include "space/tmatrix.hh"
#include "toolbox/scannerConnectors.hh"

namespace concepts {
  // forward declarations
  class Connector;
  class Cell;
}

namespace linearFEM {

  using concepts::Real;

  // *************************************************************** Element **

  class Element : public concepts::ElementWithCell<Real> {
  public:
    Element(const uint i, uint idx[]);
    Element(concepts::TColumn<Real>* T0, concepts::TColumn<Real>* T1);
    virtual const concepts::Connector& support() const = 0;
    virtual const concepts::Cell& cell() const = 0;
    virtual const concepts::TMatrixBase<Real>& T() const { return *T_; }
  protected:
    std::auto_ptr<concepts::TMatrixBase<Real> > T_;
  };

} // namespace linearFEM

namespace concepts {

  // ************************************************************** Scanners **

  template<>
  class Scan<linearFEM::Element> :
    public concepts::Scan<concepts::Element<Real> > {
  public:
    virtual linearFEM::Element& operator++(int) = 0;
  };
}

#endif // elementLin_hh