src/finiteVolume/interpolation/surfaceInterpolation/schemes/quadraticLinearFit/quadraticLinearFitPolynomial.H

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