uki(ddlmZddlmZddlmZddlZddlm Z ddlm Z ddlm Z ddlmZ dd lm Z dd lmZmZdd lmZmZdd lmZeed dGddZdZe j0ddZy)) dataclass)partial)AnyN)api)dtypes)lax)numpy)random)check_arraylikepromote_dtypes_inexact)linalgspecial)register_pytree_node_classTF)frozeninitceZdZUdZeed<eed<eed<eed<eed<ddZd Zd Ze d Z e d Z e d Z dZdZdZdZdZddZdZdZddZddZdZy) gaussian_kdeaGaussian Kernel Density Estimator JAX implementation of :class:`scipy.stats.gaussian_kde`. Parameters: dataset: arraylike, real-valued. Data from which to estimate the distribution. If 1D, shape is (n_data,). If 2D, shape is (n_dimensions, n_data). bw_method: string, scalar, or callable. Either "scott", "silverman", a scalar value, or a callable function which takes ``self`` as a parameter. weights: arraylike, optional. Weights of the same shape as the dataset. neffdatasetweights covarianceinv_covNcJtd|tj|}tjt j |tjr td|jdkDs td|j\}}|{td|t||\}}tj|}|tj|z}|j dk7r tdt#||k7r=tdt|\}tj$|d|z |j }|j'd ||j'd ||j'd dtj|d zz }|d n|}|d k(rtj(|d|dzz }nr|dk(r&tj(||d zzdz d|dzz }nGtj*|rt-|t.s|}nt1|r ||}n tdtjtj2|dd|}tj4j7|} ||d zz} | |d zz } |j'd| |j'd| y)Nrz*gaussian_kde does not support complex dataz.`dataset` input should have multiple elements.z*`weights` input should be one-dimensional.z%`weights` input should be of length n?)dtyperrrscottg silvermang@zD`bw_method` should be 'scott', 'silverman', a scalar, or a callable.F)rowvarbiasaweightsrr)r jnp atleast_2dr issubdtyperrnpcomplexfloatingNotImplementedErrorsize ValueErrorshaper atleast_1dsumndimlenfull_setattrpowerisscalar isinstancestrcallablecovr inv) selfr bw_methodrdnrfactordata_covariance data_inv_covrrs S/home/cdr/jupyterlab/.venv/lib/python3.12/site-packages/jax/_src/scipy/stats/kde.py__init__zgaussian_kde.__init__3sANG,nnW%G 7+R-?-?@  L MM <>/2L 619,JVQY&GMM, +MM)W%c4tj||||SN)object __setattr__)r:namevalues rAr2zgaussian_kde._setattrcs tT5) LrCcx|j|j|j|j|jfdfSrE)rrrrrr:s rA tree_flattenzgaussian_kde.tree_flattenis4 YY dllDOO \\   ""rCc~|j|}|jd|d|jd|d|jd|d|jd|d|jd |d |S) Nrrrrrrrrr)__new__r2)clsaux_datachildrenkdes rAtree_unflattenzgaussian_kde.tree_unflattenmsr ++c CLL!%LLHQK(LLHQK(LLx{+LLHQK( JrCc4|jjdS)Nrrr,rKs rAr<zgaussian_kde.dx <<  a  rCc4|jjdS)NrrVrKs rAr=zgaussian_kde.n|rWrCctd||j|}td|jj|j dddf|j|j }|dddfS)z.Evaluate the Gaussian KDE on the given points.evaluateFNrr _reshape_points_gaussian_kernel_evalrTrr)r:pointsresults rArZzgaussian_kde.evaluates\J'  ! !& )F "5$,,..$,,q$w:O#)88T\\;F !Q$<rCc$|j|SrErZ)r:r_s rA__call__zgaussian_kde.__call__s ==  rCctjtj|}tj|}|j|j fk7rt d|j |j|j |j fk7rt d|j tj|j|z}tjdtjz|j ztjtj|dz}d|z }t|||j |j"|S)z2Integrate the distribution weighted by a Gaussian.zmean does not have dimension z#covariance does not have dimension rrr)r$r-squeezer%r,r<r+r cho_factorrsqrtr'piproddiag_gaussian_kernel_convolverr)r:meanr8cholnorms rAintegrate_gaussianzgaussian_kde.integrate_gaussians >>#++d+ ,D .. C zzdffY 6tvvh? @@ yyTVVTVV$$ >   T__u/?/?? @D 88AI  &#((472C)D DD :D"ffuww.dE]UDMFB!SXXg7tRZZ jj*!!#-F 776BJJ& ''rCctj|\}}tj||j||j}tj |t j|j|jj|j||jjj}|jdd|f|zS)ajRandomly sample a dataset from the estimated pdf Args: key: a PRNG key used as the random key. shape: optional, a tuple of nonnegative integers specifying the result batch shape; that is, the prefix of the result shape excluding the last axis. Returns: The resampled dataset as an array with shape `(d,) + shape`. )r,p)r,rN) r splitchoicer=rmultivariate_normalr$zerosr<rrrr^)r:keyr,ind_keyeps_keyindepss rAresamplezgaussian_kde.resamples||C(GW --u EC  $ $W%(YYtvvt7L7L%M%)__+0+/<<+=+=  ?@Aq  <<3 # %%rCc$|j|S)zProbability density functionrb)r:xs rApdfzgaussian_kde.pdfs == rCctd||j|}td|jj|j dddf|j|j }|dddfS)z Log probability density functionlogpdfTNrr[)r:rr`s rArzgaussian_kde.logpdfs\Ha  QA "4ag9N#$33 6F !Q$<rCc~~~td)4This method is not implemented in the JAX interface.z>only 1D box integrations are supported; use `integrate_box_1d`r))r: low_bounds high_boundsmaxptss rA integrate_boxzgaussian_kde.integrate_boxsK H JJrCc~td)rz:dynamically changing the bandwidth method is not supportedr)r:r;s rA set_bandwidthzgaussian_kde.set_bandwidths D FFrCctjtj|tj r t dtj|}|j\}}||jk7r]|dk(r3||jk(r$tj||jdf}|Stdj||j|S)Nz1gaussian_kde does not support complex coordinatesrz2points have dimension {}, dataset has dimension {})rr&rrr'r(r)r$r%r,r<reshaper+format)r:r_r<ms rAr\zgaussian_kde._reshape_pointss 6*B,>,>?  = ?? ^^F #F <.kernel sS S3773::g.>#?@@ @C WWW  ## swws| ##rCcRtjd|dS)N)Nrrrwrr)rry)rrr_reducevaluess rAz'_gaussian_kernel_eval..s+V$JCHHV\$J%&%8)*,rC)r r,r+r choleskyr$dotr.rrjr'rhr logsumexprry) rr_rxi precisionr< whiteningreduced_kernel mapped_kernelrrrs ``` @@@rAr]r]s"8 fb)#%&&"i ll1o!XXa[A @ AA__A < ==ooit4) 7769 %& wwr9" WWSWW hhy!AgBEE (::;($!'7  CGG&,.((>*- r rC) dataclassesr functoolsrtypingrr r'jax._srcrrrr$r jax._src.numpy.utilr r jax._src.scipyr rjax._src.tree_utilrrrkjitr]rrCrArs{"!G*9 $U#OO$Od1 rC