src/errorEstimation/errorEstimate/resErrorSup.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 Description
  25     Residual error estimate for the fv source operators.
  26
  27 \*---------------------------------------------------------------------------*/
  28
  29 #include "resErrorSup.H"
  30
  31 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  32
  33 namespace Foam
  34 {
  35
  36 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  37
  38 namespace resError
  39 {
  40
  41 template<class Type>
  42 tmp<errorEstimate<Type> >
  43 Sp
  44 (
  45     const volScalarField& sp,
  46     GeometricField<Type, fvPatchField, volMesh>& vf
  47 )
  48 {
  49     return tmp<errorEstimate<Type> >
  50     (
  51         new errorEstimate<Type>
  52         (
  53             vf,
  54             sp.dimensions()*vf.dimensions()*dimVolume,
  55             sp.internalField()*vf.internalField(),
  56             scalarField(vf.internalField().size(), 0)
  57         )
  58     );
  59 }
  60
  61 template<class Type>
  62 tmp<errorEstimate<Type> >
  63 Sp
  64 (
  65     const tmp<volScalarField>& tsp,
  66     GeometricField<Type, fvPatchField, volMesh>& vf
  67 )
  68 {
  69     tmp<errorEstimate<Type> > tee = resError::Sp(tsp(), vf);
  70     tsp.clear();
  71     return tee;
  72 }
  73
  74
  75 template<class Type>
  76 tmp<errorEstimate<Type> >
  77 Sp
  78 (
  79     const dimensionedScalar& sp,
  80     GeometricField<Type, fvPatchField, volMesh>& vf
  81 )
  82 {
  83     return tmp<errorEstimate<Type> >
  84     (
  85         new errorEstimate<Type>
  86         (
  87             vf,
  88             sp.dimensions()*vf.dimensions()*dimVolume,
  89             sp.value()*vf.internalField(),
  90             scalarField(vf.internalField().size(), 0)
  91         )
  92     );
  93 }
  94
  95
  96 template<class Type>
  97 tmp<errorEstimate<Type> >
  98 SuSp
  99 (
 100     const volScalarField& sp,
 101     GeometricField<Type, fvPatchField, volMesh>& vf
 102 )
 103 {
 104     return Sp(sp, vf);
 105 }
 106
 107 template<class Type>
 108 tmp<errorEstimate<Type> >
 109 SuSp
 110 (
 111     const tmp<volScalarField>& tsp,
 112     GeometricField<Type, fvPatchField, volMesh>& vf
 113 )
 114 {
 115     tmp<errorEstimate<Type> > tee = resError::SuSp(tsp(), vf);
 116     tsp.clear();
 117     return tee;
 118 }
 119
 120
 121 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 122
 123 } // End namespace resError
 124
 125 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 126
 127 } // End namespace Foam
 128
 129 // ************************************************************************* //
 130
 131
 132 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 133