src/sampling/sampledSet/cloud/cloudSet.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 "cloudSet.H"
  27 #include "sampledSet.H"
  28 #include "meshSearch.H"
  29 #include "DynamicList.H"
  30 #include "polyMesh.H"
  31 #include "addToRunTimeSelectionTable.H"
  32 #include "word.H"
  33
  34 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  35
  36 namespace Foam
  37 {
  38     defineTypeNameAndDebug(cloudSet, 0);
  39     addToRunTimeSelectionTable(sampledSet, cloudSet, word);
  40 }
  41
  42
  43 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  44
  45 void Foam::cloudSet::calcSamples
  46 (
  47     DynamicList<point>& samplingPts,
  48     DynamicList<label>& samplingCells,
  49     DynamicList<label>& samplingFaces,
  50     DynamicList<label>& samplingSegments,
  51     DynamicList<scalar>& samplingCurveDist
  52 ) const
  53 {
  54     // Force calculation of face-diagonals
  55     (void)mesh().tetBasePtIs();
  56
  57     const meshSearch& queryMesh = searchEngine();
  58
  59     forAll(sampleCoords_, sampleI)
  60     {
  61         label cellI = queryMesh.findCell(sampleCoords_[sampleI]);
  62
  63         if (cellI != -1)
  64         {
  65             samplingPts.append(sampleCoords_[sampleI]);
  66             samplingCells.append(cellI);
  67             samplingFaces.append(-1);
  68             samplingSegments.append(0);
  69             samplingCurveDist.append(1.0 * sampleI);
  70         }
  71     }
  72 }
  73
  74
  75 void Foam::cloudSet::genSamples()
  76 {
  77     // Storage for sample points
  78     DynamicList<point> samplingPts;
  79     DynamicList<label> samplingCells;
  80     DynamicList<label> samplingFaces;
  81     DynamicList<label> samplingSegments;
  82     DynamicList<scalar> samplingCurveDist;
  83
  84     calcSamples
  85     (
  86         samplingPts,
  87         samplingCells,
  88         samplingFaces,
  89         samplingSegments,
  90         samplingCurveDist
  91     );
  92
  93     samplingPts.shrink();
  94     samplingCells.shrink();
  95     samplingFaces.shrink();
  96     samplingSegments.shrink();
  97     samplingCurveDist.shrink();
  98
  99     setSamples
 100     (
 101         samplingPts,
 102         samplingCells,
 103         samplingFaces,
 104         samplingSegments,
 105         samplingCurveDist
 106     );
 107 }
 108
 109
 110 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 111
 112 Foam::cloudSet::cloudSet
 113 (
 114     const word& name,
 115     const polyMesh& mesh,
 116     meshSearch& searchEngine,
 117     const word& axis,
 118     const List<point>& sampleCoords
 119 )
 120 :
 121     sampledSet(name, mesh, searchEngine, axis),
 122     sampleCoords_(sampleCoords)
 123 {
 124     genSamples();
 125
 126     if (debug)
 127     {
 128         write(Info);
 129     }
 130 }
 131
 132
 133 Foam::cloudSet::cloudSet
 134 (
 135     const word& name,
 136     const polyMesh& mesh,
 137     meshSearch& searchEngine,
 138     const dictionary& dict
 139 )
 140 :
 141     sampledSet(name, mesh, searchEngine, dict),
 142     sampleCoords_(dict.lookup("points"))
 143 {
 144     genSamples();
 145
 146     if (debug)
 147     {
 148         write(Info);
 149     }
 150 }
 151
 152
 153 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 154
 155 Foam::cloudSet::~cloudSet()
 156 {}
 157
 158
 159 // ************************************************************************* //