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

polynomial.C

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) 2008-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 \*---------------------------------------------------------------------------*/
00025 
00026 #include "polynomial.H"
00027 
00028 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
00029 
00030 namespace Foam
00031 {
00032     defineTypeNameAndDebug(polynomial, 0);
00033 }
00034 
00035 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
00036 
00037 Foam::polynomial::polynomial(const word& entryName, Istream& is)
00038 :
00039     DataEntry<scalar>(entryName),
00040     coeffs_(is)
00041 {
00042     if (!coeffs_.size())
00043     {
00044         FatalErrorIn("Foam::polynomial::polynomial(const word&, Istream&)")
00045             << "polynomial coefficients for entry " << this->name_
00046             << " is invalid (empty)" << nl << exit(FatalError);
00047     }
00048 }
00049 
00050 
00051 Foam::polynomial::polynomial(const polynomial& poly)
00052 :
00053     DataEntry<scalar>(poly),
00054     coeffs_(poly.coeffs_)
00055 {}
00056 
00057 
00058 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
00059 
00060 Foam::polynomial::~polynomial()
00061 {}
00062 
00063 
00064 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
00065 
00066 Foam::scalar Foam::polynomial::value(const scalar x) const
00067 {
00068     scalar y = 0.0;
00069     forAll(coeffs_, i)
00070     {
00071         y += coeffs_[i].first()*pow(x, coeffs_[i].second());
00072     }
00073 
00074     return y;
00075 }
00076 
00077 
00078 Foam::scalar Foam::polynomial::integrate(const scalar x1, const scalar x2) const
00079 {
00080     scalar intx = 0.0;
00081 
00082     forAll(coeffs_, i)
00083     {
00084         intx +=
00085             coeffs_[i].first()/(coeffs_[i].second() + 1)
00086            *(
00087                 pow(x2, coeffs_[i].second() + 1)
00088               - pow(x1, coeffs_[i].second() + 1)
00089             );
00090     }
00091 
00092     return intx;
00093 }
00094 
00095 
00096 // ************************ vim: set sw=4 sts=4 et: ************************ //
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines