`i%tddlZddlZddlmZddlmZddlmZddlmZddZ dZ dd Z dd Z dd Z dd ZdS)N)_core)internal)_util)specialFctjtjg}|s|tjtjgz }|L|j|vr!tj|j|j}ntj|j|d}ndt|tur/||vrtdtj|j|}n|j|jkrtd|S)Ndtypezoutput type not supportedzoutput shape not correct) cupy complex64 complex128float32float64r emptyshapetype RuntimeError)outputinput complex_onlytypess p/home/jaya/work/projects/VOICE-AGENT/VIET/agent-env/lib/python3.11/site-packages/cupyx/scipy/ndimage/_fourier.py_get_output_fourierr s ^T_ -E . $, -- ~ ;%  Z 5;???FFZ 59===FF f     :;; ;EKv666  $ $5666 Mc`d|z|jfzd||z dz zz}||S)z>Promote a 1d array to ndim with non-singleton size along axis.)r)sizereshape)arrndimaxisnd_shapes r _reshape_ndr#s8d{ch[(44$;?+CCH ;;x  rr c|j}t||}tj||}t j||d}t j||tt||j D]\}\}} ||kr,|dkr&tj | |j j} | |z} ntj| } | |j jd} | | z} ||zdz } | dt$jzt$jz| zz} tj| | t+| || } || z}|S) aMultidimensional Gaussian shift filter. The array is multiplied with the Fourier transform of a (separable) Gaussian kernel. Args: input (cupy.ndarray): The input array. sigma (float or sequence of float): The sigma of the Gaussian kernel. If a float, `sigma` is the same for all axes. If a sequence, `sigma` has to contain one value for each axis. n (int, optional): If `n` is negative (default), then the input is assumed to be the result of a complex fft. If `n` is larger than or equal to zero, the input is assumed to be the result of a real fft, and `n` gives the length of the array before transformation along the real transform direction. axis (int, optional): The axis of the real transform (only used when ``n > -1``). output (cupy.ndarray, optional): If given, the result of shifting the input is placed in this array. Returns: output (cupy.ndarray): The filtered output. sigmarrFcopyoutr r!)r rr_normalize_axis_indexr_fix_sequence_argrelementwise_copy enumerateziprr arangerealr fftfftfreqastypenumpypiexpr#) rr%nr!rr sigmasaxsigmakax_sizerscales rfourier_gaussianr@#sV0 :D  / /F  )$ 5 5D  $UD' : :F 5&)))!*3vv|+D+D!E!E VW ::!a%%+gV[->???C 1HCC(""7++Cjj*j77 s "$ EH ux'500 ##Dr222#  MrcR|j}t||}tj||}t j||d}t j||tt||j D]\}\}}||kr,|dkr&tj ||j j} | |z} ntj|} | |j jd} | |z} tj| | t'| ||} || z}|S)aMultidimensional uniform shift filter. The array is multiplied with the Fourier transform of a box of given size. Args: input (cupy.ndarray): The input array. size (float or sequence of float): The sigma of the box used for filtering. If a float, `size` is the same for all axes. If a sequence, `size` has to contain one value for each axis. n (int, optional): If `n` is negative (default), then the input is assumed to be the result of a complex fft. If `n` is larger than or equal to zero, the input is assumed to be the result of a real fft, and `n` gives the length of the array before transformation along the real transform direction. axis (int, optional): The axis of the real transform (only used when ``n > -1``). output (cupy.ndarray, optional): If given, the result of shifting the input is placed in this array. Returns: output (cupy.ndarray): The filtered output. rrrFr&r*r,)r rrr-rr.rr/r0r1rr r2r3r r4r5r6sincr#) rrr:r!rr sizesr<r>rs rfourier_uniformrDXs+. :D  / /F  )$ 5 5D  #D$ 7 7E 5&)))(UFL)A)ABBOT7 ::!a%%+gV[->???C 1HCC(""7++Cjj*j77 t  #3#Dr222#  Mrc^|j}t||d}tj||}t j||d}t j||tt||j D]\}\}} |dkr||kr:|dkr4tj | |j } | dtjz|z|z z} n4tj| } | dtjz|zz} tj| | t'| ||} || z}|S) aFMultidimensional Fourier shift filter. The array is multiplied with the Fourier transform of a shift operation. Args: input (cupy.ndarray): The input array. This should be in the Fourier domain. shift (float or sequence of float): The size of shift. If a float, `shift` is the same for all axes. If a sequence, `shift` has to contain one value for each axis. n (int, optional): If `n` is negative (default), then the input is assumed to be the result of a complex fft. If `n` is larger than or equal to zero, the input is assumed to be the result of a real fft, and `n` gives the length of the array before transformation along the real transform direction. axis (int, optional): The axis of the real transform (only used when ``n > -1``). output (cupy.ndarray, optional): If given, the result of shifting the input is placed in this array. Returns: output (cupy.ndarray): The shifted output (in the Fourier domain). T)rshiftrryr*r,)r rrr-rr.rr/r0r1rr r2r r7r8r4r5r9r#) rrFr:r!rr shiftsr<shiftkr>rs r fourier_shiftrIs90 :D T B B BF  )$ 5 5D  $UD' : :F 5&)))!*3vv|+D+D!E!E VW Q;;  ::!a%%+gV\:::C 3>F*Q. .CC(""7++Cux&01C ##Dr222#  Mrc|j}|dkrt|||||S|dkrtdt||}t j||}t j||d}tj ||d}tt||j D]\}\}} ||kr<|dkr6tj| |jj} | t"j|z|z z} n1tj| } | t"j|zz} | |jjd} | | z} t-| || } || z}tj|| |d kr!t1j|| |d z}||z}nC|dkr=tj|| ||tj|zz}|dz}||dzz}d |d |z<||z}|S)a Multidimensional ellipsoid Fourier filter. The array is multiplied with the fourier transform of a ellipsoid of given sizes. Args: input (cupy.ndarray): The input array. size (float or sequence of float): The size of the box used for filtering. If a float, `size` is the same for all axes. If a sequence, `size` has to contain one value for each axis. n (int, optional): If `n` is negative (default), then the input is assumed to be the result of a complex fft. If `n` is larger than or equal to zero, the input is assumed to be the result of a real fft, and `n` gives the length of the array before transformation along the real transform direction. axis (int, optional): The axis of the real transform (only used when ``n > -1``). output (cupy.ndarray, optional): If given, the result of shifting the input is placed in this array. Returns: output (cupy.ndarray): The filtered output. rz'Only 1d, 2d and 3d inputs are supportedrrrFr&r,r*g?)r)r rDNotImplementedErrorrrr-rr.rr/r0r1rr r2r3r r7r8r4r5r6r#sqrtrj1sincos) rrr:r!rr rCdistancer<r>rs rfourier_ellipsoidrSs0 :D qyyudAtV<<< axx""KLLL  / /F  )$ 5 5D  #D$ 7 7E 5&)))H(UFL)A)ABB " "OT7 ::!a%%+gV[->???C 58d?Q& &CC(""7++C 58d? "Cjj*j77 s #Dr222c>IhH%%%% qyy 8((((! (  v&&&&(TXh////! (a-F4$; eOF Mr)F)r r N)r7r r cupy._corercupyx.scipy.ndimager cupyx.scipyrrr#r@rDrIrSrrrXs %%%%%%&!!! 2222j////d....bBBBBBBr