`scipy.special._spfun_stats:multigammaln`

source: build-install/usr/lib/python3.14/site-packages/scipy/special/_spfun_stats.py:45

Signature

def multigammaln(a, d)

Summary

Returns the log of multivariate gamma, also sometimes called the generalized gamma.

Parameters

a : ndarray: The multivariate gamma is computed for each item of a.
d : int: The dimension of the space of integration.

Returns

res : ndarray: The values of the log multivariate gamma at the given points a.

Notes

The formal definition of the multivariate gamma of dimension d for a real a is

Γ_{d} (a) = \int_{A > 0} e^{- t r (A)} ∣ A ∣^{a - (d + 1) /2} d A

with the condition $a > (d - 1) /2$ , and $A > 0$ being the set of all the positive definite matrices of dimension d. Note that a is a scalar: the integrand only is multivariate, the argument is not (the function is defined over a subset of the real set).

This can be proven to be equal to the much friendlier equation

Γ_{d} (a) = π^{d (d - 1) /4} i = 1 \prod d Γ (a - (i - 1) /2) .

Examples

import numpy as np
from scipy.special import multigammaln, gammaln
a = 23.5
d = 10

✓

multigammaln(a, d)

✗

np.float64(457.1766747829046)

Verify that the result agrees with the logarithm of the equation shown above:

d*(d-1)/4*np.log(np.pi) + gammaln(a - 0.5*np.arange(0, d)).sum()

✗

np.float64(457.1766747829046)

Aliases

scipy.special.multigammaln