src/thermophysicalModels/specie/equationOfState/icoPolynomial/icoPolynomial.C

   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 \*---------------------------------------------------------------------------*/
  25
  26 #include "icoPolynomial.H"
  27 #include "IOstreams.H"
  28
  29 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  30
  31 namespace Foam
  32 {
  33
  34 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  35
  36 template<int PolySize>
  37 icoPolynomial<PolySize>::icoPolynomial(Istream& is)
  38 :
  39     specie(is),
  40     rhoCoeffs_("rhoCoeffs<" + Foam::name(PolySize) + '>', is)
  41 {
  42     rhoCoeffs_ *= this->W();
  43 }
  44
  45
  46 template<int PolySize>
  47 icoPolynomial<PolySize>::icoPolynomial(const dictionary& dict)
  48 :
  49     specie(dict),
  50     rhoCoeffs_
  51 (
  52     dict.subDict("equationOfState").lookup
  53     (
  54         "rhoCoeffs<" + Foam::name(PolySize) + '>'
  55     )
  56 )
  57 {
  58     rhoCoeffs_ *= this->W();
  59 }
  60
  61
  62 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  63
  64 template<int PolySize>
  65 void icoPolynomial<PolySize>::write(Ostream& os) const
  66 {
  67     specie::write(os);
  68
  69     dictionary dict("equationOfState");
  70     dict.add
  71     (
  72         word("rhoCoeffs<" + Foam::name(PolySize) + '>'),
  73         rhoCoeffs_/this->W()
  74     );
  75
  76     os  << indent << dict.dictName() << dict;
  77 }
  78
  79
  80 // * * * * * * * * * * * * * * * Ostream Operator  * * * * * * * * * * * * * //
  81
  82 template<int PolySize>
  83 Ostream& operator<<(Ostream& os, const icoPolynomial<PolySize>& ip)
  84 {
  85     os  << static_cast<const specie&>(ip) << tab
  86         << "rhoCoeffs<" << Foam::name(PolySize) << '>' << tab
  87         << ip.rhoCoeffs_/ip.W();
  88
  89     os.check
  90     (
  91         "Ostream& operator<<(Ostream& os, const icoPolynomial<PolySize>& ip)"
  92     );
  93
  94     return os;
  95 }
  96
  97
  98 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  99
 100 } // End namespace Foam
 101
 102 // ************************************************************************* //