element.hh

Go to the documentation of this file.

 
#ifndef hpElement3d_h
#define hpElement3d_h 
 
#include <cstring>
#include "basics/vectorsMatrices.hh"
#include "geometry/connector.hh"
#include "geometry/cell.hh"
#include "space/element.hh"
#include "space/tmatrix.hh"
 
namespace hp3D {
 
  using concepts::Real;
 
  // *************************************************************** Element **
 
  template<class F>
  class Element : public concepts::ElementWithCell<F> {
  public:
    inline Element(concepts::TColumn<Real>* T, uint p) : T_(T) {
      p_[0] = p_[1] = p_[2] = p;
    }
 
    inline Element(concepts::TColumn<Real>* T, const ushort* p) : T_(T) {
      std::memcpy(p_, p, 3*sizeof(ushort));
    }
 
    inline const ushort* p() const { return p_; }
 
    virtual const concepts::Connector3& support() const = 0;
 
    virtual const concepts::TMatrix<Real>& T() const { return T_; }
 
    void appendT(concepts::TColumn<Real>* T) { T_.append(T); }
 
    virtual concepts::Real3d vertex(uint i) const = 0;
 
    virtual const concepts::Cell3& cell() const = 0;
 
    virtual bool operator<(const Element<F>& elm) const = 0;
  protected:
    virtual std::ostream& info(std::ostream& os) const;
 
    concepts::TMatrix<Real> T_;
  private:
    ushort p_[3];
  };
 
} // namespace hp3d
 
#endif // hpElement3d_h