applications/utilities/parallelProcessing/decomposePar/lagrangianFieldDecomposerDecomposeFields.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2008 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 "lagrangianFieldDecomposer.H"
  28 #include "IOobjectList.H"
  29
  30 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  31
  32 namespace Foam
  33 {
  34
  35 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  36
  37 template<class Type>
  38 void lagrangianFieldDecomposer::readFields
  39 (
  40     const label cloudI,
  41     const IOobjectList& lagrangianObjects,
  42     PtrList<PtrList<IOField<Type> > >& lagrangianFields
  43 )
  44 {
  45     // Search list of objects for lagrangian fields
  46     IOobjectList lagrangianTypeObjects
  47     (
  48         lagrangianObjects.lookupClass(IOField<Type>::typeName)
  49     );
  50
  51     lagrangianFields.set
  52     (
  53         cloudI,
  54         new PtrList<IOField<Type> >
  55         (
  56             lagrangianTypeObjects.size()
  57         )
  58     );
  59
  60     label lagrangianFieldi=0;
  61     forAllIter(IOobjectList, lagrangianTypeObjects, iter)
  62     {
  63         lagrangianFields[cloudI].set
  64         (
  65             lagrangianFieldi++,
  66             new IOField<Type>(*iter())
  67         );
  68     }
  69 }
  70
  71
  72 template<class Type>
  73 tmp<IOField<Type> > lagrangianFieldDecomposer::decomposeField
  74 (
  75     const word& cloudName,
  76     const IOField<Type>& field
  77 ) const
  78 {
  79     // Create and map the internal field values
  80     Field<Type> procField(field, particleIndices_);
  81
  82     // Create the field for the processor
  83     return tmp<IOField<Type> >
  84     (
  85         new IOField<Type>
  86         (
  87             IOobject
  88             (
  89                 field.name(),
  90                 procMesh_.time().timeName(),
  91                 "lagrangian"/cloudName,
  92                 procMesh_,
  93                 IOobject::NO_READ,
  94                 IOobject::NO_WRITE
  95             ),
  96             procField
  97         )
  98     );
  99 }
 100
 101
 102 template<class GeoField>
 103 void lagrangianFieldDecomposer::decomposeFields
 104 (
 105     const word& cloudName,
 106     const PtrList<GeoField>& fields
 107 ) const
 108 {
 109     if (particleIndices_.size())
 110     {
 111         forAll (fields, fieldI)
 112         {
 113             decomposeField(cloudName, fields[fieldI])().write();
 114         }
 115     }
 116 }
 117
 118
 119 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 120
 121 } // End namespace Foam
 122
 123 // ************************************************************************* //