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

Question

我mvn.mvnun()在scipy.stats模块中使用该函数来计算论文中给定多元正态分布的CDF。审阅者问我，由于CDF没有封闭格式，我如何估算CDF。我猜可能会使用抽样方法。为了了解它是如何工作的，我搜索了scipy源代码仓库。但是，在scipy.stats 文件夹中，我没有看到mvn.py，而是看到了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函数的接口