apfShape.h

Go to the documentation of this file.

/*
 * Copyright 2011 Scientific Computation Research Center
 *
 * This work is open source software, licensed under the terms of the
 * BSD license as described in the LICENSE file in the top-level directory.
 */
 
#ifndef APFSHAPE_H
#define APFSHAPE_H
 
#include "apf.h"
#include "apfNew.h"
#include <sstream>
 
namespace apf {
 
class Mesh;
class MeshEntity;
 
class EntityShape
{
  public:
    virtual ~EntityShape();
    virtual void getValues(
        Mesh* m,
        MeshEntity* e,
        Vector3 const& xi,
        NewArray<double>& values) const = 0;
    virtual void getLocalGradients(
        Mesh* m,
        MeshEntity* e,
        Vector3 const& xi,
        NewArray<Vector3>& grads) const = 0;
    virtual void getVectorValues(
        Mesh* m,
        MeshEntity* e,
        Vector3 const& xi,
        NewArray<apf::Vector3>& values) const = 0;
    virtual void getLocalVectorCurls(
        Mesh* m,
        MeshEntity* e,
        Vector3 const& xi,
        NewArray<apf::Vector3>& values) const;
    virtual int countNodes() const = 0;
    virtual void alignSharedNodes(Mesh* m,
        MeshEntity* elem, MeshEntity* shared, int order[]);
};
 
class FieldShape
{
  public:
    virtual ~FieldShape();
    virtual EntityShape* getEntityShape(int type) = 0;
    virtual bool hasNodesIn(int dimension) = 0;
    virtual int countNodesOn(int type) = 0;
    virtual int getOrder() = 0;
    virtual void getNodeXi(int type, int node, Vector3& xi);
    virtual void getNodeTangent(int type, int node, Vector3& t);
    virtual bool isVectorShape();
    virtual const char* getName() const = 0;
    void registerSelf(const char* name);
};
 
FieldShape* getLagrange(int order);
 
FieldShape* getSerendipity();
 
FieldShape* getConstant(int dimension);
FieldShape* getIPShape(int dimension, int order);
FieldShape* getVoronoiShape(int dimension, int order);
FieldShape* getIPFitShape(int dimension, int order);
 
FieldShape* getHierarchic(int order);
 
FieldShape* getNedelec(int order);
 
FieldShape* getL2Shape(int order, int type);
 
FieldShape* getH1Shape(int order);
 
void projectHierarchicField(Field* to, Field* from);
 
void projectNedelecField(Field* to, Field* from);
 
void projectL2Field(Field* to, Field* from);
 
FieldShape* getShapeByName(const char* name);
 
int countElementNodes(FieldShape* s, int type);
 
void getElementNodeXis(FieldShape* s, int type,
    apf::NewArray<apf::Vector3>& xis);
 
void getElementNodeXis(FieldShape* s, Mesh* m, MeshEntity* e,
    apf::NewArray<apf::Vector3>& xis);
 
Vector3 boundaryToElementXi(
    Mesh* m,
    MeshEntity* boundary,
    MeshEntity* element,
    Vector3 const& xi);
 
}
 
#endif