a
    BCCf                     @   s   d dl mZ d dlZd dlZd dlZd dlZd dlZd dlm	Z	 ddl
mZ g dZe Ze ZdddZd	d
 Zdd Zdd Zdd Zdd Zd d dddZdd Zdd Zejdd Zdd ZdS )    )NumberN)copy_if_needed   )	good_size)r   set_workersget_workersc              
   C   sd   t | tr| f} zdd | D } W n< ty^ } z$|p8d}t| d|W Y d}~n
d}~0 0 | S )a  Convert ``x`` to an iterable sequence of int

    Parameters
    ----------
    x : value, or sequence of values, convertible to int
    name : str, optional
        Name of the argument being converted, only used in the error message

    Returns
    -------
    y : ``List[int]``
    c                 S   s   g | ]}t |qS  )operatorindex.0ar   r   W/var/www/html/django/DPS/env/lib/python3.9/site-packages/scipy/fft/_pocketfft/helper.py
<listcomp>&       z$_iterable_of_int.<locals>.<listcomp>valuez) must be a scalar or iterable of integersN)
isinstancer   	TypeError
ValueError)xnameer   r   r   _iterable_of_int   s    
&r   c                    s>  |du }|du }|sjt |d} fdd|D }t fdd|D rNtdtt|t|krjtd|st |d	}|rt|t|krtd
|rt| jkrtdt jt|  j} fddt||D }n,|rt j	}t j}n fdd|D }tdd |D r.td| dt
|t|fS )z2Handles shape and axes arguments for nd transformsNaxesc                    s"   g | ]}|d k r| j  n|qS )r   ndimr   r   r   r   r   5   r   z+_init_nd_shape_and_axes.<locals>.<listcomp>c                 3   s    | ]}| j kp|d k V  qdS )r   Nr   r   r   r   r   	<genexpr>7   r   z*_init_nd_shape_and_axes.<locals>.<genexpr>z$axes exceeds dimensionality of inputzall axes must be uniqueshapezBwhen given, axes and shape arguments have to be of the same lengthz)shape requires more axes than are presentc                    s&   g | ]\}}|d kr j | n|qS )r   )r   sr   r   r   r   r   G   r   c                    s   g | ]} j | qS r   r    r   r   r   r   r   L   r   c                 s   s   | ]}|d k V  qdS )r   Nr   )r   r!   r   r   r   r   N   r   invalid number of data points () specified)r   anyr   lensetr   rangeziplistr   tuple)r   r   r   ZnoshapeZnoaxesr   r   r   _init_nd_shape_and_axes.   s6    



r+   c                 C   sv   t | dst| } | jtjkr.t| tjS | jjdvrHt| tjS | jd}| j	d sbdnt
}tj| ||dS )zl
    Convert to array with floating or complex dtype.

    float16 values are also promoted to float32.
    dtypeZfc=ZALIGNEDT)r,   copy)hasattrnpZasarrayr,   Zfloat16Zfloat32kindZfloat64Znewbyteorderflagsr   array)r   r,   r.   r   r   r   	_asfarrayU   s    

r4   c                 C   s0   | |u rdS t |tjs&t|dr&dS | jdu S )z|
    Strict check for `arr` not sharing any data with `original`,
    under the assumption that arr = asarray(original)
    FZ	__array__N)r   r0   Zndarrayr/   base)Zarroriginalr   r   r   _datacopiedi   s
    r7   c           
      C   s   d}t dg| j }t||D ]>\}}| j| |krDt d|||< qt d| j| ||< d}qt|}|sv| | dfS t| j}t||D ]\}}|||< qt|| j}	| | |	|< |	dfS )z5Internal auxiliary function for _raw_fft, _raw_fftnd.FNr   T)	slicer   r(   r   r*   r)   r0   Zzerosr,   )
r   r   r   Z	must_copyr
   naxr!   axiszr   r   r   
_fix_shapeu   s     

r=   c                 C   s(   |dk rt d| dt| |f|fS )Nr   r"   r#   )r   r=   )r   r9   r;   r   r   r   _fix_shape_1d   s
    
r>      )NZbackwardZorthoforwardc                 C   sD   zt |  }|r|nd| W S  ty>   td| ddY n0 dS )zAReturns the pypocketfft normalization mode from the norm argumentr?   zInvalid norm value z,, should be "backward", "ortho" or "forward"N)	_NORM_MAPKeyErrorr   )Znormr@   Zinormr   r   r   _normalization   s    
rC   c                 C   s`   | d u rt tddS | dk rL| t kr4| dt 7 } q\td|  dt  n| dkr\td| S )Ndefault_workersr   r   z workers value out of range; got z, must not be less than zworkers must not be zero)getattr_config
_cpu_countr   )workersr   r   r   _workers   s    


rI   c                 c   s4   t  }tt| t_zdV  W |t_n|t_0 dS )a  Context manager for the default number of workers used in `scipy.fft`

    Parameters
    ----------
    workers : int
        The default number of workers to use

    Examples
    --------
    >>> import numpy as np
    >>> from scipy import fft, signal
    >>> rng = np.random.default_rng()
    >>> x = rng.standard_normal((128, 64))
    >>> with fft.set_workers(4):
    ...     y = signal.fftconvolve(x, x)

    N)r   rI   r	   r
   rF   rD   )rH   Zold_workersr   r   r   r      s
    r   c                   C   s   t tddS )zReturns the default number of workers within the current context

    Examples
    --------
    >>> from scipy import fft
    >>> fft.get_workers()
    1
    >>> with fft.set_workers(4):
    ...     fft.get_workers()
    4
    rD   r   )rE   rF   r   r   r   r   r      s    r   )N)numbersr   r	   os	threading
contextlibnumpyr0   Zscipy._lib._utilr   Zpypocketfftr   __all__localrF   	cpu_countrG   r   r+   r4   r7   r=   r>   rA   rC   rI   contextmanagerr   r   r   r   r   r   <module>   s,   
'
