3

Numpy 文档建议使用Random Generator 包中的所有np.random.permutation新代码。np.random.default_rng()我在文档中看到 Random Generator 包已经标准化了围绕 BitGenerator 生成各种随机分布，而不是使用我隐约熟悉的 Mersenne Twister。

我看到了一个缺点，过去只需要一行代码就可以进行简单的排列：

np.random.permutation(10)

现在变成了两行代码，对于这么简单的任务感觉有点别扭：

rng = np.random.default_rng()
rng.permutation(10)
  • 为什么这种新方法比以前的方法有所改进？
  • 为什么现有的方法不能像np.random.permutation包装这种新的首选方法一样呢？
  • np.random.default_rng().permutation(10)假设它没有被大量调用，是否有充分的理由不将这种新方法用作单线？
  • 是否有将现有代码切换到此方法的论据？
4

1 回答 1

2

一些上下文：

对于您的问题，按逻辑顺序：

为什么像 np.random.permutation 这样的现有方法不包装这个新的首选方法？

可能是因为向后兼容性问题。即使“顶级”API 不会改变，它的内部结构也足以被视为兼容性中断。

为什么这种新方法比以前的方法有所改进？

“默认情况下，Generator 使用 PCG64 提供的位，它比 RandomState 中使用的传统 MT19937 具有更好的统计特性。” （来源）。PCG64文档字符串提供了更多技术细节。

np.random.default_rng().permutation(10)假设它没有被大量调用，是否有充分的理由不将这种新方法用作单线？

我非常同意如果在模块开始时添加的代码行有点尴尬。我只想指出 NumPy 文档在文档字符串示例中直接使用这种形式，例如：

n = np.random.default_rng().standard_exponential((3, 8000))

细微的差别是在模块加载/导入时实例化一个类，而在您的表单中它可能会稍后出现。但这应该是一个微小的差异（再次假设它只使用一次或几次）。如果您查看default_rng(seed)源代码，当使用 调用时None，它会Generator(PCG64(seed))在对 进行几次快速检查后返回seed。

是否有将现有代码切换到此方法的论据？

因为我没有任何接近深度的技术知识来对算法进行很好的比较，并且还因为它取决于其他一些变量，例如您是否担心使下游代码兼容使用旧版本的 NumPy，default_rng()根本不存在。

于 2020-06-17T21:50:26.963 回答