manual.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration | Website: https://openfoam.org
5  \\ / A nd | Copyright (C) 2011-2024 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 "manual.H"
28 #include "IFstream.H"
29 #include "labelIOList.H"
30 
31 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
32 
33 namespace Foam
34 {
35 namespace decompositionMethods
36 {
38 
40  (
42  manual,
43  decomposer
44  );
45 
47  (
49  manual,
50  distributor
51  );
52 }
53 }
54 
55 
56 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
57 
59 :
60  decompositionMethod(decompositionDict),
61  decompDataFile_
62  (
63  decompositionDict.optionalSubDict
64  (
65  word(decompositionDict.lookup("method")) + "Coeffs"
66  ).lookup("dataFile")
67  )
68 {}
69 
70 
71 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
72 
74 (
75  const polyMesh& mesh,
76  const pointField& points,
77  const scalarField& pointWeights
78 )
79 {
80  labelIOList finalDecomp
81  (
82  IOobject
83  (
84  decompDataFile_,
85  mesh.facesInstance(),
86  mesh,
89  false
90  )
91  );
92 
93  // check if the final decomposition is OK
94 
95  if (finalDecomp.size() != points.size())
96  {
98  << "Size of decomposition list does not correspond "
99  << "to the number of points. Size: "
100  << finalDecomp.size() << " Number of points: "
101  << points.size()
102  << ".\n" << "Manual decomposition data read from file "
103  << decompDataFile_ << "." << endl
104  << exit(FatalError);
105  }
106 
107  if (min(finalDecomp) < 0 || max(finalDecomp) > nProcessors_ - 1)
108  {
110  << "According to the decomposition, cells assigned to "
111  << "impossible processor numbers. Min processor = "
112  << min(finalDecomp) << " Max processor = " << max(finalDecomp)
113  << ".\n" << "Manual decomposition data read from file "
114  << decompDataFile_ << "." << endl
115  << exit(FatalError);
116  }
117 
118  return finalDecomp;
119 }
120 
121 
122 // ************************************************************************* //
Macros for easy insertion into run-time selection tables.
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...
Definition: IOobject.H:99
Abstract base class for decomposition.
Decomposition given a cell-to-processor association in a file.
Definition: manual.H:51
virtual labelList decompose(const polyMesh &mesh, const pointField &cellCentres, const scalarField &cellWeights)
Return for every coordinate the wanted processor number. Use the.
Definition: manual.C:74
manual(const dictionary &decompositionDict)
Construct given the decomposition dictionary.
Definition: manual.C:58
A list of keyword definitions, which are a keyword followed by any number of values (e....
Definition: dictionary.H:162
Mesh consisting of general polyhedral cells.
Definition: polyMesh.H:80
const fileName & facesInstance() const
Return the current instance directory for faces.
Definition: polyMesh.C:971
A class for handling words, derived from string.
Definition: word.H:62
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
Definition: error.H:334
const pointField & points
addToRunTimeSelectionTable(decompositionMethod, metis, decomposer)
Namespace for OpenFOAM.
errorManipArg< error, int > exit(error &err, const int errNo=1)
Definition: errorManip.H:124
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:257
layerAndWeight min(const layerAndWeight &a, const layerAndWeight &b)
layerAndWeight max(const layerAndWeight &a, const layerAndWeight &b)
error FatalError