quadraticUpwindFitPolynomial.H
Go to the documentation of this file.
1/*---------------------------------------------------------------------------*\
2 ========= |
3 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4 \\ / O peration |
5 \\ / A nd | www.openfoam.com
6 \\/ M anipulation |
7-------------------------------------------------------------------------------
8 Copyright (C) 2011-2015 OpenFOAM Foundation
9-------------------------------------------------------------------------------
10License
11 This file is part of OpenFOAM.
12
13 OpenFOAM is free software: you can redistribute it and/or modify it
14 under the terms of the GNU General Public License as published by
15 the Free Software Foundation, either version 3 of the License, or
16 (at your option) any later version.
17
18 OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
19 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
20 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
21 for more details.
22
23 You should have received a copy of the GNU General Public License
24 along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
25
26Class
27 Foam::quadraticUpwindFitPolynomial
28
29Description
30 Quadratic polynomial for upwind biased interpolation fitting.
31
32 Can be used with the UpwindFit scheme to create a quadratic surface
33 interpolation scheme
34
35\*---------------------------------------------------------------------------*/
36
37#ifndef quadraticUpwindFitPolynomial_H
38#define quadraticUpwindFitPolynomial_H
39
40#include "vector.H"
41
42// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
43
44namespace Foam
45{
46
47/*---------------------------------------------------------------------------*\
48 Class quadraticUpwindFitPolynomial Declaration
49\*---------------------------------------------------------------------------*/
52{
53public:
54
55 // Member functions
57 static label nTerms(const direction dim)
58 {
59 return
60 (
61 dim == 1 ? 3 :
62 dim == 2 ? 6 :
63 dim == 3 ? 9 : 0
64 );
65 }
67 static void addCoeffs
68 (
69 scalar* coeffs,
70 const vector& d,
71 const scalar weight,
72 const direction dim
73 )
74 {
75 label curIdx = 0;
76
77 coeffs[curIdx++] = weight;
78 coeffs[curIdx++] = weight*d.x();
79 coeffs[curIdx++] = weight*sqr(d.x());
80
81 if (dim >= 2)
82 {
83 coeffs[curIdx++] = weight*d.y();
84 coeffs[curIdx++] = weight*d.x()*d.y();
85 //coeffs[curIdx++] = weight*d.x()*sqr(d.y());
86 coeffs[curIdx++] = weight*sqr(d.y());
87 }
88 if (dim == 3)
89 {
90 coeffs[curIdx++] = weight*d.z();
91 coeffs[curIdx++] = weight*d.x()*d.z();
92 coeffs[curIdx++] = weight*sqr(d.z());
93 }
94 }
95};
96
97
98// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
99
100} // End namespace Foam
101
102// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
103
104#endif
105
106// ************************************************************************* //
const Cmpt & z() const
Access to the vector z component.
Definition: VectorI.H:85
const Cmpt & y() const
Access to the vector y component.
Definition: VectorI.H:79
const Cmpt & x() const
Access to the vector x component.
Definition: VectorI.H:73
Quadratic polynomial for upwind biased interpolation fitting.
static void addCoeffs(scalar *coeffs, const vector &d, const scalar weight, const direction dim)
static label nTerms(const direction dim)
Namespace for OpenFOAM.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
uint8_t direction
Definition: direction.H:56