DistributionLorentz Class Reference

Description

Lorentz distribution with half width hwhm.

Definition at line 132 of file Distributions.h.

Inheritance diagram for DistributionLorentz:

Collaboration diagram for DistributionLorentz:

Public Member Functions
	DistributionLorentz ()
	DistributionLorentz (double mean, double hwhm)
	DistributionLorentz (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...
DistributionLorentz *	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	hwhm () 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 Attributes
const double &	m_hwhm
const double &	m_mean

Constructor & Destructor Documentation

DistributionLorentz() [1/3]

DistributionLorentz::DistributionLorentz

(

std::vector< double >

P

)

Definition at line 172 of file Distributions.cpp.

    : IDistribution1D(P)
    , m_mean(m_P[0])
    , m_hwhm(m_P[1])
{
    checkNodeArgs();
    if (m_hwhm < 0.0)
        throw std::runtime_error("DistributionLorentz: hwhm<0");
}

References INode::checkNodeArgs(), and m_hwhm.

Here is the call graph for this function:

DistributionLorentz() [2/3]

DistributionLorentz::DistributionLorentz	(	double	mean,
		double	hwhm
	)

Definition at line 182 of file Distributions.cpp.

    : DistributionLorentz(std::vector<double>{mean, hwhm})
{
}

References hwhm(), and mean().

Here is the call graph for this function:

DistributionLorentz() [3/3]

DistributionLorentz::DistributionLorentz

(

)

Definition at line 187 of file Distributions.cpp.

    : DistributionLorentz(0., 1.)
{
}

Referenced by clone().

Member Function Documentation

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(), equidistantPoints(), DistributionGaussian::equidistantPoints(), DistributionLogNormal::equidistantPoints(), DistributionCosine::equidistantPoints(), and DistributionTrapezoid::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(), DistributionTrapezoid::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 DistributionLorentz::className ( ) const

inlinefinalvirtual

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

Implements INode.

Definition at line 139 of file Distributions.h.

{ return "DistributionLorentz"; }

Referenced by pythonConstructor().

clone()

DistributionLorentz* DistributionLorentz::clone ( ) const

inlineoverridevirtual

Implements IDistribution1D.

Definition at line 138 of file Distributions.h.

{ return new DistributionLorentz(m_mean, m_hwhm); }

References DistributionLorentz(), m_hwhm, and m_mean.

Here is the call graph for this function:

equidistantPoints()

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

overridevirtual

generate list of sample values

Implements IDistribution1D.

Definition at line 199 of file Distributions.cpp.

{
    if (sigma_factor <= 0.0)
        sigma_factor = 2.0;
    double xmin = m_mean - sigma_factor * m_hwhm;
    double xmax = m_mean + sigma_factor * m_hwhm;
    adjustMinMaxForLimits(xmin, xmax, limits);
    return equidistantPointsInRange(nbr_samples, xmin, xmax);
}

References IDistribution1D::adjustMinMaxForLimits(), IDistribution1D::equidistantPointsInRange(), m_hwhm, and m_mean.

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(), equidistantPoints(), DistributionGaussian::equidistantPoints(), DistributionLogNormal::equidistantPoints(), DistributionCosine::equidistantPoints(), DistributionTrapezoid::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:

hwhm()

double DistributionLorentz::hwhm ( ) const

inline

Definition at line 149 of file Distributions.h.

{ return m_hwhm; }

References m_hwhm.

Referenced by DistributionLorentz().

isDelta()

bool DistributionLorentz::isDelta ( ) const

overridevirtual

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

Implements IDistribution1D.

Definition at line 210 of file Distributions.cpp.

{
    return m_hwhm == 0.0;
}

References m_hwhm.

mean()

double DistributionLorentz::mean ( ) const

inlineoverridevirtual

Returns the distribution-specific mean.

Implements IDistribution1D.

Definition at line 148 of file Distributions.h.

{ return m_mean; }

References m_mean.

Referenced by DistributionLorentz().

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> DistributionLorentz::parDefs ( ) const

inlinefinalvirtual

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

Reimplemented from INode.

Definition at line 141 of file Distributions.h.

    {
        return {{"Mean", "", "para_tooltip", -INF, +INF, 0},
                {"HWHM", "", "para_tooltip", -INF, +INF, 0}};
    }

References INF.

probabilityDensity()

double DistributionLorentz::probabilityDensity ( double  x ) const

overridevirtual

Returns the distribution-specific probability density for value x.

Implements IDistribution1D.

Definition at line 192 of file Distributions.cpp.

{
    if (m_hwhm == 0.0)
        return DoubleEqual(x, m_mean) ? 1.0 : 0.0;
    return m_hwhm / (m_hwhm * m_hwhm + (x - m_mean) * (x - m_mean)) / M_PI;
}

References m_hwhm, m_mean, and M_PI.

pythonConstructor()

std::string DistributionLorentz::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 215 of file Distributions.cpp.

{
    return Py::Fmt::printFunction(className(), m_mean, units, m_hwhm, units);
}

References className(), m_hwhm, m_mean, 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_hwhm

const double& DistributionLorentz::m_hwhm

private

Definition at line 163 of file Distributions.h.

Referenced by DistributionLorentz(), clone(), equidistantPoints(), hwhm(), isDelta(), probabilityDensity(), and pythonConstructor().

m_mean

const double& DistributionLorentz::m_mean

private

Definition at line 162 of file Distributions.h.

Referenced by clone(), equidistantPoints(), mean(), 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().

---

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