Vpf UddlmZmZddlZddlZddlZddlmZddlmZddl m Z ddl m Z ddlmZmZddlmZdeed efd Zdeed efd Zd ed ed efdZd ed ed efdZdddeedZeeeeegeffed<ded efdZded eeeeffdZ ded eeeeffdZ!ej"ej#ddedeededed ed ef d!Z$ d%dedeededed ef d$Z%dS)&)CallableSequenceN)api)util)lax) ArrayLikeArray)safe_ziparrsreturnc@tjtj|SN) functoolsreduceoperatormulr s V/var/www/html/nettyfy-visnx/env/lib/python3.11/site-packages/jax/_src/scipy/ndimage.py_nonempty_prodr  (, - --c@tjtj|Sr)rrraddrs r _nonempty_sumrrrindexsizecL|dz }tj||zd|zz|z S)N)jnpabs)rrss r_mirror_index_fixerr#"s- qA 7EAI!a%(1, - --rcftjtd|zdzd|zdzdz dS)Nrr)r floor_divider#rrs r_reflect_index_fixerr''s5  /% 1T6!8DDqH! L LLrc|Srr&s rr*+sErc4tj|d|dz S)Nrr)r clipr&s rr*r*,s38E1dQh#?#?rc ||zSrr)r&s rr*r*-s  r)constantnearestwrapmirrorreflect _INDEX_FIXERSacvtj|jtjr|nt j|Sr)r issubdtypedtypeintegerrround)r4s r_round_half_away_from_zeror:3s)nQWck22 D ! Dr coordinatect|tj}|jd}||fgSNr)r:astyper int32r7type)r;rweights r_nearest_indices_and_weightsrB7sB $Z 0 0 7 7 B B%    # #& &/ rctj|}||z }d|z }|tj}||f|dz|fgSr=)r floorr>r?)r;lower upper_weight lower_weightrs r_linear_indices_and_weightsrH=sQ )J  %e#,\!, ,,sy ! !% , %!)\!: ;;r)r)static_argnumsinput coordinatesordermodecvalc"tj|}d|D}tj||j}t||jkr5t dt||jt|}|5td|tt|dkrd}nd}|dkrt} n|dkrt} ntd g} t||jD]_\} } | | } g}| D]5\}}||| }||| }||||f6| |`g}t!j| D]}t%j|\}}}t)d |Dr ||}n;t+jt.j|}tj||||}|t5||zt7|}tj|jtjrt=|}||jS) Nc6g|]}tj|Sr))r asarray).0cs r z$_map_coordinates..Is 999S[^^999rzAcoordinates must be a sequence of length input.ndim, but {} != {}zajax.scipy.ndimage.map_coordinates does not yet support mode {}. Currently supported modes are {}.r.cd|k||kzS)Nrr)r&s rr*z"_map_coordinates..WsAJ54<#@rcdS)NTr)r&s rr*z"_map_coordinates..Ys4rrrz=jax.scipy.ndimage.map_coordinates currently requires order<=1c3K|]}|duV dS)TNr))rTvalids r z#_map_coordinates..ps& 1 1U5D= 1 1 1 1 1 1r) r rSr7lenndim ValueErrorformatr3getNotImplementedErrorsetrBrHzipshapeappend itertoolsproductrunzip3allrrrand_whererrr6r8r:r>)rLrMrNrOrP input_arrcoordinate_arrs index_fixeris_valid interp_funvalid_1d_interpolationsr;r interp_nodes valid_interprrA fixed_indexrZoutputsitemsindices validitiesweights contribution all_validresults r_map_coordinatesr}Esk%  )99[999/ T9? + +$''  &s;'7'7 H H J JJ!!$''+  ,,2F4]9K9K,L,L N NN Z@@HH''H aZZ-JJ zz,JJ G I IIoy??11j$:j))LL%88 vKt,,khud##e;v67777""<0000 ' "9:;;e#';u#5#5 GZ 1 1j 1 1 111Dw'll"8=*==iYy)G*'**\9::::  ! !&^IOS[110 ' / /F y ' ''rr.c(t|||||S)ag Map the input array to new coordinates using interpolation. JAX implementation of :func:`scipy.ndimage.map_coordinates` Given an input array and a set of coordinates, this function returns the interpolated values of the input array at those coordinates. Args: input: N-dimensional input array from which values are interpolated. coordinates: length-N sequence of arrays specifying the coordinates at which to evaluate the interpolated values order: The order of interpolation. JAX supports the following: * 0: Nearest-neighbor * 1: Linear mode: Points outside the boundaries of the input are filled according to the given mode. JAX supports one of ``('constant', 'nearest', 'mirror', 'wrap', 'reflect')``. Note the ``'wrap'`` mode in JAX behaves as ``'grid-wrap'`` mode in SciPy, and ``'constant'`` mode in JAX behaves as ``'grid-constant'`` mode in SciPy. This discrepancy was caused by a former bug in those modes in SciPy (https://github.com/scipy/scipy/issues/2640), which was first fixed in JAX by changing the behavior of the existing modes, and later on fixed in SciPy, by adding modes with new names, rather than fixing the existing ones, for backwards compatibility reasons. Default is 'constant'. cval: Value used for points outside the boundaries of the input if ``mode='constant'`` Default is 0.0. Returns: The interpolated values at the specified coordinates. Examples: >>> input = jnp.arange(12.0).reshape(3, 4) >>> input Array([[ 0., 1., 2., 3.], [ 4., 5., 6., 7.], [ 8., 9., 10., 11.]], dtype=float32) >>> coordinates = [jnp.array([0.5, 1.5]), ... jnp.array([1.5, 2.5])] >>> jax.scipy.ndimage.map_coordinates(input, coordinates, order=1) Array([3.5, 8.5], dtype=float32) Note: Interpolation near boundaries differs from the scipy function, because JAX fixed an outstanding bug; see https://github.com/google/jax/issues/11097. This function interprets the ``mode`` argument as documented by SciPy, but not as implemented by SciPy. )r})rLrMrNrOrPs rmap_coordinatesrsh %eT4 @ @@r)r.r~)&collections.abcrrrrfrjax._srcrrjaxr jax.numpynumpyr jax._src.typingrr jax._src.utilr rcrrintr#r'r3dictstr__annotations__r:listtuplerBrHpartialjitr}rr)rrrs/........,,,,,,,,)))))).%.U......E.....u.C.E.... MMSMUMMMM*)?? , ,!# ;; tC5#,"5667E%EEEEEEUtE%BR