a MSic@sddlZddlZddlZddlmZddlmZddlmZm Z ddZ gdZ gdZ gd Z gd Ze e gZe egZdd d Zejjd dZGdddeZdS)N) constraints) Distribution) broadcast_all lazy_propertycCs*t|}|}|r&|||}q|SN)listpop)ycoefresultr Y/var/www/html/django/DPS/env/lib/python3.9/site-packages/torch/distributions/von_mises.py _eval_poly s r)g?g$ @g03@g,?N?g2t?gIx?gtHZr?) e3E?g-5?gՒ+Hub?gJNYgTPÂ?g'gZ?gUL+ߐg;^p?)?gY?g(z?g*O?gZ9?g.h?gӰ٩=5?) rg.kg?VmgtZOZ?g<Q g'8`?gP⥝gqJ:N?g;PJ4qcCs|dks|dksJ|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)r _COEF_SMALLabslog _COEF_LARGEtorchwhere)xorderr Zsmalllarger 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)rzerosshapeboolrallrandrunbindcosmathpiranyrsignacos) loc concentrationZ proposal_rrdoneuu1u2u3zfcacceptr r r _rejection_sample4s , r8cseZdZdZejejdZejZdZ 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_rsuperr9__init__)selfr-r. validate_args batch_shape event_shapetaurho __class__r r r?XszVonMises.__init__cCsL|jr|||jt||j}|tdtjt |jdd}|S)Nr rr) _validate_args_validate_sampler.rr'r-r(rr)r)r@valuelog_probr r r rLds  "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) _extended_shaperemptyr-rrr8r.r=)r@ sample_shaper"rr r r sampleks zVonMises.samplecs`ztt||WStyZ|jd}|j|}|j|}t||||dYS0dS)NrI)rA) r>r9expandNotImplementedError__dict__getr-r.type)r@rBrAr-r.rFr r rQvs    zVonMises.expandcCs|jS)z8 The provided mean is the circular one. r-r@r r r meansz VonMises.meancCs|jSrrVrWr r r modesz VonMises.modecCs$dt|jddt|jddS)z< The provided variance is the circular one. rrHr)rr.exprWr r r variances zVonMises.variance)N)__name__ __module__ __qualname____doc__rrealpositivearg_constraintssupport has_rsampler?rLrno_gradr;rPrQpropertyrXrYrr[ __classcell__r r rFr r9Ds    r9)r)r(r torch.jitZtorch.distributionsr torch.distributions.distributionrtorch.distributions.utilsrrrZ_I0_COEF_SMALLZ_I0_COEF_LARGEZ_I1_COEF_SMALLZ_I1_COEF_LARGErrrjitscript_if_tracingr8r9r r r r s