U
    mh]                     @   s  d dl mZmZ d dlZd dlmZmZ d dlZd dlm	Z
 d dlmZmZmZmZmZmZmZmZ eeeedZdd ed	D Zd
d ed	D ZeeeedZeeeedZed Zdd Zdd Zdd Zdd Z dd Z!dd Z"dddZ#dddZ$dddZ%dd d!Z&G d"d# d#Z'G d$d% d%Z(G d&d' d'e(Z)G d(d) d)e(Z*G d*d+ d+e(Z+G d,d- d-e(Z,G d.d/ d/e)Z-G d0d1 d1e)Z.G d2d3 d3e)Z/G d4d5 d5e*Z0G d6d7 d7e*Z1G d8d9 d9e*Z2G d:d; d;e+Z3G d<d= d=e+Z4G d>d? d?e+Z5G d@dA dAe,Z6G dBdC dCe,Z7G dDdE dEe,Z8G dFdG dGZ9G dHdI dIe9Z:G dJdK dKe9Z;G dLdM dMe9Z<G dNdO dOe9Z=G dPdQ dQe9Z>G dRdS dSe9Z?G dTdU dUe9Z@G dVdW dWe9ZAG dXdY dYe9ZBG dZd[ d[e9ZCG d\d] d]e9ZDG d^d_ d_e9ZEG d`da daZFG dbdc dceFZGG ddde deeFZHG dfdg dgeGZIG dhdi dieGZJG djdk dkeGZKG dldm dmeFZLG dndo doeFZMG dpdq dqeFZNG drds dseFZOG dtdu dueFZPG dvdw dweFZQG dxdy dyeHZRG dzd{ d{eHZSG d|d} d}eHZTG d~d dZUG dd deUZVG dd deUZWG dd deUZXG dd deUZYG dd deUZZG dd deUZ[G dd deUZ\G dd deUZ]G dd deUZ^G dd deUZ_G dd deUZ`G dd deUZaG dd dZbG dd dZcdS )    )joindirnameN)assert_array_almost_equalassert_equal)raises)dctidctdstidstdctnidctndstnidstnztest.npzc                 C   s   g | ]}t d |  qS )zx%dMDATA.0i r   L/tmp/pip-unpacked-wheel-9oz00b9h/scipy/fftpack/tests/test_real_transforms.py
<listcomp>   s     r      c                 C   s   g | ]}t d |  qS )zy%dr   r   r   r   r   r      s     zfftw_double_ref.npzzfftw_single_ref.npzZsizesc                 C   sl   t d|d ||}t t j|}|t jkr6t}n|t jkrFt}nt |d| |f  |}|||fS )Nr      z	dct_%d_%d	npZlinspaceastyperesult_typefloat32doubleFFTWDATA_DOUBLEFFTWDATA_SINGLE
ValueErrortypesizedtxdatayr   r   r   fftw_dct_ref   s    

r)   c                 C   sl   t d|d ||}t t j|}|t jkr6t}n|t jkrFt}nt |d| |f  |}|||fS )Nr   r   z	dst_%d_%dr   r"   r   r   r   fftw_dst_ref'   s    

r*   c                 K   s   t j| dd} t| jd D ](}t| |ddf f|| |ddf< qt| jd D ](}t| dd|f f|| dd|f< qT| S )z,Calculate reference values for testing dct2.Tcopyr   Nr   )r   arrayrangeshaper   r&   kwargsrowcolr   r   r   
dct_2d_ref4   s    &&r4   c                 K   s   t j| dd} t| jd D ](}t| |ddf f|| |ddf< qt| jd D ](}t| dd|f f|| dd|f< qT| S )z-Calculate reference values for testing idct2.Tr+   r   Nr   )r   r-   r.   r/   r   r0   r   r   r   idct_2d_ref>   s    &&r5   c                 K   s   t j| dd} t| jd D ](}t| |ddf f|| |ddf< qt| jd D ](}t| dd|f f|| dd|f< qT| S )z,Calculate reference values for testing dst2.Tr+   r   Nr   )r   r-   r.   r/   r	   r0   r   r   r   
dst_2d_refH   s    &&r6   c                 K   s   t j| dd} t| jd D ](}t| |ddf f|| |ddf< qt| jd D ](}t| dd|f f|| dd|f< qT| S )z-Calculate reference values for testing idst2.Tr+   r   Nr   )r   r-   r.   r/   r
   r0   r   r   r   idst_2d_refR   s    &&r7   c           	   
   C   s4  t j| dd} t| }|d }t |}d\}}|dkrTt d| }t d| }t|D ]}td|d D ]4}||  || |  t t j| | |  7  < qn||  || d  7  < ||  || |d   |d	 dkrdnd
 7  < q\|dkr0|d  dt d	 9  < ||d   dt d	 9  < |S )z/Calculate textbook definition version of DCT-I.Tr+   r   )r      ortho      ?       @r   r8   )r   r-   lenzerossqrtr.   cospi)	r&   normNMr(   Zm0mknr   r   r   
naive_dct1\   s"    
22
rH   c                 C   s   t j| dd} t| }|d }t |}t|D ]J}t|D ]<}||  d| |  t t j|d  |d  |  7  < q<q0|dkr|t d| 9 }|S )z0Calculate textbook definition version  of DST-I.Tr+   r   r8   r:   r9         ?r   r-   r=   r>   r.   sinrA   r?   )r&   rB   rC   rD   r(   rF   rG   r   r   r   
naive_dst1q   s    
<rL   c                 C   s   t j| dd} t| }t |}t|D ]F}t|D ]8}||  | | t t j|d  |d  |  7  < q4q(|dkr|t d| 9 }n|d9 }|S )z0Calculate textbook definition version of DCT-IV.Tr+   rI   r9   r;   r8   )r   r-   r=   r>   r.   r@   rA   r?   r&   rB   rC   r(   rF   rG   r   r   r   
naive_dct4   s    
8rN   c                 C   s   t j| dd} t| }t |}t|D ]F}t|D ]8}||  | | t t j|d  |d  |  7  < q4q(|dkr|t d| 9 }n|d9 }|S )z0Calculate textbook definition version of DST-IV.Tr+   rI   r9   r;   r8   rJ   rM   r   r   r   
naive_dst4   s    
8rO   c                   @   s<   e Zd Zdd Zdd Zdd Zdd Zd	d
 Zdd ZdS )TestComplexc                 C   s8   t dtjdtjd }dt td }t|| d S )N              ?   dtype)r   r   arange	complex64r   selfr(   r&   r   r   r   test_dct_complex64   s    zTestComplex.test_dct_complex64c                 C   s2   t tdd }dt td }t|| d S NrR   rQ   )r   r   rU   r   rW   r   r   r   test_dct_complex   s    zTestComplex.test_dct_complexc                 C   s2   t tdd }dt td }t|| d S rZ   )r   r   rU   r   rW   r   r   r   test_idct_complex   s    zTestComplex.test_idct_complexc                 C   s8   t tjdtjdd }dt td }t|| d S )NrR   rS   rQ   )r	   r   rU   rV   r   rW   r   r   r   test_dst_complex64   s    zTestComplex.test_dst_complex64c                 C   s2   t tdd }dt td }t|| d S rZ   )r	   r   rU   r   rW   r   r   r   test_dst_complex   s    zTestComplex.test_dst_complexc                 C   s2   t tdd }dt td }t|| d S rZ   )r
   r   rU   r   rW   r   r   r   test_idst_complex   s    zTestComplex.test_idst_complexN)	__name__
__module____qualname__rY   r[   r\   r]   r^   r_   r   r   r   r   rP      s   rP   c                   @   s$   e Zd Zdd Zdd Zdd ZdS )_TestDCTBasec                 C   s   d | _ d| _d | _d S N   rdtdecr#   rX   r   r   r   setup_method   s    z_TestDCTBase.setup_methodc                 C   sh   t D ]^}t| j|| j\}}}t|| jd}t|j| t|t	| |t	| | j
d| d qd S Nr#   Size %d faileddecimalerr_msg)FFTWDATA_SIZESr)   r#   rg   r   r   rT   r   r   maxrh   )rX   r   r&   yrr%   r(   r   r   r   test_definition   s    z_TestDCTBase.test_definitionc                 C   s   d}dD ]}t j||}t|| jd}t|D ]&}t|| t|| | jd| jd q0|j}t|d| jd}t|D ]6}t|d d |f t|d d |f | jd| jd qvqd S )Nr8   )   r   	          @   rl   ro   r   )axisr#   )	r   randomrandnr   r#   r.   r   rh   T)rX   ntr   r&   r(   jr   r   r   	test_axis   s    (z_TestDCTBase.test_axisN)r`   ra   rb   rj   rt   r   r   r   r   r   rc      s   rc   c                   @   s   e Zd Zdd ZdS )_TestDCTIBasec                 C   sx   t t j| j}tD ]^}t j|| jd}t|ddd}t|dd}t|j	| t
|t | |t | | jd qd S NrS   r9   r   rB   r#   rB   rz   )r   r   r   rg   Xr-   r   rH   r   rT   r   rr   rh   rX   r%   xrr&   r(   y2r   r   r   test_definition_ortho   s    z#_TestDCTIBase.test_definition_orthoNr`   ra   rb   r   r   r   r   r   r      s   r   c                   @   s   e Zd Zdd ZdS )_TestDCTIIBasec                 C   s`   t t j| j}tttD ]@\}}t j||d}t|ddd}t	|j
| t||| jd qd S )NrS   r9   r8   r   rz   )r   r   r   rg   zipr   Yr-   r   r   rT   r   rh   )rX   r%   r   rs   r&   r(   r   r   r   test_definition_matlab   s    z%_TestDCTIIBase.test_definition_matlabN)r`   ra   rb   r   r   r   r   r   r      s   r   c                   @   s   e Zd Zdd ZdS )_TestDCTIIIBasec                 C   sf   t t j| j}tD ]L}t j|| jd}t|ddd}t|ddd}t|j| t	||| j
d qd S )NrS   r9   r8   r      rz   )r   r   r   rg   r   r-   r   r   rT   r   rh   )rX   r%   r   r&   r(   xir   r   r   r      s    z%_TestDCTIIIBase.test_definition_orthoNr   r   r   r   r   r      s   r   c                   @   s   e Zd Zdd ZdS )_TestDCTIVBasec                 C   sx   t t j| j}tD ]^}t j|| jd}t|ddd}t|dd}t|j	| t
|t | |t | | jd qd S NrS   r9      r   r   rz   )r   r   r   rg   r   r-   r   rN   r   rT   r   rr   rh   r   r   r   r   r     s    z$_TestDCTIVBase.test_definition_orthoNr   r   r   r   r   r      s   r   c                   @   s   e Zd Zdd ZdS )TestDCTIDoublec                 C   s   t j| _d| _d| _d S N
   r   r   r   rg   rh   r#   ri   r   r   r   rj     s    zTestDCTIDouble.setup_methodNr`   ra   rb   rj   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )TestDCTIFloatc                 C   s   t j| _d| _d| _d S Nr   r   r   r   rg   rh   r#   ri   r   r   r   rj     s    zTestDCTIFloat.setup_methodNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )TestDCTIIntc                 C   s   t | _d| _d| _d S NrR   r   intrg   rh   r#   ri   r   r   r   rj     s    zTestDCTIInt.setup_methodNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )TestDCTIIDoublec                 C   s   t j| _d| _d| _d S Nr   r8   r   ri   r   r   r   rj   "  s    zTestDCTIIDouble.setup_methodNr   r   r   r   r   r   !  s   r   c                   @   s   e Zd Zdd ZdS )TestDCTIIFloatc                 C   s   t j| _d| _d| _d S NrR   r8   r   ri   r   r   r   rj   )  s    zTestDCTIIFloat.setup_methodNr   r   r   r   r   r   (  s   r   c                   @   s   e Zd Zdd ZdS )TestDCTIIIntc                 C   s   t | _d| _d| _d S r   r   ri   r   r   r   rj   0  s    zTestDCTIIInt.setup_methodNr   r   r   r   r   r   /  s   r   c                   @   s   e Zd Zdd ZdS )TestDCTIIIDoublec                 C   s   t j| _d| _d| _d S Nre   r   r   ri   r   r   r   rj   7  s    zTestDCTIIIDouble.setup_methodNr   r   r   r   r   r   6  s   r   c                   @   s   e Zd Zdd ZdS )TestDCTIIIFloatc                 C   s   t j| _d| _d| _d S NrR   r   r   ri   r   r   r   rj   >  s    zTestDCTIIIFloat.setup_methodNr   r   r   r   r   r   =  s   r   c                   @   s   e Zd Zdd ZdS )TestDCTIIIIntc                 C   s   t | _d| _d| _d S r   r   ri   r   r   r   rj   E  s    zTestDCTIIIInt.setup_methodNr   r   r   r   r   r   D  s   r   c                   @   s   e Zd Zdd ZdS )TestDCTIVDoublec                 C   s   t j| _d| _d| _d S )N   r   r   ri   r   r   r   rj   L  s    zTestDCTIVDouble.setup_methodNr   r   r   r   r   r   K  s   r   c                   @   s   e Zd Zdd ZdS )TestDCTIVFloatc                 C   s   t j| _d| _d| _d S r   r   ri   r   r   r   rj   S  s    zTestDCTIVFloat.setup_methodNr   r   r   r   r   r   R  s   r   c                   @   s   e Zd Zdd ZdS )TestDCTIVIntc                 C   s   t | _d| _d| _d S r   r   ri   r   r   r   rj   Z  s    zTestDCTIVInt.setup_methodNr   r   r   r   r   r   Y  s   r   c                   @   s   e Zd Zdd Zdd ZdS )_TestIDCTBasec                 C   s   d | _ d| _d | _d S rd   rf   ri   r   r   r   rj   a  s    z_TestIDCTBase.setup_methodc                 C   s   t D ]}t| j|| j\}}}t|| jd}| jdkrH|d|d   }n|d|  }t|j| t|t	| |t	| | j
d| d qd S Nrl   r   r8   rm   rn   )rq   r)   r#   rg   r   r   rT   r   r   rr   rh   rX   r   r   rs   r%   r&   r   r   r   rt   f  s    
z_TestIDCTBase.test_definitionNr`   ra   rb   rj   rt   r   r   r   r   r   `  s   r   c                   @   s   e Zd Zdd ZdS )TestIDCTIDoublec                 C   s   t j| _d| _d| _d S r   r   ri   r   r   r   rj   x  s    zTestIDCTIDouble.setup_methodNr   r   r   r   r   r   w  s   r   c                   @   s   e Zd Zdd ZdS )TestIDCTIFloatc                 C   s   t j| _d| _d| _d S r   r   ri   r   r   r   rj     s    zTestIDCTIFloat.setup_methodNr   r   r   r   r   r   ~  s   r   c                   @   s   e Zd Zdd ZdS )TestIDCTIIntc                 C   s   t | _d| _d| _d S r   r   ri   r   r   r   rj     s    zTestIDCTIInt.setup_methodNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )TestIDCTIIDoublec                 C   s   t j| _d| _d| _d S r   r   ri   r   r   r   rj     s    zTestIDCTIIDouble.setup_methodNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )TestIDCTIIFloatc                 C   s   t j| _d| _d| _d S r   r   ri   r   r   r   rj     s    zTestIDCTIIFloat.setup_methodNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )TestIDCTIIIntc                 C   s   t | _d| _d| _d S r   r   ri   r   r   r   rj     s    zTestIDCTIIInt.setup_methodNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )TestIDCTIIIDoublec                 C   s   t j| _d| _d| _d S r   r   ri   r   r   r   rj     s    zTestIDCTIIIDouble.setup_methodNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )TestIDCTIIIFloatc                 C   s   t j| _d| _d| _d S r   r   ri   r   r   r   rj     s    zTestIDCTIIIFloat.setup_methodNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )TestIDCTIIIIntc                 C   s   t | _d| _d| _d S r   r   ri   r   r   r   rj     s    zTestIDCTIIIInt.setup_methodNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )TestIDCTIVDoublec                 C   s   t j| _d| _d| _d S Nr   r   r   ri   r   r   r   rj     s    zTestIDCTIVDouble.setup_methodNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )TestIDCTIVFloatc                 C   s   t j| _d| _d| _d S NrR   r   r   ri   r   r   r   rj     s    zTestIDCTIVFloat.setup_methodNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )TestIDCTIVIntc                 C   s   t | _d| _d| _d S r   r   ri   r   r   r   rj     s    zTestIDCTIVInt.setup_methodNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd Zdd ZdS )_TestDSTBasec                 C   s   d | _ d | _d | _d S Nrf   ri   r   r   r   rj     s    z_TestDSTBase.setup_methodc                 C   sh   t D ]^}t| j|| j\}}}t|| jd}t|j| t|t	| |t	| | j
d| d qd S rk   )rq   r*   r#   rg   r	   r   rT   r   r   rr   rh   )rX   r   r   rs   r%   r(   r   r   r   rt     s    z_TestDSTBase.test_definitionNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )_TestDSTIBasec                 C   sx   t t j| j}tD ]^}t j|| jd}t|ddd}t|dd}t|j	| t
|t | |t | | jd qd S r   )r   r   r   rg   r   r-   r	   rL   r   rT   r   rr   rh   r   r   r   r   r     s    z#_TestDSTIBase.test_definition_orthoNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )_TestDSTIVBasec                 C   sd   t t j| j}tD ]J}t j|| jd}t|ddd}t|dd}t|j	| t
||| jd qd S r   )r   r   r   rg   r   r-   r	   rO   r   rT   r   rh   r   r   r   r   r     s    z$_TestDSTIVBase.test_definition_orthoNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )TestDSTIDoublec                 C   s   t j| _d| _d| _d S Nr   r   r   ri   r   r   r   rj     s    zTestDSTIDouble.setup_methodNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )TestDSTIFloatc                 C   s   t j| _d| _d| _d S r   r   ri   r   r   r   rj     s    zTestDSTIFloat.setup_methodNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )TestDSTIIntc                 C   s   t | _d| _d| _d S r   r   ri   r   r   r   rj     s    zTestDSTIInt.setup_methodNr   r   r   r   r   r      s   r   c                   @   s   e Zd Zdd ZdS )TestDSTIIDoublec                 C   s   t j| _d| _d| _d S Nre   r8   r   ri   r   r   r   rj     s    zTestDSTIIDouble.setup_methodNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )TestDSTIIFloatc                 C   s   t j| _d| _d| _d S N   r8   r   ri   r   r   r   rj     s    zTestDSTIIFloat.setup_methodNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )TestDSTIIIntc                 C   s   t | _d| _d| _d S r   r   ri   r   r   r   rj     s    zTestDSTIIInt.setup_methodNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )TestDSTIIIDoublec                 C   s   t j| _d| _d| _d S r   r   ri   r   r   r   rj     s    zTestDSTIIIDouble.setup_methodNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )TestDSTIIIFloatc                 C   s   t j| _d| _d| _d S Nru   r   r   ri   r   r   r   rj   $  s    zTestDSTIIIFloat.setup_methodNr   r   r   r   r   r   #  s   r   c                   @   s   e Zd Zdd ZdS )TestDSTIIIIntc                 C   s   t | _d| _d| _d S r   r   ri   r   r   r   rj   +  s    zTestDSTIIIInt.setup_methodNr   r   r   r   r   r   *  s   r   c                   @   s   e Zd Zdd ZdS )TestDSTIVDoublec                 C   s   t j| _d| _d| _d S r   r   ri   r   r   r   rj   2  s    zTestDSTIVDouble.setup_methodNr   r   r   r   r   r   1  s   r   c                   @   s   e Zd Zdd ZdS )TestDSTIVFloatc                 C   s   t j| _d| _d| _d S )Nr   r   ri   r   r   r   rj   9  s    zTestDSTIVFloat.setup_methodNr   r   r   r   r   r   8  s   r   c                   @   s   e Zd Zdd ZdS )TestDSTIVIntc                 C   s   t | _d| _d| _d S r   r   ri   r   r   r   rj   @  s    zTestDSTIVInt.setup_methodNr   r   r   r   r   r   ?  s   r   c                   @   s   e Zd Zdd Zdd ZdS )_TestIDSTBasec                 C   s   d | _ d | _d | _d S r   rf   ri   r   r   r   rj   G  s    z_TestIDSTBase.setup_methodc                 C   s   t D ]}t| j|| j\}}}t|| jd}| jdkrH|d|d   }n|d|  }t|j| t|t	| |t	| | j
d| d qd S r   )rq   r*   r#   rg   r
   r   rT   r   r   rr   rh   r   r   r   r   rt   L  s    
z_TestIDSTBase.test_definitionNr   r   r   r   r   r   F  s   r   c                   @   s   e Zd Zdd ZdS )TestIDSTIDoublec                 C   s   t j| _d| _d| _d S r   r   ri   r   r   r   rj   ^  s    zTestIDSTIDouble.setup_methodNr   r   r   r   r   r   ]  s   r   c                   @   s   e Zd Zdd ZdS )TestIDSTIFloatc                 C   s   t j| _d| _d| _d S r   r   ri   r   r   r   rj   e  s    zTestIDSTIFloat.setup_methodNr   r   r   r   r   r   d  s   r   c                   @   s   e Zd Zdd ZdS )TestIDSTIIntc                 C   s   t | _d| _d| _d S r   r   ri   r   r   r   rj   l  s    zTestIDSTIInt.setup_methodNr   r   r   r   r   r   k  s   r   c                   @   s   e Zd Zdd ZdS )TestIDSTIIDoublec                 C   s   t j| _d| _d| _d S r   r   ri   r   r   r   rj   s  s    zTestIDSTIIDouble.setup_methodNr   r   r   r   r   r   r  s   r   c                   @   s   e Zd Zdd ZdS )TestIDSTIIFloatc                 C   s   t j| _d| _d| _d S r   r   ri   r   r   r   rj   z  s    zTestIDSTIIFloat.setup_methodNr   r   r   r   r   r   y  s   r   c                   @   s   e Zd Zdd ZdS )TestIDSTIIIntc                 C   s   t | _d| _d| _d S r   r   ri   r   r   r   rj     s    zTestIDSTIIInt.setup_methodNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )TestIDSTIIIDoublec                 C   s   t j| _d| _d| _d S r   r   ri   r   r   r   rj     s    zTestIDSTIIIDouble.setup_methodNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )TestIDSTIIIFloatc                 C   s   t j| _d| _d| _d S Nr   r   r   ri   r   r   r   rj     s    zTestIDSTIIIFloat.setup_methodNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )TestIDSTIIIIntc                 C   s   t | _d| _d| _d S r   r   ri   r   r   r   rj     s    zTestIDSTIIIInt.setup_methodNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )TestIDSTIVDoublec                 C   s   t j| _d| _d| _d S r   r   ri   r   r   r   rj     s    zTestIDSTIVDouble.setup_methodNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )TestIDSTIVFloatc                 C   s   t j| _d| _d| _d S Nr   r   r   ri   r   r   r   rj     s    zTestIDSTIVFloat.setup_methodNr   r   r   r   r   r     s   r   c                   @   s   e Zd Zdd ZdS )TestIDSTIVntc                 C   s   t | _d| _d| _d S r   r   ri   r   r   r   rj     s    zTestIDSTIVnt.setup_methodNr   r   r   r   r   r     s   r   c                   @   sL   e Zd ZdZejejgZdd Zdd Z	dd Z
dd	 Zd
d Zdd ZdS )TestOverwritezCheck input overwrite behavior.c                 K   sP   |  }	||	|||||d d|j|j|j|||f }
|sLt|	|d|
 d d S )N)overwrite_xz%%s(%s%r, %r, axis=%r, overwrite_x=%r)zspurious overwrite in %s)rp   )r,   r`   rT   r/   r   )rX   r&   routiner#   Zfftsizer{   rB   r   kwZx2sigr   r   r   _check  s         zTestOverwrite._checkc           	      C   s   t jd t |t jr8t jj| dt jj|   }nt jj| }||}dD ].}dD ]$}dD ]}| |||d ||| qbqZqRd S )N  rQ   )r   r8   r   r   )TF)Nr9   )r   r|   seedZ
issubdtypeZcomplexfloatingr}   r   r   )	rX   r   rT   r/   r{   r'   r#   r   rB   r   r   r   	_check_1d  s    
zTestOverwrite._check_1dc                 C   s@   | j D ]4}| t|dd | t|dd | t|dd qd S N)rw   r<   )rw   r8   r   )r8   rw   r   )real_dtypesr   r   rX   rT   r   r   r   test_dct  s    
zTestOverwrite.test_dctc                 C   s@   | j D ]4}| t|dd | t|dd | t|dd qd S r   )r   r   r   r   r   r   r   	test_idct  s    
zTestOverwrite.test_idctc                 C   s@   | j D ]4}| t|dd | t|dd | t|dd qd S r   )r   r   r	   r   r   r   r   test_dst  s    
zTestOverwrite.test_dstc                 C   s@   | j D ]4}| t|dd | t|dd | t|dd qd S r   )r   r   r
   r   r   r   r   	test_idst  s    
zTestOverwrite.test_idstN)r`   ra   rb   __doc__r   r   Zfloat64r   r   r   r   r   r   r   r   r   r   r   r     s   	r   c                   @   s  e Zd ZdZddddgZddgZejdZ	d	Z
e	je
 Zejd
eefeefgejdddddgdddgdddgdddggejdeejddgdd Zejdeefeefgejdeejdedd Zejdeefeefgejdeejdddgdd Zejd
eefeefgdd Zejdeegd d! Zejd
eefeefgejddddgdddggd"d# ZdS )$Test_DCTN_IDCTNre   r   r8   r   r   Nr9   r   )rx   rw   zfforward,finverseaxes)r   r   )r   r   r   )r<   r   r<   dct_typerB   c                 C   s6   || j |||d}|||||d}t| j |dd d S )Nr#   r   rB   r   rz   r'   r   )rX   fforwardfinverser   r   rB   tmpr   r   r   test_axes_round_trip  s    
z$Test_DCTN_IDCTN.test_axes_round_tripzfforward,fforward_refc                 C   s4   || j |d |d}|| j ||d}t||dd d S )Nr   r#   rB      rz   r   )rX   r   Zfforward_refr   rB   y1r   r   r   r   test_dctn_vs_2d_reference  s    z)Test_DCTN_IDCTN.test_dctn_vs_2d_referencezfinverse,finverse_refc                 C   s>   t | j||d}||||d}||||d}t||dd d S )Nr  r  rz   )r   r'   r   )rX   r   Zfinverse_refr   rB   fdatar  r   r   r   r   test_idctn_vs_2d_reference  s    z*Test_DCTN_IDCTN.test_idctn_vs_2d_referencec              	   C   s   t tdd || j| jjd dd W 5 Q R X t tdd || j| jjd d d W 5 Q R X t tdd || j| jjdd W 5 Q R X d S )NzBwhen given, axes and shape arguments have to be of the same length)matchr   r   r/   r   )assert_raisesr!   r'   r/   )rX   r   r   r   r   r   test_axes_and_shape  s    ""z#Test_DCTN_IDCTN.test_axes_and_shaper   c                 C   s    || j dd d}t|jd d S )N)   r  r	  )r'   r   r/   )rX   r   r   r   r   r   
test_shape#  s    zTest_DCTN_IDCTN.test_shapec                 C   s8   || j d |dd}||d |dd}t| j || jd d S )Nr9   )r/   r   rB   rz   )r'   r   rh   )rX   r   r   r   r   r   r   r   test_shape_is_none_with_axes(  s    z,Test_DCTN_IDCTN.test_shape_is_none_with_axes)r`   ra   rb   rh   r   Znormsr   r|   ZRandomStateZrstater/   r}   r'   pytestmarkZparametrizer   r   r   r   r  r4   r6   r  r5   r7   r  r  r  r  r   r   r   r   r     sf   

      

  r   )N)N)N)N)dos.pathr   r   Znumpyr   Znumpy.testingr   r   r  r   r
  Zscipy.fftpack._realtransformsr   r   r	   r
   r   r   r   r   load__file__r   r.   r   r   r   r    rq   r)   r*   r4   r5   r6   r7   rH   rL   rN   rO   rP   rc   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   <module>   s   (







 "5