Home Research Papers Posters Manuals Talks Download People Help Internal

libs/ck-libs/datatransfer/GenericElement.h

Go to the documentation of this file.

#ifndef _CSAR_GenericElement_H_
#define _CSAR_GenericElement_H_

#include "ckvector3d.h"
typedef CkVector3d CVector;
typedef CkVector3d CPoint;

class ConcreteElement {
public:
    virtual CPoint getNodeLocation(int i) const =0;
};

class ConcreteElementNodeData : public ConcreteElement {
public:
    virtual const double *getNodeData(int i) const =0;
};

class GenericElement {
protected:
  enum {maxSize=20}; // Maximum number of nodes per element we can handle
  unsigned int _size; 
public:
  GenericElement(unsigned int s = 4)
    : _size(s)
  {};
  unsigned int size() const 
  {
    return(_size);
  };
  unsigned int nedges() const
  {
    switch(_size){
    case 4:
    case 10:
      return (6);
    case 8:
    case 20:
      return (12);
    default:
      return(0);
    }
    return(0);
  };
  unsigned int nfaces() const
  {
    switch(_size){
    case 4:
    case 10:
      return(4);
    case 8:
    case 20:
      return(6);
    default:
      return(0);
    }
    return(0);
  };
  void shape_func(const CVector &natc,
          double []) const;
  void dshape_func(const CVector &natc,
           double [][3]) const;
  void jacobian(const CPoint p[],
        const CVector &natc,
        CVector J[]) const;
  
  void interpolate_natural(int nValuesPerNode,
           const ConcreteElementNodeData &src, // Source element
           const CVector &nc,
           double *dest) const;
  
  bool interpolate(int nValuesPerNode,
           const ConcreteElementNodeData &src, // Source element
           const CPoint &p,
           double *dest) const
  {
    CVector natc;
    if (!element_contains_point(p,src,natc)) return false;
    interpolate_natural(nValuesPerNode,src,natc,dest);
    return true;
  }

  bool element_contains_point(const CPoint &p, //    Target Mesh point
                          const ConcreteElement &e, // Source element
                  CVector &natc) const; // Returns Target point natural coords

  void shapef_jacobian_at(const CPoint &p,
              CVector &natc,
                      const ConcreteElement &e, // Source element
              CVector &fvec,CVector fjac[]) const;
  
};

#endif

---