src/foam/primitives/FadOne/FadOneField.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | foam-extend: Open Source CFD
   4    \\    /   O peration     | Version:     3.2
   5     \\  /    A nd           | Web:         http://www.foam-extend.org
   6      \\/     M anipulation  | For copyright notice see file Copyright
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of foam-extend.
  10
  11     foam-extend is free software: you can redistribute it and/or modify it
  12     under the terms of the GNU General Public License as published by the
  13     Free Software Foundation, either version 3 of the License, or (at your
  14     option) any later version.
  15
  16     foam-extend is distributed in the hope that it will be useful, but
  17     WITHOUT ANY WARRANTY; without even the implied warranty of
  18     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  19     General Public License for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 Class
  25     FadOneFields
  26
  27 \*----------------------------------------------------------------------------*/
  28
  29 #include "FadOneField.H"
  30
  31 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  32
  33 template<int nVars>
  34 Foam::tmp<Foam::Field<Foam::FadOne<nVars> > >
  35 Foam::ValueFadOneField(const UList<scalar>& u)
  36 {
  37     tmp<Foam::Field<FadOne<nVars> > > tr
  38     (
  39         new Foam::Field<FadOne<nVars> > (u.size())
  40     );
  41
  42     Field<FadOne<nVars> >& r = tr();
  43
  44     forAll (r, i)
  45     {
  46         r[i].value() = u[i];
  47     }
  48
  49     return tr;
  50 }
  51
  52
  53 template<int nVars>
  54 Foam::tmp<Foam::scalarField> Foam::FadOneValue(const Field<FadOne<nVars> >& u)
  55 {
  56     tmp<scalarField> tr(new scalarField(u.size()));
  57     scalarField& r = tr();
  58
  59     forAll (r, i)
  60     {
  61         r[i] = u[i].value();
  62     }
  63
  64     return tr;
  65 }
  66
  67
  68 template<int nVars>
  69 void Foam::FadOneSetValue
  70 (
  71     Field<FadOne<nVars> >& u,
  72     const scalarField& val
  73 )
  74 {
  75     forAll (u, i)
  76     {
  77         u[i].value() = val[i];
  78     }
  79 }
  80
  81
  82 template<int nVars>
  83 Foam::tmp<Foam::scalarField> Foam::FadOneDeriv
  84 (
  85     const Field<FadOne<nVars> >& u,
  86     const direction d
  87 )
  88 {
  89     tmp<scalarField> tr(new scalarField(u.size()));
  90     scalarField& r = tr();
  91
  92     forAll (r, i)
  93     {
  94         r[i] = u[i].deriv(d);
  95     }
  96
  97     return tr;
  98 }
  99
 100
 101 template<int nVars>
 102 void Foam::FadOneSetDeriv
 103 (
 104     Field<FadOne<nVars> >& u,
 105     const direction d,
 106     const scalarField& der
 107 )
 108 {
 109     forAll (u, i)
 110     {
 111         u[i].deriv(d) = der[i];
 112     }
 113 }
 114
 115
 116 // ************************************************************************* //