FreeFOAM The Cross-Platform CFD Toolkit
Hosted by SourceForge:
Get FreeFOAM at SourceForge.net.
            Fast, secure and Free Open Source software downloads

cubicUpwindFitPolynomial.H

Go to the documentation of this file.
00001 /*---------------------------------------------------------------------------*\
00002   =========                 |
00003   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
00004    \\    /   O peration     |
00005     \\  /    A nd           | Copyright (C) 1991-2010 OpenCFD Ltd.
00006      \\/     M anipulation  |
00007 -------------------------------------------------------------------------------
00008 License
00009     This file is part of OpenFOAM.
00010 
00011     OpenFOAM is free software: you can redistribute it and/or modify it
00012     under the terms of the GNU General Public License as published by
00013     the Free Software Foundation, either version 3 of the License, or
00014     (at your option) any later version.
00015 
00016     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
00017     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
00018     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
00019     for more details.
00020 
00021     You should have received a copy of the GNU General Public License
00022     along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
00023 
00024 Class
00025     Foam::cubicUpwindFitPolynomial
00026 
00027 Description
00028     Cubic polynomial for upwind biased interpolation fitting.
00029 
00030     Can be used with the UpwindFit scheme to create a cubic surface
00031     interpolation scheme
00032 
00033 \*---------------------------------------------------------------------------*/
00034 
00035 #ifndef cubicUpwindFitPolynomial_H
00036 #define cubicUpwindFitPolynomial_H
00037 
00038 #include <OpenFOAM/vector.H>
00039 
00040 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00041 
00042 namespace Foam
00043 {
00044 
00045 /*---------------------------------------------------------------------------*\
00046                  Class cubicUpwindFitPolynomial Declaration
00047 \*---------------------------------------------------------------------------*/
00048 
00049 class cubicUpwindFitPolynomial
00050 {
00051 public:
00052 
00053     // Member functions
00054 
00055         static label nTerms(const direction dim)
00056         {
00057             return
00058             (
00059                 dim == 1 ? 4 :
00060                 dim == 2 ? 8 :
00061                 dim == 3 ? 14 : 0
00062             );
00063         }
00064 
00065         static void addCoeffs
00066         (
00067             scalar* coeffs,
00068             const vector& d,
00069             const scalar weight,
00070             const direction dim
00071         )
00072         {
00073             register label curIdx = 0;
00074 
00075             coeffs[curIdx++] = weight;
00076             coeffs[curIdx++] = weight*d.x();
00077             coeffs[curIdx++] = weight*sqr(d.x());
00078             coeffs[curIdx++] = weight*pow(d.x(),3);
00079 
00080             if (dim >= 2)
00081             {
00082                 coeffs[curIdx++] = weight*d.y();
00083                 coeffs[curIdx++] = weight*d.x()*d.y();
00084                 coeffs[curIdx++] = weight*sqr(d.y());
00085                 coeffs[curIdx++] = weight*d.x()*sqr(d.y());
00086             }
00087             if (dim == 3)
00088             {
00089                 coeffs[curIdx++] = weight*d.z();
00090                 coeffs[curIdx++] = weight*d.x()*d.z();
00091                 coeffs[curIdx++] = weight*d.y()*d.z();
00092                 coeffs[curIdx++] = weight*sqr(d.z());
00093                 coeffs[curIdx++] = weight*d.x()*d.y()*d.z();
00094                 coeffs[curIdx++] = weight*d.x()*sqr(d.z());
00095             }
00096         }
00097 };
00098 
00099 
00100 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00101 
00102 } // End namespace Foam
00103 
00104 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00105 
00106 #endif
00107 
00108 // ************************ vim: set sw=4 sts=4 et: ************************ //
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines