如何查看numpy.random.exponential的源代码？

Question

我想看看numpy.random.exponential是否用F^{-1} (U)方法实现，其中F是指数分布的c.d.f，U是均匀分布。

我尝试了numpy.source(random.exponential)，但是返回了“这个对象不可用”。这是否意味着这个函数不是用Python编写的？

我也尝试过inspect.getsource(random.exponential)，但是返回了一个错误，说它不是模块、函数等等。

Answer 1

numpy的源代码是在github，所以您可以使用github的源搜索。

通常，库的这些部分不是用纯python实现的。

python-部件(关于您的问题)是这里。

更相关的代码部分来自distributions.c。

double rk_standard_exponential(rk_state *state)
{
    /* We use -log(1-U) since U is [0, 1) */
    return -log(1.0 - rk_double(state));
}

double rk_exponential(rk_state *state, double scale)
{
    return scale * rk_standard_exponential(state);
}

Answer 2

大量的numpy函数是用C/C++和Fortran编写的。numpy.source()只返回用Python编写的对象的源代码。它是写在NumPy网站上的。

您可以在其NumPy页面上找到所有的GitHub函数。您需要的一个是用C编写的。下面是指向文件的链接。