37#ifndef cubicUpwindFitPolynomial_H
38#define cubicUpwindFitPolynomial_H
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);
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());
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());
const Cmpt & z() const
Access to the vector z component.
const Cmpt & y() const
Access to the vector y component.
const Cmpt & x() const
Access to the vector x component.
Cubic 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)
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)