applications/test/DimensionedField/DimensionedFieldTest.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2009 OpenCFD Ltd.
   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 the
  13     Free Software Foundation; either version 2 of the License, or (at your
  14     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, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 \*---------------------------------------------------------------------------*/
  26
  27 #include "argList.H"
  28 #include "Time.H"
  29 #include "DimensionedFields.H"
  30 #include "DimensionedSphericalTensorField.H"
  31 #include "vector.H"
  32 #include "tensor.H"
  33 #include "GeoMesh.H"
  34
  35 using namespace Foam;
  36
  37 namespace Foam
  38 {
  39
  40 class vMesh
  41 {
  42
  43 public:
  44
  45     vMesh()
  46     {}
  47
  48     label size() const
  49     {
  50         return 10;
  51     }
  52 };
  53
  54 };
  55
  56 template<>
  57 const word Foam::DimensionedField<scalar, GeoMesh<vMesh> >::typeName
  58 (
  59     "dimenionedScalarField"
  60 );
  61
  62 template<>
  63 const word Foam::DimensionedField<vector, GeoMesh<vMesh> >::typeName
  64 (
  65     "dimenionedVectorField"
  66 );
  67
  68 template<>
  69 const word Foam::DimensionedField<tensor, GeoMesh<vMesh> >::typeName
  70 (
  71     "dimenionedTensorField"
  72 );
  73
  74 template<>
  75 const word Foam::DimensionedField<sphericalTensor, GeoMesh<vMesh> >::typeName
  76 (
  77     "dimenionedSphericalTensorField"
  78 );
  79
  80
  81 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  82 // Main program:
  83
  84 int main(int argc, char *argv[])
  85 {
  86
  87 #   include "setRootCase.H"
  88 #   include "createTime.H"
  89
  90     vMesh vm;
  91
  92     DimensionedField<scalar, GeoMesh<vMesh> > dsf
  93     (
  94         IOobject
  95         (
  96             "dsf",
  97             runTime.timeName(),
  98             runTime,
  99             IOobject::MUST_READ,
 100             IOobject::NO_WRITE
 101         ),
 102         vm
 103     );
 104
 105     Info<< dsf << endl;
 106     dsf += dsf;
 107     dsf -= dimensionedScalar("5", dsf.dimensions(), 5.0);
 108     Info<< dsf << endl;
 109
 110     Info<< sqr(dsf + dsf) - sqr(dsf + dsf) << endl;
 111
 112     DimensionedField<vector, GeoMesh<vMesh> > dvf
 113     (
 114         IOobject
 115         (
 116             "dvf",
 117             runTime.timeName(),
 118             runTime,
 119             IOobject::MUST_READ,
 120             IOobject::NO_WRITE
 121         ),
 122         vm
 123     );
 124
 125     Info<< (dvf ^ (dvf ^ dvf)) << endl;
 126
 127     Info << "End\n" << endl;
 128
 129     return 0;
 130 }
 131
 132
 133 // ************************************************************************* //