-
Notifications
You must be signed in to change notification settings - Fork 2
/
polynomialBasis.h
44 lines (36 loc) · 1.35 KB
/
polynomialBasis.h
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
// Gmsh - Copyright (C) 1997-2018 C. Geuzaine, J.-F. Remacle
//
// See the LICENSE.txt file for license information. Please report all
// issues on https://gitlab.onelab.info/gmsh/gmsh/issues
#ifndef _POLYNOMIAL_BASIS_H_
#define _POLYNOMIAL_BASIS_H_
#include <math.h>
#include <map>
#include <vector>
#include "fullMatrix.h"
#include "nodalBasis.h"
#include <iostream>
class polynomialBasis : public nodalBasis {
public:
// for now the only implemented polynomial basis are nodal poly basis, we use
// the type of the corresponding gmsh element as type
fullMatrix<double> monomials;
fullMatrix<double> coefficients;
polynomialBasis(){};
polynomialBasis(int tag);
~polynomialBasis();
virtual inline int getNumShapeFunctions() const
{
return coefficients.size1();
}
virtual void f(double u, double v, double w, double *sf) const;
virtual void f(const fullMatrix<double> &coord, fullMatrix<double> &sf) const;
virtual void df(const fullMatrix<double> &coord,
fullMatrix<double> &dfm) const;
virtual void df(double u, double v, double w, double grads[][3]) const;
virtual void ddf(double u, double v, double w, double hess[][3][3]) const;
virtual void dddf(double u, double v, double w,
double third[][3][3][3]) const;
void evaluateMonomials(double u, double v, double w, double p[]) const;
};
#endif