`scipy.stats._qmvnt:_qmvn`

Signature

def _qmvn(m, covar, low, high, rng, lattice = cbc, n_batches = 10)

Multivariate normal integration over box bounds.

m : int > n_batches: The number of points to sample. This number will be divided into n_batches batches that apply random offsets of the sampling lattice for each batch in order to estimate the error.
covar : (n, n) float array: Possibly singular, positive semidefinite symmetric covariance matrix.
low, high : (n,) float array: The low and high integration bounds.
rng : Generator, optional: default_rng(), yada, yada
lattice : 'cbc' or callable: The type of lattice rule to use to construct the integration points.
n_batches : int > 0, optional: The number of QMC batches to apply.