uki ddlZddlmZddlmZddlmZddlm Z ddl m Z m Z ddl mZmZded efd Zded ed efd Zded ed efd Zded ed efdZy)N)lax)numpy)_const)promote_dtypes_inexact)gammalnxlogy)Array ArrayLikexreturnctj|d}tj|dkDdt|dz dkzS)Nraxisgư>)jnpsumallabs)r x_sums Y/home/cdr/jupyterlab/.venv/lib/python3.12/site-packages/jax/_src/scipy/stats/dirichlet.py _is_simplexrs9 ''!! % QQ 3uqy>D#8 99alphac&tt||S)aVDirichlet log probability distribution function. JAX implementation of :obj:`scipy.stats.dirichlet` ``logpdf``. The Dirichlet probability density function is .. math:: f(\mathbf{x}) = \frac{1}{B(\mathbf{\alpha})} \prod_{i=1}^K x_i^{\alpha_i - 1} where :math:`B(\mathbf{\alpha})` is the :func:`~jax.scipy.special.beta` function in a :math:`K`-dimensional vector space. Args: x: arraylike, value at which to evaluate the PDF alpha: arraylike, distribution shape parameter Returns: array of logpdf values. See Also: :func:`jax.scipy.stats.dirichlet.pdf` )_logpdfrr rs rlogpdfrs0 (E2 33rc |jdk7rtd|j|jd|jd|jddz fvr%td|jd|jt|d}|jd|jdk7r>t j |t j||jddgd}t jt|tt j|z }|jdkDr4t j||jd |jdz zzd }t jt jtt j|||d|}t jt||tj S) Nrz1`alpha` must be one-dimensional; got alpha.shape=rz[`x` must have either the same number of entries as `alpha` or one entry fewer; got x.shape=z, alpha.shape=T)keepdimsr)r)r)ndim ValueErrorshape _lax_constr concatenatersubrrbroadcast_in_dimrwherernpinf)r ronenormalize_term log_probss rrr8su ZZ1_ 9%++G WWQZ A A(:;; ))*  }N  1a#WWQZ5;;q>! CGGCq4)@ABKA7775>*WSWWU^-DD.VVaZ   daffqj6I(I4 PEggcggeCGGE3$7;!DnU) ;q>9rvvg 66rc@tjt||S)aODirichlet probability distribution function. JAX implementation of :obj:`scipy.stats.dirichlet` ``pdf``. The Dirichlet probability density function is .. math:: f(\mathbf{x}) = \frac{1}{B(\mathbf{\alpha})} \prod_{i=1}^K x_i^{\alpha_i - 1} where :math:`B(\mathbf{\alpha})` is the :func:`~jax.scipy.special.beta` function in a :math:`K`-dimensional vector space. Args: x: arraylike, value at which to evaluate the PDF alpha: arraylike, distribution shape parameter Returns: array of pdf values. See Also: :func:`jax.scipy.stats.dirichlet.logpdf` )rexprrs rpdfr/Ls0 5! ""r)rr(jax._srcrrjax._src.lax.laxrr#jax._src.numpy.utilrjax._src.scipy.specialrrjax._src.typingr r rrrr/rrr6s}!161,:5:U: 4i4 4e447u7U7u7(#9#Y#5#r