Element.h

Go to the documentation of this file.

#ifndef ELEMENT
#define ELEMENT

#include <vector>
#include "AuxDefs.h"
#include "ElementGeometry.h"
#include "BasisFunctions.h"


// TODO: add P1nDBoundaryTrace and P1nDMap

class Element
{
 public:
  
  inline Element(){}
  inline virtual ~Element(){}
  inline Element(const Element &){}
  virtual Element *  clone() const = 0;
  
  // Accessors/Mutators:

  virtual const size_t getNumFields() const = 0;  

  virtual const size_t getDof(size_t field) const = 0; 

  virtual const size_t getNumDerivatives(size_t field) const = 0;  

  virtual const std::vector <double> &getShapes(size_t field) const = 0;  

  virtual const std::vector <double> &getDShapes(size_t field) const = 0;  

  virtual const std::vector <double> &getIntegrationWeights(size_t field) const = 0; 

  virtual const std::vector <double> &getIntegrationPtCoords(size_t field) const = 0;  

  virtual const double getShape(size_t field, size_t quad, size_t shapenumber) const = 0;

  virtual const double getDShape(size_t field, size_t quad, size_t shapenumber, size_t dir)
    const = 0;

  virtual const ElementGeometry & getGeometry() const = 0; 

};




class Element_: public Element {
private:
  void copy (const Element_& that) {
    for (size_t i = 0; i < that.LocalShapes.size (); i++) {
          LocalShapes.push_back ((that.LocalShapes[i])->clone ());
    }
  }

  void destroy () {
    for (size_t i = 0; i < LocalShapes.size (); i++) {
          delete LocalShapes[i];
          LocalShapes[i] = NULL;
    }
  }

public:
  inline Element_ () : Element() {
  }

  inline virtual ~Element_ () {
    destroy ();
  }
  
  Element_ (const Element_ &OldElement_) :
    Element (OldElement_) {
    copy (OldElement_);
  }

  Element_& operator = (const Element_& that) {
    if (this != &that) {
      destroy ();
      copy (that);
    }
    return (*this);
  }

  virtual Element_ * clone () const = 0;

  inline const size_t getDof (size_t field) const {
    return LocalShapes[getFieldShapes (field)]->getBasisDimension ();
  }
  inline const size_t getNumDerivatives (size_t field) const {
    return LocalShapes[getFieldShapes (field)]-> getNumberOfDerivativesPerFunction ();
  }
  inline const std::vector<double> &getShapes (size_t field) const {
    return LocalShapes[getFieldShapes (field)]->getShapes ();
  }
  inline const std::vector<double> &getDShapes (size_t field) const {
    return LocalShapes[getFieldShapes (field)]->getDShapes ();
  }
  inline const std::vector<double> &getIntegrationWeights (size_t field) const {
    return LocalShapes[getFieldShapes (field)]->getIntegrationWeights ();
  }
  inline const std::vector<double> &getIntegrationPtCoords (size_t field) const {
    return LocalShapes[getFieldShapes (field)]->getQuadraturePointCoordinates ();
  }

  inline const double getShape (size_t field, size_t quad, size_t shapenumber) const {
    return getShapes (field)[quad * getDof (field) + shapenumber];
  }

  inline const double getDShape (size_t field, size_t quad, size_t shapenumber, size_t dir) const {
    return getDShapes (field)[quad * getDof (field) * getNumDerivatives (field) + shapenumber
        * getNumDerivatives (field) + dir];
  }

protected:
  inline void addBasisFunctions (const BasisFunctions &BasisFunctionsPointer) {
    LocalShapes.push_back (BasisFunctionsPointer.clone ());
    // XXX: amber: clone has no effect (see ShapesEvaluated__ copy constructor)  and it seems
    // that clone is not needed
    // LocalShapes.push_back( (const_cast<BasisFunctions*> (&BasisFunctionsPointer)) );
  }

  virtual size_t getFieldShapes (size_t Field) const = 0;

  inline size_t getNumShapes () const {
    return LocalShapes.size ();
  }

private:
  std::vector<BasisFunctions *> LocalShapes;
};

class SpecificElementFamily {
};

class LocalToGlobalMap {
public:
  inline LocalToGlobalMap () {
  }
  inline virtual ~LocalToGlobalMap () {
  }
  inline LocalToGlobalMap (const LocalToGlobalMap &) {
  }
  virtual LocalToGlobalMap * clone () const = 0;

  // //! @param ElementMapped: GlobalElementIndex of the Element to be mapped. 
  // //! This sets
  // //! the default Element  object whose field degrees of freedom are
  // //! mapped.
  // // XXX: commented out (amber)
  // // virtual void set (const GlobalElementIndex &ElementMapped) = 0;
 
  // //! @param field field number in the element, 0\f$ \le \f$ field \f$\le\f$ Nfields-1\n
  // //! @param dof   number of degree of freedom in that field, 
  // //! \f$ 0 \le \f$ dof \f$\le\f$ Ndofs-1 \n
  // //! map returns the GlobalDofIndex associated to degree of freedom "dof" of field "field" 
  // //! in the default Element object. The latter is  
  // //! set with the function set(const Element &).
  // // XXX: commented out (amber)
  // // virtual const GlobalDofIndex map (size_t field, size_t dof) const = 0;

  virtual const GlobalDofIndex map (size_t field, size_t dof, const GlobalElementIndex & ElementMapped) const = 0;

  virtual const size_t getNumElements () const = 0;

  virtual const size_t getNumFields (const GlobalElementIndex & ElementMapped) const = 0;

  virtual const size_t getNumDof (const GlobalElementIndex & ElementMapped, size_t field) const = 0;

  virtual const size_t getTotalNumDof () const = 0;
};

#endif

---