o i%@s~ddlZddlZddlmZddlmZddlmZddlmZdddZ d d Z dd d Z dddZ dddZ dddZdS)N)_core)internal)_util)specialFcCstjtjg}|s|tjtjg7}|dur0|j|vr$tj|j|jd}|Stj|j|dd}|St|turH||vr>t dtj|j|d}|S|j|jkrRt d|S)Ndtypezoutput type not supportedzoutput shape not correct) cupy complex64 complex128float32float64remptyshapetype RuntimeError)outputinput complex_onlytypesrY/home/ubuntu/veenaModal/venv/lib/python3.10/site-packages/cupyx/scipy/ndimage/_fourier.py_get_output_fourier s"     rcCs*d||jfd||d}||S)z>Promote a 1d array to ndim with non-singleton size along axis.)r)sizereshape)arrndimaxisnd_shaperrr _reshape_nds  r rc Cs|j}t||}t||}t||d}t||tt ||j D]S\}\}} ||kr?|dkr?t j | |j jd} | |} nt j| } | j|j jdd} | | 9} ||d} | dtjtj| 9} t j| | dt| ||d } || 9}q#|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. sigmarrFcopyoutrr)rrr_normalize_axis_indexr_fix_sequence_argrelementwise_copy enumerateziprr arangerealrfftfftfreqastypenumpypiexpr ) rr!nrrrsigmasaxsigmakax_sizerscalerrrfourier_gaussian#s$       r<c Cs|j}t||}t||}t||d}t||tt ||j D]A\}\}}||kr?|dkr?t j ||j jd} | |} nt j|} | j|j jdd} | |9} t j| | dt| ||d} || 9}q#|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()rrrr)rr*rr+r,r-rr r.r/rr0r1r2sincr ) rrr6rrrsizesr8r:rrrrfourier_uniformXs       r?c Cs|j}t||dd}t||}t||d}t||tt ||j D]H\}\}} |dkr0q%||krL|dkrLt j | |j d} | dtj||9} nt j| } | dtj|} t j| | dt| ||d} || 9}q%|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()rrrr)rr*rr+r,r-rr r.rr3r4r0r1r5r ) rr@r6rrrshiftsr8shiftkr:rrrr fourier_shifts"    rCc Cs|j}|dkrt|||||S|dkrtdt||}t||}t||d}t ||d}t t ||j D]F\}\}} ||krZ|dkrZt j| |jjd} | tj||9} n t j| } | tj|9} | j|jjdd} | | 9} t| ||d } || }q9t j||d |d krtj||d |d 9}||}n|dkrt j||d ||t |8}|d9}||d}d |d |<||9}|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)rr?NotImplementedErrorrrr)rr*rr+r,r-rr r.r/rr3r4r0r1r2r sqrtrj1sincos) rrr6rrrr>distancer8r:rrrrfourier_ellipsoidsB        rL)F)rrN)r3r r cupy._corercupyx.scipy.ndimager cupyx.scipyrrr r<r?rCrLrrrrs       5 21