o
    Zh{                     @   sx  U d dl Z d dlZd dlmZ d dlmZ d dlZd dlmZmZ ddl	m
Z
 ddlmZ ddlmZ dd	lmZ dd
lmZ ddlmZ ddlmZ ddlmZ ddlmZ ddlmZ ddlmZ ddlm Z  ddl!m"Z" ddl#m$Z$ ddl%m&Z& ddl'm(Z( ddl)m*Z*m+Z+m,Z, ddl-m.Z.m/Z/ ddl0m1Z1 ddl2m3Z3 ddl4m5Z5 ddl6m7Z7 ddl8m9Z9 ddl:m;Z; ddl<m=Z=m>Z? i Z@eAeBeCeCf ef eDd< i ZEeAeBeCeCf ef eDd < d!d"gZFd#d! ZGeG d$d% d%ZHd&d' ZId(d) ZJd*d+ ZKd,d- ZLd.ed/ed0efd1d"ZMeGe
e
d2d3 ZNeGeed4d5 ZOeGeed6d7 ZPeGeed8d9 ZQeGeed:d; ZReGeed<d= ZSeGeed>d? ZTeGeed@dA ZUeGeedBdC ZVeGe"e"dDdE ZWeGe e dFdG ZXeGe$e$dHdI ZYeGe(e(dJdK ZZeGe,e,dLdM Z[eGe/e,dNdO Z\eGe,e/dPdQ Z]eGe/e/dRdS Z^eGe1e1dTdU Z_eGe3e3dVdW Z`eGe5e5dXdY ZaeGe7e7dZd[ ZbeGe9e9d\d] ZceGe;e;d^d_ ZdeGe
e7d`da ZeeGeedbdc ZfeGee5ddde ZgeGeedfdg ZheGeedhdi ZieGee1djdk ZjeGee;dldm ZkeGee5dndo ZleGeedpdq ZmeGee1drds ZneGee;dtdu ZoeGeeeGeeeGee5eGee;dvdw ZpeGeedxdy ZqeGee"dzd{ ZreGee1d|d} ZseGeeeGeeeGee5eGee;d~d ZteGeedd ZueGee"dd ZveGee1dd ZweGe"eeGe"eeGe"eeGe"eeGe"e5eGe"e;dd ZxeGe"e1dd ZyeGe(eeGe(eeGe(eeGe(eeGe(e5eGe(e;dd ZzeGe(e1dd Z{eGe1eeGe1eeGe1eeGe1eeGe1e5eGe1e;dd Z|eGe1e"dd Z}eGe1e(dd Z~eGe5eeGe5eeGe5e;dd ZeGe5edd ZeGe5edd ZeGe5e1dd ZeGe7e
eGe7edd ZeGe;edd ZeGe;edd ZeGe;edd ZeGe;edd ZeGe;e"dd ZeGe;e1dd ZeGe;e5dd ZeGe&e&dd ZeGeedd Zdd ZdS )    N)total_ordering)Callable)infTensor   )	Bernoulli)Beta)Binomial)Categorical)Cauchy)ContinuousBernoulli)	Dirichlet)Distribution)ExponentialFamily)Exponential)Gamma)	Geometric)Gumbel)
HalfNormal)Independent)Laplace)_batch_lowrank_logdet_batch_lowrank_mahalanobisLowRankMultivariateNormal)_batch_mahalanobisMultivariateNormal)Normal)OneHotCategorical)Pareto)Poisson)TransformedDistribution)Uniform)_sum_rightmosteuler_constant_KL_REGISTRY_KL_MEMOIZEregister_klkl_divergencec                    sV   t  tst trtd  t ts"ttr"td  fdd}|S )a[  
    Decorator to register a pairwise function with :meth:`kl_divergence`.
    Usage::

        @register_kl(Normal, Normal)
        def kl_normal_normal(p, q):
            # insert implementation here

    Lookup returns the most specific (type,type) match ordered by subclass. If
    the match is ambiguous, a `RuntimeWarning` is raised. For example to
    resolve the ambiguous situation::

        @register_kl(BaseP, DerivedQ)
        def kl_version1(p, q): ...
        @register_kl(DerivedP, BaseQ)
        def kl_version2(p, q): ...

    you should register a third most-specific implementation, e.g.::

        register_kl(DerivedP, DerivedQ)(kl_version1)  # Break the tie.

    Args:
        type_p (type): A subclass of :class:`~torch.distributions.Distribution`.
        type_q (type): A subclass of :class:`~torch.distributions.Distribution`.
    z6Expected type_p to be a Distribution subclass but got z6Expected type_q to be a Distribution subclass but got c                    s   | t  f< t  | S N)r$   r%   clear)funtype_ptype_q E/var/www/auris/lib/python3.10/site-packages/torch/distributions/kl.py	decoratorV   s   zregister_kl.<locals>.decorator)
isinstancetype
issubclassr   	TypeError)r,   r-   r0   r.   r+   r/   r&   3   s   c                   @   s*   e Zd ZdgZdd Zdd Zdd ZdS )	_Matchtypesc                 G   s
   || _ d S r(   r6   )selfr6   r.   r.   r/   __init__b   s   
z_Match.__init__c                 C   s   | j |j kS r(   r7   )r8   otherr.   r.   r/   __eq__e   s   z_Match.__eq__c                 C   s:   t | j|jD ]\}}t||s dS ||ur dS qdS )NFT)zipr6   r3   )r8   r:   xyr.   r.   r/   __le__h   s   
z_Match.__le__N)__name__
__module____qualname__	__slots__r9   r;   r?   r.   r.   r.   r/   r5   ^   s
    r5   c           	         s    fddt D }|stS tdd |D j\}}tdd |D j\}}t ||f }t ||f }||urMtd j dj d|j d|j d		t |S )
zP
    Find the most specific approximate match, assuming single inheritance.
    c                    s,   g | ]\}}t  |rt |r||fqS r.   )r3   ).0Zsuper_pZsuper_qr+   r.   r/   
<listcomp>u   s    z _dispatch_kl.<locals>.<listcomp>c                 s   s    | ]}t | V  qd S r(   )r5   rD   mr.   r.   r/   	<genexpr>   s    z_dispatch_kl.<locals>.<genexpr>c                 s   s    | ]	}t t| V  qd S r(   )r5   reversedrF   r.   r.   r/   rH      s    zAmbiguous kl_divergence(z, z). Please register_kl())r$   NotImplementedminr6   warningswarnr@   RuntimeWarning)	r,   r-   matchesZleft_pZleft_qZright_qZright_pZleft_funZ	right_funr.   r+   r/   _dispatch_klq   s&   rQ   c                 C   s   t | tS )zI
    Helper function for obtaining infinite KL Divergence throughout
    )torchZ	full_liker   Ztensorr.   r.   r/   _infinite_like   s   rT   c                 C   s   t j| | S )z2
    Utility function for calculating x log x
    )rR   ZspecialZxlogyrS   r.   r.   r/   _x_log_x   s   rU   c                 C   sD   |  d}|  d}| d|| dd}|| jdd S )zp
    Utility function for calculating the trace of XX^{T} with X having arbitrary trailing batch dimensions
       N)sizeZreshapepowsumshape)ZbmatnrG   Z
flat_tracer.   r.   r/   _batch_trace_XXT   s   

r^   pqreturnc                 C   s   zt t| t|f }W n ty(   tt| t|}|t t| t|f< Y nw |tu r;td| jj d|jj || |S )a"  
    Compute Kullback-Leibler divergence :math:`KL(p \| q)` between two distributions.

    .. math::

        KL(p \| q) = \int p(x) \log\frac {p(x)} {q(x)} \,dx

    Args:
        p (Distribution): A :class:`~torch.distributions.Distribution` object.
        q (Distribution): A :class:`~torch.distributions.Distribution` object.

    Returns:
        Tensor: A batch of KL divergences of shape `batch_shape`.

    Raises:
        NotImplementedError: If the distribution types have not been registered via
            :meth:`register_kl`.
    z(No KL(p || q) is implemented for p type z and q type )r%   r2   KeyErrorrQ   rK   NotImplementedError	__class__r@   )r_   r`   r*   r.   r.   r/   r'      s   
c                 C   s   | j tjj|j tjj| j   }t||j dk< d|| j dk< d| j  tjj|jtjj| j  }t||j dk< d|| j dk< || S Nr   r   )probsrR   nnZ
functionalZsoftpluslogitsr   r_   r`   t1t2r.   r.   r/   _kl_bernoulli_bernoulli   s   rl   c           	      C   s   | j | j }|j |j }|j  |j  |  }| j  | j  |  }| j |j  t| j  }| j|j t| j }|| t| }|| | | | S r(   )concentration1concentration0lgammarR   digamma)	r_   r`   Zsum_params_pZsum_params_qrj   rk   t3t4t5r.   r.   r/   _kl_beta_beta   s   rt   c                 C   sh   | j |j k  rtd| j | j| j|j  | j   |j    }| j |j k}t|| ||< |S )NzKKL between Binomials where q.total_count > p.total_count is not implemented)Ztotal_countanyrc   rf   rh   log1prT   )r_   r`   klZinf_idxsr.   r.   r/   _kl_binomial_binomial   s   (rx   c                 C   sD   | j | j|j  }t||j dk|< d|| j dk|< |dS )Nr   rV   )rf   rh   r   Z	expand_asr[   )r_   r`   tr.   r.   r/   _kl_categorical_categorical   s   
rz   c                 C   sL   | j | j|j  }|  t| j  }|  t|j  }|| | S r(   )meanrh   _cont_bern_log_normrR   rv   rf   r_   r`   rj   rk   rq   r.   r.   r/   -_kl_continuous_bernoulli_continuous_bernoulli   s   r~   c                 C   s|   | j d}|j d}| |  }| j  |j   d}| j |j  }| j  | d }|| || d S )NrV   )concentrationr[   ro   rp   	unsqueeze)r_   r`   Zsum_p_concentrationZsum_q_concentrationrj   rk   rq   rr   r.   r.   r/   _kl_dirichlet_dirichlet  s   r   c                 C   s"   |j | j  }|  }|| d S Nr   ratelog)r_   r`   Z
rate_ratiorj   r.   r.   r/   _kl_exponential_exponential  s   
r   c                 C   s   t | t |kstddd | jD }|j}| j| }tjj| |dd}|j| | }t|||D ]\}}}	|| |	 }
|t	|
t
|j8 }q4|S )NzThe cross KL-divergence between different exponential families cannot                             be computed using Bregman divergencesc                 S   s   g | ]}|   qS r.   )detachZrequires_grad_)rD   npr.   r.   r/   rE      s    z+_kl_expfamily_expfamily.<locals>.<listcomp>T)Zcreate_graph)r2   rc   Z_natural_paramsZ_log_normalizerrR   ZautogradZgradr[   r<   r"   lenevent_shape)r_   r`   Z	p_nparamsZ	q_nparamsZ	lg_normalZ	gradientsresultZpnpZqnpgtermr.   r.   r/   _kl_expfamily_expfamily  s   
r   c                 C   sn   |j | j|j   }t|j t| j  }| j |j  t| j  }|j| j | j | j  }|| | | S r(   )r   r   r   rR   ro   rp   r_   r`   rj   rk   rq   rr   r.   r.   r/   _kl_gamma_gamma+  s
   r   c                 C   sl   | j |j  }|j|j  }| j|j  }|  | | }|t }t|d|   | }|| | dt  S r   )scalelocr   _euler_gammarR   expro   )r_   r`   Zct1Zct2Zct3rj   rk   rq   r.   r.   r/   _kl_gumbel_gumbel4  s   r   c                 C   s$   |    t|j | j  |j S r(   )entropyrR   rv   rf   rh   r_   r`   r.   r.   r/   _kl_geometric_geometric?  s   $r   c                 C      t | j|jS r(   )_kl_normal_normal	base_distr   r.   r.   r/   _kl_halfnormal_halfnormalD     r   c                 C   sV   | j |j  }| j|j  }|  }||j  }|t| | j   }|| | d S r   )r   r   absr   rR   r   )r_   r`   scale_ratioZloc_abs_diffrj   rk   rq   r.   r.   r/   _kl_laplace_laplaceI  s   

r   c                 C   s   | j |j kr
tdt|j|j|jt| j| j| j }t|j|j|j| j |j}|jj|j	d }t
jj|j|dd}| j|j d}t| j|j 	d }t|| j 	d }t|| j}	|| | |	 }
d||
 | | j d   S )NzKL-divergence between two Low Rank Multivariate Normals with                          different event shapes cannot be computedrW   FupperrV         ?r   )r   
ValueErrorr   _unbroadcasted_cov_factor_unbroadcasted_cov_diag_capacitance_trilr   r   mTr   rR   linalgsolve_triangularr[   r^   rsqrtsqrtmatmul)r_   r`   term1term3	qWt_qDinvAterm21term22Zterm23Zterm24term2r.   r.   r/   7_kl_lowrankmultivariatenormal_lowrankmultivariatenormalT  s6   
	r   c           	      C   s   | j |j kr
tdt|j|j|jd| jjddd 	d  }t
|j|j|j| j |j}|jj|jd }tjj|j|dd}t| j|j d }t|| j}|| }d|| | | j d	   S )
NKL-divergence between two (Low Rank) Multivariate Normals with                          different event shapes cannot be computedrX   rW   rV   Zdim1Zdim2Fr   r   r   )r   r   r   r   r   r   _unbroadcasted_scale_trildiagonalr   r[   r   r   r   r   rR   r   r   r^   r   r   )	r_   r`   r   r   r   r   r   r   r   r.   r.   r/   0_kl_multivariatenormal_lowrankmultivariatenormalv  s.   
	r   c                 C   s$  | j |j kr
tdd|jjddd d t| j| j| j	 }t
|j|j| j }tj|jjd d | jjd d }| j d }|j|||f }| j||| jdf }t| j |||f }ttjj||dd}	ttjj||dd}
|	|
 }d	|| | | j d   S )
Nr   rX   rW   rV   r   r   Fr   r   )r   r   r   r   r   r[   r   r   r   r   r   r   rR   _C_infer_sizer\   expandZ
cov_factorrY   Z
diag_embedr   r^   r   r   )r_   r`   r   r   combined_batch_shaper]   q_scale_trilZp_cov_factorZ
p_cov_diagr   r   r   r.   r.   r/   0_kl_lowrankmultivariatenormal_multivariatenormal  s>   

r   c           	      C   s   | j |j kr
td|jjddd d| jjddd d }tj|jj	d d | jj	d d }| j d }|j
|||f }| j
|||f }ttjj||dd}t|j|j| j }|d|| |   S )	NzvKL-divergence between two Multivariate Normals with                          different event shapes cannot be computedrW   rV   r   r   Fr   r   )r   r   r   r   r   r[   rR   r   r   r\   r   r^   r   r   r   r   )	r_   r`   Z
half_term1r   r]   r   Zp_scale_trilr   r   r.   r.   r/   )_kl_multivariatenormal_multivariatenormal  s(   
r   c                 C   sB   | j |j  d}| j|j |j  d}d|| d |   S NrX   r   r   r   rZ   r   r   )r_   r`   Z	var_ratiorj   r.   r.   r/   r     s   r   c                 C   r   r(   )rz   Z_categoricalr   r.   r.   r/   '_kl_onehotcategorical_onehotcategorical  r   r   c                 C   sX   | j |j  }|j| j }|j|  }|  }|| | d }t|| jj|jjk < |S r   )r   alphar   r   supportlower_bound)r_   r`   r   Zalpha_ratiorj   rk   r   r.   r.   r/   _kl_pareto_pareto  s   
r   c                 C   s&   | j | j  |j    | j |j   S r(   r   r   r.   r.   r/   _kl_poisson_poisson     &r   c                 C   s.   | j |j krt| j|jkrtt| j|jS r(   )Z
transformsrc   r   r'   r   r   r.   r.   r/   _kl_transformed_transformed  s
   r   c                 C   s<   |j |j | j | j   }t||j| jk|j | j k B < |S r(   )highlowr   r   r_   r`   r   r.   r.   r/   _kl_uniform_uniform  s   r   c                 C   s    |    | j|j  |j  S r(   )r   rf   r   r   r   r.   r.   r/   _kl_bernoulli_poisson   s    r   c                 C   s,   |    | j|j  t|j  |  S r(   )r   r{   rh   rR   rv   rf   r|   r   r.   r.   r/   _kl_beta_continuous_bernoulli  s   
r   c                 C   
   t | jS r(   )rT   rm   r   r.   r.   r/   _kl_beta_infinity     
r   c                 C   s,   |    |j  |j| j| j| j    S r(   )r   r   r   rm   rn   r   r.   r.   r/   _kl_beta_exponential  s   r   c                 C   sp   |    }|j |j|j   }|jd | j | j| j    }|j| j | j| j  }|| | | S r   )r   r   ro   r   r   rm   rp   rn   r   r.   r.   r/   _kl_beta_gamma  s   
r   c           	      C   s   | j | j | j  }|jd}|   }d|d tj   }|d|  | j | j d  |d d }|j| }|jdd }|| || | |  S r   )	rm   rn   r   rZ   r   mathpir   r   )	r_   r`   ZE_beta
var_normalrj   rk   rq   rr   rs   r.   r.   r/   _kl_beta_normal+  s   

r   c                 C   s>   |    |j|j   }t||j| jjk|j| jjk B < |S r(   )r   r   r   r   r   r   r   upper_boundr   r.   r.   r/   _kl_beta_uniform:  s    r   c                 C   r   r(   )rT   rf   r   r.   r.   r/   !_kl_continuous_bernoulli_infinityD  r   r   c                 C   s"   |    t|j |j| j  S r(   )r   rR   r   r   r{   r   r.   r.   r/   $_kl_continuous_bernoulli_exponentialI  s   "r   c                 C   sz   |    }dtdtj t|j|j   t|j }| jt| j	 d|j | j	  dt|j  }|| | S )Nr   g       @)
r   r   r   r   rR   Zsquarer   r   Zvariancer{   r}   r.   r.   r/   _kl_continuous_bernoulli_normalR  s   
( r   c              	   C   sV   |    |j|j   }ttt|j| jj	t
|j| jjt|t |S r(   )r   r   r   r   rR   wheremaxger   r   ler   	ones_liker   r   r.   r.   r/    _kl_continuous_bernoulli_uniform^  s   r   c                 C   r   r(   rT   r   r   r.   r.   r/   _kl_exponential_infinityk     
r   c                 C   sB   |j | j  }|j t| }|| |j  |jt  dt  S r   )r   r   rR   r   ro   r   )r_   r`   ratiorj   r.   r.   r/   _kl_exponential_gammas  s   r   c                 C   sR   | j |j }|j|j }| d }t|| |d  }| }|| | | S r   )r   r   r   r   rR   r   
reciprocal)r_   r`   scale_rate_prodloc_scale_ratiorj   rk   rq   r.   r.   r/   _kl_exponential_gumbel  s   r   c                 C   sp   |j d}| jd}dt|| d tj  }| }|j| j }|jdd }|d || | |  S r   )	r   rZ   r   rR   r   r   r   r   r   )r_   r`   r   Zrate_sqrrj   rk   rq   rr   r.   r.   r/   _kl_exponential_normal  s   r   c                 C   r   r(   )rT   r   r   r.   r.   r/   _kl_gamma_infinity  r   r   c                 C   s&   |    |j  |j| j | j  S r(   )r   r   r   r   r   r.   r.   r/   _kl_gamma_exponential  r   r   c                 C   s~   | j |j }|j|j }| jd | j  | j  | j }| | j|  }t|d|	  
| j  | }|| | S r   )r   r   r   r   rp   ro   r   rR   r   r   rZ   )r_   r`   Zbeta_scale_prodr   rj   rk   rq   r.   r.   r/   _kl_gamma_gumbel  s    r   c                 C   s   |j d}| jd}dt|| d tj  | j | j  }d| jd| j  | }|j	| j | j }d|j	d }|| jd | j
   || | |  S r   )r   rZ   r   rR   r   r   r   r   ro   r   rp   )r_   r`   r   Zbeta_sqrrj   rk   rq   rr   r.   r.   r/   _kl_gamma_normal  s"   r   c                 C   r   r(   rT   r   r   r.   r.   r/   _kl_gumbel_infinity     
r   c                 C   sx   | j |j  }|tdtj   }tj| d dd }| j| j t  |j |j  dd }| | | td  S )NrX   r      r   )r   r   r   r   r   rZ   r   r   )r_   r`   Zparam_ratiorj   rk   rq   r.   r.   r/   _kl_gumbel_normal  s
   &r   c                 C   r   r(   r   r   r.   r.   r/   _kl_laplace_infinity  r   r   c                 C   s~   |j d}| j d| }dtd| tj  }d| jd }| j|j }d|jd }| | || | |  d S r   )r   rZ   rR   r   r   r   r   )r_   r`   r   Zscale_sqr_var_ratiorj   rk   rq   rr   r.   r.   r/   _kl_laplace_normal  s   r   c                 C   r   r(   r   r   r.   r.   r/   _kl_normal_infinity  r   r   c                 C   s|   | j |j }| j|j d}|j |j }| d }|| }t| d|  | }| | | ddtdtj    S r   )r   r   rZ   r   rR   r   r   r   )r_   r`   Zmean_scale_ratioZvar_scale_sqr_ratior   rj   rk   rq   r.   r.   r/   _kl_normal_gumbel  s   &r   c                 C   s   | j |j  }| j|j }|| j }t|}tdtj | j td|d  }|t	td|  }| || |j  ddtdtj    S )NrX   g      r   r   )
r   r   rR   r   r   r   r   r   rZ   erf)r_   r`   Zloc_diffr   Zloc_diff_scale_ratiorj   rk   rq   r.   r.   r/   _kl_normal_laplace  s   

(,r   c                 C   r   r(   )rT   r   r   r.   r.   r/   _kl_pareto_infinity  s   
r   c                 C   sZ   | j |j }| j|  }| j }| j| | jd  }|| | d }t|| jdk< |S r   )r   r   r   r   r   r   )r_   r`   r   rj   rk   rq   r   r.   r.   r/   _kl_pareto_exponential!  s   
r   c                 C   s   | j  | j  }| j | }|j |j|j   }d|j | }|j| j | j  | jd  }|| | | d }t|| jdk< |S r   )r   r   r   r   r   ro   r   r   r_   r`   common_termrj   rk   rq   rr   r   r.   r.   r/   _kl_pareto_gamma,  s   r  c           	      C   s   d|j d }| j | jd  }tdtj |j  | j | j   }| j }| j|d | jd  }| j| |j d}|| || |  d }t	|| jdk< |S )NrX   r   )
r   rZ   r   r   r   r   r   r   r   r   )	r_   r`   r   r  rj   rk   rq   rr   r   r.   r.   r/   _kl_pareto_normal;  s   &
r  c                 C   r   r(   r   r   r.   r.   r/   _kl_poisson_infinityH  s   
r  c                 C   s   | j | j }t|}|jd t| j t| j |  | }|jd td| j  td| j  |  | }|j |j  |j|j   }|| | | }t|| j |j	j
k| j|j	jk B < |S r   )r   r   rR   r   rm   rU   rn   ro   r   r   r   r   r  r.   r.   r/   _kl_uniform_betaN  s.   
 r  c              	   C   sh   |    | j|j  t|j  |  }ttt	| j
|jjt| j|jjt|t |S r(   )r   r{   rh   rR   rv   rf   r|   r   r   r   r   r   r   r   r   r   r   r   r   r.   r.   r/    _kl_uniform_continuous_bernoullif  s    
r  c                 C   sB   |j | j| j  d | j| j |j    }t|| j|jjk < |S )NrX   )r   r   r   r   r   r   r   r   r.   r.   r/   _kl_uniform_exponetialx  s   ,r  c                 C   s   | j | j }| }|j |j|j   }d|j t| j t| j |  | }|j| j | j  d }| | | | }t|| j|jj	k < |S )Nr   rX   )
r   r   r   r   ro   r   rU   r   r   r   r  r.   r.   r/   _kl_uniform_gamma  s   r	  c                 C   sn   |j | j| j  }| j|j |j  }| j|j |j  }| d||   }|t| t|   }|| S )Nr   )r   r   r   r   r   rR   r   )r_   r`   r  Zhigh_loc_diffZlow_loc_diffrj   rk   r.   r.   r/   _kl_uniform_gumbel  s   r
  c                 C   st   | j | j }ttjd |j |  }|dd }| j | j d|j  d d}|d||  |jd  S )NrX      r   )	r   r   r   r   r   r   r   rZ   r   )r_   r`   r  rj   rk   rq   r.   r.   r/   _kl_uniform_normal  s
    r  c                 C   sl   | j | j }|j|j|j |  }t| j t| j | | }||jd  | }t|| j|jj	k < |S r   )
r   r   r   r   rZ   r   rU   r   r   r   )r_   r`   Zsupport_uniformrj   rk   r   r.   r.   r/   _kl_uniform_pareto  s   r  c                 C   s*   | j |j krtt| j|j}t|| j S r(   )Zreinterpreted_batch_ndimsrc   r'   r   r"   r   r.   r.   r/   _kl_independent_independent  s   r  c                 C   sD   | j |j  d| j|j d  }d| j  |j   }|| S )NrX      r   ri   r.   r.   r/   _kl_cauchy_cauchy  s   (r  c                  C   sb   dg} t tdd dD ]\}}| d|j d|j d qd| }tjr/t j|7  _d	S d	S )
zHAppends a list of implemented KL functions to the doc for kl_divergence.zLKL divergence is currently implemented for the following distribution pairs:c                 S   s   | d j | d j fS re   )r@   )Zp_qr.   r.   r/   <lambda>  s    z_add_kl_info.<locals>.<lambda>)keyz* :class:`~torch.distributions.z#` and :class:`~torch.distributions.`z
	N)sortedr$   appendr@   joinr'   __doc__)rowsr_   r`   Zkl_infor.   r.   r/   _add_kl_info  s   
r  )r   rM   	functoolsr   typingr   rR   r   r   Z	bernoullir   betar   Zbinomialr	   Zcategoricalr
   Zcauchyr   Zcontinuous_bernoullir   Z	dirichletr   distributionr   Z
exp_familyr   Zexponentialr   gammar   Z	geometricr   Zgumbelr   Zhalf_normalr   Zindependentr   Zlaplacer   Zlowrank_multivariate_normalr   r   r   Zmultivariate_normalr   r   normalr   Zone_hot_categoricalr   Zparetor   Zpoissonr   Ztransformed_distributionr    uniformr!   utilsr"   r#   r   r$   dicttupler2   __annotations__r%   __all__r&   r5   rQ   rT   rU   r^   r'   rl   rt   rx   rz   r~   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r  r  r  r  r  r  r	  r
  r  r  r  r  r  r.   r.   r.   r/   <module>   s  
+
&















!

$








	




	



























	

