src/decompositionMethods/decompositionMethods/decompositionMethod/decompositionMethod.C

   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 InClass
  26     decompositionMethod
  27
  28 Description
  29
  30 \*---------------------------------------------------------------------------*/
  31
  32 #include "decompositionMethod.H"
  33
  34 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  35
  36 namespace Foam
  37 {
  38     defineTypeNameAndDebug(decompositionMethod, 0);
  39     defineRunTimeSelectionTable(decompositionMethod, dictionary);
  40     defineRunTimeSelectionTable(decompositionMethod, dictionaryMesh);
  41 }
  42
  43 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  44
  45 Foam::autoPtr<Foam::decompositionMethod> Foam::decompositionMethod::New
  46 (
  47     const dictionary& decompositionDict
  48 )
  49 {
  50     word decompositionMethodTypeName(decompositionDict.lookup("method"));
  51
  52     Info<< "Selecting decompositionMethod "
  53         << decompositionMethodTypeName << endl;
  54
  55     dictionaryConstructorTable::iterator cstrIter =
  56         dictionaryConstructorTablePtr_->find(decompositionMethodTypeName);
  57
  58     if (cstrIter == dictionaryConstructorTablePtr_->end())
  59     {
  60         FatalErrorIn
  61         (
  62             "decompositionMethod::New"
  63             "(const dictionary& decompositionDict)"
  64         )   << "Unknown decompositionMethod "
  65             << decompositionMethodTypeName << endl << endl
  66             << "Valid decompositionMethods are : " << endl
  67             << dictionaryConstructorTablePtr_->toc()
  68             << exit(FatalError);
  69     }
  70
  71     return autoPtr<decompositionMethod>(cstrIter()(decompositionDict));
  72 }
  73
  74
  75 Foam::autoPtr<Foam::decompositionMethod> Foam::decompositionMethod::New
  76 (
  77     const dictionary& decompositionDict,
  78     const polyMesh& mesh
  79 )
  80 {
  81     word decompositionMethodTypeName(decompositionDict.lookup("method"));
  82
  83     Info<< "Selecting decompositionMethod "
  84         << decompositionMethodTypeName << endl;
  85
  86     dictionaryMeshConstructorTable::iterator cstrIter =
  87         dictionaryMeshConstructorTablePtr_->find(decompositionMethodTypeName);
  88
  89     if (cstrIter == dictionaryMeshConstructorTablePtr_->end())
  90     {
  91         FatalErrorIn
  92         (
  93             "decompositionMethod::New"
  94             "(const dictionary& decompositionDict, "
  95             "const polyMesh& mesh)"
  96         )   << "Unknown decompositionMethod "
  97             << decompositionMethodTypeName << endl << endl
  98             << "Valid decompositionMethods are : " << endl
  99             << dictionaryMeshConstructorTablePtr_->toc()
 100             << exit(FatalError);
 101     }
 102
 103     return autoPtr<decompositionMethod>(cstrIter()(decompositionDict, mesh));
 104 }
 105
 106
 107 Foam::labelList Foam::decompositionMethod::decompose
 108 (
 109     const labelList& fineToCoarse,
 110     const pointField& coarsePoints
 111 )
 112 {
 113     // Decompose based on agglomerated points
 114     labelList coarseDistribution(decompose(coarsePoints));
 115
 116     // Rework back into decomposition for original mesh_
 117     labelList fineDistribution(fineToCoarse.size());
 118
 119     forAll(fineDistribution, i)
 120     {
 121         fineDistribution[i] = coarseDistribution[fineToCoarse[i]];
 122     }
 123
 124     return fineDistribution;
 125 }
 126
 127
 128 void Foam::decompositionMethod::calcCellCells
 129 (
 130     const polyMesh& mesh,
 131     const labelList& fineToCoarse,
 132     const label nCoarse,
 133     labelListList& cellCells
 134 )
 135 {
 136     if (fineToCoarse.size() != mesh.nCells())
 137     {
 138         FatalErrorIn
 139         (
 140             "decompositionMethod::calcCellCells"
 141             "(const labelList&, labelListList&) const"
 142         )   << "Only valid for mesh agglomeration." << exit(FatalError);
 143     }
 144
 145     List<DynamicList<label> > dynCellCells(nCoarse);
 146
 147     forAll(mesh.faceNeighbour(), faceI)
 148     {
 149         label own = fineToCoarse[mesh.faceOwner()[faceI]];
 150         label nei = fineToCoarse[mesh.faceNeighbour()[faceI]];
 151
 152         if (own != nei)
 153         {
 154             if (findIndex(dynCellCells[own], nei) == -1)
 155             {
 156                 dynCellCells[own].append(nei);
 157             }
 158             if (findIndex(dynCellCells[nei], own) == -1)
 159             {
 160                 dynCellCells[nei].append(own);
 161             }
 162         }
 163     }
 164
 165     cellCells.setSize(dynCellCells.size());
 166     forAll(dynCellCells, coarseI)
 167     {
 168         cellCells[coarseI].transfer(dynCellCells[coarseI]);
 169     }
 170 }
 171
 172
 173 // ************************************************************************* //