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
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
|
<center><a href="https://gitlab.com/petsc/petsc/-/blob/966382dc56242773704ef5f5cee7aa2db3ebc577/include/petscfe.h">Actual source code: petscfe.h</a></center><br>
<html>
<head>
<title></title>
<meta name="generator" content="c2html 0.9.6">
<meta name="date" content="2025-04-30T18:14:50+00:00">
</head>
<body bgcolor="#FFFFFF">
<pre width=80>
<a name="line1"> 1: </a><font color="#B22222">/*</font>
<a name="line2"> 2: </a><font color="#B22222"> Objects which encapsulate finite element spaces and operations</font>
<a name="line3"> 3: </a><font color="#B22222">*/</font>
<a name="line4"> 4: </a><font color="#A020F0">#pragma once</font>
<a name="line5"> 5: </a>#include <A href="../include/petscdm.h.html"><petscdm.h></A>
<a name="line6"> 6: </a>#include <A href="../include/petscdt.h.html"><petscdt.h></A>
<a name="line7"> 7: </a>#include <A href="../include/petscfetypes.h.html"><petscfetypes.h></A>
<a name="line8"> 8: </a>#include <A href="../include/petscdstypes.h.html"><petscdstypes.h></A>
<a name="line9"> 9: </a>#include <A href="../include/petscspace.h.html"><petscspace.h></A>
<a name="line10"> 10: </a>#include <A href="../include/petscdualspace.h.html"><petscdualspace.h></A>
<a name="line12"> 12: </a><font color="#B22222">/* MANSEC = <a href="../manualpages/DM/DM.html">DM</a> */</font>
<a name="line13"> 13: </a><font color="#B22222">/* SUBMANSEC = FE */</font>
<a name="line15"> 15: </a><font color="#B22222">/*MC</font>
<a name="line16"> 16: </a><font color="#B22222"> <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> - Structure for geometric information for `<a href="../manualpages/FE/PetscFE.html">PetscFE</a>`</font>
<a name="line18"> 18: </a><font color="#B22222"> Level: intermediate</font>
<a name="line20"> 20: </a><font color="#B22222"> Note:</font>
<a name="line21"> 21: </a><font color="#B22222"> This is a struct, not a `<a href="../manualpages/Sys/PetscObject.html">PetscObject</a>`</font>
<a name="line23"> 23: </a><font color="#B22222">.seealso: `<a href="../manualpages/FE/PetscFE.html">PetscFE</a>`, `<a href="../manualpages/FE/PetscFEGeomCreate.html">PetscFEGeomCreate</a>()`, `<a href="../manualpages/FE/PetscFEGeomDestroy.html">PetscFEGeomDestroy</a>()`, `<a href="../manualpages/FE/PetscFEGeomGetChunk.html">PetscFEGeomGetChunk</a>()`, `<a href="../manualpages/FE/PetscFEGeomRestoreChunk.html">PetscFEGeomRestoreChunk</a>()`, `<a href="../manualpages/FE/PetscFEGeomGetPoint.html">PetscFEGeomGetPoint</a>()`, `<a href="../manualpages/FE/PetscFEGeomGetCellPoint.html">PetscFEGeomGetCellPoint</a>()`,</font>
<a name="line24"> 24: </a><font color="#B22222"> `<a href="../manualpages/FE/PetscFEGeomComplete.html">PetscFEGeomComplete</a>()`, `<a href="../manualpages/SPACE/PetscSpace.html">PetscSpace</a>`, `<a href="../manualpages/DUALSPACE/PetscDualSpace.html">PetscDualSpace</a>`</font>
<a name="line25"> 25: </a><font color="#B22222">M*/</font>
<a name="line26"> 26: </a><font color="#4169E1">typedef</font> <font color="#4169E1">struct</font> {
<a name="line27"> 27: </a> <font color="#B22222">// We can represent several different types of geometry, which we call modes:</font>
<a name="line28"> 28: </a> <font color="#B22222">// basic: dim == dE, only bulk data</font>
<a name="line29"> 29: </a> <font color="#B22222">// These are normal dim-cells</font>
<a name="line30"> 30: </a> <font color="#B22222">// embedded: dim < dE, only bulk data</font>
<a name="line31"> 31: </a> <font color="#B22222">// These are dim-cells embedded in a higher dimension, as an embedded manifold</font>
<a name="line32"> 32: </a> <font color="#B22222">// boundary: dim < dE, bulk and face data</font>
<a name="line33"> 33: </a> <font color="#B22222">// These are dim-cells on the boundary of a dE-mesh</font>
<a name="line34"> 34: </a> <font color="#B22222">// cohesive: dim < dE, bulk and face data</font>
<a name="line35"> 35: </a> <font color="#B22222">// These are dim-cells in the interior of a dE-mesh</font>
<a name="line36"> 36: </a> <font color="#B22222">// affine:</font>
<a name="line37"> 37: </a> <font color="#B22222">// For all modes, the transforms between real and reference are affine</font>
<a name="line38"> 38: </a> <a href="../manualpages/FE/PetscFEGeomMode.html">PetscFEGeomMode</a> mode; <font color="#B22222">// The type of geometric data stored</font>
<a name="line39"> 39: </a> <a href="../manualpages/Sys/PetscBool.html">PetscBool</a> isAffine; <font color="#B22222">// Flag for affine transforms</font>
<a name="line40"> 40: </a> <font color="#B22222">// Sizes</font>
<a name="line41"> 41: </a> <a href="../manualpages/Sys/PetscInt.html">PetscInt</a> dim; <font color="#B22222">// dim: topological dimension and reference coordinate dimension</font>
<a name="line42"> 42: </a> <a href="../manualpages/Sys/PetscInt.html">PetscInt</a> dimEmbed; <font color="#B22222">// dE: real coordinate dimension</font>
<a name="line43"> 43: </a> <a href="../manualpages/Sys/PetscInt.html">PetscInt</a> numCells; <font color="#B22222">// Nc: Number of mesh points represented in the arrays (points are assumed to be the same <a href="../manualpages/DM/DMPolytopeType.html">DMPolytopeType</a>)</font>
<a name="line44"> 44: </a> <a href="../manualpages/Sys/PetscInt.html">PetscInt</a> numPoints; <font color="#B22222">// Np: Number of evaluation points represented in the arrays</font>
<a name="line45"> 45: </a> <font color="#B22222">// Bulk data</font>
<a name="line46"> 46: </a> const <a href="../manualpages/Sys/PetscReal.html">PetscReal</a> *xi; <font color="#B22222">// xi[dim] The first point in each cell in reference coordinates</font>
<a name="line47"> 47: </a> <a href="../manualpages/Sys/PetscReal.html">PetscReal</a> *v; <font color="#B22222">// v[Nc*Np*dE]: The first point in each cell in real coordinates</font>
<a name="line48"> 48: </a> <a href="../manualpages/Sys/PetscReal.html">PetscReal</a> *J; <font color="#B22222">// J[Nc*Np*dE*dE]: The Jacobian of the map from reference to real coordinates (if nonsquare it is completed with orthogonal columns)</font>
<a name="line49"> 49: </a> <a href="../manualpages/Sys/PetscReal.html">PetscReal</a> *invJ; <font color="#B22222">// invJ[Nc*Np*dE*dE]: The inverse of the Jacobian of the map from reference to real coordinates (if nonsquare it is completed with orthogonal columns)</font>
<a name="line50"> 50: </a> <a href="../manualpages/Sys/PetscReal.html">PetscReal</a> *detJ; <font color="#B22222">// detJ[Nc*Np]: The determinant of J, and if J is non-square it is the volume change</font>
<a name="line51"> 51: </a> <font color="#B22222">// Face data</font>
<a name="line52"> 52: </a> <a href="../manualpages/Sys/PetscReal.html">PetscReal</a> *n; <font color="#B22222">// n[Nc*Np*dE]: For faces, the normal to the face in real coordinates, outward for the first supporting cell</font>
<a name="line53"> 53: </a> <a href="../manualpages/Sys/PetscInt.html">PetscInt</a> (*face)[4]; <font color="#B22222">// face[Nc][s*2]: For faces, the local face number (cone index) and orientation for this face in each supporting cell</font>
<a name="line54"> 54: </a> <a href="../manualpages/Sys/PetscReal.html">PetscReal</a> *suppJ[2]; <font color="#B22222">// sJ[s][Nc*Np*dE*dE]: For faces, the Jacobian for each supporting cell</font>
<a name="line55"> 55: </a> <a href="../manualpages/Sys/PetscReal.html">PetscReal</a> *suppInvJ[2]; <font color="#B22222">// sInvJ[s][Nc*Np*dE*dE]: For faces, the inverse Jacobian for each supporting cell</font>
<a name="line56"> 56: </a> <a href="../manualpages/Sys/PetscReal.html">PetscReal</a> *suppDetJ[2]; <font color="#B22222">// sdetJ[s][Nc*Np]: For faces, the Jacobian determinant for each supporting cell</font>
<a name="line57"> 57: </a>} <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a>;
<a name="line59"> 59: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/DM/PetscFEInitializePackage.html">PetscFEInitializePackage</a>(void)</font></strong>;
<a name="line61"> 61: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEGeomCreate.html">PetscFEGeomCreate</a>(<a href="../manualpages/DT/PetscQuadrature.html">PetscQuadrature</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/FE/PetscFEGeomMode.html">PetscFEGeomMode</a>, <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> **)</font></strong>;
<a name="line62"> 62: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEGeomGetChunk.html">PetscFEGeomGetChunk</a>(<a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> **)</font></strong>;
<a name="line63"> 63: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEGeomRestoreChunk.html">PetscFEGeomRestoreChunk</a>(<a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> **)</font></strong>;
<a name="line64"> 64: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEGeomGetPoint.html">PetscFEGeomGetPoint</a>(<a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, const <a href="../manualpages/Sys/PetscReal.html">PetscReal</a>[], <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *)</font></strong>;
<a name="line65"> 65: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEGeomGetCellPoint.html">PetscFEGeomGetCellPoint</a>(<a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *)</font></strong>;
<a name="line66"> 66: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEGeomComplete.html">PetscFEGeomComplete</a>(<a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *)</font></strong>;
<a name="line67"> 67: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEGeomDestroy.html">PetscFEGeomDestroy</a>(<a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> **)</font></strong>;
<a name="line69"> 69: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/DUALSPACE/PetscDualSpaceApply.html">PetscDualSpaceApply</a>(<a href="../manualpages/DUALSPACE/PetscDualSpace.html">PetscDualSpace</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscReal.html">PetscReal</a>, <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> (*)(<a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscReal.html">PetscReal</a>, const <a href="../manualpages/Sys/PetscReal.html">PetscReal</a>[], <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a> *, void *), void *, <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a> *)</font></strong>;
<a name="line70"> 70: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/DUALSPACE/PetscDualSpaceApplyDefault.html">PetscDualSpaceApplyDefault</a>(<a href="../manualpages/DUALSPACE/PetscDualSpace.html">PetscDualSpace</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscReal.html">PetscReal</a>, <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> (*)(<a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscReal.html">PetscReal</a>, const <a href="../manualpages/Sys/PetscReal.html">PetscReal</a>[], <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a> *, void *), void *, <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a> *)</font></strong>;
<a name="line72"> 72: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/DUALSPACE/PetscDualSpaceTransform.html">PetscDualSpaceTransform</a>(<a href="../manualpages/DUALSPACE/PetscDualSpace.html">PetscDualSpace</a>, <a href="../manualpages/DUALSPACE/PetscDualSpaceTransformType.html">PetscDualSpaceTransformType</a>, <a href="../manualpages/Sys/PetscBool.html">PetscBool</a>, <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[])</font></strong>;
<a name="line73"> 73: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/DUALSPACE/PetscDualSpaceTransformGradient.html">PetscDualSpaceTransformGradient</a>(<a href="../manualpages/DUALSPACE/PetscDualSpace.html">PetscDualSpace</a>, <a href="../manualpages/DUALSPACE/PetscDualSpaceTransformType.html">PetscDualSpaceTransformType</a>, <a href="../manualpages/Sys/PetscBool.html">PetscBool</a>, <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[])</font></strong>;
<a name="line74"> 74: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/DUALSPACE/PetscDualSpaceTransformHessian.html">PetscDualSpaceTransformHessian</a>(<a href="../manualpages/DUALSPACE/PetscDualSpace.html">PetscDualSpace</a>, <a href="../manualpages/DUALSPACE/PetscDualSpaceTransformType.html">PetscDualSpaceTransformType</a>, <a href="../manualpages/Sys/PetscBool.html">PetscBool</a>, <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[])</font></strong>;
<a name="line75"> 75: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/DUALSPACE/PetscDualSpacePullback.html">PetscDualSpacePullback</a>(<a href="../manualpages/DUALSPACE/PetscDualSpace.html">PetscDualSpace</a>, <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[])</font></strong>;
<a name="line76"> 76: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/DUALSPACE/PetscDualSpacePushforward.html">PetscDualSpacePushforward</a>(<a href="../manualpages/DUALSPACE/PetscDualSpace.html">PetscDualSpace</a>, <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[])</font></strong>;
<a name="line77"> 77: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/DUALSPACE/PetscDualSpacePushforwardGradient.html">PetscDualSpacePushforwardGradient</a>(<a href="../manualpages/DUALSPACE/PetscDualSpace.html">PetscDualSpace</a>, <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[])</font></strong>;
<a name="line78"> 78: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/DUALSPACE/PetscDualSpacePushforwardHessian.html">PetscDualSpacePushforwardHessian</a>(<a href="../manualpages/DUALSPACE/PetscDualSpace.html">PetscDualSpace</a>, <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[])</font></strong>;
<a name="line80"> 80: </a>PETSC_EXTERN <a href="../manualpages/Sys/PetscClassId.html">PetscClassId</a> PETSCFE_CLASSID;
<a name="line82"> 82: </a><font color="#B22222">/*J</font>
<a name="line83"> 83: </a><font color="#B22222"> <a href="../manualpages/FE/PetscFEType.html">PetscFEType</a> - String with the name of a PETSc finite element space</font>
<a name="line85"> 85: </a><font color="#B22222"> Level: beginner</font>
<a name="line87"> 87: </a><font color="#B22222"> Note:</font>
<a name="line88"> 88: </a><font color="#B22222"> Currently, the classes are concerned with the implementation of element integration</font>
<a name="line90"> 90: </a><font color="#B22222">.seealso: `<a href="../manualpages/FE/PetscFESetType.html">PetscFESetType</a>()`, `<a href="../manualpages/FE/PetscFE.html">PetscFE</a>`</font>
<a name="line91"> 91: </a><font color="#B22222">J*/</font>
<a name="line92"> 92: </a><font color="#4169E1">typedef const char *<a href="../manualpages/FE/PetscFEType.html">PetscFEType</a>;</font>
<a name="line93"> 93: </a><strong><font color="#228B22">#define <a href="../manualpages/FE/PETSCFEBASIC.html">PETSCFEBASIC</a> </font><font color="#666666">"basic"</font><font color="#228B22"></font></strong>
<a name="line94"> 94: </a><strong><font color="#228B22">#define <a href="../manualpages/FE/PETSCFEOPENCL.html">PETSCFEOPENCL</a> </font><font color="#666666">"opencl"</font><font color="#228B22"></font></strong>
<a name="line95"> 95: </a><strong><font color="#228B22">#define <a href="../manualpages/FE/PETSCFECOMPOSITE.html">PETSCFECOMPOSITE</a> </font><font color="#666666">"composite"</font><font color="#228B22"></font></strong>
<a name="line96"> 96: </a><strong><font color="#228B22">#define <a href="../manualpages/FE/PETSCFEVECTOR.html">PETSCFEVECTOR</a> </font><font color="#666666">"vector"</font><font color="#228B22"></font></strong>
<a name="line98"> 98: </a>PETSC_EXTERN <a href="../manualpages/Sys/PetscFunctionList.html">PetscFunctionList</a> PetscFEList;
<a name="line99"> 99: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFECreate.html">PetscFECreate</a>(<a href="../manualpages/Sys/MPI_Comm.html">MPI_Comm</a>, <a href="../manualpages/FE/PetscFE.html">PetscFE</a> *)</font></strong>;
<a name="line100">100: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEDestroy.html">PetscFEDestroy</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a> *)</font></strong>;
<a name="line101">101: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFESetType.html">PetscFESetType</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/FE/PetscFEType.html">PetscFEType</a>)</font></strong>;
<a name="line102">102: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEGetType.html">PetscFEGetType</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/FE/PetscFEType.html">PetscFEType</a> *)</font></strong>;
<a name="line103">103: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFESetUp.html">PetscFESetUp</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>)</font></strong>;
<a name="line104">104: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFESetFromOptions.html">PetscFESetFromOptions</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>)</font></strong>;
<a name="line105">105: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEViewFromOptions.html">PetscFEViewFromOptions</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/Sys/PetscObject.html">PetscObject</a>, const char[])</font></strong>;
<a name="line106">106: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFESetName.html">PetscFESetName</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, const char[])</font></strong>;
<a name="line107">107: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFECreateVector.html">PetscFECreateVector</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscBool.html">PetscBool</a>, <a href="../manualpages/Sys/PetscBool.html">PetscBool</a>, <a href="../manualpages/FE/PetscFE.html">PetscFE</a> *)</font></strong>;
<a name="line109">109: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEView.html">PetscFEView</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/Viewer/PetscViewer.html">PetscViewer</a>)</font></strong>;
<a name="line110">110: </a><strong><font color="#4169E1"><a name="PetscFERegister"></a>PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFERegister.html">PetscFERegister</a>(const char[], <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> (*)(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>)</font></strong>);
<a name="line111">111: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> PetscFERegisterDestroy(void)</font></strong>;
<a name="line112">112: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFECreateDefault.html">PetscFECreateDefault</a>(<a href="../manualpages/Sys/MPI_Comm.html">MPI_Comm</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscBool.html">PetscBool</a>, const char[], <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/FE/PetscFE.html">PetscFE</a> *)</font></strong>;
<a name="line113">113: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFECreateByCell.html">PetscFECreateByCell</a>(<a href="../manualpages/Sys/MPI_Comm.html">MPI_Comm</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/DM/DMPolytopeType.html">DMPolytopeType</a>, const char[], <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/FE/PetscFE.html">PetscFE</a> *)</font></strong>;
<a name="line114">114: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFECreateLagrange.html">PetscFECreateLagrange</a>(<a href="../manualpages/Sys/MPI_Comm.html">MPI_Comm</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscBool.html">PetscBool</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/FE/PetscFE.html">PetscFE</a> *)</font></strong>;
<a name="line115">115: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFECreateLagrangeByCell.html">PetscFECreateLagrangeByCell</a>(<a href="../manualpages/Sys/MPI_Comm.html">MPI_Comm</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/DM/DMPolytopeType.html">DMPolytopeType</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/FE/PetscFE.html">PetscFE</a> *)</font></strong>;
<a name="line116">116: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFECreateFromSpaces.html">PetscFECreateFromSpaces</a>(<a href="../manualpages/SPACE/PetscSpace.html">PetscSpace</a>, <a href="../manualpages/DUALSPACE/PetscDualSpace.html">PetscDualSpace</a>, <a href="../manualpages/DT/PetscQuadrature.html">PetscQuadrature</a>, <a href="../manualpages/DT/PetscQuadrature.html">PetscQuadrature</a>, <a href="../manualpages/FE/PetscFE.html">PetscFE</a> *)</font></strong>;
<a name="line117">117: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFELimitDegree.html">PetscFELimitDegree</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/FE/PetscFE.html">PetscFE</a> *)</font></strong>;
<a name="line119">119: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEGetDimension.html">PetscFEGetDimension</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a> *)</font></strong>;
<a name="line120">120: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEGetSpatialDimension.html">PetscFEGetSpatialDimension</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a> *)</font></strong>;
<a name="line121">121: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFESetNumComponents.html">PetscFESetNumComponents</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>)</font></strong>;
<a name="line122">122: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEGetNumComponents.html">PetscFEGetNumComponents</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a> *)</font></strong>;
<a name="line123">123: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEGetTileSizes.html">PetscFEGetTileSizes</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a> *, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a> *, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a> *, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a> *)</font></strong>;
<a name="line124">124: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFESetTileSizes.html">PetscFESetTileSizes</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>)</font></strong>;
<a name="line125">125: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFESetBasisSpace.html">PetscFESetBasisSpace</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/SPACE/PetscSpace.html">PetscSpace</a>)</font></strong>;
<a name="line126">126: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEGetBasisSpace.html">PetscFEGetBasisSpace</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/SPACE/PetscSpace.html">PetscSpace</a> *)</font></strong>;
<a name="line127">127: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFESetDualSpace.html">PetscFESetDualSpace</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/DUALSPACE/PetscDualSpace.html">PetscDualSpace</a>)</font></strong>;
<a name="line128">128: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEGetDualSpace.html">PetscFEGetDualSpace</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/DUALSPACE/PetscDualSpace.html">PetscDualSpace</a> *)</font></strong>;
<a name="line129">129: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFESetQuadrature.html">PetscFESetQuadrature</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/DT/PetscQuadrature.html">PetscQuadrature</a>)</font></strong>;
<a name="line130">130: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEGetQuadrature.html">PetscFEGetQuadrature</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/DT/PetscQuadrature.html">PetscQuadrature</a> *)</font></strong>;
<a name="line131">131: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFESetFaceQuadrature.html">PetscFESetFaceQuadrature</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/DT/PetscQuadrature.html">PetscQuadrature</a>)</font></strong>;
<a name="line132">132: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEGetFaceQuadrature.html">PetscFEGetFaceQuadrature</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/DT/PetscQuadrature.html">PetscQuadrature</a> *)</font></strong>;
<a name="line133">133: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFECopyQuadrature.html">PetscFECopyQuadrature</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/FE/PetscFE.html">PetscFE</a>)</font></strong>;
<a name="line134">134: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEGetNumDof.html">PetscFEGetNumDof</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, const <a href="../manualpages/Sys/PetscInt.html">PetscInt</a> **)</font></strong>;
<a name="line136">136: </a><font color="#B22222">/* TODO: Need a function to reuse the memory when retabulating the same FE at different points */</font>
<a name="line137">137: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEGetCellTabulation.html">PetscFEGetCellTabulation</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/DT/PetscTabulation.html">PetscTabulation</a> *)</font></strong>;
<a name="line138">138: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEGetFaceTabulation.html">PetscFEGetFaceTabulation</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/DT/PetscTabulation.html">PetscTabulation</a> *)</font></strong>;
<a name="line139">139: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEGetFaceCentroidTabulation.html">PetscFEGetFaceCentroidTabulation</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/DT/PetscTabulation.html">PetscTabulation</a> *)</font></strong>;
<a name="line140">140: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFECreateTabulation.html">PetscFECreateTabulation</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, const <a href="../manualpages/Sys/PetscReal.html">PetscReal</a>[], <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/DT/PetscTabulation.html">PetscTabulation</a> *)</font></strong>;
<a name="line141">141: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEComputeTabulation.html">PetscFEComputeTabulation</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, const <a href="../manualpages/Sys/PetscReal.html">PetscReal</a>[], <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/DT/PetscTabulation.html">PetscTabulation</a>)</font></strong>;
<a name="line142">142: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscTabulationDestroy.html">PetscTabulationDestroy</a>(<a href="../manualpages/DT/PetscTabulation.html">PetscTabulation</a> *)</font></strong>;
<a name="line144">144: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFERefine.html">PetscFERefine</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/FE/PetscFE.html">PetscFE</a> *)</font></strong>;
<a name="line145">145: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEGetHeightSubspace.html">PetscFEGetHeightSubspace</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/FE/PetscFE.html">PetscFE</a> *)</font></strong>;
<a name="line147">147: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> PetscFECreateCellGeometry(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/DT/PetscQuadrature.html">PetscQuadrature</a>, <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *)</font></strong>;
<a name="line148">148: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> PetscFEDestroyCellGeometry(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *)</font></strong>;
<a name="line149">149: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEPushforward.html">PetscFEPushforward</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[])</font></strong>;
<a name="line150">150: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEPushforwardGradient.html">PetscFEPushforwardGradient</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[])</font></strong>;
<a name="line151">151: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEPushforwardHessian.html">PetscFEPushforwardHessian</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[])</font></strong>;
<a name="line153">153: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEIntegrate.html">PetscFEIntegrate</a>(<a href="../manualpages/DT/PetscDS.html">PetscDS</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *, const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], <a href="../manualpages/DT/PetscDS.html">PetscDS</a>, const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[])</font></strong>;
<a name="line154">154: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEIntegrateBd.html">PetscFEIntegrateBd</a>(<a href="../manualpages/DT/PetscDS.html">PetscDS</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, void (*)(<a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, const <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>[], const <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>[], const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], const <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>[], const <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>[], const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], <a href="../manualpages/Sys/PetscReal.html">PetscReal</a>, const <a href="../manualpages/Sys/PetscReal.html">PetscReal</a>[], const <a href="../manualpages/Sys/PetscReal.html">PetscReal</a>[], <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[]), <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *, const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], <a href="../manualpages/DT/PetscDS.html">PetscDS</a>, const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[])</font></strong>;
<a name="line155">155: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEIntegrateResidual.html">PetscFEIntegrateResidual</a>(<a href="../manualpages/DT/PetscDS.html">PetscDS</a>, <a href="../manualpages/DT/PetscFormKey.html">PetscFormKey</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *, const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], <a href="../manualpages/DT/PetscDS.html">PetscDS</a>, const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], <a href="../manualpages/Sys/PetscReal.html">PetscReal</a>, <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[])</font></strong>;
<a name="line156">156: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEIntegrateBdResidual.html">PetscFEIntegrateBdResidual</a>(<a href="../manualpages/DT/PetscDS.html">PetscDS</a>, <a href="../manualpages/DT/PetscWeakForm.html">PetscWeakForm</a>, <a href="../manualpages/DT/PetscFormKey.html">PetscFormKey</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *, const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], <a href="../manualpages/DT/PetscDS.html">PetscDS</a>, const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], <a href="../manualpages/Sys/PetscReal.html">PetscReal</a>, <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[])</font></strong>;
<a name="line157">157: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEIntegrateHybridResidual.html">PetscFEIntegrateHybridResidual</a>(<a href="../manualpages/DT/PetscDS.html">PetscDS</a>, <a href="../manualpages/DT/PetscDS.html">PetscDS</a>, <a href="../manualpages/DT/PetscFormKey.html">PetscFormKey</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *, const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], <a href="../manualpages/DT/PetscDS.html">PetscDS</a>, const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], <a href="../manualpages/Sys/PetscReal.html">PetscReal</a>, <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[])</font></strong>;
<a name="line158">158: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEIntegrateJacobian.html">PetscFEIntegrateJacobian</a>(<a href="../manualpages/DT/PetscDS.html">PetscDS</a>, <a href="../manualpages/FE/PetscFEJacobianType.html">PetscFEJacobianType</a>, <a href="../manualpages/DT/PetscFormKey.html">PetscFormKey</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *, const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], <a href="../manualpages/DT/PetscDS.html">PetscDS</a>, const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], <a href="../manualpages/Sys/PetscReal.html">PetscReal</a>, <a href="../manualpages/Sys/PetscReal.html">PetscReal</a>, <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[])</font></strong>;
<a name="line159">159: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEIntegrateBdJacobian.html">PetscFEIntegrateBdJacobian</a>(<a href="../manualpages/DT/PetscDS.html">PetscDS</a>, <a href="../manualpages/DT/PetscWeakForm.html">PetscWeakForm</a>, <a href="../manualpages/FE/PetscFEJacobianType.html">PetscFEJacobianType</a>, <a href="../manualpages/DT/PetscFormKey.html">PetscFormKey</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *, const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], <a href="../manualpages/DT/PetscDS.html">PetscDS</a>, const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], <a href="../manualpages/Sys/PetscReal.html">PetscReal</a>, <a href="../manualpages/Sys/PetscReal.html">PetscReal</a>, <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[])</font></strong>;
<a name="line160">160: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEIntegrateHybridJacobian.html">PetscFEIntegrateHybridJacobian</a>(<a href="../manualpages/DT/PetscDS.html">PetscDS</a>, <a href="../manualpages/DT/PetscDS.html">PetscDS</a>, <a href="../manualpages/FE/PetscFEJacobianType.html">PetscFEJacobianType</a>, <a href="../manualpages/DT/PetscFormKey.html">PetscFormKey</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/FE/PetscFEGeom.html">PetscFEGeom</a> *, const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], <a href="../manualpages/DT/PetscDS.html">PetscDS</a>, const <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[], <a href="../manualpages/Sys/PetscReal.html">PetscReal</a>, <a href="../manualpages/Sys/PetscReal.html">PetscReal</a>, <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a>[])</font></strong>;
<a name="line162">162: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFECompositeGetMapping.html">PetscFECompositeGetMapping</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a> *, const <a href="../manualpages/Sys/PetscReal.html">PetscReal</a> *[], const <a href="../manualpages/Sys/PetscReal.html">PetscReal</a> *[], const <a href="../manualpages/Sys/PetscReal.html">PetscReal</a> *[])</font></strong>;
<a name="line164">164: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> PetscFECreateHeightTrace(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/FE/PetscFE.html">PetscFE</a> *)</font></strong>;
<a name="line165">165: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> PetscFECreatePointTrace(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/Sys/PetscInt.html">PetscInt</a>, <a href="../manualpages/FE/PetscFE.html">PetscFE</a> *)</font></strong>;
<a name="line167">167: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEOpenCLSetRealType.html">PetscFEOpenCLSetRealType</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/Sys/PetscDataType.html">PetscDataType</a>)</font></strong>;
<a name="line168">168: </a><strong><font color="#4169E1">PETSC_EXTERN <a href="../manualpages/Sys/PetscErrorCode.html">PetscErrorCode</a> <a href="../manualpages/FE/PetscFEOpenCLGetRealType.html">PetscFEOpenCLGetRealType</a>(<a href="../manualpages/FE/PetscFE.html">PetscFE</a>, <a href="../manualpages/Sys/PetscDataType.html">PetscDataType</a> *)</font></strong>;
<a name="line170">170: </a><font color="#A020F0">#ifdef PETSC_HAVE_LIBCEED</font>
<a name="line172">172: </a><font color="#A020F0"> #ifndef PLEXFE_QFUNCTION</font>
<a name="line173">173: </a><strong><font color="#228B22"> #define PLEXFE_QFUNCTION(fname, f0_name, f1_name) \</font></strong>
<a name="line174">174: </a><strong><font color="#228B22"> CEED_QFUNCTION(PlexQFunction##fname)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *const *out) \</font></strong>
<a name="line175">175: </a><strong><font color="#228B22"> { \</font></strong>
<a name="line176">176: </a><strong><font color="#228B22"> const CeedScalar *u = in[0], *du = in[1], *qdata = in[2]; \</font></strong>
<a name="line177">177: </a><strong><font color="#228B22"> CeedScalar *v = out[0], *dv = out[1]; \</font></strong>
<a name="line178">178: </a><strong><font color="#228B22"> const <a href="../manualpages/Sys/PetscInt.html">PetscInt</a> Nc = 1; \</font></strong>
<a name="line179">179: </a><strong><font color="#228B22"> const <a href="../manualpages/Sys/PetscInt.html">PetscInt</a> cdim = 2; \</font></strong>
<a name="line180">180: </a><strong><font color="#228B22">\</font></strong>
<a name="line181">181: </a><strong><font color="#228B22"> CeedPragmaSIMD for (CeedInt i = 0; i < Q; ++i) \</font></strong>
<a name="line182">182: </a><strong><font color="#228B22"> { \</font></strong>
<a name="line183">183: </a><strong><font color="#228B22"> const <a href="../manualpages/Sys/PetscInt.html">PetscInt</a> uOff[2] = {0, Nc}; \</font></strong>
<a name="line184">184: </a><strong><font color="#228B22"> const <a href="../manualpages/Sys/PetscInt.html">PetscInt</a> uOff_x[2] = {0, Nc * cdim}; \</font></strong>
<a name="line185">185: </a><strong><font color="#228B22"> const CeedScalar x[2] = {qdata[i + Q * 1], qdata[i + Q * 2]}; \</font></strong>
<a name="line186">186: </a><strong><font color="#228B22"> const CeedScalar invJ[2][2] = { \</font></strong>
<a name="line187">187: </a><strong><font color="#228B22"> {qdata[i + Q * 3], qdata[i + Q * 5]}, \</font></strong>
<a name="line188">188: </a><strong><font color="#228B22"> {qdata[i + Q * 4], qdata[i + Q * 6]} \</font></strong>
<a name="line189">189: </a><strong><font color="#228B22"> }; \</font></strong>
<a name="line190">190: </a><strong><font color="#228B22"> const CeedScalar u_x[2] = {invJ[0][0] * du[i + Q * 0] + invJ[1][0] * du[i + Q * 1], invJ[0][1] * du[i + Q * 0] + invJ[1][1] * du[i + Q * 1]}; \</font></strong>
<a name="line191">191: </a><strong><font color="#228B22"> <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a> f0[Nc]; \</font></strong>
<a name="line192">192: </a><strong><font color="#228B22"> <a href="../manualpages/Sys/PetscScalar.html">PetscScalar</a> f1[Nc * cdim]; \</font></strong>
<a name="line193">193: </a><strong><font color="#228B22">\</font></strong>
<a name="line194">194: </a><strong><font color="#228B22"> for (<a href="../manualpages/Sys/PetscInt.html">PetscInt</a> k = 0; k < Nc; ++k) f0[k] = 0; \</font></strong>
<a name="line195">195: </a><strong><font color="#228B22"> for (<a href="../manualpages/Sys/PetscInt.html">PetscInt</a> k = 0; k < Nc * cdim; ++k) f1[k] = 0; \</font></strong>
<a name="line196">196: </a><strong><font color="#228B22"> f0_name(2, 1, 0, uOff, uOff_x, u, NULL, u_x, NULL, NULL, NULL, NULL, NULL, 0.0, x, 0, NULL, f0); \</font></strong>
<a name="line197">197: </a><strong><font color="#228B22"> f1_name(2, 1, 0, uOff, uOff_x, u, NULL, u_x, NULL, NULL, NULL, NULL, NULL, 0.0, x, 0, NULL, f1); \</font></strong>
<a name="line198">198: </a><strong><font color="#228B22">\</font></strong>
<a name="line199">199: </a><strong><font color="#228B22"> dv[i + Q * 0] = qdata[i + Q * 0] * (invJ[0][0] * f1[0] + invJ[0][1] * f1[1]); \</font></strong>
<a name="line200">200: </a><strong><font color="#228B22"> dv[i + Q * 1] = qdata[i + Q * 0] * (invJ[1][0] * f1[0] + invJ[1][1] * f1[1]); \</font></strong>
<a name="line201">201: </a><strong><font color="#228B22"> v[i] = qdata[i + Q * 0] * f0[0]; \</font></strong>
<a name="line202">202: </a><strong><font color="#228B22"> } \</font></strong>
<a name="line203">203: </a><strong><font color="#228B22"> return CEED_ERROR_SUCCESS; \</font></strong>
<a name="line204">204: </a><strong><font color="#228B22"> }</font></strong>
<a name="line205">205: </a><font color="#A020F0"> #endif</font>
<a name="line207">207: </a><font color="#A020F0">#else</font>
<a name="line209">209: </a><font color="#A020F0"> #ifndef PLEXFE_QFUNCTION</font>
<a name="line210">210: </a><strong><font color="#228B22"> #define PLEXFE_QFUNCTION(fname, f0_name, f1_name)</font></strong>
<a name="line211">211: </a><font color="#A020F0"> #endif</font>
<a name="line213">213: </a><font color="#A020F0">#endif</font>
</pre>
</body>
</html>
|
