The open source CFD toolbox
cubicUpwindFitPolynomial.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-------------------------------------------------------------------------------
11 This file is part of OpenFOAM.
12
13 OpenFOAM is free software: you can redistribute it and/or modify it
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::cubicUpwindFitPolynomial
28
29Description
30 Cubic polynomial for upwind biased interpolation fitting.
31
32 Can be used with the UpwindFit scheme to create a cubic surface
33 interpolation scheme
34
35\*---------------------------------------------------------------------------*/
36
37#ifndef cubicUpwindFitPolynomial_H
38#define cubicUpwindFitPolynomial_H
39
40#include "vector.H"
41
42// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
43
44namespace Foam
45{
46
47/*---------------------------------------------------------------------------*\
48 Class cubicUpwindFitPolynomial Declaration
49\*---------------------------------------------------------------------------*/
52{
53public:
54
55 // Member functions
57 static label nTerms(const direction dim)
58 {
59 return
60 (
61 dim == 1 ? 4 :
62 dim == 2 ? 8 :
63 dim == 3 ? 14 : 0
64 );
65 }
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 coeffs[curIdx++] = weight*pow(d.x(),3);
81
82 if (dim >= 2)
83 {
84 coeffs[curIdx++] = weight*d.y();
85 coeffs[curIdx++] = weight*d.x()*d.y();
86 coeffs[curIdx++] = weight*sqr(d.y());
87 coeffs[curIdx++] = weight*d.x()*sqr(d.y());
88 }
89 if (dim == 3)
90 {
91 coeffs[curIdx++] = weight*d.z();
92 coeffs[curIdx++] = weight*d.x()*d.z();
93 coeffs[curIdx++] = weight*d.y()*d.z();
94 coeffs[curIdx++] = weight*sqr(d.z());
95 coeffs[curIdx++] = weight*d.x()*d.y()*d.z();
96 coeffs[curIdx++] = weight*d.x()*sqr(d.z());
97 }
98 }
99};
100
101
102// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
103
104} // End namespace Foam
105
106// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
107
108#endif
109
110// ************************************************************************* //
const Cmpt & z() const
Definition: VectorI.H:85
const Cmpt & y() const