applications/utilities/preProcessing/mapFields/MapLagrangianFields.H

   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 InNamespace
  26     Foam
  27
  28 Description
  29     Gets the indices of (source)particles that have been appended to the
  30     target cloud and maps the lagrangian fields accordingly.
  31
  32 \*---------------------------------------------------------------------------*/
  33
  34 #ifndef MapLagrangianFields_H
  35 #define MapLagrangianFields_H
  36
  37 #include <OpenFOAM/cloud.H>
  38 #include <OpenFOAM/GeometricField.H>
  39 #include <sampling/meshToMesh.H>
  40 #include <OpenFOAM/IOobjectList.H>
  41
  42 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  43
  44 namespace Foam
  45 {
  46
  47 //- Gets the indices of (source)particles that have been appended to the
  48 //  target cloud and maps the lagrangian fields accordingly.
  49 template<class Type>
  50 void MapLagrangianFields
  51 (
  52     const string& cloudName,
  53     const IOobjectList& objects,
  54     const meshToMesh& meshToMeshInterp,
  55     const labelList& addParticles
  56 )
  57 {
  58     const fvMesh& meshTarget = meshToMeshInterp.toMesh();
  59
  60     IOobjectList fields = objects.lookupClass(IOField<Type>::typeName);
  61
  62     forAllIter(IOobjectList, fields, fieldIter)
  63     {
  64         Info<< "    mapping lagrangian field " << fieldIter()->name() << endl;
  65
  66         // Read field (does not need mesh)
  67         IOField<Type> fieldSource(*fieldIter());
  68
  69         // Map
  70         IOField<Type> fieldTarget
  71         (
  72             IOobject
  73             (
  74                 fieldIter()->name(),
  75                 meshTarget.time().timeName(),
  76                 cloud::prefix/cloudName,
  77                 meshTarget,
  78                 IOobject::NO_READ,
  79                 IOobject::NO_WRITE,
  80                 false
  81             ),
  82             addParticles.size()
  83         );
  84         forAll(addParticles, i)
  85         {
  86             fieldTarget[i] = fieldSource[addParticles[i]];
  87         }
  88
  89         // Write field
  90         fieldTarget.write();
  91     }
  92 }
  93
  94
  95 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  96
  97 } // End namespace Foam
  98
  99 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 100
 101 #endif
 102
 103 // ************************ vim: set sw=4 sts=4 et: ************************ //