src/finiteVolume/fvMesh/extendedStencil/cellToFace/extendedCellToFaceStencil.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2008-2010 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 "extendedCellToFaceStencil.H"
  27 #include "globalIndex.H"
  28 #include "syncTools.H"
  29 #include "SortableList.H"
  30
  31 /* * * * * * * * * * * * * * * Static Member Data  * * * * * * * * * * * * * */
  32
  33 defineTypeNameAndDebug(Foam::extendedCellToFaceStencil, 0);
  34
  35
  36 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  37
  38 void Foam::extendedCellToFaceStencil::writeStencilStats
  39 (
  40     Ostream& os,
  41     const labelListList& stencil,
  42     const mapDistribute& map
  43 )
  44 {
  45     label sumSize = 0;
  46     label nSum = 0;
  47     label minSize = labelMax;
  48     label maxSize = labelMin;
  49
  50     forAll(stencil, i)
  51     {
  52         const labelList& sCells = stencil[i];
  53
  54         if (sCells.size() > 0)
  55         {
  56             sumSize += sCells.size();
  57             nSum++;
  58             minSize = min(minSize, sCells.size());
  59             maxSize = max(maxSize, sCells.size());
  60         }
  61     }
  62     reduce(sumSize, sumOp<label>());
  63     reduce(nSum, sumOp<label>());
  64
  65     reduce(minSize, minOp<label>());
  66     reduce(maxSize, maxOp<label>());
  67
  68     os  << "Stencil size :" << nl
  69         << "    average : " << scalar(sumSize)/nSum << nl
  70         << "    min     : " << minSize << nl
  71         << "    max     : " << maxSize << nl
  72         << endl;
  73
  74     // Sum all sent data
  75     label nSent = 0;
  76     label nLocal = 0;
  77     forAll(map.subMap(), procI)
  78     {
  79         if (procI != Pstream::myProcNo())
  80         {
  81             nSent += map.subMap()[procI].size();
  82         }
  83         else
  84         {
  85             nLocal += map.subMap()[procI].size();
  86         }
  87     }
  88
  89     os  << "Local data size : " << returnReduce(nLocal, sumOp<label>()) << nl
  90         << "Sent data size  : " << returnReduce(nSent, sumOp<label>()) << nl
  91         << endl;
  92 }
  93
  94
  95 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  96
  97 Foam::extendedCellToFaceStencil::extendedCellToFaceStencil(const polyMesh& mesh)
  98 :
  99     mesh_(mesh)
 100 {
 101     // Check for transformation - not supported.
 102     const polyBoundaryMesh& patches = mesh.boundaryMesh();
 103
 104     forAll(patches, patchI)
 105     {
 106         if (isA<coupledPolyPatch>(patches[patchI]))
 107         {
 108             const coupledPolyPatch& cpp =
 109                 refCast<const coupledPolyPatch>(patches[patchI]);
 110
 111             if (!cpp.parallel() || cpp.separated())
 112             {
 113                 FatalErrorIn
 114                 (
 115                     "extendedCellToFaceStencil::extendedCellToFaceStencil"
 116                     "(const polyMesh&)"
 117                 )   << "Coupled patches with transformations not supported."
 118                     << endl
 119                     << "Problematic patch " << cpp.name() << exit(FatalError);
 120             }
 121         }
 122     }
 123 }
 124
 125
 126 // ************************************************************************* //