
    sgR$                        d dl Zd dlmZmZ d dlZd dlZd dlmZm	Z	m
Z
mZmZmZmZmZ d dlmZ d dlmZ d dlmZ d dlmZmZ eej0                  j3                  d      gZej0                  j6                  Z ej8                  d      Z ed	
      ej0                  j=                  dee	feefg      ej0                  j=                  dg d      ej0                  j=                  dg d      ej0                  j=                  dd dg      ej0                  j=                  dg d      ej0                  j=                  ddd	g      d                                                  Z ed	dg      ej0                  j=                  dee	feefg      ej0                  j=                  dg d      ej0                  j=                  dej@                  ejB                  ejD                  ejF                  ejH                  g      ej0                  j=                  dd dg      ej0                  j=                  dg d      ej0                  j=                  dd	dg      d                                                  Z% ed	
      ej0                  j=                  de
efeefg      ej0                  j=                  dg d      ej0                  j=                  dg d      ej0                  j=                  dg d      ej0                  j=                  ddd	g      d                                           Z& ed	dg      ej0                  j=                  de
efeefg      ej0                  j=                  dg d      ej0                  j=                  dg d      ej0                  j=                  dej@                  ejB                  ejD                  ejF                  ejH                  g      ej0                  j=                  dg d      ej0                  j=                  ddd	g      d                                                   Z' ed	
      ej0                  j=                  d!g d"      ej0                  j=                  dg d      ej0                  j=                  dg d      d#                             Z( ed	
      ej0                  j=                  d!eee
eg      ej0                  j=                  dg d      d$                      Z) ed	
      ej0                  j=                  dg d%      ej0                  j=                  d&ed'fedfed'fg      d(                      Z* ed)d*gd	+      ej0                  j=                  dg d%      d,               Z+ ed)d*gd	+      ej0                  j=                  dg d%      ej0                  j=                  d!eeg      d-                      Z, ed)d*gd	+      ej0                  j=                  dg d%      ej0                  j=                  d!eeg      d.                      Z- ed	d/g      ej0                  j=                  d!ee	e
eeeeeg      d0               Z.y)1    N)assert_allcloseassert_array_equal)dctidctdctnidctndstidstdstnidstn)fftpack)array_api_compatible)copyxp_assert_closeskip_xp_backends   T)cpu_onlyzforward, backwardtype)   r         n)r   r   r      
      axisr   norm)NbackwardorthoforwardorthogonalizeFc                 l   |j                  t        j                  j                  ||            } | |||||      }	 ||	||||      }
t	        |
|       dgdz  }d||<   |j                  t        j
                  t        j                   |	      |d            } |||||||      }t	        ||       y )N)r   r   r!   r   r   r   r   r   edgemoder!   )asarraynprandomrandr   pad)r    r   r   r   r   r   r!   xpxyzr-   y2z2s                 W/var/www/html/venv/lib/python3.12/site-packages/scipy/fft/tests/test_real_transforms.pytest_identity_1dr5      s     	

299>>!Q'(A4d]KADt$mLAAq(Q,CCI	BFF2::a=#F;	<B	"dAtT	GBB    z/`overwrite_x` only supported for NumPy backend.)np_onlyreasonsdtypeoverwrite_xc                 H   t         j                  j                  dd      j                  |      }|j	                         } | |||||      }	|	j	                         }
 ||	||||      }|s(t        ||dd       t        ||       t        |	|
       y t        ||dd       y )N      )r   r   r:   ư>rtolatol)r*   r+   r,   astyper   r   r   )r    r   r   r9   r   r   r:   r/   x_origr0   y_origr1   s               r4   test_identity_1d_overwriterE   -   s     			q!##E*AVVXF4d;GAVVXFDt$KHA14d31f%1f%648r6   zshape, axes)	)r   r   r   )rF   r   )rF   N)rF   r   r   )r      N)rH   rG   )r   r      N)rJ   r   )rJ   )r   r   c                 .   |j                  t        j                  j                  |            }|t        j                  ||      } | |||||      }	 ||	||||      }
t	        |
|       |dg|j
                  z  }nBt        |t              rdg|j
                  z  }d||<   ndg|j
                  z  }|D ]  }d||<   	 |j                  t        j                  t        j                   |	      |d            } |||||||      }t	        ||       y )N)axesr   r!   r$   r#   r%   r&   r(   )	r)   r*   r+   taker   ndim
isinstanceintr-   )r    r   r   shaperM   r   r!   r.   r/   r0   r1   r-   ar2   r3   s                  r4   test_identity_ndrT   G   s   * 	

299##E*+At$4d]KADt$mLAAq|h	D#	hD	h 	ACF	 
BFF2::a=#F;	<B	"dE4]	KBBr6   ))r   r   r   )rU   r   )rU   Nc                 r   t         j                  j                  |      j                  |      }|j                         }	|t        j                  ||      } | ||||      }
|
j                         } ||
|||      }|rt        ||	dd       y t        ||dd       t        ||	       t        |
|       y )N)rM   r   r>   r?   )r*   r+   rB   r   rN   r   r   )r    r   r   rR   rM   r9   r   r:   r/   rC   r0   rD   r1   s                r4   test_identity_nd_overwriterW   v   s    $ 			&&u-AVVXFt$4d.AVVXFDt$/A64814d31f%1f%r6   func)r   r	   r   r   c                     t         j                  j                  dd      }|j                   t	        t
        |       |||            }|j                  |      } t	        t        |       |||      }t        ||       y )Nr=   r   )r   )r*   r+   r,   r)   getattrr   fftr   )rX   r   r   r.   r/   fftpack_resfft_ress          r4   test_fftpack_equivaliencer^      sg    
 			q"A**3WWd3At$GHK


1A gc4 Dt4GG[)r6   c                     |j                  t        j                  j                  d            }dD ](  \  }} | ||||      } | |||      }t	        ||       * y )Nd   ))r    F)r   F)r   Tr   r   r!   )r   r   r)   r*   r+   r,   r   )rX   r   r.   r/   r   r   rS   bs           r4   test_orthogonalize_defaultrd      s^    
 	

299>>#&'A e
 D>D)1r6   )r   r   r    z
func, typer   c                     |j                  t        j                  j                  d            } | |||d      } | |||d      }t	        ||       y )Nr`   Tra   Frb   )rX   r   r   r.   r/   y1r2   s          r4   test_orthogonalize_nooprg      sI     	

299>>#&'A	adT	:B	adU	;BBr6   z	jax.numpyz)jax arrays do not support item assignment)r8   r   c                 V   |j                  t        j                  j                  d            }t	        ||      }|dxx   t
        z  cc<   |dxx   t
        z  cc<   t        |d| d      }t        |d| d      }|dxx   t
        z  cc<   |dxx   t
        z  cc<   t        ||       y )	Nr`   r.   r   r   Tra   F)r)   r*   r+   r,   r   SQRT_2r   r   )r   r.   r/   x2rf   r2   s         r4   test_orthogonalize_dct1rm      s    
 	

299>>#&'A	aBBqEVOErFfF	QQT	6B	Rad%	8BqEVOErFfFBr6   c                     |j                  t        j                  j                  d            } | |d|d      } | |d|d      }|| t        k(  rdndxx   t
        z  cc<   t        ||       y )Nr`   r   Tra   Fr   rj   )r)   r*   r+   r,   r   rk   r   )rX   r   r.   r/   rf   r2   s         r4   test_orthogonalize_dcst2ro      sa     	

299>>#&'A	aad$	7B	aad%	8BDCKqR F* Br6   c                     |j                  t        j                  j                  d            }t	        ||      }|| t
        k(  rdndxx   t        z  cc<    | |d|d      } | |d|d      }t        ||       y )	Nr`   ri   r   rj   r   Tra   F)r)   r*   r+   r,   r   r   rk   r   )rX   r   r.   r/   rl   rf   r2   s          r4   test_orthogonalize_dcst3rq      sl     	

299>>#&'A	aBBDCKqR F* 	aad$	7B	bqt5	9BBr6   z,array-likes only supported for NumPy backendc                     ddgddggddgddggddgddggg}t         ||       || j                  |                   y )Ng      ?)r   r)   )r.   rX   r/   s      r4   test_array_likers      sX     *sCj	!*sCj	!*sCj	!	#A DGT"**Q-01r6   )/numpyr*   numpy.testingr   r   pytestmath	scipy.fftr   r   r   r   r	   r
   r   r   r[   scipyr   scipy.conftestr   scipy._lib._array_apir   r   markusefixtures
pytestmarkr   sqrtrk   parametrizer5   float16float32float64	complex64
complex128rE   rT   rW   r^   rd   rg   rm   ro   rq   rs    r6   r4   <module>r      s;    =   D D D   / 7"FKK$;$;<N$OP
;;// 	1 4 ,T{S$K.HI.23!Q(!GH5$-8 9 I ) 4 / J !  $MNP,T{S$K.HI.2::rzz2::#%<<#@ A!Q(!GHu69 7 I )A / JP9" 4 ,e}tUm.LM.
 !GH5$-8 9 I / N !"< $MNP,e}tUm.LM. 2::rzz2::#%<<#@ A!GH6& 7 IA / NP&* 4 !?@.!GH* I / A !* 4 #sD$!78. / 9 ! 4 !AB!HsAha(" ## C ! +FG! !AB C! +FG! !AB#s, - C!
 +FG! !AB#s, - C!
 $IJL#tT5#tT5!QR2 SL2r6   