1 /*---------------------------------------------------------------------------*\
3 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
5 \\ / A nd | Copyright (C) 1991-2009 OpenCFD Ltd.
7 -------------------------------------------------------------------------------
9 This file is part of OpenFOAM.
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.
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
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
29 Metis domain decomposition
34 \*---------------------------------------------------------------------------*/
39 #include "decompositionMethod.H"
44 /*---------------------------------------------------------------------------*\
45 Class metisDecomp Declaration
46 \*---------------------------------------------------------------------------*/
50 public decompositionMethod
54 const polyMesh& mesh_;
57 // Private Member Functions
61 const List<int>& adjncy,
62 const List<int>& xadj,
63 List<int>& finalDecomp
66 //- Disallow default bitwise copy construct and assignment
67 void operator=(const metisDecomp&);
68 metisDecomp(const metisDecomp&);
73 //- Runtime type information
79 //- Construct given the decomposition dictionary and mesh
82 const dictionary& decompositionDict,
89 virtual ~metisDecomp()
95 virtual bool parallelAware() const
97 // Metis does not know about proc boundaries
101 //- Return for every coordinate the wanted processor number. Use the
102 // mesh connectivity (if needed)
103 virtual labelList decompose(const pointField&);
105 //- Return for every coordinate the wanted processor number. Gets
106 // passed agglomeration map (from fine to coarse cells) and coarse cell
107 // location. Can be overridden by decomposers that provide this
108 // functionality natively.
109 virtual labelList decompose
111 const labelList& agglom,
115 //- Return for every coordinate the wanted processor number. Explicitly
116 // provided mesh connectivity.
117 // The connectivity is equal to mesh.cellCells() except for
118 // - in parallel the cell numbers are global cell numbers (starting
119 // from 0 at processor0 and then incrementing all through the
121 // - the connections are across coupled patches
122 virtual labelList decompose
124 const labelListList& globalCellCells,
128 //- Helper to convert cellcells into Metis storage
129 static void calcMetisCSR
131 const labelListList& globalCellCells,
139 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
141 } // End namespace Foam
143 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
147 // ************************************************************************* //