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

Question

在本文中，我使用scipy.stats模块中的函数mvn.mvnun()来计算给定多元正态分布的CDF。评审员问我如何估计CDF，因为CDF没有封闭表单。我想可以用抽样的方法。为了了解它是如何工作的，我搜索了scipy源代码库。但是，在scipy.statsfolder中，我没有看到mvn.py，而是看到了{a2}。在

在文件中，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函数的接口