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 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 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);

void projectHierarchicField(Field* to, Field* from);

FieldShape* getShapeByName(const char* name);

int countElementNodes(FieldShape* s, int type);

Vector3 boundaryToElementXi(
    Mesh* m,
    MeshEntity* boundary,
    MeshEntity* element,
    Vector3 const& xi);

}

#endif