crvBezierShapes.h

Go to the documentation of this file.

/*
 * Copyright 2015 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 CRVBEZIERSHAPES_H
#define CRVBEZIERSHAPES_H

namespace crv {

void BlendedTriangleGetValues(apf::Mesh* m, apf::MeshEntity* e,
    apf::Vector3 const& xi, apf::NewArray<double>& values);
void BlendedTriangleGetLocalGradients(apf::Mesh* m, apf::MeshEntity* e,
    apf::Vector3 const& xi, apf::NewArray<apf::Vector3>& grads);
void BlendedTetGetValues(apf::Mesh* m, apf::MeshEntity* e,
    apf::Vector3 const& xi, apf::NewArray<double>& values);
void BlendedTetGetLocalGradients(apf::Mesh* m, apf::MeshEntity* e,
    apf::Vector3 const& xi, apf::NewArray<apf::Vector3>& grads);

typedef void (*bezierShape)(int P,
    apf::Vector3 const& xi,
    apf::NewArray<double>& values);

typedef void (*bezierShapeGrads)(int P,
    apf::Vector3 const& xi,
    apf::NewArray<apf::Vector3>& grads);

extern const bezierShape bezier[apf::Mesh::TYPES];
extern const bezierShapeGrads bezierGrads[apf::Mesh::TYPES];

void getBezierTransformationMatrix(int type, int P,
    mth::Matrix<double>& A,
    const apf::Vector3 *range);
void getBezierTransformationMatrix(int parentType,
    int childType, int P,
    mth::Matrix<double>& A,
    const apf::Vector3* childRange);

}

#endif