a
    Qe                     @   s^   d dl Zd dlmZmZmZ d dlmZmZm	Z	 G dd dZ
G dd dZG dd	 d	ZdS )
    N)assert_equalassert_almost_equalassert_allclose)logitexpit	log_expitc                   @   s,   e Zd Zdd Zdd Zdd Zdd Zd	S )
	TestLogitc                 C   sp   t ddd}t j||d}t jdd t|}W d    n1 sF0    Y  t|| t|jt | d S )Nr      
   dtypeignore)divide)nplinspacearrayerrstater   r   r   r   selfr   expectedaactual r   R/var/www/sunrise/env/lib/python3.9/site-packages/scipy/special/tests/test_logit.pycheck_logit_out   s    &
zTestLogit.check_logit_outc                 C   s:   t jt j ddddddddt jg
t jd	}| d
| d S )NgcY@ gwQDB.gԙ̿gc?gC?C.?gN7?Q?gcY@ @r   f4)r   r   inffloat32r   r   r   r   r   r   test_float32   s    zTestLogit.test_float32c                 C   s4   t t j ddddddddt jg
}| d	| d S )
Ngu? g$/Qr   g̿g?gDB.?g$/Q?gu? @f8)r   r   r   r   r   r   r   r   test_float64   s    zTestLogit.test_float64c                 C   s^   t t jgd }t jdd" tt g d}W d    n1 sF0    Y  t|| d S )N   r   )invalid)g      g       g       @g      @)r   r   nanr   r   r   )r   r   r   r   r   r   test_nan"   s    0zTestLogit.test_nanN)__name__
__module____qualname__r   r    r"   r&   r   r   r   r   r      s   
r   c                   @   s,   e Zd Zdd Zdd Zdd Zdd Zd	S )
	TestExpitc                 C   sD   t ddd}t j||d}t|}t|| t|jt | d S )Nr#   r
   r   )r   r   r   r   r   r   r   r   r   r   r   check_expit_out+   s
    
zTestExpit.check_expit_outc                 C   s$   t jg dt jd}| d| d S )N)
{j?\wb֥?g~ܟ?g5?s4 ?ǻe?g+?>S?w?gݙ?gl?r   r   )r   r   r   r,   r   r   r   r   r    2   s    
zTestExpit.test_float32c                 C   s   t g d}| d| d S )N)
r-   r.   gIq?g!/?r/   r0   gw?4S?r1   g4ؙ?g!/l?r!   )r   r   r,   r   r   r   r   r"   :   s    zTestExpit.test_float64c                 C   sz   t jt jt jfD ]d}dD ]Z}t j||d}tt|ddd tt| ddd tt|j| tt| j| qqd S )N)X   Y   i    i\,  i],  r   g      ?g#B;)Zatolg        )	r   r   float64Z
longdoubler   r   r   r   r   )r   r   nr   r   r   
test_largeB   s    zTestExpit.test_largeN)r'   r(   r)   r,   r    r"   r7   r   r   r   r   r*   *   s   r*   c                   @   s,   e Zd Zdd Zdd Zdd Zdd Zd	S )
TestLogExpitc                 C   s$   t g d}t|}t|| d S )N)g     g     pg     @g     Ar   r   r   r   r   xyr   r   r   test_large_negativeN   s    z TestLogExpit.test_large_negativec                 C   s.   t g d}t|}t|t g d d S )N)g     p@g     @@g     @)       r>   r>   r9   r:   r   r   r   test_large_positiveS   s    z TestLogExpit.test_large_positivec                 C   s0   t g d}t|}g d}t||dd d S )N)皙&.r   &.>皙?r	   r
   d   i  r4   i  i  )g     @gC   4g(U $g\ÁcgIgOg>B.g9B.gB.g]g$zԿgLVg]gCx^&ْgó(8g   g'1     gV瞯<Zrtol)r   r   r   r   r   r;   r<   r   r   r   r   test_basic_float64\   s    zTestLogExpit.test_basic_float64c                 C   sB   t jg dt jd}t|}t jg dt jd}t||dd d S )N)r@   rA   rB   rC   rD   rE   rF   r   rG   rH   r	   r
   rI   r   )g      @g      4g $gcgm{qg8j	C.rM   rM   g@g %zԿgcAVgα!gƠ>rJ   )r   r   r   r   r   rK   r   r   r   test_basic_float32{   s    

zTestLogExpit.test_basic_float32N)r'   r(   r)   r=   r?   rL   rN   r   r   r   r   r8   L   s   	r8   )numpyr   Znumpy.testingr   r   r   Zscipy.specialr   r   r   r   r*   r8   r   r   r   r   <module>   s
   #"