35 namespace distributions
65 if (i && values[i].
first() - values[i - 1].
first() < 0)
68 << typeName <<
": The probability density function is not "
72 if (values[i].
second() < 0)
75 << typeName <<
": The probability density function is not "
84 x_[i] = values[i].
first();
91 PDF_[i] =
integerPow(x_[i],
q())*values[i].second();
120 if (
q() == d.
q())
return;
171 const scalar x0 =
min(), x1 =
max(), d = 0.1*(x1 - x0);
176 result[0] =
Foam::max(x0 - d, q() < 0 ? x0/2 : rootVSmall);
181 result[x_.size() + 2] = x1;
182 result[x_.size() + 3] = x1 + d;
#define forAll(list, i)
Loop across all elements in list.
Macros for easy insertion into run-time selection tables.
A 1D array of objects of type <T>, where the size of the vector is known and used for subscript bound...
void resize(const label)
Alias for setSize(const label)
void size(const label)
Override size to be inconsistent with allocated storage.
Pre-declare related SubField type.
T & first()
Return the first element of the list.
T & last()
Return the last element of the list.
A list of keyword definitions, which are a keyword followed by any number of values (e....
Accumulating histogram of values. Specified bin resolution automatic generation of bins.
void report() const
Report.
virtual void validatePositive(const dictionary &dict) const
Validate that the lower bound is positive.
label q() const
Return the effective distribution size exponent.
virtual void validateBounds(const dictionary &dict) const
Validate that the bounds are monotonic.
Distribution in which the probability density function is given as a table of values.
virtual scalar min() const
Return the minimum value.
virtual scalar sample() const
Sample the distribution.
virtual tmp< scalarField > PDF(const scalarField &x) const
Return the distribution probability density function.
virtual tmp< scalarField > x(const label n) const
Return coordinates to plot across the range of the distribution.
tabulatedDensity(const dictionary &dict, Random &rndGen, const label sampleQ)
Construct from a dictionary.
virtual ~tabulatedDensity()
Destructor.
virtual scalar max() const
Return the maximum value.
virtual scalar mean() const
Return the mean value.
scalar sample() const
Sample the distribution.
static tmp< scalarField > integrate(const scalarField &x, const scalarField &y)
Integrate the values y with respect to the coordinates x.
static tmp< scalarField > integrateX(const scalarField &x, const scalarField &y)
Integrate the values x*y with respect to the coordinates x.
A class for managing temporary objects.
T & ref() const
Return non-const reference or generate a fatal error.
#define FatalIOErrorInFunction(ios)
Report an error message using Foam::FatalIOError.
defineTypeNameAndDebug(exponential, 0)
addToRunTimeSelectionTable(distribution, exponential, dictionary)
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
errorManip< error > abort(error &err)
scalar integerPow(const scalar x, const label e)
Compute the power of the number x to the integer e.
labelList second(const UList< labelPair > &p)
Field< scalar > scalarField
Specialisation of Field<T> for scalar.
labelList first(const UList< labelPair > &p)
layerAndWeight min(const layerAndWeight &a, const layerAndWeight &b)
layerAndWeight max(const layerAndWeight &a, const layerAndWeight &b)