Vpfd e7de7de6fdZ?de7de7de6fdZ@de7de6fdZAde7de6fdZBde7de6fdZCde7de6fdZDde7de6fdZEd ZFd!ZGd"ZHd#ZId$ZJd%ZKd&ZLd'ZMd(ZNd)ZOGd*d+eZPd,e6de6fd-ZQd.ZRd/ZSd0ZTd1ZUd2ZVd3ZWd4ZXd5ZYd6ZZd7Z[d8Z\e"e'e'e'gd9Z]e2j^e]e2j_eXeOd:;e1j`e]eGeGeFe#e'd<Zae1j`ead=e2j^eaee(e4j9e#e'd>Zbe2j^ebee(e4j:e1j`ebd?e"e'e'gd@Zce2j^ecee(e4j;e1j`eceLeMe"e'e'gdAZde2j^ede2j_eXeSd:;e"e'e'gdBZee2j^eee2j_eXeVd:;e1j`edeIeHe"e'e'gdCZfe2j^efe2j_eXeUd:;e1j`efeKeJe"e'e'gdDZge2j^ege2j_eXeWd:;e"e'e'gdEZhe2j^ehee(e4j@e#e'dFZie2j^eie2j_e\d:;e1jjeidGe#e'dHZke2j^ekee(e4jBdIZle1jjekele#e'dJZme1j`emdKe2j^emee(e4jCe#e'dLZne1j`endMe2j^enee(e4jDe#e'dNZoe1jjeodOe2j^eoee(e4jEdS)Pz` Special functions LAX decompositions for special functions into their StableHLO counterparts. )EnumN)partial)$add bitwise_and bitwise_not bitwise_orbroadcast_in_dimbroadcast_shapesconvert_element_typediveqexp full_likegegtleloglog1pltmulneneg reciprocalreduceselectsignsqrtsquarestandard_naryop standard_unopsub_const_dtype_float_nary_lower_hlo_ones_isnan_reduce) while_loop)dtypes)ad)mlir)chlo)Array ArrayLikeabxreturnc:t|||S)z1Elementwise regularized incomplete beta integral.)regularized_incomplete_beta_pbind)r0r1r2s T/var/www/html/nettyfy-visnx/env/lib/python3.11/site-packages/jax/_src/lax/special.pybetaincr8(s & + +Aq! 4 44c6t|S)z7Elementwise log gamma: :math:`\mathrm{log}(\Gamma(x))`.)lgamma_pr6r2s r7lgammar=,s q  r9c6t|S)z%Elementwise digamma: :math:`\psi(x)`.) digamma_pr6r<s r7digammar@0   r9mc8t||S)z-Elementwise polygamma: :math:`\psi^{(m)}(x)`.) polygamma_pr6)rBr2s r7 polygammarE4s   !Q  r9c8t||S)z2Elementwise regularized incomplete gamma function.)igamma_pr6r0r2s r7igammarI8s q!  r9c8t||S)z@Elementwise complementary regularized incomplete gamma function.) igammac_pr6rHs r7igammacrL<s 1  r9c8t||S)zDElementwise derivative of the regularized incomplete gamma function.)igamma_grad_a_pr6rHs r7 igamma_grad_arO@s   a # ##r9c8t||S)z5Elementwise derivative of samples from `Gamma(a, 1)`.)random_gamma_grad_pr6rHs r7random_gamma_gradrRDs  ! !!Q ' ''r9qc8t||S)z6Elementwise Hurwitz zeta function: :math:`\zeta(x, q)`)zeta_pr6)r2rSs r7zetarVHs Q  r9c6t|S)zpExponentially scaled modified Bessel function of order 0: :math:`\mathrm{i0e}(x) = e^{-|x|} \mathrm{i0}(x)` ) bessel_i0e_pr6r<s r7 bessel_i0erYL   1  r9c6t|S)zpExponentially scaled modified Bessel function of order 1: :math:`\mathrm{i1e}(x) = e^{-|x|} \mathrm{i1}(x)` ) bessel_i1e_pr6r<s r7 bessel_i1er]RrZr9c6t|S)z4Elementwise error function: :math:`\mathrm{erf}(x)`.)erf_pr6r<s r7erfr`Xs Ar9c6t|S)z]Elementwise complementary error function: :math:`\mathrm{erfc}(x) = 1 - \mathrm{erf}(x)`.)erfc_pr6r<s r7erfcrc\s Qr9c6t|S)zAElementwise inverse error function: :math:`\mathrm{erf}^{-1}(x)`.) erf_inv_pr6r<s r7erf_invrfarAr9ct|t|zt||zz }t|dz t| z|dz t|zz|z }||zSN)r=rrr)gr0r1r2lbeta partial_xs r7 betainc_gradxrmesr ))fQii &Q-- /%1q5E1"II%q5CFF"#%*+,,) Qr9c td)Nz7Betainc gradient with respect to a and b not supported. ValueError)rjr0r1r2s r7betainc_grad_not_implementedrqksLMMMr9c|t| |t|z t|zzt|z zSN)rr&rr=rjr0r2s r7 igamma_gradxruns< S!q588|s1vv--q 9 : : ::r9c(|t||zSrs)rOrts r7 igamma_gradarwqs ]1a  r9c&t||| Srs)rurts r7 igammac_gradxryt q!Q   r9c&t||| Srs)rwrts r7 igammac_gradar|wrzr9c td)Nz5polygamma gradient with respect to m is not supportedrorjrBr2s r7polygamma_gradmrzsJKKKr9c `|tt|t|d|zSrh)rErr"r~s r7polygamma_gradxr}s) Ys1fQll++Q / / //r9c D ddd d d fd} f d}dgR}ttt|t|jd|} dd | t | d| g} t ||| } | S) Nrricd|}t|}|}t||Srs)rr)values iterationiterations_remain_condvalues_unconverged_cond kIterationIdxkValuesUnconvergedIdxnum_iterationss r7 while_cond_fnz7lentz_thompson_barnett_algorithm..while_cond_fns;}%I >::$%:; -/F G GGr9c  |}|g R}|g R}t|t|| }t|}tt t ||||}t|t || }tt t ||||}t|}t ||}t | |}|dz|<|| <|| <|| <tt t|t|d} t| |<|S)Nri?) rr rrrabsrrrr!r"_any)rrpartial_numeratorpartial_denominatorcsmall_constantddeltahtolerance_comparisoninputskCIdxkDIdxkHIdxrrnth_partial_denominatornth_partial_numeratorsmall thresholds r7 while_body_fnz7lentz_thompson_barnett_algorithm..while_body_fnsi}%I--iA&AAA11)EfEEE %6u !F!FGGAq%((Nr#a&&.))>1==A %6u !F!FGGAr#a&&.))>1==A1 A 1IIE F5M5!!A&MF=F5MF5MF5Mc#eVE3-?-?"@"@AA9MM$()=$>$>F ! Mr9T)rrrr shaperr))rrrrrrrrrrrrrrrrs`````` @@@@@r7 lentz_thompson_barnett_algorithmrs!- % % %HHHHHHH 00/;F;;; R'((% 0 0 e%8%> C C ""! d1Yq^^A && m]F ; ;& r9c  |j  fd} fd}tttt|t|dt|t|dt |t|dt |t|d}t ||t|dz||zt|dzz }|}t |||}t |||}t ||tt|d|}ttj krdndtj j dz tj j dz |||||g } t|t|zt||zz } | t!t#||zt%| |zz| z z|z } t |t|t'd | } t || tt| d| S) Nc^t|g}t|t|dzt|d}t|t|d}|t|dz }|t|dz}t|}t|d} t|d} ||z ||z|zz|z|| |zz|| |zz| zzz } |||z z|z|| |zz| z || |zzzz } t|d} t || | }t || |S)z The partial numerator for the incomplete beta function is given here: http://dlmf.nist.gov/8.17.E23 Note that there is a special case: the partial numerator for the first iteration is one. rrri@r)r r rr r)rr0r1r2iteration_bcastiteration_is_eveniteration_is_oneiteration_minus_onerBonetwoeven_numerator odd_numerator one_numerator numeratordtypers r7nth_partial_betainc_numeratorzGregularized_incomplete_beta_impl..nth_partial_betainc_numeratorsV 'y%<.nth_partial_betainc_denominatorsQ&y%<>>>>>.44444 Z q)Aq//Bq)Aq//22)4)4q)Aq//9Q??3355-Q1a0QUYq#=N=N5NOPP & 1%%! 6**! 3yA#:#:;;!72:--333 <   "Q & . .u 5 5|E""&*225997; q!9 AYY "VAE]] 2( CFFQJrQ$>$I J J JQ N& -1eEll!;!;V D D& !63y/C/CV+L+L M MMr9ceZdZdZdZdZdS) IgammaModerirrN)__name__ __module__ __qualname__VALUE DERIVATIVESAMPLE_DERIVATIVErr9r7rrs %*r9r predicatesct|d}|j}ttt |}t ||t |S)NF)r"rtuplerangelenrr)rfpredicates_shapeall_dimensionss r7rrsL Z!%s#3445566.  Az> : ::r9c d}fd}||t|dt|d|t|dt|df}t|||} | d} | d} tjkr| |z|z St |t |t |dzz } tjkr|| | z| zz|z Stjkr| | | zz |z|z Std)Nc,t|dS)Nr)r)valss r7cond_fnz_igamma_series..cond_fns Q==r9c |\}}}}}}}|t|dz}|||z z||z||zz z }||z}|||z z}||z} tjkr/t|||z t j jk}n;t|t||z t j jk}|t|||dt|||dt|||dt|||dt|||dt|||dfS)Nrrirrr) r"rrrr*rrrr) renabledrransr2dc_dadans_da conditionalrmodes r7body_fnz_igamma_series..body_fnsQ,0)GQ35' F1bMMA QUOq1uQ/ /EoG QU A 'C zS6<3F3F3J)JKKkk #EGO 4 4 U8K8K8O OQQk  Waa!! Waa!! Wc47## Waa!! WeT!W%% WgtAw'' r9rirrrzInvalid IgammaMode) rr)rrrr@r"rrrp) axr2r0rrrrr init_valsrrr dlogax_das `` r7_igamma_seriesrs06 Q !Q1a!Yq!__ aOO) GWi 0 0$ Q# G' Z  "H>!ffwq6!Q<</000) Z """ y7* +a // z+++ sY & '! +a // ) * **r9ctt|t|}t|t|d}t|t|t d}tt |t|dt |t|d}tt|t|dt||}|t|z|z t|z }t |ttj |j  } t|}tt!t||| |g} t#|t|dt%|||t| ||t&jz t+|||t| t||t&j} t#|t-|d| } t#|t-|d| } t#t||t-|t d| } | S)Nrinfrir)rr'r r"rrrrrrr=r*rmaxrrr(r_igammac_continued_fractionrrrr) r0r2ris_nan x_is_zero x_is_infinity domain_error use_igammacr underflowroutputs r7 igamma_implr0s  fQii + +&F1aLL!!)Qq%,,//00-Bq&A,,//Ava||1D1DEE,Bq&A,,//Aq::+3q66zA~q !"c&,u--122233) 2ww"  j)\9fEFF H H'  1aLL!"aK,M,M"')9;;;2q![+k2J2JKK**,,   & )Yq!__f 5 5& -1a& 9 9& *\622AuU||,,f 6 6& -r9ctj|jd}fd}t|d|z }||zt|dz} t|d} t |d} |} |t|dz} | |z}| |z }t |d}t |d}t |d}t |d}| }|||zz |z }||||| | | || | |||||f}t |||}|d}t jkr||zS|d}t|t|z }t j kr,t|tt|||St j kr.tt|t|||zStd)Ncr|^}}}}}}}t|t|dkt|S)Ni)rr"r)rr_ans_t_y_xr_s r7rz,_igammac_continued_fraction..cond_fnMs9'+$GT2r2q1 q6!T??*DMM : ::r9c|\}}}}}}}}} } } } } }}|t|dz}|t|dz}|t|dz}||z}||z| |zz }||z| |zz }t|t|d}||z }t|tt t |||t |d}t|||}| |z|z | |zz | |zz}||z|z | |zz | |zz}t|t |||zz ||}t|t||z t |d}|} |}|} |}| } |} |} |}tt|tt|}t|t| t| | } t|t|t||}t|t| t| | } t|t|t||}t|t| t| | } t|t| t| | } t|t| t| | } t|t|t||}tj krt||k}n"t||t|k}|t|||dt|||dt|||dt|||d|t|||dt|||dt|| |dt|| |d t|| |d t|| |d t|| |d t|||d t|||dfS)Nrirrrrr ) r"rrrr r!rrrrrrr)rrrtyzrpkm1qkm1pkm2qkm2dpkm2_dadqkm2_dadpkm1_dadqkm1_darycpkqk qk_is_nonzerordpk_dadqk_da dans_da_newgrad_conditionalrescalerrrs r7rz,_igammac_continued_fraction..body_fnQsK:>8Wc1aAtT48Xw F1aLLA F1aLLA F1aLLA QB D2I B D2I Br6"a==))M RA}c#c#qkk1"5"566 !QHHA 3 ' 'C \D 8b= 04!8 ;F \D 8b= 04!8 ;FFS6\,A2(F(FPPKm!+"788'3355 D D D DHHHHR*VB__5566G '3tVD#%6%677 > >D '3tVD#%6%677 > >D '3tVD#%6%677 > >D '3tVD#%6%677 > >Dgs8VHc-B-BCCXNNHgs8VHc-B-BCCXNNHgs8VHc-B-BCCXNNHgs8VHc-B-BCCXNNH zS11kk6"2C888::k  d1g & & DG $ $ DG $ $ DG $ $ tAw ' ' tAw ' ' tAw ' ' tAw ' ' 48 , , 48 , , 48 , , 48 , , d2h / / 11r9rirrzInvalid mode: )r*rrr"rr)rrrr@rrrrrrp)rr2r0rrrrrrrrr r r r rrr rrrrrrrrs ` @r7rrJs U#;;;@1@1@1@1@1@1D QllQ!!efQll! Qll! 1a$ $ VAq\\ $ Q$ t #1oo! q!__( q!__( q!__(R( h &$ .'qQTTTT8Xw@) GWi 0 0$ Q# Z  "9 H'!ff ") Z """ r3s3 **G44 5 55 z+++ s7CY//0014 5 55 ,d,, - --r9c tt|t|dt|t|d}tt|t|dt||}|t |z|z t |z }t|t t j|j }tt||}t|}t|||t|||tj}t|||t|t||tj} t!|t|d|z | } t#|t|t%d} t!| t'| d| } t!|t'|d| S)Nrrir)rrr"rrr=r*rrrrrrrrrrr rr) r0r2r out_of_range use_igammarrr igamma_calligammac_cf_callrrs r7 igammac_implrsBq&A,,//Ava||1D1DEE,"Qq! --r!Qxx88*3q66zA~q !"c&,u--122233)  <;; < <' 2ww"r1aWj)I)I$j&688+/AqZ00115*:JLL/ *fQll[8/ J J&Qq%,,//00- -61!5!5v > >&  i1oov 6 66r9cvtt|t|}t|t|d}tt |t|dt |t|d}t t|t|dt||}|t|z|z t|z }t |ttj |j  }t|}tttt||||} t|t!|||t | ||t"j t'|||t | t||t"j} t|t| d| } tt||t|t)d| } | SNrrir)rr'r rrrrrrr=r*rrrrrrrrrr r0r2rrrrrrrrrs r7igamma_grad_a_implr#s fQii + +&IaNN##)Bq)Aq//22Bq)Aq//4J4JKK,Bq)Aa..112a88<<+3q66zA~q !"c&,u--122233) 2ww"  :j// )$+$+,244 5 5' + Q;w +L+L!& (=???2q![+k2J2JKK*/11 2 2& )Yva00& 9 9& *\622AuU||,,f 6 6& -r9ctt|t|}t|t|d}tt |t|dt |t|d}t t|t|dt||}|t|z|z t|z }t |ttj |j j  }t|}tttt||||} t!|t#|||t | ||t$j t)|||t | t||t$j} t!|t| d| } t!t||t|t+d| } | Sr!)rr'r rrrrrrr=r*rrrrrrrrrrrr"s r7random_gamma_grad_implr%s fQii + +&IaNN##)Bq)Aa..112a1Q3H3HII,Bq)Aa..112a88<<+3q66zA~q !"c&,qw//344455) 2ww"  :j //($*$*+133 4 4' + Q;w +L+L!& (DFFF2q![+k2J2JKK). 0LNN O O& )Yva00& 9 9& *\622AuU||,,f 6 6& -r9cfd}|S)Nctd|Dfd|D}|dj}|tjkp|tjk}|rd|D}tj}n|}|d|i}|rt||}|S)Nc3$K|] }|jV dSrs)r.0r0s r7 z>_up_and_broadcast..up_and_broadcast..s$*A*Aq17*A*A*A*A*A*Ar9c pg|]2}t|tt|j3Sr)r listrndim)r*r0broadcasted_shapes r7 z?_up_and_broadcast..up_and_broadcast..s7 V V VA Q 14af 3F3F G G V V Vr9rcBg|]}t|tjSr)r rrr)s r7r0z?_up_and_broadcast..up_and_broadcast..s% @ @ @a"1bj11 @ @ @r9r)r rr*bfloat16rfloat16rr )argsa_dtype needs_upcasta_x_typerr/doits @r7up_and_broadcastz+_up_and_broadcast..up_and_broadcasts(*A*AD*A*A*AB V V V VQU V V VD1gmGfo-FBJ1FL @ @4 @ @ @dhhh T4 (x ( (F5#FG44f Mr9r)r8r9s` r7_up_and_broadcastr:s$ r9ct|d}t|d}t|d}|D]}|}|}||z|z t||z} d||z zS)Nr?)r)r2 coefficientsb0b1b2rs r7evaluate_chebyshev_polynomialrAsl1~~"1~~"1~~" ''a B B R"yA &BB Rr9ctjgd}tjgd}t|}t|d}t|d}t|d}t ||z|z |}t t ||z |z |t |}t|dk||S)z Computes an approximation to the modified Bessel function of the first kind, zeroth order. The following implementation follows Cephes' F32 and F64 implementation of i0e. )0 K5dMv;p>"c쑾$>'doҾY(X?>ZY&+|t(?RBuZ?I ^qa?!N-Ί>?-4pKw?Wӿ*5N?)["d,->mրVX>na>+A>Rx?I墌k? b?r<r@@ @rarrayrrrAr rrr2 i0e_coeffs_a i0e_coeffs_bhalfr thirty_two result_le_8 result_gt_8s r7 _i0e_impl32rg s <<<  ,   , !ff! 1c  $!S#D!!*-dQhnlKK+1*q.32F2>@@AEaJJ+ S+{ 3 33r9ctjgd}tjgd}t|}t|d}t|d}t|d}t ||z|z |}t t ||z |z |t |}t|dk||S)N)g4!\Tg}b3rCrDrErFrGrHrIrJrKrLrMrNrOrPrQrRrSrT)gT`g0fFVg!rUrVrWrXrYrZr[r<rr\r]r^r`s r7 _i0e_impl64ri-s>>>," ! ! !,  !ff! 1c  $!S#D!!*-dQhnlKK+1*q.32F2>@@AEaJJ+ S+{ 3 33r9c|jtjkrt|S|jtjkrt |S|j}|tj}tt ||Srs)rrfloat64rirrgrr )r2x_dtypes r7bessel_i0e_implrmZsjW  q>>w"* q>>gG A  A 8 88r9regularized_incomplete_betaF)multiple_resultsr=c<t|t|Srs)rr@rjr2s r7rrssQ !3!3r9r@c Zt|tt|d|Srh)rrEr"rqs r7rrrrxs"#a6!Q<<)C)C"D"Dr9rErIrOrLrRrVrYcL|t|t||zz zSrs)r]r)rjrr2s r7rrrrs jmmd1ggk.I)Jr9r]c^tjt|j}t ||k}t ||t ||}t||t|t|zzz }t ||t |d}||zS)Nr<) r*rr#rrrrrYrr)rjrr2r x_is_not_tinysafe_xdy_dxs r7_bessel_i1e_jvprys VAYY##a&&3,- -Ia$5$5 6 6& V  qDLL:f3E3E$EF F% yC'8'8 9 9% Ur9r`c tt|dtjtjz t|t t t|SNrrr"rrpirrrrqs r7rrrrsLc&B,?"@"@"%aS^^)<)<"="=??r9rcc tt|dtjtjz t|t t t|S)Ngr|rqs r7rrrrsNs6!S2725>>-A#B#B#&q#c&))nn*=*=#>#> @ @r9rfc tt|tjtjdz t|t t |Sr{)rr"rrr}rr)rjrr2s r7rrrrsHF1bgbennr6I,J,J,/3vc{{3C3C,D,D)F)Fr9)p__doc__enumrnumpyr functoolsrjax._src.lax.laxrrrrr r r r r rrrrrrrrrrrrrrrrrrr r!r"r#r$r%r&r'r(jax._src.lax.control_flowr)jax._srcr*jax._src.interpretersr+r,jax._src.lib.mlir.dialectsr-jax._src.typingr.r/r8r=r@rErIrLrOrRrVrYr]r`rcrfrmrqrurwryr|rrrrrrrrrrr#r%r:rArgrirmr5register_lowering lower_fundefjvpr;r?rDrGrNrKrQrUrXdefjvp2r\ryr_rbrerr9r7rs@  OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO100000$$$$$$&&&&&&++++++,,,,,,,,5y5Y5955555iEyU  y U    iI%yY5$Y$9$$$$$((y(U((((I)) ) 9I% yU NNN;;;!!!      LLL000 000f:N:N:Nx ;U;u;;;; 1+1+1+f4g.g.g.R777$((&!4!4!4F+4+4+4Z 9 9 9!0 VV;!=!=4$.""#CDD"')))***  ' = * *  (33444x$+!F!FGGG M&) , , y''/4<"H"HIII  )DDEEEovv. << {GGOT^$L$LMMM  +888 ?FF+X 6 6x0A0A+0N0NAF"H"H"HIII"/66"2ODD%t~&7&78J&K&K7<>>>???  (L,/// OVV,i 8 8 y%t~&7&7 &E&E7<>>>???  )]M222%ovv&68KLL*%t~&7&78N&O&O7<>>>??? &&)6 2 2vww BBCCC}V\22 |%t~o7<>>>???  <JJKKK}V\22 |AACCC  <))) fe$$  %??@@@uggotx@@AAA vv & &  &@@AAAvww BBCCC M&) , ,  9FFGGGy''/4<"H"HIIIIIr9