numpy random.normal 使用什么算法？它只是正态分布的高斯公式吗？它使用梅森捻线机吗？

Question

我是初学者，数学算法不是我的菜，但我很好奇 numpy random.normal 是否使用 Mersenne twister，因为那样它就不是正态分布？

Answer 1

Numpy 使用 Ziggurat 方法。至少 1.19 版是这样。你可以在这里查看源代码：https://github.com/numpy/numpy/blob/91d9bbebc159ada6cccd0e1fcec1926e26b19c6f/numpy/random/src/distributions/distributions.c#L137-L170。