ElementBoundaryTrace.h

Go to the documentation of this file.

#ifndef ELEMENTBOUNDARYTRACES
#define ELEMENTBOUNDARYTRACES

#include <vector>
#include "AuxDefs.h"
#include "Element.h"

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

  virtual size_t getNumTraceFaces() const = 0;

  virtual const std::vector<size_t> & getTraceFaceIds() const = 0;


  virtual size_t getTraceNumberOfFace(size_t FaceIndex) const = 0;
  
  virtual const Element & getTrace(size_t FaceIndex) const = 0;
  
  inline const Element & operator[](size_t FaceIndex) 
    { return getTrace(FaceIndex); }

  virtual const std::vector<double> & getNormal(size_t FaceIndex) const = 0;

  virtual size_t dofMap(size_t FaceIndex, size_t field, size_t dof) const = 0;
};




class ElementBoundaryTraces_: public ElementBoundaryTraces {
private:
  void copy (const ElementBoundaryTraces_& that) {
    FaceNumbers = that.FaceNumbers;

    for (size_t i = 0; i < that.FaceElements.size (); i++) {
          FaceElements.push_back (that.FaceElements[i]->clone ());
    }
  }

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

public:
  ElementBoundaryTraces_ () {
  }

  virtual ~ElementBoundaryTraces_ () {
    destroy();
  }

  ElementBoundaryTraces_ (const ElementBoundaryTraces_ & that) :
    ElementBoundaryTraces (that) {

    copy (that);
  }

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

  virtual ElementBoundaryTraces_ * clone () const = 0;

  size_t getNumTraceFaces () const {
    return FaceElements.size ();
  }

  const std::vector<size_t> & getTraceFaceIds () const {
    return FaceNumbers;
  }

  inline size_t getTraceNumberOfFace (size_t FaceIndex) const {
    for (size_t i = 0; i < FaceNumbers.size (); i++) {
      if (FaceNumbers[i] == FaceIndex) {
        return i;
      }
    }
    return -1;
  }

  virtual const Element & getTrace (size_t FaceIndex) const {
    return *FaceElements[FaceIndex];
  }

protected:
  void addFace (const Element_ * NewFace, const size_t FaceNumber) {
    FaceElements.push_back (NewFace->clone ());
    FaceNumbers.push_back (FaceNumber);
  }

private:
  std::vector<const Element *> FaceElements;
  std::vector<size_t> FaceNumbers;
};

#endif

---