crvBezier.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 CRVBEZIER_H
#define CRVBEZIER_H

#include "crv.h"
#include "mth.h"

namespace crv {

double Bij(const int i, const int j,const double u, const double v);
double Bijk(const int i, const int j, const int k, const double u,
    const double v, const double w);
double Bijkl(const int i, const int j, const int k, const int l,
    const double u, const double v, const double w, const double t);

double Bij(const int ij[], const double xi[]);
double Bijk(const int ijk[], const double xi[]);
double Bijkl(const int ijkl[], const double xi[]);

int computeTriNodeIndex(int P, int i, int j);
int computeTetNodeIndex(int P, int i, int j, int k);

int getNumControlPoints(int type, int order);
int getNumInternalControlPoints(int type, int order);

void getTriNodesFromTetNodes(int f, int P,
    apf::NewArray<apf::Vector3>& tetNodes,
    apf::NewArray<apf::Vector3>& triNodes);
void getTriDetJacNodesFromTetDetJacNodes(int f, int P,
    apf::NewArray<double>& tetNodes,
    apf::NewArray<double>& triNodes);
void getFullRepFromBlended(int type,
    apf::NewArray<double>& transformCoefficients,
    apf::NewArray<apf::Vector3>& elemNodes);
double computeTriJacobianDetFromBezierFormulation(apf::Mesh* m,
    apf::MeshEntity* e, apf::Vector3& xi);
double computeTetJacobianDetFromBezierFormulation(apf::Mesh* m,
    apf::MeshEntity* e, apf::Vector3& xi);

void getBezierNodeXi(int type, int P, int node, apf::Vector3& xi);
void collectNodeXi(int parentType, int childType, int P,
    const apf::Vector3* range, apf::NewArray<apf::Vector3>& xi);
void elevateBezierEdge(int P, int r, apf::NewArray<apf::Vector3>& nodes,
    apf::NewArray<apf::Vector3>& elevatedNodes);
void elevateBezierTriangle(int P, int r, apf::NewArray<apf::Vector3>& nodes,
    apf::NewArray<apf::Vector3>& elevatedNodes);
void elevateBezierTet(int P, int r, apf::NewArray<apf::Vector3>& nodes,
    apf::NewArray<apf::Vector3>& elevatedNodes);
void elevateBezier(int type, int P, int r,
    apf::NewArray<apf::Vector3>& nodes,
    apf::NewArray<apf::Vector3>& elevatedNodes);

void elevateBezierCurve(apf::Mesh2* m, apf::MeshEntity* edge, int n, int r);

void subdivideBezierEdge(int P, double t, apf::NewArray<apf::Vector3>& nodes,
    apf::NewArray<apf::Vector3> (&subNodes)[2]);
void subdivideBezierTriangle(int P, apf::Vector3& p,
    apf::NewArray<apf::Vector3>& nodes,
    apf::NewArray<apf::Vector3> (&subNodes)[3]);
void subdivideBezierTriangle(int P, apf::NewArray<apf::Vector3>& nodes,
    apf::NewArray<apf::Vector3> (&subNodes)[4]);
void subdivideBezierTet(int P, apf::Vector3& p,
    apf::NewArray<apf::Vector3>& nodes,
    apf::NewArray<apf::Vector3> (&subNodes)[4]);

void convertInterpolationPoints(int n, int ne,
    apf::NewArray<apf::Vector3>& nodes,
    apf::NewArray<double>& c,
    apf::NewArray<apf::Vector3>& newNodes);

void convertInterpolationPoints(apf::Mesh2* m, apf::MeshEntity* e,
    int n, int ne, apf::NewArray<double>& c);

void getBezierTransformationCoefficients(int P, int type,
    apf::NewArray<double>& c);
void getInternalBezierTransformationCoefficients(apf::Mesh* m, int P, int blend,
    int type, apf::NewArray<double>& c);
void getGregoryTransformationCoefficients(int type, apf::NewArray<double>& c);
void getGregoryBlendedTransformationCoefficients(int blend, int type,
    apf::NewArray<double>& c);

void snapToInterpolate(apf::Mesh2* m, apf::MeshEntity* e, bool isNew = false);

void getTransformationMatrix(apf::Mesh* m, apf::MeshEntity* e,
    mth::Matrix<double>& A, const apf::Vector3 *range);

}

#endif