均匀真实分布中的c ++错误？ [复制]

Question

我正在尝试使用std::uniform_real_distribution<float>(a, b) 生成随机浮点数，我发现输出等于上限b 的情况。根据： http://www.cplusplus.com/reference/random/uniform_real_distribution http://en.cppreference.com/w/cpp/numeric/random/uniform_real_distribution 这不应该发生。

对于gcc-4.9.2 和clang-3.5.0，我会出现以下崩溃：

#include <iomanip>
#include <iostream>
#include <limits>
#include <random>

int main() {
    float a = 1.0f;
    float b = 1.001f;
    size_t seed = 293846;
    size_t n = 9830;
    std::mt19937 rg(seed);
    std::uniform_real_distribution< float > u(a, b);
    for (size_t i = 0; i < n; ++i) {
        float v = u(rg);
        if (not (v < b)) {
            std::cerr << "error: i=" << i
                      << " v=" << std::scientific << std::setprecision(std::numeric_limits< double >::max_digits10) << v
                      << " b=" << b << std::endl;
            abort();
        }
    }
}

我看到的输出是：

error: i=9829 v=1.00100004673004150e+00 b=1.00100004673004150e+00
Aborted (core dumped)

这是标准库的错误吗？我错过了什么吗？

编辑。虽然不是完全重复，但此问题中提出的问题确实已包含在此处引用的其他问题的答案中。

有关此问题状态的最详细信息是此处的 LLVM 错误线程：https://llvm.org/bugs/show_bug.cgi?id=18767

Answer 1

虽然[rand.dist.uni.real] 要求在实数的半开区间上均匀分布，但它不包含任何关于如何将这些四舍五入为浮点表示的语言。

所以从我能找到的信息来看，一个符合要求的实现产生比较等于上限的浮点数很可能是合法的（并且要求）。