c++ - 将a==b作为参数传递给uniform_real_distribution ok

Question

可以将 std::uniform_real_distribution 与 a==b 一起使用吗？这种情况可以解释为没有随机性。

在 GNU C++ 库中:

explicit
param_type(_RealType __a = _RealType(0),
           _RealType __b = _RealType(1))
: _M_a(__a), _M_b(__b)
{
   _GLIBCXX_DEBUG_ASSERT(_M_a <= _M_b);
}

断言接受a==b，尽管它可能需要除以零。这是吗

• 标准的扩展？
• 一个无论如何都行得通的巧妙设计？
• 在某些 CPU 上发生错误？

Answer 1

相关规范段落[26.5.8.2.2, rand.dist.uni.real]/1说:

A uniform_real_distribution random number distribution produces random numbers x, a ≤ x < b, distributed according to the constant probability density function

p(x | a, b) = 1/(b − a) .

[Note: This implies that p(x | a, b) is undefined when a == b. — end note]

请注意，a <b 是对分布的约束。

有趣的是，构造函数(第 2 段)的要求更宽松:

explicit uniform_real_distribution(RealType a = 0.0, RealType b = 1.0);

Requires: a ≤ b [...]

这意味着您可以构造一个具有a == b的分布对象，但它描述的分布是未定义的。也就是说，您可能不会尝试使用这些参数通过此类分布对象生成值(这实际上可能会导致除以零或其他任何值)。

另请参阅N3926为什么这两个约束并不不一致。该问题最初是作为 LWG 2168 提出的。并已通过添加第一段中引用的注释得到解决。该论文提请注意采用不同分布参数的 operator() 重载；结果是约束仅在生成值时发挥作用，而不是在构造分布对象时发挥作用。