DistributionTrapezoid Class Reference

Description

Trapezoidal distribution.

Definition at line 304 of file Distributions.h.

Inheritance diagram for DistributionTrapezoid:

Collaboration diagram for DistributionTrapezoid:

Public Member Functions
	DistributionTrapezoid ()
	DistributionTrapezoid (double center, double left, double middle, double right)
	DistributionTrapezoid (std::vector< double > P)
void	checkNodeArgs () const
	Raises exception if a parameter value is invalid. More...
std::string	className () const final
	Returns the class name, to be hard-coded in each leaf class that inherits from INode. More...
DistributionTrapezoid *	clone () const override
std::vector< double >	equidistantPoints (size_t nbr_samples, double sigma_factor, const RealLimits &limits=RealLimits()) const override
	generate list of sample values More...
virtual std::vector< double >	equidistantPointsInRange (size_t nbr_samples, double xmin, double xmax) const
	Returns equidistant interpolation points from xmin to xmax. More...
std::vector< ParameterSample >	equidistantSamples (size_t nbr_samples, double sigma_factor=0., const RealLimits &limits=RealLimits()) const
	Returns equidistant samples, using intrinsic parameters, weighted with probabilityDensity(). More...
std::vector< ParameterSample >	equidistantSamplesInRange (size_t nbr_samples, double xmin, double xmax) const
	Returns equidistant samples from xmin to xmax, weighted with probabilityDensity(). More...
double	getLeftWidth () const
double	getMiddleWidth () const
double	getRightWidth () const
bool	isDelta () const override
	Returns true if the distribution is in the limit case of a Dirac delta distribution. More...
double	mean () const override
	Returns the distribution-specific mean. More...
virtual std::vector< const INode * >	nodeChildren () const
	Returns all children. More...
std::vector< const INode * >	nodeOffspring () const
	Returns all descendants. More...
std::vector< ParaMeta >	parDefs () const final
	Returns the parameter definitions, to be hard-coded in each leaf class. More...
double	probabilityDensity (double x) const override
	Returns the distribution-specific probability density for value x. More...
std::string	pythonConstructor (const std::string &units) const override
	Prints distribution with constructor parameters in given units. ba.DistributionGaussian(2.0*deg, 0.02*deg) More...
virtual void	transferToCPP ()
	Used for Python overriding of clone (see swig/tweaks.py) More...

Protected Member Functions
void	adjustMinMaxForLimits (double &xmin, double &xmax, const RealLimits &limits) const
	modifies xmin and xmax if they are outside of limits More...
std::vector< ParameterSample >	generateSamplesFromValues (const std::vector< double > &sample_values) const
	Returns weighted samples from given interpolation points and probabilityDensity(). More...

Protected Attributes
std::vector< double >	m_P

Private Member Functions
void	adjustLimitsToNonZeroSamples (double &min, double &max, size_t nbr_samples) const

Private Attributes
const double &	m_center
const double &	m_left
const double &	m_middle
const double &	m_right

Constructor & Destructor Documentation

DistributionTrapezoid() [1/3]

DistributionTrapezoid::DistributionTrapezoid

(

std::vector< double >

P

)

Definition at line 394 of file Distributions.cpp.

    : IDistribution1D(P)
    , m_center(m_P[0])
    , m_left(m_P[1])
    , m_middle(m_P[2])
    , m_right(m_P[3])
{
    checkNodeArgs();
    if (m_left < 0.0)
        throw std::runtime_error("DistributionTrapezoid: leftWidth < 0");
    if (m_middle < 0.0)
        throw std::runtime_error("DistributionTrapezoid: middleWidth < 0");
    if (m_right < 0.0)
        throw std::runtime_error("DistributionTrapezoid: rightWidth < 0");
}

References INode::checkNodeArgs(), m_left, m_middle, and m_right.

Here is the call graph for this function:

DistributionTrapezoid() [2/3]

DistributionTrapezoid::DistributionTrapezoid	(	double	center,
		double	left,
		double	middle,
		double	right
	)

Definition at line 410 of file Distributions.cpp.

    : DistributionTrapezoid(std::vector<double>{center, left, middle, right})
{
}

DistributionTrapezoid() [3/3]

DistributionTrapezoid::DistributionTrapezoid

(

)

Definition at line 416 of file Distributions.cpp.

    : DistributionTrapezoid(0., 0., 1., 0.)
{
}

Referenced by clone().

Member Function Documentation

adjustLimitsToNonZeroSamples()

void DistributionTrapezoid::adjustLimitsToNonZeroSamples ( double &  min, double &  max, size_t  nbr_samples  ) const

private

Definition at line 457 of file Distributions.cpp.

{
    if (nbr_samples <= 1)
        return;
    size_t N = nbr_samples;
    if (m_left > 0.0)
        ++N;
    if (m_right > 0.0)
        ++N;
    if (N == nbr_samples)
        return;
    double step = (max - min) / (N - 1);
    if (m_left > 0.0)
        min += step;
    if (m_right > 0.0)
        max -= step;
}

References m_left, m_right, and N.

Referenced by equidistantPoints().

adjustMinMaxForLimits()

void IDistribution1D::adjustMinMaxForLimits ( double &  xmin, double &  xmax, const RealLimits &  limits  ) const

protectedinherited

modifies xmin and xmax if they are outside of limits

Definition at line 81 of file Distributions.cpp.

{
    if (limits.hasLowerLimit() && xmin < limits.lowerLimit())
        xmin = limits.lowerLimit();
    if (limits.hasUpperLimit() && xmax > limits.upperLimit())
        xmax = limits.upperLimit();
    if (xmin > xmax) {
        std::ostringstream ostr;
        ostr << "IDistribution1D::adjustMinMaxForLimits() -> Error. Can't' adjust ";
        ostr << "xmin:" << xmin << " xmax:" << xmax << " for given limits " << limits << std::endl;
        throw std::runtime_error(ostr.str());
    }
}

References RealLimits::hasLowerLimit(), RealLimits::hasUpperLimit(), RealLimits::lowerLimit(), and RealLimits::upperLimit().

Referenced by DistributionGate::equidistantPoints(), DistributionLorentz::equidistantPoints(), DistributionGaussian::equidistantPoints(), DistributionLogNormal::equidistantPoints(), DistributionCosine::equidistantPoints(), and equidistantPoints().

Here is the call graph for this function:

checkNodeArgs()

void INode::checkNodeArgs ( ) const

inherited

Raises exception if a parameter value is invalid.

Definition at line 27 of file INode.cpp.

{
    size_t nP = m_P.size();
    if (parDefs().size() != nP) {
        std::cerr << "BUG in class " << className() << std::endl;
        std::cerr << "#m_P = " << nP << std::endl;
        std::cerr << "#PDf = " << parDefs().size() << std::endl;
        for (const ParaMeta& pm : parDefs())
            std::cerr << "  PDf: " << pm.name << std::endl;
        ASSERT(0);
    }
    ASSERT(parDefs().size() == nP);
    for (size_t i = 0; i < nP; ++i) {
        const ParaMeta pm = parDefs()[i];
 
        RealLimits limits = RealLimits::limitless();
        if (pm.vMin == -INF) {
            ASSERT(pm.vMax == +INF);
            // nothing to do
        } else if (pm.vMax == +INF) {
            ASSERT(pm.vMin == 0);
            limits = RealLimits::nonnegative();
        } else {
            limits = RealLimits::limited(pm.vMin, pm.vMax);
        }
        limits.check(pm.name, m_P[i]);
    }
}

References ASSERT, RealLimits::check(), INode::className(), INF, RealLimits::limited(), RealLimits::limitless(), INode::m_P, ParaMeta::name, RealLimits::nonnegative(), INode::parDefs(), ParaMeta::vMax, and ParaMeta::vMin.

Referenced by BarGauss::BarGauss(), BarLorentz::BarLorentz(), Bipyramid4::Bipyramid4(), Box::Box(), CantellatedCube::CantellatedCube(), Cone::Cone(), ConstantBackground::ConstantBackground(), CosineRippleBox::CosineRippleBox(), CosineRippleGauss::CosineRippleGauss(), CosineRippleLorentz::CosineRippleLorentz(), Cylinder::Cylinder(), DistributionCosine::DistributionCosine(), DistributionGate::DistributionGate(), DistributionGaussian::DistributionGaussian(), DistributionLogNormal::DistributionLogNormal(), DistributionLorentz::DistributionLorentz(), DistributionTrapezoid(), Dodecahedron::Dodecahedron(), EllipsoidalCylinder::EllipsoidalCylinder(), FootprintGauss::FootprintGauss(), FootprintSquare::FootprintSquare(), FuzzySphere::FuzzySphere(), GaussSphere::GaussSphere(), HemiEllipsoid::HemiEllipsoid(), HollowSphere::HollowSphere(), HorizontalCylinder::HorizontalCylinder(), Icosahedron::Icosahedron(), LongBoxGauss::LongBoxGauss(), LongBoxLorentz::LongBoxLorentz(), PlatonicOctahedron::PlatonicOctahedron(), PlatonicTetrahedron::PlatonicTetrahedron(), Prism3::Prism3(), Prism6::Prism6(), Profile1DCauchy::Profile1DCauchy(), Profile1DCosine::Profile1DCosine(), Profile1DGate::Profile1DGate(), Profile1DGauss::Profile1DGauss(), Profile1DTriangle::Profile1DTriangle(), Profile1DVoigt::Profile1DVoigt(), Profile2DCauchy::Profile2DCauchy(), Profile2DCone::Profile2DCone(), Profile2DGate::Profile2DGate(), Profile2DGauss::Profile2DGauss(), Profile2DVoigt::Profile2DVoigt(), Pyramid2::Pyramid2(), Pyramid3::Pyramid3(), Pyramid4::Pyramid4(), Pyramid6::Pyramid6(), RotationEuler::RotationEuler(), RotationX::RotationX(), RotationY::RotationY(), RotationZ::RotationZ(), SawtoothRippleBox::SawtoothRippleBox(), SawtoothRippleGauss::SawtoothRippleGauss(), SawtoothRippleLorentz::SawtoothRippleLorentz(), Sphere::Sphere(), Spheroid::Spheroid(), TruncatedCube::TruncatedCube(), TruncatedSphere::TruncatedSphere(), and TruncatedSpheroid::TruncatedSpheroid().

Here is the call graph for this function:

className()

std::string DistributionTrapezoid::className ( ) const

inlinefinalvirtual

Returns the class name, to be hard-coded in each leaf class that inherits from INode.

Implements INode.

Definition at line 314 of file Distributions.h.

{ return "DistributionTrapezoid"; }

Referenced by pythonConstructor().

clone()

DistributionTrapezoid* DistributionTrapezoid::clone ( ) const

inlineoverridevirtual

Implements IDistribution1D.

Definition at line 310 of file Distributions.h.

    {
        return new DistributionTrapezoid(m_center, m_left, m_middle, m_right);
    }

References DistributionTrapezoid(), m_center, m_left, m_middle, and m_right.

Here is the call graph for this function:

equidistantPoints()

std::vector< double > DistributionTrapezoid::equidistantPoints ( size_t  nbr_samples, double  sigma_factor, const RealLimits &  limits = RealLimits()  ) const

overridevirtual

generate list of sample values

Implements IDistribution1D.

Definition at line 436 of file Distributions.cpp.

{
    double xmin = m_center - m_middle / 2.0 - m_left;
    double xmax = xmin + m_left + m_middle + m_right;
    adjustLimitsToNonZeroSamples(xmin, xmax, nbr_samples);
    adjustMinMaxForLimits(xmin, xmax, limits);
    return equidistantPointsInRange(nbr_samples, xmin, xmax);
}

References adjustLimitsToNonZeroSamples(), IDistribution1D::adjustMinMaxForLimits(), IDistribution1D::equidistantPointsInRange(), m_center, m_left, m_middle, and m_right.

Here is the call graph for this function:

equidistantPointsInRange()

std::vector< double > IDistribution1D::equidistantPointsInRange ( size_t  nbr_samples, double  xmin, double  xmax  ) const

virtualinherited

Returns equidistant interpolation points from xmin to xmax.

Definition at line 70 of file Distributions.cpp.

{
    if (nbr_samples < 2 || DoubleEqual(xmin, xmax))
        return {mean()};
    std::vector<double> result(nbr_samples);
    for (size_t i = 0; i < nbr_samples; ++i)
        result[i] = xmin + i * (xmax - xmin) / (nbr_samples - 1.0);
    return result;
}

References IDistribution1D::mean().

Referenced by DistributionGate::equidistantPoints(), DistributionLorentz::equidistantPoints(), DistributionGaussian::equidistantPoints(), DistributionLogNormal::equidistantPoints(), DistributionCosine::equidistantPoints(), equidistantPoints(), and IDistribution1D::equidistantSamplesInRange().

Here is the call graph for this function:

equidistantSamples()

std::vector< ParameterSample > IDistribution1D::equidistantSamples ( size_t  nbr_samples, double  sigma_factor = 0., const RealLimits &  limits = RealLimits()  ) const

inherited

Returns equidistant samples, using intrinsic parameters, weighted with probabilityDensity().

Definition at line 43 of file Distributions.cpp.

{
    if (nbr_samples == 0)
        throw std::runtime_error("IDistribution1D::generateSamples: "
                                 "number of generated samples must be bigger than zero");
    if (isDelta())
        return {ParameterSample(mean())};
    return generateSamplesFromValues(equidistantPoints(nbr_samples, sigma_factor, limits));
}

References IDistribution1D::equidistantPoints(), IDistribution1D::generateSamplesFromValues(), IDistribution1D::isDelta(), and IDistribution1D::mean().

Here is the call graph for this function:

equidistantSamplesInRange()

std::vector< ParameterSample > IDistribution1D::equidistantSamplesInRange ( size_t  nbr_samples, double  xmin, double  xmax  ) const

inherited

Returns equidistant samples from xmin to xmax, weighted with probabilityDensity().

Definition at line 58 of file Distributions.cpp.

{
    if (nbr_samples == 0)
        throw std::runtime_error("IDistribution1D::generateSamples: "
                                 "number of generated samples must be bigger than zero");
    if (isDelta())
        return {ParameterSample(mean())};
    return generateSamplesFromValues(equidistantPointsInRange(nbr_samples, xmin, xmax));
}

References IDistribution1D::equidistantPointsInRange(), IDistribution1D::generateSamplesFromValues(), IDistribution1D::isDelta(), and IDistribution1D::mean().

Here is the call graph for this function:

generateSamplesFromValues()

std::vector< ParameterSample > IDistribution1D::generateSamplesFromValues ( const std::vector< double > &  sample_values ) const

protectedinherited

Returns weighted samples from given interpolation points and probabilityDensity().

Definition at line 99 of file Distributions.cpp.

{
    std::vector<ParameterSample> result;
    double norm_factor = 0.0;
    for (double value : sample_values) {
        double pdf = probabilityDensity(value);
        result.emplace_back(value, pdf);
        norm_factor += pdf;
    }
    if (norm_factor <= 0.0)
        throw std::runtime_error("IDistribution1D::generateSamples: "
                                 "total probability must be bigger than zero");
    for (ParameterSample& sample : result)
        sample.weight /= norm_factor;
    return result;
}

References IDistribution1D::probabilityDensity().

Referenced by IDistribution1D::equidistantSamples(), and IDistribution1D::equidistantSamplesInRange().

Here is the call graph for this function:

getLeftWidth()

double DistributionTrapezoid::getLeftWidth ( ) const

inline

Definition at line 326 of file Distributions.h.

{ return m_left; }

References m_left.

getMiddleWidth()

double DistributionTrapezoid::getMiddleWidth ( ) const

inline

Definition at line 327 of file Distributions.h.

{ return m_middle; }

References m_middle.

getRightWidth()

double DistributionTrapezoid::getRightWidth ( ) const

inline

Definition at line 328 of file Distributions.h.

{ return m_right; }

References m_right.

isDelta()

bool DistributionTrapezoid::isDelta ( ) const

overridevirtual

Returns true if the distribution is in the limit case of a Dirac delta distribution.

Implements IDistribution1D.

Definition at line 446 of file Distributions.cpp.

{
    return (m_left + m_middle + m_right) == 0.0;
}

References m_left, m_middle, and m_right.

mean()

double DistributionTrapezoid::mean ( ) const

inlineoverridevirtual

Returns the distribution-specific mean.

Implements IDistribution1D.

Definition at line 325 of file Distributions.h.

{ return m_center; }

References m_center.

nodeChildren()

std::vector< const INode * > INode::nodeChildren ( ) const

virtualinherited

Returns all children.

Reimplemented in ISimulation2D, ISimulation, ParticleCoreShell, ParticleComposition, Particle, MesoCrystal, IParticle, Crystal, MultiLayer, Layer, LayerInterface, ParticleLayout, InterferenceRadialParaCrystal, InterferenceFinite3DLattice, InterferenceFinite2DLattice, Interference3DLattice, Interference2DSuperLattice, Interference2DParaCrystal, Interference2DLattice, Interference1DLattice, ConvolutionDetectorResolution, IDetector, and Beam.

Definition at line 56 of file INode.cpp.

{
    return {};
}

Referenced by NodeUtils::AllDescendantsOfType(), NodeUtils::ChildNodesOfType(), ISampleNode::containedMaterials(), and INode::nodeOffspring().

nodeOffspring()

std::vector< const INode * > INode::nodeOffspring ( ) const

inherited

Returns all descendants.

Definition at line 61 of file INode.cpp.

{
    std::vector<const INode*> result;
    result.push_back(this);
    for (const auto* child : nodeChildren()) {
        for (const auto* p : child->nodeOffspring())
            result.push_back(p);
    }
    return result;
}

References INode::nodeChildren().

Here is the call graph for this function:

parDefs()

std::vector<ParaMeta> DistributionTrapezoid::parDefs ( ) const

inlinefinalvirtual

Returns the parameter definitions, to be hard-coded in each leaf class.

Reimplemented from INode.

Definition at line 316 of file Distributions.h.

    {
        return {{"Center", "", "para_tooltip", -INF, +INF, 0},
                {"LeftWidth", "", "para_tooltip", -INF, +INF, 0},
                {"MiddleWidth", "", "para_tooltip", -INF, +INF, 0},
                {"RightWidth", "", "para_tooltip", -INF, +INF, 0}};
    }

References INF.

probabilityDensity()

double DistributionTrapezoid::probabilityDensity ( double  x ) const

overridevirtual

Returns the distribution-specific probability density for value x.

Implements IDistribution1D.

Definition at line 421 of file Distributions.cpp.

{
    double height = 2.0 / (m_left + 2.0 * m_middle + m_right);
    double min = m_center - m_middle / 2.0 - m_left;
    if (x < min)
        return 0.0;
    if (x < min + m_left)
        return (x - min) * height / m_left;
    if (x < min + m_left + m_middle)
        return height;
    if (x < min + m_left + m_middle + m_right)
        return height - (x - min - m_left - m_middle) * height / m_right;
    return 0.0;
}

References m_center, m_left, m_middle, and m_right.

pythonConstructor()

std::string DistributionTrapezoid::pythonConstructor ( const std::string &  units ) const

overridevirtual

Prints distribution with constructor parameters in given units. ba.DistributionGaussian(2.0*deg, 0.02*deg)

Implements IDistribution1D.

Definition at line 451 of file Distributions.cpp.

{
    return Py::Fmt::printFunction(
        className(), {{m_center, units}, {m_left, units}, {m_middle, units}, {m_right, units}});
}

References className(), m_center, m_left, m_middle, m_right, and Py::Fmt::printFunction().

Here is the call graph for this function:

transferToCPP()

virtual void ICloneable::transferToCPP ( )

inlinevirtualinherited

Used for Python overriding of clone (see swig/tweaks.py)

Definition at line 32 of file ICloneable.h.

Member Data Documentation

m_center

const double& DistributionTrapezoid::m_center

private

Definition at line 342 of file Distributions.h.

Referenced by clone(), equidistantPoints(), mean(), probabilityDensity(), and pythonConstructor().

m_left

const double& DistributionTrapezoid::m_left

private

Definition at line 343 of file Distributions.h.

Referenced by DistributionTrapezoid(), adjustLimitsToNonZeroSamples(), clone(), equidistantPoints(), getLeftWidth(), isDelta(), probabilityDensity(), and pythonConstructor().

m_middle

const double& DistributionTrapezoid::m_middle

private

Definition at line 344 of file Distributions.h.

Referenced by DistributionTrapezoid(), clone(), equidistantPoints(), getMiddleWidth(), isDelta(), probabilityDensity(), and pythonConstructor().

m_P

std::vector<double> INode::m_P

protectedinherited

Definition at line 63 of file INode.h.

Referenced by IFootprintFactor::IFootprintFactor(), INode::checkNodeArgs(), IProfile1D::pythonConstructor(), IProfile2D::pythonConstructor(), IFormFactor::pythonConstructor(), Profile1DVoigt::pythonConstructor(), and Profile2DVoigt::pythonConstructor().

m_right

const double& DistributionTrapezoid::m_right

private

Definition at line 345 of file Distributions.h.

Referenced by DistributionTrapezoid(), adjustLimitsToNonZeroSamples(), clone(), equidistantPoints(), getRightWidth(), isDelta(), probabilityDensity(), and pythonConstructor().

---

The documentation for this class was generated from the following files:

• /home/build/builds/o5h8MZZm/0/mlz/bornagain/Param/Distrib/Distributions.h
• /home/build/builds/o5h8MZZm/0/mlz/bornagain/Param/Distrib/Distributions.cpp