U d @sddlZddlZddlZddlmZddlmZddlmZm Z ddZ ddd d d d d gZ dddddddddg Z dddddddgZ dddd d!d"d#d$d%g Ze e gZe egZd,d&d'Zejjd(d)ZGd*d+d+eZdS)-N) constraints) Distribution) broadcast_all lazy_propertycCs*t|}|}|r&|||}q|SN)listpop)yZcoefresultr A/tmp/pip-unpacked-wheel-ua33x9lu/torch/distributions/von_mises.py _eval_poly s r g?g$ @g03@g,?N?g2t?gIx?gtHZr?ge3E?g-5?gՒ+Hub?gJNYgTPÂ?g'gZ?gUL+ߐg;^p??gY?g(z?g*O?gZ9?g.h?gӰ٩=5?g.kg?VmgtZOZ?g<Q g'8`?gP⥝gqJ:N?g;PJ4qcCs|dks|dkst|d}||}t|t|}|dkrF||}|}d|}|d|t|t|}t|dk||}|S)zX Returns ``log(I_order(x))`` for ``x > 0``, where `order` is either 0 or 1. rg@r)AssertionErrorr _COEF_SMALLabslog _COEF_LARGEtorchwhere)xorderr smallZlarger r r r _log_modified_bessel_fns "rc Cstj|jtj|jd}|stjd|j|j|jd}|\}}}t t j |} d|| || } ||| } | d| |dk| | d| dkB} | rt| |d| |}|| B}q|t j |dt j t j S)Ndtypedevice)rrr)rzerosshapeboolrallZrandrZunbindcosmathpiranyrsignacos) loc concentrationZ proposal_rrdoneuu1u2u3zfcacceptr r r _rejection_sample4s , r5cseZdZdZejejdZejZdZ dfdd Z ddZ e e fd d Zfd d Zed dZeddZeddZZS)VonMisesa" A circular von Mises distribution. This implementation uses polar coordinates. The ``loc`` and ``value`` args can be any real number (to facilitate unconstrained optimization), but are interpreted as angles modulo 2 pi. Example:: >>> m = dist.VonMises(torch.tensor([1.0]), torch.tensor([1.0])) >>> m.sample() # von Mises distributed with loc=1 and concentration=1 tensor([1.9777]) :param torch.Tensor loc: an angle in radians. :param torch.Tensor concentration: concentration parameter )r*r+FNcst||\|_|_|jj}t}ddd|jd}|d|d|j}d|dd||_tt | |||dS)Nrr) rr*r+r!rSizesqrt _proposal_rsuperr6__init__)selfr*r+ validate_args batch_shapeZ event_shapetaurho __class__r r r<XszVonMises.__init__cCsL|jr|||jt||j}|tdtjt |jdd}|S)Nrrr) _validate_argsZ_validate_sampler+rr$r*r%rr&r)r=valuelog_probr r r rGds  "zVonMises.log_probcCs6||}tj||jj|jjd}t|j|j|j|S)z The sampling algorithm for the von Mises distribution is based on the following paper: Best, D. J., and Nicholas I. Fisher. "Efficient simulation of the von Mises distribution." Applied Statistics (1979): 152-157. r) Z_extended_shaperemptyr*rrr5r+r:)r=Z sample_shaper!rr r r sampleks zVonMises.samplecsbztt||WStk r\|jd}|j|}|j|}t||||dYSXdS)NrE)r>) r;r6expandNotImplementedError__dict__getr*r+type)r=r?r>r*r+rBr r rJvs   zVonMises.expandcCs|jS)z8 The provided mean is the circular one. r*r=r r r meansz VonMises.meancCs|jSrrOrPr r r modesz VonMises.modecCs$dt|jddt|jddS)z< The provided variance is the circular one. rrDr)rr+exprPr r r variances zVonMises.variance)N)__name__ __module__ __qualname____doc__rrealZpositiveZarg_constraintsZsupportZ has_rsampler<rGrZno_gradr8rIrJpropertyrQrRrrT __classcell__r r rBr r6Ds    r6)r)r%rZ torch.jitZtorch.distributionsrZ torch.distributions.distributionrZtorch.distributions.utilsrrr Z_I0_COEF_SMALLZ_I0_COEF_LARGEZ_I1_COEF_SMALLZ_I1_COEF_LARGErrrZjitZscript_if_tracingr5r6r r r r s4