如何找到 scipy.stats.mvn 的源代码？

Question

我在论文中使用 scipy.stats 模块中的函数 mvn.mvnun() 来计算给定多元正态分布的 CDF。审稿人问我如何估计 CDF，因为 CDF 没有封闭形式。我想可能会使用抽样方法。为了了解它是如何工作的，我搜索了 scipy 源代码库。但是，在 scipy.stats 中取而代之的是 folder, I didn't see mvn.py, I saw mvn.pyf。

在文件中，mvn.mvnun好像是这样定义的，

!    -*- f90 -*-
! Note: the context of this file is case sensitive.

python module mvn ! in 
    interface  ! in :mvn
        subroutine mvnun(d,n,lower,upper,means,covar,maxpts,abseps,releps,value,inform) ! in :mvn:mvndst.f
            integer intent(hide) :: d=shape(means,0)
            integer intent(hide) :: n=shape(means,1)
            double precision dimension(d) :: lower
            double precision dimension(d) :: upper
            double precision dimension(d,n) :: means
            double precision dimension(d,d) :: covar
            integer intent(optional) :: maxpts=d*1000
            double precision intent(optional) :: abseps=1e-6
            double precision intent(optional) :: releps=1e-6
            double precision intent(out) :: value
            integer intent(out) :: inform
        end subroutine mvnun

但是，没有那个函数的详细定义。我想知道在哪里可以找到 mvnun 的源代码，告诉我它是如何计算 CDF 的？

Answer 1

实现是用 Fortran 语言编写的。 pyf 文件定义了 mvndst.f

中 Fortran 函数的接口

如何找到 scipy.stats.mvn 的源代码？

How can I find the source code for scipy.stats.mvn?

python

statistics

scipy

probability-theory