src/sampling/sampledSurface/sampledTriSurfaceMesh/sampledTriSurfaceMeshTemplates.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2010-2011 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
  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 "sampledTriSurfaceMesh.H"
  27
  28 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  29
  30 template <class Type>
  31 Foam::tmp<Foam::Field<Type> >
  32 Foam::sampledTriSurfaceMesh::sampleField
  33 (
  34     const GeometricField<Type, fvPatchField, volMesh>& vField
  35 ) const
  36 {
  37     // One value per face
  38     tmp<Field<Type> > tvalues(new Field<Type>(sampleElements_.size()));
  39     Field<Type>& values = tvalues();
  40
  41     if (sampleSource_ == cells)
  42     {
  43         // Sample cells
  44
  45         forAll(sampleElements_, triI)
  46         {
  47             values[triI] = vField[sampleElements_[triI]];
  48         }
  49     }
  50     else
  51     {
  52         // Sample boundary faces
  53
  54         const polyBoundaryMesh& pbm = mesh().boundaryMesh();
  55         label nBnd = mesh().nFaces()-mesh().nInternalFaces();
  56
  57         // Create flat boundary field
  58
  59         Field<Type> bVals(nBnd, pTraits<Type>::zero);
  60
  61         forAll(vField.boundaryField(), patchI)
  62         {
  63             label bFaceI = pbm[patchI].start() - mesh().nInternalFaces();
  64
  65             SubList<Type>
  66             (
  67                 bVals,
  68                 vField.boundaryField()[patchI].size(),
  69                 bFaceI
  70             ).assign(vField.boundaryField()[patchI]);
  71         }
  72
  73         // Sample in flat boundary field
  74
  75         forAll(sampleElements_, triI)
  76         {
  77             label faceI = sampleElements_[triI];
  78             values[triI] = bVals[faceI-mesh().nInternalFaces()];
  79         }
  80     }
  81
  82     return tvalues;
  83 }
  84
  85
  86 template <class Type>
  87 Foam::tmp<Foam::Field<Type> >
  88 Foam::sampledTriSurfaceMesh::interpolateField
  89 (
  90     const interpolation<Type>& interpolator
  91 ) const
  92 {
  93     // One value per vertex
  94     tmp<Field<Type> > tvalues(new Field<Type>(sampleElements_.size()));
  95     Field<Type>& values = tvalues();
  96
  97     if (sampleSource_ == cells)
  98     {
  99         // Sample cells.
 100
 101         forAll(sampleElements_, pointI)
 102         {
 103             values[pointI] = interpolator.interpolate
 104             (
 105                 samplePoints_[pointI],
 106                 sampleElements_[pointI]
 107             );
 108         }
 109     }
 110     else
 111     {
 112         // Sample boundary faces.
 113
 114         forAll(samplePoints_, pointI)
 115         {
 116             label faceI = sampleElements_[pointI];
 117
 118             values[pointI] = interpolator.interpolate
 119             (
 120                 samplePoints_[pointI],
 121                 mesh().faceOwner()[faceI],
 122                 faceI
 123             );
 124         }
 125     }
 126
 127     return tvalues;
 128 }
 129
 130
 131 // ************************************************************************* //