U Kc:C@s.ddlZddlmZmZmZmZmZmZmZddl m Z ddl Z ddl m ZddlmZddlmZddlmZmZmZmZddlmZddd d d d d dddddddddddddgZede de de dfZddddhZeeeeeddd Ze j ee j d!d"d#Z!d`eeed%d&d'Z"eeed(feed(fed)d*d+Z#e$eeeeeeed,d-d.Z%e$eeeeeeeed/d0d1Z&e$eeeeeeed,d2d3Z'ee j(j)j*edaeeeeeed5d6dZ+ee j(j)j,edbeeeeeed5d7dZ-ee j(j)j.edceeeeeed5d8d Z/ee j(j)j0eddeeeeeed5d9dZ1ee j(j)j2edeeeeeeed5d:d Z3ee j(j)j4edfeeeeeed5d;dZ5Gdd?d@Z7eeedAdBdCZ8e$eeed(feed(feeedDdEdFZ9ee j(j)j:edgeeeeeeedGdHd Z;ee j(j)j<edheeeeeeedGdIdZ=ee j(j)j>edieeeeeeedGdJdZ?ee j(j)j@edjeeeeeeedGdKdZAGdLdMdMeZBe$eeeeeeBdNdOdPZCee j(j)jDedkeeeeeeedGdQdZEee j(j)jFedleeeeeeedGdRd ZGee j(j)jHedmeeeeeeedGdTdZIee j(j)jJedneeeeeeedGdUdZKee j(j)jLedoeeeeeeedGdVdZMee j(j)jNedpeeeeeeedGdWdZOee j(j)jPedqeeeeeeedGdXd ZQee j(j)jRedreeeeeeedGdYdZSeeeeedZd[d\ZTee j(j)jUdseeeed]d^dZVee j(j)jWdteeeed]d_dZXdS)uN)IterableList NamedTupleOptionalSequenceTupleUnion)Literal)register_decomposition)checkDimsType ShapeTypeTensorLikeType) out_wrapperfftfft2fftnhffthfft2hfftnrfftrfft2rfftnifftifft2ifftnihfftihfft2ihfftnirfftirfft2irfftnfftshift ifftshiftforwardbackwardortho)xnorm signal_numelr$returncsfttkfdddkr0|dt|S| rFdkpPdkpP|oPdk}|rb|d|S|S)z3Apply normalization to the un-normalized FFT resultcs dS)NzInvalid normalization mode: r+r(r+3/tmp/pip-unpacked-wheel-gikjz4vx/torch/_refs/fft.py/z_apply_norm..r&Nr%r$)r _NORM_VALUESmathsqrt)r'r(r)r$ normalizer+r,r- _apply_norm+s r5)dtyperequire_complexr*cCs*|jr |S|jst}|r&t|}|S)z@Helper to promote a dtype to one supported by the FFT primitives) is_complexZis_floating_pointtorchZget_default_dtypeutilsZcorresponding_complex_dtype)r6r7r+r+r-_promote_type_fft:s r;F)tr7r*cCs(|j}t||}||kr|St||S)zEHelper to promote a tensor to a dtype supported by the FFT primitives)r6r;primsZconvert_element_type)r<r7Zcur_typenew_typer+r+r-_maybe_promote_tensor_fftIs  r?.)r'dimssizesr*cCst|t|kstd}|j}dgt|d}tt|D]}||dkrNq<|||||krd}t|d||d}|||||||<|||||kr<|||d||}q<|rt||S|S)z Fixes the shape of x such that x.size(dims[i]) == sizes[i], either by zero-padding, or by slicing x starting from 0. FrTr0)lenAssertionErrorshaperangeZnarrowr9Zconstant_pad_nd)r'r@rAZ must_copyZx_sizesZ pad_amountiZpad_idxr+r+r-_resize_fft_inputTs rI) func_nameinputndimr(r$r*c st|dd}t|j|f}dk r(nd|j|d}t|dkfdddk rpt|||ddfd}|r~t|}t j |||d }t ||||d S) zBCommon code for performing any complex to real FFT (irfft or hfft)Tr7NrBr0cs ddSNzInvalid number of data points (z ) specifiedr+r+rLr+r-r.zr/z_fft_c2r..)r@rArM last_dim_sizer(r)r$) r?r:canonicalize_dimndimrFr rIr9conjr=fft_c2rr5) rJrKrLrMr(r$r@rRoutputr+rPr-_fft_c2rns  rY)rJrKrLrMr(r$onesidedr*c s~tjj fddttj|f}|dk rHt||ftj ||d}t ||j ||}|rt|St |S)zBCommon code for performing any real to complex FFT (rfft or ihfft)csdjS)Nz0 expects a floating point input tensor, but got r6r+rJrKr+r-r.r/z_fft_r2c..NrMrZ)r r6r8r?r:rTrUrIr=fft_r2cr5rFr9rV) rJrKrLrMr(r$rZr@retr+r\r-_fft_r2cs  r`csbtjjfddtj|f}|dk r>t||ftj||d}t ||j ||S)zCCommon code for performing any complex to complex FFT (fft or ifft)csdjSNz) expects a complex input tensor, but got r[r+r\r+r-r.r/z_fft_c2c..NrMr$) r r6r8r:rTrUrIr=fft_c2cr5rF)rJrKrLrMr(r$r@r_r+r\r-_fft_c2cs  rdrC)rKrLrMr(r*c Cs6|jjrtd||||ddStd||||dddSdS)NrTr$Fr$rZr6r8rdr`rKrLrMr(r+r+r-rsc Cs6|jjrtd||||ddStd||||dddSdS)NrFrerfrgrhr+r+r-rsc Cstd||||dddS)NrTrfr`rhr+r+r-rscCstd||||ddS)NrFrerYrhr+r+r-rscCstd||||ddS)NrTrerjrhr+r+r-rsc Cstd||||dddS)NrFTrfrirhr+r+r-rsc@s.eZdZUeedfed<eedfed<dS) _ShapeAndDims.rFr@N__name__ __module__ __qualname__rint__annotations__r+r+r+r-rks rk)rKrFrMr*csB|j|j|dk rNt|ts$|f}t|}ttt|t|kdd|dk rt|tsf|f}t|dkp~t|t|kddt|tkfdd|dkrt t }t fddt ||D}n4|dkrt t }t }nt fdd|D}|D]td kfd dqt ||d S) zTConvert the shape and dim arguments into a canonical form where neither are optionalNcSsdS)NzFFT dims must be uniquer+r+r+r+r-r.r/z6_canonicalize_fft_shape_and_dim_args..cSsdS)Nz=When given, dim and shape arguments must have the same lengthr+r+r+r+r-r.r/csdddS)NzGot shape with z" values but input tensor only has z dimensions.r+r+) input_dimtransform_ndimr+r-r.r/c3s&|]\}}|dkr|n|VqdS)rCNr+).0sd input_sizesr+r- 'sz7_canonicalize_fft_shape_and_dim_args..c3s|]}|VqdSNr+rtrvrwr+r-ry0srcs ddSrOr+r+rPr+r-r.3r/)rFr@) rUrF isinstancerr:Zcanonicalize_dimsr rDsettuplerGziprk)rKrFrMZret_dimsZ ret_shaper+)rrrxrLrsr-$_canonicalize_fft_shape_and_dim_argss>        r)xsr*cCsd}|D] }||9}q|S)zCompute product of a listr0r+)rprodr'r+r+r-_prod8s r) function_namerKrFrMr(r$r*csHtjjfddt||}tj|||d}t||t||dS)zECommon code for n-dimensional complex to complex FFTs (fftn or ifftn)csdjSrar[r+rrKr+r-r.Kr/z_fftn_c2c..rbrS)r r6r8rIr=rcr5r)rrKrFrMr(r$r'rXr+rr- _fftn_c2c@s   r)rKrurMr(r*cCs0t|||\}}t|dd}td||||ddS)NTrNrrerr?rrKrurMr(rFr'r+r+r-rSs cCs0t|||\}}t|dd}td||||ddS)NTrNrFrerrr+r+r-r`s csdtjj fddt||\}}tddt||tj|dd}t||t |ddS)Ncs djS)Nz2rfftn expects a real-valued input tensor, but got r[r+rKr+r-r.wr/zrfftn..FrNTr]rS) r r6r8rr?rIr=r^r5r)rKrurMr(rFoutr+rr-rms   cstjj fddt||\}}tt|dkddtddt||tj|dddd }t|d krt |||ddd }t |St |}tj ||dddd }t ||t |dd S) Ncs djS)Nz3ihfftn expects a real-valued input tensor, but got r[r+rr+r-r.r/zihfftn..rcSsdS)Nz'ihfftn must transform at least one axisr+r+r+r+r-r.r/FrNrCTr]r0rSrb)r r6r8rrDr?rIr=r^r5rV conj_physicalrcr)rKrurMr(rFtmpr+rr-rs      c@s6eZdZUeedfed<eedfed<eed<dS)_CanonicalizeC2rReturn.rFrMrRNrlr+r+r+r-rs r)fnamerKrurMr*cst|||\}}tt|dkfdd|dks>|ddkrVd|j|ddn|dtdkfddt|}dd|d<tt||d S) zCanonicalize shape and dim arguments for n-dimensional c2r transforms, as well as calculating the last_dim_size which is shape[dim[-1]] for the outputrcs dS)Nz! must transform at least one axisr+r+)rr+r-r.r/z:_canonicalize_fft_c2r_shape_and_dim_args..NrCrBr0cs ddSrOr+r+)rRr+r-r.r/)rFrMrR)rr rDrFlistrr~)rrKrurMrFZ shape_listr+)rrRr-(_canonicalize_fft_c2r_shape_and_dim_argss  rcs^td|||\}}}t|dd}t|||}tj|||dt|tfdd|DddS) Nr!TrNrQc3s|]}j|VqdSrzrFr{rr+r-ryszirfftn..Fre)rr?rIr=rWr5r)rKrurMr(rFrRr+rr-r!s   cCstd|||\}}}t|dd}t|||}t|dkrPtj||ddddn|}t||t|dddd}t|}tj ||dd|d}t|||ddS) NrTrNr0rCrbrerQ) rr?rIrDr=rcr5rrrW)rKrurMr(rFrRrrr+r+r-rs   ( rCcCstjj||||dSN)rurMr()r9rrrKrurMr(r+r+r-rscCstjj||||dSr)r9rrrr+r+r-rscCstjj||||dSr)r9rrrr+r+r-rscCstjj||||dSr)r9rr!rr+r+r-r scCstjj||||dSr)r9rrrr+r+r-rscCstjj||||dSr)r9rrrr+r+r-rs)rMr'r*cCs2|dkrtt|jSt|ts&|gSt|SdS)zIConvert Optional[DimsType] to a simple list, defaulting to all dimensionsN)rrGrUr|r)rMr'r+r+r-_default_alldims&s  r)rKrMr*cs*t|}fdd|D}t||S)Ncsg|]}j|dqS)rBrr{rr+r- 3szfftshift..rr9ZrollrKrMr@shiftr+rr-r"0s cs*t|}fdd|D}t||S)Ncsg|]}j|ddqS)r0rBrr{rr+r-r:szifftshift..rrr+rr-r#7s )F)NrCN)NrCN)NrCN)NrCN)NrCN)NrCN)NNN)NNN)NNN)NNN)NNN)NNN)NrN)NrN)NrN)NrN)NrN)NrN)N)N)Yr2typingrrrrrrrZtyping_extensionsr r9Z torch._primsZ_primsr=Ztorch._prims_commonZ _prims_commonr:Z torch._decompr r r r rZtorch._prims_common.wrappersr__all__ZNormTyper1rpboolr5r6r;r?rIstrrYr`rdopsZatenZfft_fftrZfft_ifftrZfft_rfftrZ fft_irfftrZfft_hfftrZ fft_ihfftrrkrrrZfft_fftnrZ fft_ifftnrZ fft_rfftnrZ fft_ihfftnrrrZ fft_irfftnr!Z fft_hfftnrZfft_fft2rZ fft_ifft2rZ fft_rfft2rZ fft_irfft2r Z fft_hfft2rZ fft_ihfft2rrZ fft_fftshiftr"Z fft_ifftshiftr#r+r+r+r-sz$                     6