U
    Hvf                     @   s\   d dl Zzd dlZW n ek
r(   Y nX dd Zdd Zdd Zdd	 Zed
krXe  dS )    Nc                 C   s   t t | t dS )Nz-0.5)mpmathZlambertwexpZmpf)x r   I/tmp/pip-unpacked-wheel-96ln3f52/scipy/special/_precompute/wrightomega.pympmath_wrightomega	   s    r   c                 C   s   | }t | }t|| | S N)r   abs)r   Zseriesdesiredr   r   r   wrightomega_series_error   s    r   c                 C   s"   t | }t| }t|| | S r   )r   r   r   r	   )r   Zexponential_approxr
   r   r   r   wrightomega_exp_error   s    
r   c               
   C   s   dt tj } td dD ]2}td t|}W 5 Q R X t|||| k  qtd dD ]2}td t|}W 5 Q R X t|||| k  q\d S )N   zSeries Error)g     j@g    _Bg  4&kCg@xDd   z	Exp error)iiiii8ipiDi)	npZfinfofloatZepsprintr   Zworkdpsr   r   )Zdesired_errorr   errorr   r   r   main   s    r   __main__)	Znumpyr   r   ImportErrorr   r   r   r   __name__r   r   r   r   <module>   s   