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
|
/*****************************************************************************
*
* Copyright (c) 2003-2020 by The University of Queensland
* http://www.uq.edu.au
*
* Primary Business: Queensland, Australia
* Licensed under the Apache License, version 2.0
* http://www.apache.org/licenses/LICENSE-2.0
*
* Development until 2012 by Earth Systems Science Computational Center (ESSCC)
* Development 2012-2013 by School of Earth Sciences
* Development from 2014-2017 by Centre for Geoscience Computing (GeoComp)
* Development from 2019 by School of Earth and Environmental Sciences
**
*****************************************************************************/
/****************************************************************************
Finley: integration schemes for element shapes Tri, Quad, Hex, Tet, Line,
Point
*****************************************************************************/
#ifndef __FINLEY_QUADRATURE_H__
#define __FINLEY_QUADRATURE_H__
#include "Finley.h"
#define MAX_numQuadNodesLine 10
namespace finley {
typedef enum {
PointQuad,
LineQuad,
TriQuad,
RecQuad,
TetQuad,
HexQuad,
NoQuad // marks end of list
} QuadTypeId;
typedef void (Quad_getNodes) (int, std::vector<double>&, std::vector<double>&);
typedef int (Quad_getNumNodes) (int);
typedef int (Quad_getMacro) (int numSubElements, int numQuadNodes,
const double* quadNodes,
const double* quadWeights,
int numF, const double* dFdv,
int new_len, double* new_quadNodes,
double* new_quadWeights, double* new_dFfv);
struct QuadInfo {
/// quadrature type id
QuadTypeId TypeId;
/// the name in text form e.g. "Line", "Rec", ...
const char* Name;
/// number of spatial dimensions
int numDim;
/// number of vertices of the element
int numVertices;
/// function that returns the quadrature points for a given order
Quad_getNodes* getQuadNodes;
/// function that returns the number of quadrature nodes for a given
/// accuracy order
Quad_getNumNodes* getNumQuadNodes;
/// transfers a given quadrature scheme to a macro element structure
Quad_getMacro *getMacro;
};
/****** Interfaces ******/
Quad_getMacro Quad_MacroPoint;
Quad_getMacro Quad_MacroLine;
Quad_getMacro Quad_MacroTri;
Quad_getMacro Quad_MacroRec;
Quad_getMacro Quad_MacroTet;
Quad_getMacro Quad_MacroHex;
Quad_getNodes Quad_getNodesTri;
Quad_getNodes Quad_getNodesTet;
Quad_getNodes Quad_getNodesRec;
Quad_getNodes Quad_getNodesHex;
Quad_getNodes Quad_getNodesLine;
Quad_getNodes Quad_getNodesPoint;
Quad_getNodes Quad_getNodesTriOnFace;
Quad_getNodes Quad_getNodesRecOnFace;
Quad_getNodes Quad_getNodesLineOnFace;
Quad_getNodes Quad_getNodesPointOnFace;
Quad_getNodes Quad_getNodesTriMacro;
Quad_getNodes Quad_getNodesTetMacro;
Quad_getNodes Quad_getNodesRecMacro;
Quad_getNodes Quad_getNodesHexMacro;
Quad_getNodes Quad_getNodesLineMacro;
Quad_getNumNodes Quad_getNumNodesPoint;
Quad_getNumNodes Quad_getNumNodesLine;
Quad_getNumNodes Quad_getNumNodesTri;
Quad_getNumNodes Quad_getNumNodesRec;
Quad_getNumNodes Quad_getNumNodesTet;
Quad_getNumNodes Quad_getNumNodesHex;
void Quad_makeNodesOnFace(int, int, double*, double*, Quad_getNodes);
const QuadInfo* QuadInfo_getInfo(QuadTypeId id);
} // namespace finley
#endif // __FINLEY_QUADRATURE_H__
|
