/usr/include/gmsh/bezierBasis.h is in libgmsh-dev 2.15.0+dfsg1-3.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 | // Gmsh - Copyright (C) 1997-2016 C. Geuzaine, J.-F. Remacle
//
// See the LICENSE.txt file for license information. Please report all
// bugs and problems to the public mailing list <gmsh@onelab.info>.
#ifndef _BEZIER_BASIS_H_
#define _BEZIER_BASIS_H_
#include <vector>
#include "fullMatrix.h"
#include "FuncSpaceData.h"
class MElement;
class bezierBasisRaiser;
class bezierBasis {
private :
// the 'numLagCoeff' first exponents are related to 'real' values
int _numLagCoeff;
int _numDivisions, _dimSimplex;
const FuncSpaceData _data;
bezierBasisRaiser *_raiser;
friend class bezierBasisRaiser;
fullMatrix<double> _exponents;
public :
fullMatrix<double> matrixLag2Bez;
fullMatrix<double> matrixBez2Lag;
fullMatrix<double> subDivisor;
// Constructors
bezierBasis(FuncSpaceData data);
~bezierBasis();
// get methods
inline int getDim() const {return _exponents.size2();}
inline int getType() const {return _data.elementType();}
inline int getOrder() const {return _data.spaceOrder();}
inline int getDimSimplex() const {return _dimSimplex;}
inline int getNumLagCoeff() const {return _numLagCoeff;}
inline int getNumDivision() const {return _numDivisions;}
inline int getNumSubNodes() const {return subDivisor.size1();}
inline FuncSpaceData getFuncSpaceData() const {return _data;}
bezierBasisRaiser* getRaiser() const;
// generate Bezier points
void generateBezierPoints(fullMatrix<double>&) const;
// transform coeff Lagrange into Bezier coeff
void lag2Bez(const fullMatrix<double> &lag, fullMatrix<double> &bez) const;
// Subdivide Bezier coefficients
void subdivideBezCoeff(const fullMatrix<double> &coeff,
fullMatrix<double> &subCoeff) const;
void subdivideBezCoeff(const fullVector<double> &coeff,
fullVector<double> &subCoeff) const;
// Interpolation of n functions on N points :
// coeffs(numCoeff, n) and uvw(N, dim)
// => result(N, n)
void interpolate(const fullMatrix<double> &coeffs,
const fullMatrix<double> &uvw,
fullMatrix<double> &result,
bool bezCoord = false) const;
void interpolate(const fullVector<double> &coeffs,
const fullMatrix<double> &uvw,
fullVector<double> &result,
bool bezCoord = false) const {
int size = uvw.size1();
result.resize(size);
fullMatrix<double> c(const_cast<double*>(coeffs.getDataPtr()), size, 1);
fullMatrix<double> r(const_cast<double*>(result.getDataPtr()), size, 1);
interpolate(c, uvw, r, bezCoord);
}
private :
void _construct();
void _constructPyr();
void _FEpoints2BezPoints(fullMatrix<double>&) const;
};
class bezierBasisRaiser {
// Let f, g, h be three function whose Bezier coefficients are given.
// This class allows to compute the Bezier coefficients of f*g and f*g*h.
private :
class _Data {
friend class bezierBasisRaiser;
private:
int i, j, k;
double val;
public:
_Data(double vv, int ii, int jj = -1, int kk = -1) :
i(ii), j(jj), k(kk), val(vv) {}
};
std::vector<std::vector<_Data> > _raiser2, _raiser3;
const bezierBasis *_bfs;
public:
bezierBasisRaiser(const bezierBasis *bezier) : _bfs(bezier) {
_fillRaiserData();
};
const bezierBasis* getRaisedBezierBasis(int multipliedBy2or3);
// const bezierBasis* getRaisedBezierBasis(int raised) const;
void computeCoeff(const fullVector<double> &coeffA,
const fullVector<double> &coeffB,
fullVector<double> &coeffSquare);
void computeCoeff(const fullMatrix<double> &coeffA,
const fullMatrix<double> &coeffB,
fullMatrix<double> &coeffSquare);
void computeCoeff(const fullVector<double> &coeffA,
const fullVector<double> &coeffB,
const fullVector<double> &coeffC,
fullVector<double> &coeffCubic);
void computeCoeff(const fullVector<double> &coeffA,
const fullMatrix<double> &coeffB,
const fullMatrix<double> &coeffC,
fullMatrix<double> &coeffCubic);
private:
void _fillRaiserData();
void _fillRaiserDataPyr();
};
#endif
|