Loading...
Searching...
No Matches
reconstruction.hpp
1/*---------------------------------------------------------------------------*\
2 *
3 * bitpit
4 *
5 * Copyright (C) 2015-2021 OPTIMAD engineering Srl
6 *
7 * -------------------------------------------------------------------------
8 * License
9 * This file is part of bitpit.
10 *
11 * bitpit is free software: you can redistribute it and/or modify it
12 * under the terms of the GNU Lesser General Public License v3 (LGPL)
13 * as published by the Free Software Foundation.
14 *
15 * bitpit is distributed in the hope that it will be useful, but WITHOUT
16 * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
17 * FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
18 * License for more details.
19 *
20 * You should have received a copy of the GNU Lesser General Public License
21 * along with bitpit. If not, see <http://www.gnu.org/licenses/>.
22 *
23\*---------------------------------------------------------------------------*/
24
25#ifndef __BTPIT_RECONSTRUCTION_HPP__
26#define __BTPIT_RECONSTRUCTION_HPP__
27
28#include <array>
29#include <iostream>
30#include <memory>
31#include <vector>
32
33#include "bitpit_patchkernel.hpp"
34
35namespace bitpit {
36
37class ReconstructionPolynomial {
38
39friend class ReconstructionKernel;
40friend class ReconstructionAssembler;
41
42public:
43 ReconstructionPolynomial();
44 ReconstructionPolynomial(uint8_t degree, uint8_t dimensions, const std::array<double, 3> &origin, int nFields=1);
45
46 ReconstructionPolynomial(const ReconstructionPolynomial &other);
47 ReconstructionPolynomial(ReconstructionPolynomial &&other) = default;
48 ReconstructionPolynomial & operator = (const ReconstructionPolynomial &other);
49 ReconstructionPolynomial & operator=(ReconstructionPolynomial &&other) = default;
50
51 void swap(ReconstructionPolynomial &other) noexcept;
52
53 void initialize(uint8_t degree, uint8_t dimensions, const std::array<double, 3> &origin, int nFields=1, bool release = true);
54 void clear(bool release = true);
55
56 uint8_t getDegree() const;
57 uint8_t getDimensions() const;
58 const std::array<double, 3> & getOrigin() const;
59
60 uint16_t getCoefficientCount() const;
61
62 int getFieldCount() const;
63
64 const double * getCoefficients() const;
65 double * getCoefficients();
66 const double * getCoefficients(int field) const;
67 double * getCoefficients(int field);
68 const double * getDegreeCoefficients(uint8_t degree, int field = 0) const;
69 double * getDegreeCoefficients(uint8_t degree, int field = 0);
70
71 void computeValue(const std::array<double, 3> &point, int field, double *values) const;
72 void computeValues(const std::array<double, 3> &point, double *values) const;
73 void computeValues(const std::array<double, 3> &point, int nFields, double *values) const;
74 void computeValues(const std::array<double, 3> &point, int nFields, int offset, double *values) const;
75 void computeValue(int degree, const std::array<double, 3> &point, int field, double *values) const;
76 void computeValues(int degree, const std::array<double, 3> &point, double *values) const;
77 void computeValues(int degree, const std::array<double, 3> &point, int nFields, double *values) const;
78 void computeValues(int degree, const std::array<double, 3> &point, int nFields, int offset, double *values) const;
79
80 void computeValueLimited(const std::array<double, 3> &point, const double *limiters, int field, double *value) const;
81 void computeValuesLimited(const std::array<double, 3> &point, const double *limiters, double *values) const;
82 void computeValuesLimited(const std::array<double, 3> &point, const double *limiters, int nFields, double *values) const;
83 void computeValuesLimited(const std::array<double, 3> &point, const double *limiters, int nFields, int offset, double *values) const;
84 void computeValueLimited(int degree, const std::array<double, 3> &point, const double *limiters, int field, double *values) const;
85 void computeValuesLimited(int degree, const std::array<double, 3> &point, const double *limiters, double *values) const;
86 void computeValuesLimited(int degree, const std::array<double, 3> &point, const double *limiters, int nFields, double *values) const;
87 void computeValuesLimited(int degree, const std::array<double, 3> &point, const double *limiters, int nFields, int offset, double *values) const;
88
89 void computeDerivative(const std::array<double, 3> &point, const std::array<double, 3> &direction, int field, double *derivative) const;
90 void computeDerivatives(const std::array<double, 3> &point, const std::array<double, 3> &direction, double *derivatives) const;
91 void computeDerivatives(const std::array<double, 3> &point, const std::array<double, 3> &direction, int nFields, double *derivatives) const;
92 void computeDerivatives(const std::array<double, 3> &point, const std::array<double, 3> &direction, int nFields, int offset, double *derivatives) const;
93 void computeDerivative(int degree, const std::array<double, 3> &point, const std::array<double, 3> &direction, int field, double *derivative) const;
94 void computeDerivatives(int degree, const std::array<double, 3> &point, const std::array<double, 3> &direction, double *derivatives) const;
95 void computeDerivatives(int degree, const std::array<double, 3> &point, const std::array<double, 3> &direction, int nFields, double *derivatives) const;
96 void computeDerivatives(int degree, const std::array<double, 3> &point, const std::array<double, 3> &direction, int nFields, int offset, double *derivatives) const;
97
98 void computeDerivativeLimited(const std::array<double, 3> &point, const std::array<double, 3> &direction, const double *limiters, int field, double *derivative) const;
99 void computeDerivativesLimited(const std::array<double, 3> &point, const std::array<double, 3> &direction, const double *limiters, double *derivatives) const;
100 void computeDerivativesLimited(const std::array<double, 3> &point, const std::array<double, 3> &direction, const double *limiters, int nFields, double *derivatives) const;
101 void computeDerivativesLimited(const std::array<double, 3> &point, const std::array<double, 3> &direction, const double *limiters, int nFields, int offset, double *derivatives) const;
102 void computeDerivativeLimited(int degree, const std::array<double, 3> &point, const std::array<double, 3> &direction, const double *limiters, int field, double *derivative) const;
103 void computeDerivativesLimited(int degree, const std::array<double, 3> &point, const std::array<double, 3> &direction, const double *limiters, double *derivatives) const;
104 void computeDerivativesLimited(int degree, const std::array<double, 3> &point, const std::array<double, 3> &direction, const double *limiters, int nFields, double *derivatives) const;
105 void computeDerivativesLimited(int degree, const std::array<double, 3> &point, const std::array<double, 3> &direction, const double *limiters, int nFields, int offset, double *derivatives) const;
106
107 void computeGradient(const std::array<double, 3> &point, int field, std::array<double, 3> *gradient) const;
108 void computeGradients(const std::array<double, 3> &point, std::array<double, 3> *gradients) const;
109 void computeGradients(const std::array<double, 3> &point, int nFields, std::array<double, 3> *gradients) const;
110 void computeGradients(const std::array<double, 3> &point, int nFields, int offset, std::array<double, 3> *gradients) const;
111 void computeGradient(int degree, const std::array<double, 3> &point, int field, std::array<double, 3> *gradient) const;
112 void computeGradients(int degree, const std::array<double, 3> &point, std::array<double, 3> *gradients) const;
113 void computeGradients(int degree, const std::array<double, 3> &point, int nFields, std::array<double, 3> *gradients) const;
114 void computeGradients(int degree, const std::array<double, 3> &point, int nFields, int offset, std::array<double, 3> *gradients) const;
115
116 void computeGradientLimited(const std::array<double, 3> &point, const double *limiters, int field, std::array<double, 3> *gradient) const;
117 void computeGradientsLimited(const std::array<double, 3> &point, const double *limiters, std::array<double, 3> *gradients) const;
118 void computeGradientsLimited(const std::array<double, 3> &point, const double *limiters, int nFields, std::array<double, 3> *gradients) const;
119 void computeGradientsLimited(const std::array<double, 3> &point, const double *limiters, int nFields, int offset, std::array<double, 3> *gradients) const;
120 void computeGradientLimited(int degree, const std::array<double, 3> &point, const double *limiters, int field, std::array<double, 3> *gradient) const;
121 void computeGradientsLimited(int degree, const std::array<double, 3> &point, const double *limiters, std::array<double, 3> *gradients) const;
122 void computeGradientsLimited(int degree, const std::array<double, 3> &point, const double *limiters, int nFields, std::array<double, 3> *gradients) const;
123 void computeGradientsLimited(int degree, const std::array<double, 3> &point, const double *limiters, int nFields, int offset, std::array<double, 3> *gradients) const;
124
125 void display(std::ostream &out) const;
126
127protected:
128 static const uint8_t MAX_DEGREE;
129 static const uint8_t MAX_DIMENSIONS;
130
131 static uint16_t getCoefficientCount(uint8_t degree, uint8_t dimensions);
132 static const std::vector<uint16_t> & getCoefficientsCount(uint8_t dimensions);
133 static uint16_t countCoefficients(uint8_t degree, uint8_t dimensions);
134 static const std::vector<uint16_t> & getDegreeCoefficientsCount(uint8_t dimensions);
135
136 static uint16_t getDegreeCoefficientCount(uint8_t degree, uint8_t dimensions);
137 static uint16_t countDegreeCoefficients(uint8_t degree, uint8_t dimensions);
138
139 static void evalPointBasisValues(uint8_t degree, uint8_t dimensions, const std::array<double, 3> &origin, const std::array<double, 3> &point, double *csi);
140 static void evalPointBasisDerivatives(uint8_t degree, uint8_t dimensions, const std::array<double, 3> &origin, const std::array<double, 3> &point, const std::array<double, 3> &direction, double *dcsi);
141
142 static void evalCellBasisValues(uint8_t degree, uint8_t dimensions, const std::array<double, 3> &origin, const Cell &cell, const std::array<double, 3> *vertexCoords, double *csi);
143
144private:
145 static const bool ENABLE_FAST_PATH_OPTIMIZATIONS;
146 static const int MAX_STACK_WORKSPACE_SIZE;
147
148 const static std::vector<std::vector<uint16_t>> m_countCoefficientCache;
149 const static std::vector<std::vector<uint16_t>> m_countDegreeCoefficientCache;
150
151 static std::vector<std::vector<uint16_t>> generateCountCoefficientCache();
152 static std::vector<std::vector<uint16_t>> generateCountDegreeCoefficientCache();
153
154 uint8_t m_degree;
155 uint8_t m_dimensions;
156
157 int m_nFields;
158
159 uint16_t m_nCoeffs;
160 std::unique_ptr<double[]> m_coeffs;
161
162 std::array<double, 3> m_origin;
163
164 std::size_t computeFieldCoefficientsOffset(uint8_t degree, int field) const;
165 std::size_t getFieldCoefficientsStride() const;
166
167};
168
169class ReconstructionKernel {
170
171public:
172 ReconstructionKernel();
173 ReconstructionKernel(uint8_t degree, uint8_t dimensions, int nEquations);
174
175 ReconstructionKernel(const ReconstructionKernel &other);
176 ReconstructionKernel(ReconstructionKernel &&other) = default;
177 ReconstructionKernel & operator = (const ReconstructionKernel &other);
178 ReconstructionKernel & operator=(ReconstructionKernel &&other) = default;
179
180 void swap(ReconstructionKernel &other) noexcept;
181
182 void initialize(uint8_t degree, uint8_t dimensions, int nEquations, bool release = true);
183
184 uint8_t getDegree() const;
185 uint8_t getDimensions() const;
186
187 uint16_t getCoefficientCount() const;
188
189 int getEquationCount() const;
190
191 const double * getPolynomialWeights() const;
192 double * getPolynomialWeights();
193
194 void assemblePolynomial(const std::array<double, 3> &origin, const double *values, ReconstructionPolynomial *polynomial) const;
195 void assemblePolynomial(const std::array<double, 3> &origin, int nFields, const double **values, ReconstructionPolynomial *polynomial) const;
196 void assemblePolynomial(uint8_t degree, const std::array<double, 3> &origin, const double *values, ReconstructionPolynomial *polynomial) const;
197 void assemblePolynomial(uint8_t degree, const std::array<double, 3> &origin, int nFields, const double **values, ReconstructionPolynomial *polynomial) const;
198
199 void updatePolynomial(const double *values, ReconstructionPolynomial *polynomial) const;
200 void updatePolynomial(int nFields, const double **values, ReconstructionPolynomial *polynomial) const;
201 void updatePolynomial(uint8_t degree, const double *values, ReconstructionPolynomial *polynomial) const;
202 void updatePolynomial(uint8_t degree, int nFields, const double **values, ReconstructionPolynomial *polynomial) const;
203
204 void computeValueWeights(const std::array<double, 3> &origin, const std::array<double, 3> &point, double *valueWeights) const;
205 void computeValueWeights(uint8_t degree, const std::array<double, 3> &origin, const std::array<double, 3> &point, double *valueWeights) const;
206 void computeValueWeights(int nEquations, uint8_t degree, const std::array<double, 3> &origin, const std::array<double, 3> &point, double *valueWeights) const;
207 void computeValueLimitedWeights(const std::array<double, 3> &origin, const std::array<double, 3> &point, const double *limiters, double *valueWeights) const;
208 void computeValueLimitedWeights(uint8_t degree, const std::array<double, 3> &origin, const std::array<double, 3> &point, const double *limiters, double *valueWeights) const;
209 void computeValueLimitedWeights(int nEquations, uint8_t degree, const std::array<double, 3> &origin, const std::array<double, 3> &point, const double *limiters, double *valueWeights) const;
210
211 void computeDerivativeWeights(const std::array<double, 3> &origin, const std::array<double, 3> &point, const std::array<double, 3> &direction, double *derivativeWeights) const;
212 void computeDerivativeWeights(uint8_t degree, const std::array<double, 3> &origin, const std::array<double, 3> &point, const std::array<double, 3> &direction, double *derivativeWeights) const;
213 void computeDerivativeWeights(int nEquations, uint8_t degree, const std::array<double, 3> &origin, const std::array<double, 3> &point, const std::array<double, 3> &direction, double *derivativeWeights) const;
214 void computeDerivativeLimitedWeights(const std::array<double, 3> &origin, const std::array<double, 3> &point, const std::array<double, 3> &direction, const double *limiters, double *derivativeWeights) const;
215 void computeDerivativeLimitedWeights(uint8_t degree, const std::array<double, 3> &origin, const std::array<double, 3> &point, const std::array<double, 3> &direction, const double *limiters, double *derivativeWeights) const;
216 void computeDerivativeLimitedWeights(int nEquations, uint8_t degree, const std::array<double, 3> &origin, const std::array<double, 3> &point, const std::array<double, 3> &direction, const double *limiters, double *derivativeWeights) const;
217
218 void computeGradientWeights(const std::array<double, 3> &origin, const std::array<double, 3> &point, std::array<double, 3> *gradientWeights) const;
219 void computeGradientWeights(uint8_t degree, const std::array<double, 3> &origin, const std::array<double, 3> &point, std::array<double, 3> *gradientWeights) const;
220 void computeGradientWeights(int nEquations, uint8_t degree, const std::array<double, 3> &origin, const std::array<double, 3> &point, std::array<double, 3> *gradientWeights) const;
221 void computeGradientLimitedWeights(const std::array<double, 3> &origin, const std::array<double, 3> &point, const double *limiters, std::array<double, 3> *gradientWeights) const;
222 void computeGradientLimitedWeights(uint8_t degree, const std::array<double, 3> &origin, const std::array<double, 3> &point, const double *limiters, std::array<double, 3> *gradientWeights) const;
223 void computeGradientLimitedWeights(int nEquations, uint8_t degree, const std::array<double, 3> &origin, const std::array<double, 3> &point, const double *limiters, std::array<double, 3> *gradientWeights) const;
224
225 void display(std::ostream &out, double tolerance = 1.e-10) const;
226
227protected:
228 void applyLimiter(uint8_t degree, const double *limiters, double *coeffs) const;
229
230private:
231 static const bool ENABLE_FAST_PATH_OPTIMIZATIONS;
232 static const int MAX_STACK_WORKSPACE_SIZE;
233
234 std::unique_ptr<double[]> m_weights;
235
236 int m_nEquations;
237
238 uint16_t m_nCoeffs;
239
240 uint8_t m_degree;
241 uint8_t m_dimensions;
242
243};
244
245class ReconstructionAssembler {
246
247public:
248 enum ReconstructionType {
249 TYPE_CONSTRAINT,
250 TYPE_LEAST_SQUARE
251 };
252
253 ReconstructionAssembler();
254 ReconstructionAssembler(uint8_t degree, uint8_t dimensions);
255
256 void swap(ReconstructionAssembler &other) noexcept;
257
258 void initialize(uint8_t degree, uint8_t dimensions, bool release = true);
259 void clear(bool release = true);
260
261 uint8_t getDegree() const;
262 uint8_t getDimensions() const;
263
264 uint16_t getCoefficientCount() const;
265
266 int countConstraints() const;
267 int countLeastSquares() const;
268 int countEquations() const;
269
270 void addPointValueEquation(ReconstructionType type, const std::array<double, 3> &origin, const std::array<double, 3> &point, double scaleFactor = 1.);
271 void addPointDerivativeEquation(ReconstructionType type, const std::array<double, 3> &origin, const std::array<double, 3> &point, const std::array<double, 3> &direction, double scaleFactor = 1.);
272 void addCellAverageEquation(ReconstructionType type, const Cell &cell, const std::array<double, 3> &origin, const std::array<double, 3> *vertexCoords, double scaleFactor = 1.);
273
274 void assembleKernel(ReconstructionKernel *kernel) const;
275
276 void updateKernel(ReconstructionKernel *kernel) const;
277
278private:
279 static const double SVD_ZERO_THRESHOLD;
280
281 uint8_t m_degree;
282 uint8_t m_dimensions;
283
284 uint16_t m_nCoeffs;
285
286 std::vector<int> m_constraintsOrder;
287 std::vector<int> m_leastSquaresOrder;
288
289 std::vector<double> m_leastSquaresScaleFactors;
290
291 std::vector<double> m_A;
292 std::vector<double> m_C;
293
294 mutable std::vector<double> m_sigma;
295 mutable std::vector<double> m_U;
296 mutable std::vector<double> m_S;
297 mutable std::vector<double> m_Vt;
298 mutable std::vector<double> m_SVDWorkspace;
299 mutable std::vector<double> m_w;
300
301 double * _addEquation(ReconstructionType type, double scaleFactor);
302
303 void computePseudoInverse(int m, int n, double tolerance, double *A) const;
304
305};
306
307class Reconstruction : public ReconstructionKernel, public ReconstructionAssembler {
308
309public:
310 Reconstruction(uint8_t degree, uint8_t dimensions);
311
312 void swap(Reconstruction &other) noexcept;
313
314 void initialize(uint8_t degree, uint8_t dimensions, bool release = true);
315 void clear(bool release = true);
316
317 using ReconstructionAssembler::getDegree;
318 using ReconstructionAssembler::getDimensions;
319
320 using ReconstructionAssembler::getCoefficientCount;
321
322 void assemble();
323
324};
325
326}
327
328#endif
bitpit::Cell
The Cell class defines the cells.
Definition cell.hpp:42
bitpit::ReconstructionAssembler::swap
void swap(ReconstructionAssembler &other) noexcept
Definition reconstruction.cpp:3153
bitpit::ReconstructionAssembler::addPointValueEquation
void addPointValueEquation(ReconstructionType type, const std::array< double, 3 > &origin, const std::array< double, 3 > &point, double scaleFactor=1.)
Definition reconstruction.cpp:3304
bitpit::ReconstructionAssembler::assembleKernel
void assembleKernel(ReconstructionKernel *kernel) const
Definition reconstruction.cpp:3403
bitpit::ReconstructionAssembler::initialize
void initialize(uint8_t degree, uint8_t dimensions, bool release=true)
Definition reconstruction.cpp:3180
bitpit::ReconstructionAssembler::addPointDerivativeEquation
void addPointDerivativeEquation(ReconstructionType type, const std::array< double, 3 > &origin, const std::array< double, 3 > &point, const std::array< double, 3 > &direction, double scaleFactor=1.)
Definition reconstruction.cpp:3322
bitpit::ReconstructionAssembler::addCellAverageEquation
void addCellAverageEquation(ReconstructionType type, const Cell &cell, const std::array< double, 3 > &origin, const std::array< double, 3 > *vertexCoords, double scaleFactor=1.)
Definition reconstruction.cpp:3341
bitpit::ReconstructionAssembler::updateKernel
void updateKernel(ReconstructionKernel *kernel) const
Definition reconstruction.cpp:3424
bitpit::ReconstructionKernel
The ReconstructionKernel class allows to evaluate the weight of a reconstruction polynomial previousl...
Definition reconstruction.hpp:169
bitpit::ReconstructionKernel::computeGradientWeights
void computeGradientWeights(const std::array< double, 3 > &origin, const std::array< double, 3 > &point, std::array< double, 3 > *gradientWeights) const
Definition reconstruction.cpp:2864
bitpit::ReconstructionKernel::computeValueWeights
void computeValueWeights(const std::array< double, 3 > &origin, const std::array< double, 3 > &point, double *valueWeights) const
Definition reconstruction.cpp:2469
bitpit::ReconstructionKernel::computeGradientLimitedWeights
void computeGradientLimitedWeights(const std::array< double, 3 > &origin, const std::array< double, 3 > &point, const double *limiters, std::array< double, 3 > *gradientWeights) const
Definition reconstruction.cpp:2942
bitpit::ReconstructionKernel::updatePolynomial
void updatePolynomial(const double *values, ReconstructionPolynomial *polynomial) const
Definition reconstruction.cpp:2312
bitpit::ReconstructionKernel::computeDerivativeLimitedWeights
void computeDerivativeLimitedWeights(const std::array< double, 3 > &origin, const std::array< double, 3 > &point, const std::array< double, 3 > &direction, const double *limiters, double *derivativeWeights) const
Definition reconstruction.cpp:2745
bitpit::ReconstructionKernel::swap
void swap(ReconstructionKernel &other) noexcept
Definition reconstruction.cpp:2071
bitpit::ReconstructionKernel::initialize
void initialize(uint8_t degree, uint8_t dimensions, int nEquations, bool release=true)
Definition reconstruction.cpp:2090
bitpit::ReconstructionKernel::display
void display(std::ostream &out, double tolerance=1.e-10) const
Definition reconstruction.cpp:3095
bitpit::ReconstructionKernel::computeValueLimitedWeights
void computeValueLimitedWeights(const std::array< double, 3 > &origin, const std::array< double, 3 > &point, const double *limiters, double *valueWeights) const
Definition reconstruction.cpp:2546
bitpit::ReconstructionKernel::computeDerivativeWeights
void computeDerivativeWeights(const std::array< double, 3 > &origin, const std::array< double, 3 > &point, const std::array< double, 3 > &direction, double *derivativeWeights) const
Definition reconstruction.cpp:2660
bitpit::ReconstructionKernel::operator=
ReconstructionKernel & operator=(const ReconstructionKernel &other)
Definition reconstruction.cpp:2056
bitpit::ReconstructionKernel::assemblePolynomial
void assemblePolynomial(const std::array< double, 3 > &origin, const double *values, ReconstructionPolynomial *polynomial) const
Definition reconstruction.cpp:2193
bitpit::ReconstructionKernel::getPolynomialWeights
const double * getPolynomialWeights() const
Definition reconstruction.cpp:2164
bitpit::ReconstructionKernel::applyLimiter
void applyLimiter(uint8_t degree, const double *limiters, double *coeffs) const
Definition reconstruction.cpp:3067
bitpit::ReconstructionPolynomial
The ReconstructionPolynomial class allows to apply a reconstruction polynomial previously assembled.
Definition reconstruction.hpp:37
bitpit::ReconstructionPolynomial::countCoefficients
static uint16_t countCoefficients(uint8_t degree, uint8_t dimensions)
Definition reconstruction.cpp:148
bitpit::ReconstructionPolynomial::evalPointBasisDerivatives
static void evalPointBasisDerivatives(uint8_t degree, uint8_t dimensions, const std::array< double, 3 > &origin, const std::array< double, 3 > &point, const std::array< double, 3 > &direction, double *dcsi)
Definition reconstruction.cpp:267
bitpit::ReconstructionPolynomial::computeDerivativeLimited
void computeDerivativeLimited(const std::array< double, 3 > &point, const std::array< double, 3 > &direction, const double *limiters, int field, double *derivative) const
Definition reconstruction.cpp:1268
bitpit::ReconstructionPolynomial::computeValue
void computeValue(const std::array< double, 3 > &point, int field, double *values) const
Definition reconstruction.cpp:697
bitpit::ReconstructionPolynomial::getDegreeCoefficients
const double * getDegreeCoefficients(uint8_t degree, int field=0) const
Definition reconstruction.cpp:641
bitpit::ReconstructionPolynomial::computeDerivativesLimited
void computeDerivativesLimited(const std::array< double, 3 > &point, const std::array< double, 3 > &direction, const double *limiters, double *derivatives) const
Definition reconstruction.cpp:1289
bitpit::ReconstructionPolynomial::display
void display(std::ostream &out) const
Definition reconstruction.cpp:1972
bitpit::ReconstructionPolynomial::computeValuesLimited
void computeValuesLimited(const std::array< double, 3 > &point, const double *limiters, double *values) const
Definition reconstruction.cpp:875
bitpit::ReconstructionPolynomial::initialize
void initialize(uint8_t degree, uint8_t dimensions, const std::array< double, 3 > &origin, int nFields=1, bool release=true)
Definition reconstruction.cpp:485
bitpit::ReconstructionPolynomial::swap
void swap(ReconstructionPolynomial &other) noexcept
Definition reconstruction.cpp:464
bitpit::ReconstructionPolynomial::getCoefficients
const double * getCoefficients() const
Definition reconstruction.cpp:591
bitpit::ReconstructionPolynomial::evalPointBasisValues
static void evalPointBasisValues(uint8_t degree, uint8_t dimensions, const std::array< double, 3 > &origin, const std::array< double, 3 > &point, double *csi)
Definition reconstruction.cpp:218
bitpit::ReconstructionPolynomial::computeGradientsLimited
void computeGradientsLimited(const std::array< double, 3 > &point, const double *limiters, std::array< double, 3 > *gradients) const
Definition reconstruction.cpp:1735
bitpit::ReconstructionPolynomial::getCoefficientsCount
static const std::vector< uint16_t > & getCoefficientsCount(uint8_t dimensions)
Definition reconstruction.cpp:132
bitpit::ReconstructionPolynomial::computeGradients
void computeGradients(const std::array< double, 3 > &point, std::array< double, 3 > *gradients) const
Definition reconstruction.cpp:1532
bitpit::ReconstructionPolynomial::computeDerivative
void computeDerivative(const std::array< double, 3 > &point, const std::array< double, 3 > &direction, int field, double *derivative) const
Definition reconstruction.cpp:1079
bitpit::ReconstructionPolynomial::computeGradientLimited
void computeGradientLimited(const std::array< double, 3 > &point, const double *limiters, int field, std::array< double, 3 > *gradient) const
Definition reconstruction.cpp:1716
bitpit::ReconstructionPolynomial::evalCellBasisValues
static void evalCellBasisValues(uint8_t degree, uint8_t dimensions, const std::array< double, 3 > &origin, const Cell &cell, const std::array< double, 3 > *vertexCoords, double *csi)
Definition reconstruction.cpp:338
bitpit::ReconstructionPolynomial::getDegreeCoefficientCount
static uint16_t getDegreeCoefficientCount(uint8_t degree, uint8_t dimensions)
Definition reconstruction.cpp:167
bitpit::ReconstructionPolynomial::computeValueLimited
void computeValueLimited(const std::array< double, 3 > &point, const double *limiters, int field, double *value) const
Definition reconstruction.cpp:856
bitpit::ReconstructionPolynomial::computeDerivatives
void computeDerivatives(const std::array< double, 3 > &point, const std::array< double, 3 > &direction, double *derivatives) const
Definition reconstruction.cpp:1093
bitpit::ReconstructionPolynomial::computeValues
void computeValues(const std::array< double, 3 > &point, double *values) const
Definition reconstruction.cpp:709
bitpit::ReconstructionPolynomial::operator=
ReconstructionPolynomial & operator=(const ReconstructionPolynomial &other)
Definition reconstruction.cpp:449
bitpit::ReconstructionPolynomial::getOrigin
const std::array< double, 3 > & getOrigin() const
Definition reconstruction.cpp:561
bitpit::ReconstructionPolynomial::countDegreeCoefficients
static uint16_t countDegreeCoefficients(uint8_t degree, uint8_t dimensions)
Definition reconstruction.cpp:199
bitpit::ReconstructionPolynomial::getDegreeCoefficientsCount
static const std::vector< uint16_t > & getDegreeCoefficientsCount(uint8_t dimensions)
Definition reconstruction.cpp:183
bitpit::ReconstructionPolynomial::computeGradient
void computeGradient(const std::array< double, 3 > &point, int field, std::array< double, 3 > *gradient) const
Definition reconstruction.cpp:1520
bitpit::Reconstruction::clear
void clear(bool release=true)
Definition reconstruction.cpp:3659
bitpit::Reconstruction::initialize
void initialize(uint8_t degree, uint8_t dimensions, bool release=true)
Definition reconstruction.cpp:3647
bitpit::Reconstruction::Reconstruction
Reconstruction(uint8_t degree, uint8_t dimensions)
Definition reconstruction.cpp:3620
bitpit::Reconstruction::swap
void swap(Reconstruction &other) noexcept
Definition reconstruction.cpp:3632
Generated on Mon Apr 28 2025 20:53:28 for bitpit by   doxygen 1.13.2