BornAgain  1.19.79
Open-source research software to simulate and fit neutron and x-ray reflectometry and grazing-incidence small-angle scattering
DistributionGate Class Reference

Description

Uniform distribution function with half width hwhm.

Definition at line 90 of file Distributions.h.

Inheritance diagram for DistributionGate:
[legend]
Collaboration diagram for DistributionGate:
[legend]

Public Member Functions

 DistributionGate ()
 
 DistributionGate (double min, double max)
 
 DistributionGate (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...
 
DistributionGateclone () const override
 
std::vector< double > equidistantPoints (size_t nbr_samples, double sigma_factor, const RealLimits &limits=RealLimits()) const override
 Returns 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< ParameterSampleequidistantSamples (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< ParameterSampleequidistantSamplesInRange (size_t nbr_samples, double xmin, double xmax) const
 Returns equidistant samples from xmin to xmax, weighted with probabilityDensity(). More...
 
bool isDelta () const override
 Returns true if the distribution is in the limit case of a Dirac delta distribution. More...
 
double max () const
 
double mean () const override
 Returns the distribution-specific mean. More...
 
double min () const
 
virtual std::vector< const INode * > nodeChildren () const
 Returns all children. More...
 
std::vector< const INode * > nodeOffspring () const
 Returns all descendants. More...
 
std::vector< ParaMetaparDefs () 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< ParameterSamplegenerateSamplesFromValues (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_max
 
const double & m_min
 

Constructor & Destructor Documentation

◆ DistributionGate() [1/3]

DistributionGate::DistributionGate ( std::vector< double >  P)

Definition at line 120 of file Distributions.cpp.

121  : IDistribution1D(P)
122  , m_min(m_P[0])
123  , m_max(m_P[1])
124 {
125  checkNodeArgs();
126  if (m_max < m_min)
127  throw std::runtime_error("DistributionGate: max<min");
128 }
const double & m_max
const double & m_min
IDistribution1D(const std::vector< double > &PValues)
void checkNodeArgs() const
Raises exception if a parameter value is invalid.
Definition: INode.cpp:27
std::vector< double > m_P
Definition: INode.h:63

References INode::checkNodeArgs(), m_max, and m_min.

Here is the call graph for this function:

◆ DistributionGate() [2/3]

DistributionGate::DistributionGate ( double  min,
double  max 
)

Definition at line 130 of file Distributions.cpp.

131  : DistributionGate(std::vector<double>{min, max})
132 {
133 }
double max() const
double min() const

References max(), and min().

Here is the call graph for this function:

◆ DistributionGate() [3/3]

DistributionGate::DistributionGate ( )

Definition at line 135 of file Distributions.cpp.

136  : DistributionGate(0., 1.)
137 {
138 }

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.

83 {
84  if (limits.hasLowerLimit() && xmin < limits.lowerLimit())
85  xmin = limits.lowerLimit();
86  if (limits.hasUpperLimit() && xmax > limits.upperLimit())
87  xmax = limits.upperLimit();
88  if (xmin > xmax) {
89  std::ostringstream ostr;
90  ostr << "IDistribution1D::adjustMinMaxForLimits() -> Error. Can't' adjust ";
91  ostr << "xmin:" << xmin << " xmax:" << xmax << " for given limits " << limits << std::endl;
92  throw std::runtime_error(ostr.str());
93  }
94 }
bool hasUpperLimit() const
if has upper limit
Definition: RealLimits.cpp:66
double upperLimit() const
Returns upper limit.
Definition: RealLimits.cpp:71
double lowerLimit() const
Returns lower limit.
Definition: RealLimits.cpp:49
bool hasLowerLimit() const
if has lower limit
Definition: RealLimits.cpp:44

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

Referenced by equidistantPoints(), DistributionLorentz::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.

28 {
29  size_t nP = m_P.size();
30  if (parDefs().size() != nP) {
31  std::cerr << "BUG in class " << className() << std::endl;
32  std::cerr << "#m_P = " << nP << std::endl;
33  std::cerr << "#PDf = " << parDefs().size() << std::endl;
34  for (const ParaMeta& pm : parDefs())
35  std::cerr << " PDf: " << pm.name << std::endl;
36  ASSERT(0);
37  }
38  ASSERT(parDefs().size() == nP);
39  for (size_t i = 0; i < nP; ++i) {
40  const ParaMeta pm = parDefs()[i];
41 
43  if (pm.vMin == -INF) {
44  ASSERT(pm.vMax == +INF);
45  // nothing to do
46  } else if (pm.vMax == +INF) {
47  ASSERT(pm.vMin == 0);
48  limits = RealLimits::nonnegative();
49  } else {
50  limits = RealLimits::limited(pm.vMin, pm.vMax);
51  }
52  limits.check(pm.name, m_P[i]);
53  }
54 }
#define ASSERT(condition)
Definition: Assert.h:45
const double INF
Definition: INode.h:26
virtual std::vector< ParaMeta > parDefs() const
Returns the parameter definitions, to be hard-coded in each leaf class.
Definition: INode.h:51
virtual std::string className() const =0
Returns the class name, to be hard-coded in each leaf class that inherits from INode.
Limits for a real fit parameter.
Definition: RealLimits.h:24
static RealLimits limitless()
Creates an object without bounds (default)
Definition: RealLimits.cpp:139
void check(const std::string &name, double value) const
Throws if value is outside limits. Parameter 'name' is for exception message.
Definition: RealLimits.cpp:170
static RealLimits nonnegative()
Creates an object which can have only positive values with 0. included.
Definition: RealLimits.cpp:124
static RealLimits limited(double left_bound_value, double right_bound_value)
Creates an object bounded from the left and right.
Definition: RealLimits.cpp:134
Metadata of one model parameter.
Definition: INode.h:29
double vMin
Definition: INode.h:33
double vMax
Definition: INode.h:34
std::string name
Definition: INode.h:30

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(), DistributionGaussian::DistributionGaussian(), DistributionLogNormal::DistributionLogNormal(), DistributionLorentz::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 DistributionGate::className ( ) const
inlinefinalvirtual

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

Implements INode.

Definition at line 97 of file Distributions.h.

97 { return "DistributionGate"; }

Referenced by pythonConstructor().

◆ clone()

DistributionGate* DistributionGate::clone ( ) const
inlineoverridevirtual

Implements IDistribution1D.

Definition at line 96 of file Distributions.h.

96 { return new DistributionGate(m_min, m_max); }

References DistributionGate(), m_max, and m_min.

Here is the call graph for this function:

◆ equidistantPoints()

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

Returns list of sample values.

Implements IDistribution1D.

Definition at line 149 of file Distributions.cpp.

151 {
152  double xmin = m_min;
153  double xmax = m_max;
154  adjustMinMaxForLimits(xmin, xmax, limits);
155  return equidistantPointsInRange(nbr_samples, xmin, xmax);
156 }
void adjustMinMaxForLimits(double &xmin, double &xmax, const RealLimits &limits) const
modifies xmin and xmax if they are outside of limits
virtual std::vector< double > equidistantPointsInRange(size_t nbr_samples, double xmin, double xmax) const
Returns equidistant interpolation points from xmin to xmax.

References IDistribution1D::adjustMinMaxForLimits(), IDistribution1D::equidistantPointsInRange(), m_max, and m_min.

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.

72 {
73  if (nbr_samples < 2 || DoubleEqual(xmin, xmax))
74  return {mean()};
75  std::vector<double> result(nbr_samples);
76  for (size_t i = 0; i < nbr_samples; ++i)
77  result[i] = xmin + i * (xmax - xmin) / (nbr_samples - 1.0);
78  return result;
79 }
virtual double mean() const =0
Returns the distribution-specific mean.

References IDistribution1D::mean().

Referenced by equidistantPoints(), DistributionLorentz::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.

46 {
47  if (nbr_samples == 0)
48  throw std::runtime_error("IDistribution1D::generateSamples: "
49  "number of generated samples must be bigger than zero");
50  if (isDelta())
51  return {ParameterSample(mean())};
52  return generateSamplesFromValues(equidistantPoints(nbr_samples, sigma_factor, limits));
53 }
std::vector< ParameterSample > generateSamplesFromValues(const std::vector< double > &sample_values) const
Returns weighted samples from given interpolation points and probabilityDensity().
virtual bool isDelta() const =0
Returns true if the distribution is in the limit case of a Dirac delta distribution.
virtual std::vector< double > equidistantPoints(size_t nbr_samples, double sigma_factor, const RealLimits &limits=RealLimits()) const =0
Returns equidistant interpolation points, with range computed in distribution-specific way from mean ...
A parameter value with a weight, as obtained when sampling from a distribution.

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.

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

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.

100 {
101  std::vector<ParameterSample> result;
102  double norm_factor = 0.0;
103  for (double value : sample_values) {
104  double pdf = probabilityDensity(value);
105  result.emplace_back(value, pdf);
106  norm_factor += pdf;
107  }
108  if (norm_factor <= 0.0)
109  throw std::runtime_error("IDistribution1D::generateSamples: "
110  "total probability must be bigger than zero");
111  for (ParameterSample& sample : result)
112  sample.weight /= norm_factor;
113  return result;
114 }
virtual double probabilityDensity(double x) const =0
Returns the distribution-specific probability density for value x.

References IDistribution1D::probabilityDensity().

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

Here is the call graph for this function:

◆ isDelta()

bool DistributionGate::isDelta ( ) const
overridevirtual

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

Implements IDistribution1D.

Definition at line 158 of file Distributions.cpp.

159 {
160  return DoubleEqual(m_min, m_max);
161 }

References m_max, and m_min.

◆ max()

double DistributionGate::max ( ) const
inline

Definition at line 108 of file Distributions.h.

108 { return m_max; }

References m_max.

Referenced by DistributionGate().

◆ mean()

double DistributionGate::mean ( ) const
inlineoverridevirtual

Returns the distribution-specific mean.

Implements IDistribution1D.

Definition at line 106 of file Distributions.h.

106 { return (m_min + m_max) / 2.0; }

References m_max, and m_min.

◆ min()

double DistributionGate::min ( ) const
inline

Definition at line 107 of file Distributions.h.

107 { return m_min; }

References m_min.

Referenced by DistributionGate().

◆ nodeChildren()

◆ nodeOffspring()

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

Returns all descendants.

Definition at line 61 of file INode.cpp.

62 {
63  std::vector<const INode*> result;
64  result.push_back(this);
65  for (const auto* child : nodeChildren()) {
66  for (const auto* p : child->nodeOffspring())
67  result.push_back(p);
68  }
69  return result;
70 }
virtual std::vector< const INode * > nodeChildren() const
Returns all children.
Definition: INode.cpp:56

References INode::nodeChildren().

Here is the call graph for this function:

◆ parDefs()

std::vector<ParaMeta> DistributionGate::parDefs ( ) const
inlinefinalvirtual

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

Reimplemented from INode.

Definition at line 99 of file Distributions.h.

100  {
101  return {{"Min", "", "para_tooltip", -INF, +INF, 0},
102  {"Max", "", "para_tooltip", -INF, +INF, 0}};
103  }

References INF.

◆ probabilityDensity()

double DistributionGate::probabilityDensity ( double  x) const
overridevirtual

Returns the distribution-specific probability density for value x.

Implements IDistribution1D.

Definition at line 140 of file Distributions.cpp.

141 {
142  if (x < m_min || x > m_max)
143  return 0.0;
144  if (DoubleEqual(m_min, m_max))
145  return 1.0;
146  return 1.0 / (m_max - m_min);
147 }

References m_max, and m_min.

◆ pythonConstructor()

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

164 {
165  return Py::Fmt::printFunction(className(), m_min, units, m_max, units);
166 }
std::string className() const final
Returns the class name, to be hard-coded in each leaf class that inherits from INode.
Definition: Distributions.h:97
std::string printFunction(const std::string &name, const std::vector< std::pair< double, std::string >> &arguments)
Print a function in the form "<name>(<arguments>)". arguments will be processed by printArguments(),...
Definition: PyFmt.cpp:168

References className(), m_max, m_min, 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_max

const double& DistributionGate::m_max
private

◆ m_min

const double& DistributionGate::m_min
private

◆ m_P


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