---

BasisFamilyBase.h

Go to the documentation of this file.

#ifndef BASISFAMILYBASE_H
#define BASISFAMILYBASE_H

#include "SundanceDefs.h"
#include "Cell.h"
#include "CellType.h"
#include "TSFArray.h"
#include "Point.h"
#include "DenseSerialVector.h"
#include "XMLObject.h"
#include "MultiIndex.h"
#include "TemporaryMapper.h"

namespace Sundance
{

  using namespace TSF;
  using std::string;

  using std::ostream;

  class Expr;

  /** \ingroup LowLevelFE
   * Base class for basis functions.
   */

  class BasisFamilyBase
    {
    public:
      /** Construct a basis family of a given order. */
      BasisFamilyBase(const string& name, int order);
      /** TUVD */
      virtual ~BasisFamilyBase() {;}

      /** return the order of the basis family */
      int order() const {return order_;}

      /** return the number of DOFs associated with this basis on
       * the given cell type */
      int nNodes(const CellTopologyCode& cellType) const ;

      /** return an identifying name for the basis family */
      const string& name() const {return name_;}

      /** return a hash code */
      int hashCode() const {return hashCode_;}

      /** return true if the subtype is a
       * vector basis (e.g., Raviart-Thomas) */
      virtual bool isVectorBasis() const {return false;}

      /** return true if the subtype is a scalar basis */
      virtual bool isScalarBasis() const {return false;}

      /** evaluate on a set of points given in reference coordinates */
      virtual void refEval(const CellTopologyCode& cellType,
                           const TSFArray<Point>& pts,
                           const MultiIndex& deriv,
                           TSFArray<DenseSerialVector>& result) const ;

      /**
       * assign DOF numbers for a set of function IDs on a specified cell
       */
      virtual void setNodeNumbers(const Cell& cell,
                                  const TSFArray<int>& funcID,
                                  TemporaryMapper& tmpMap) const = 0 ;

      /** get the DOF numbers for a set of function IDs on a specified
       * cell
       */
      virtual void getNodeNumbers(const Cell& cell,
                                  const TSFArray<int>& funcID,
                                  TemporaryMapper& tmpMap,
                                  TSFArray<TSFArray<int> >& indices) const = 0 ;

      /**
       * find the DOF list for one of the facets of a specified cell
       */
      virtual void getFacetDOFs(const Cell& c, int facetDim, int facetIndex,
                                const TSFArray<int>& maximalDOFs,
                                TSFArray<int>& facetDOFs) const = 0 ;

      /**
       * Find the points at the nodes. Some points may have been processed already
       * and are marked as inactive with the bool array. Inactive points are not
       * returned.
       */
      virtual void getNodalPoints(const Cell& cell,
                                  const TSFArray<bool>& activeNodes,
                                  TSFArray<Point>& values) const ;

      /**
       * Write to XML
       */
      virtual XMLObject toXML() const ;

    protected:
      /** check that a point's dimension is appropriate for
       * a given cell type */
      void pointCheck(const Point& pt, int expectedDim,
                      const string& cellType) const ;

      /** get number of nodes for a point cell */
      virtual int nNodesOnPoint() const ;

      /** get number of nodes for a line cell */
      virtual int nNodesOnLine() const ;

      /** get number of nodes for a triangle cell */
      virtual int nNodesOnTriangle() const ;

      /** get number of nodes for a quad cell  */
      virtual int nNodesOnQuad() const ;

      /** get number of nodes for a tet cell */
      virtual int nNodesOnTet() const ;

      /** get number of nodes for a hex cell */
      virtual int nNodesOnBrick() const ;

      /** evaluate on a point cell */
      virtual void evalOnPoint(const Point& pt,
                               const MultiIndex& deriv,
                               DenseSerialVector& result) const ;
      /** evaluate on a line cell  */
      virtual void evalOnLine(const Point& pt,
                              const MultiIndex& deriv,
                              DenseSerialVector& result) const ;

      /** evaluate on a triangle cell  */
      virtual void evalOnTriangle(const Point& pt,
                                  const MultiIndex& deriv,
                                  DenseSerialVector& result) const ;

      /** evaluate on a quad cell  */
      virtual void evalOnQuad(const Point& pt,
                              const MultiIndex& deriv,
                              DenseSerialVector& result) const ;

      /** evaluate on a tet cell  */
      virtual void evalOnTet(const Point& pt,
                             const MultiIndex& deriv,
                             DenseSerialVector& result) const ;
      /** evaluate on a brick cell  */
      virtual void evalOnBrick(const Point& pt,
                               const MultiIndex& deriv,
                               DenseSerialVector& result) const ;

      /** order of the basis function */
      int order_;
      /** an identifying name for the basis function type */
      string name_;
      /** hash code */
      int hashCode_;
    };
}

#endif

---

Contact:

Kevin Long (krlong@ca.sandia.gov)


Documentation generated by