distributionNew.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-2026 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 "distribution.H"
27 
28 // * * * * * * * * * * * * * * * * Selectors * * * * * * * * * * * * * * * * //
29 
31 (
32  const unitSet& units,
33  const dictionary& dict,
34  const label sampleQ,
36  const bool report
37 )
38 {
39  const word distributionType(dict.lookup("type"));
40 
41  if (report)
42  {
43  Info<< indentOrNl << "Selecting " << typeName
44  << " type " << distributionType << endl;
45  }
46 
47  dictionaryConstructorTable::iterator cstrIter =
48  dictionaryConstructorTablePtr_->find(distributionType);
49 
50  if (cstrIter == dictionaryConstructorTablePtr_->end())
51  {
53  << "Unknown " << typeName << " type " << distributionType
54  << nl << nl << "Valid " << typeName << " types are:" << nl
55  << dictionaryConstructorTablePtr_->sortedToc()
56  << exit(FatalIOError);
57  }
58 
59  const dictionary& distributionDict =
60  dict.optionalSubDict(distributionType + "Distribution");
61 
62  printDictionary print(distributionDict);
63 
64  autoPtr<distribution> distributionPtr
65  (
66  cstrIter()
67  (
68  units,
69  distributionDict,
70  sampleQ,
71  std::move(rndGen)
72  )
73  );
74 
75  if (report)
76  {
77  Info<< indent << "min/average/max value = "
78  << distributionPtr->min() << '/'
79  << distributionPtr->mean() << '/'
80  << distributionPtr->max() << endl;
81  }
82 
83  return distributionPtr;
84 }
85 
86 
88 (
89  const unitSet& units,
90  const dictionary& dict,
91  const label sampleQ,
92  const randomGenerator::seed& s,
93  const bool global,
94  const bool report
95 )
96 {
97  return New(units, dict, sampleQ, randomGenerator(s, global), report);
98 }
99 
100 
102 (
103  autoPtr<distribution>& dPtr,
104  const label sampleQ
105 )
106 {
107  if (dPtr->sampleQ_ == sampleQ)
108  {
109  return autoPtr<distribution>(dPtr.ptr());
110  }
111  else
112  {
113  autoPtr<distribution> result(dPtr->clone(sampleQ));
114  dPtr.clear();
115  return result;
116  }
117 }
118 
119 
120 // ************************************************************************* //
T * ptr()
Return object pointer for reuse.
Definition: autoPtrI.H:90
void clear()
Delete object (if the pointer is valid) and set pointer to.
Definition: autoPtrI.H:126
A list of keywords followed by any number of values (e.g. words and numbers) or sub-dictionaries.
Definition: dictionary.H:162
label sampleQ() const
Access the sample size exponent.
Definition: distribution.C:134
static autoPtr< distribution > New(const unitSet &units, const dictionary &dict, const label sampleQ, randomGenerator &&rndGen, const bool report=true)
Select from dictionary and a random generator.
Enables the printing of a dictionary and subsequently looked-up defaulted entries.
Random number generator.
Template function which returns the un-mangled name of a given type. Useful for types which do not ha...
Unit conversion structure. Contains the associated dimensions and the multiplier with which to conver...
Definition: unitSet.H:68
A class for handling words, derived from string.
Definition: word.H:63
#define FatalIOErrorInFunction(ios)
Report an error message using Foam::FatalIOError.
Definition: error.H:346
gmvFile<< "tracers "<< particles.size()<< nl;forAllConstIter(lagrangian::Cloud< passiveParticle >, particles, iter){ gmvFile<< iter().position().x()<< " ";}gmvFile<< nl;forAllConstIter(lagrangian::Cloud< passiveParticle >, particles, iter){ gmvFile<< iter().position().y()<< " ";}gmvFile<< nl;forAllConstIter(lagrangian::Cloud< passiveParticle >, particles, iter){ gmvFile<< iter().position().z()<< " ";}gmvFile<< nl;forAll(lagrangianScalarNames, i){ word name=lagrangianScalarNames[i];IOField< scalar > s(IOobject(name, runTime.name(), lagrangian::cloud::prefix, mesh, IOobject::MUST_READ, IOobject::NO_WRITE))
errorManipArg< error, int > exit(error &err, const int errNo=1)
Definition: errorManip.H:124
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
Definition: label.H:59
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:288
messageStream Info
IOerror FatalIOError
Ostream & indentOrNl(Ostream &os)
Indent stream or add newline if indent level == 0.
Definition: Ostream.H:250
Ostream & indent(Ostream &os)
Indent stream.
Definition: Ostream.H:243
tmp< DimensionedField< TypeR, GeoMesh, Field > > New(const tmp< DimensionedField< TypeR, GeoMesh, Field >> &tdf1, const word &name, const dimensionSet &dimensions)
static const char nl
Definition: Ostream.H:297
dictionary dict
randomGenerator rndGen(653213)