src/topoChangerFvMesh/rawTopoChangerFvMesh/rawTopoChangerFvMeshTemplates.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
   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
  13     the Free Software Foundation, either version 3 of the License, or
  14     (at your 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, see <http://www.gnu.org/licenses/>.
  23
  24 \*---------------------------------------------------------------------------*/
  25
  26 #include "rawTopoChangerFvMesh.H"
  27 #include "Time.H"
  28
  29 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  30
  31 template<class Type, template<class> class PatchField, class GeoMesh>
  32 void Foam::rawTopoChangerFvMesh::setUnmappedValues
  33 (
  34     GeometricField<Type, PatchField, GeoMesh>& fld,
  35     const PackedBoolList& mappedFace,
  36     const GeometricField<Type, PatchField, GeoMesh>& baseFld
  37 )
  38 {
  39     //Pout<< "Checking field " << fld.name() << endl;
  40
  41     forAll(fld.boundaryField(), patchI)
  42     {
  43         PatchField<Type>& fvp = const_cast<PatchField<Type>&>
  44         (
  45             fld.boundaryField()[patchI]
  46         );
  47
  48         const label start = fvp.patch().start();
  49         forAll(fvp, i)
  50         {
  51             if (!mappedFace[start+i])
  52             {
  53                 //Pout<< "** Resetting unassigned value on patch "
  54                 //    << fvp.patch().name()
  55                 //    << " localface:" << i
  56                 //    << " to:" << baseFld.boundaryField()[patchI][i] << endl;
  57                 fvp[i] = baseFld.boundaryField()[patchI][i];
  58             }
  59         }
  60     }
  61 }
  62
  63
  64 template<class Type, template<class> class PatchField, class GeoMesh>
  65 void Foam::rawTopoChangerFvMesh::zeroUnmappedValues
  66 (
  67     const PackedBoolList& mappedFace
  68 ) const
  69 {
  70     typedef GeometricField<Type, PatchField, GeoMesh> FieldType;
  71
  72     const wordList fldNames(names(FieldType::typeName));
  73
  74     forAll(fldNames, i)
  75     {
  76         //Pout<< "Checking field " << fldNames[i] << endl;
  77
  78         FieldType& fld = const_cast<FieldType&>
  79         (
  80             lookupObject<FieldType>(fldNames[i])
  81         );
  82
  83         setUnmappedValues
  84         (
  85             fld,
  86             mappedFace,
  87             FieldType
  88             (
  89                 IOobject
  90                 (
  91                     "zero",
  92                     time().timeName(),
  93                     *this,
  94                     IOobject::NO_READ,
  95                     IOobject::NO_WRITE,
  96                     false
  97                 ),
  98                 *this,
  99                 dimensioned<Type>("0", fld.dimensions(), pTraits<Type>::zero)
 100             )
 101         );
 102     }
 103 }
 104
 105
 106 // ************************************************************************* //