| Server IP : 118.27.122.248 / Your IP : 216.73.217.130 Web Server : Apache System : Linux web0264.sh.tyo1 4.18.0-553.79.1.lve.el7h.x86_64 #1 SMP Wed Oct 15 16:34:46 UTC 2025 x86_64 User : c9415830 ( 11735) PHP Version : 8.4.17 Disable Function : NONE MySQL : OFF | cURL : ON | WGET : ON | Perl : ON | Python : ON | Sudo : OFF | Pkexec : OFF Directory : /opt/alt/python36/lib64/python3.6/test/__pycache__/ |
Upload File : |
3
�w2_[� � @ sp d dl mZmZmZ d dlmZ d dlZd dlZd dlZd dl Z d dl
Z
d dl Z d dl Z dZ
ed�Zed�Zed�Ze jjZd*\ZZee d+kZed kr�e jd ZneZej j e�p�ej Z ej j e d �Z!ej j e d
�Z"dd� Z#dd� Z$dd� Z%dd� Z&dd� Z'dd� Z(dd� Z)d d � Z*d,d d!�Z+G d"d#� d#e,�Z-G d$d%� d%ej.�Z/G d&d'� d'ej.�Z0d(d)� Z1ed k�rle1� dS )-� )� run_unittest�verbose�requires_IEEE_754)�supportNg�h㈵��>�nan�infz-inf� ��7y�AC��):��@� �__main__zmath_testcases.txtzcmath_testcases.txtc C s. t jdt jd| ��d }|dk r*|d }|S )a� Convert a non-NaN float x to an integer, in such a way that
adjacent floats are converted to adjacent integers. Then
abs(ulps(x) - ulps(y)) gives the difference in ulps between two
floats.
The results from this function will only make sense on platforms
where native doubles are represented in IEEE 754 binary64 format.
Note: 0.0 and -0.0 are converted to 0 and -1, respectively.
z<qz<dr � �? l )�struct�unpack�pack)�x�n� r �3/opt/alt/python36/lib64/python3.6/test/test_math.py�to_ulps# s
r c C s� t t| ��} tj| �s tj| �r$| S tjdtjd| ��d }tjdtjd|d ��d }tj|�r�tjdtjd|d ��d }| | S || S dS )a� Return the value of the least significant bit of a
float x, such that the first float bigger than x is x+ulp(x).
Then, given an expected result x and a tolerance of n ulps,
the result y should be such that abs(y-x) <= n * ulp(x).
The results from this function will only make sense on platforms
where native doubles are represented in IEEE 754 binary64 format.
z<qz<dr � N)�abs�float�math�isnan�isinfr r r )r r Zx_nextZx_prevr r r �ulp4 s
r c C s | rdt | | d @ � S dS )zANumber of '1' bits in binary expansion of a nonnnegative integer.r r )�count_set_bits)r r r r r a s r c C sD || d? }|sdS |dkr | S | | dB }t | |�t ||� S dS )z�Product of integers in range(start, stop, 2), computed recursively.
start and stop should both be odd, with start <= stop.
r N)�partial_product)�start�stopZ
numfactorsZmidr r r r e s r c C sd d }}xJt t| j� ��D ]6}|t| |d ? d dB | |? d dB �9 }||9 }qW || t| � > S )z�Factorial of nonnegative integer n, via "Binary Split Factorial Formula"
described at http://www.luschny.de/math/factorial/binarysplitfact.html
r )�reversed�range�
bit_lengthr r )r �innerZouter�ir r r � py_factorials s
* r& c C sL t t| �t|� �}t | | �}||ks0||kr4dS d}|j||||�S dS )z�Given finite floats `expected` and `got`, check that they're
approximately equal to within the given number of ulps or the
given absolute tolerance, whichever is bigger.
Returns None on success and an error message on failure.
NzAerror = {:.3g} ({:d} ulps); permitted error = {:.3g} or {:d} ulps)r r �format)�expected�got�ulp_tol�abs_tolZ ulp_errorZ abs_errorZfmtr r r �
ulp_abs_check~ s r, c
c s� t | ���}x�|D ]x}d|kr.|d|jd�� }|j� s8q|jd�\}}|j� \}}}|j� }|d } |dd� }
||t|�t| �|
fV qW W dQ R X dS )z�Parse a file with test values
-- starts a comment
blank lines, or lines containing only a comment, are ignored
other lines are expected to have the form
id fn arg -> expected [flag]*
z--Nz->r r )�open�index�strip�splitr ) �fname�fp�line�lhs�rhs�id�fn�arg�
rhs_pieces�exp�flagsr r r �parse_mtestfile� s
r< c
c s� t | ���}x�|D ]�}|jd�s|j� r*q|jd�\}}|j� \}}}}|j� } | d | d }
} | dd� } ||t|�t|�t|
�t| �| fV qW W dQ R X dS )z�Parse a file with test values
Empty lines or lines starting with -- are ignored
yields id, fn, arg_real, arg_imag, exp_real, exp_imag
z--z->r r r N)r- �
startswithr/ r0 r )
r1 r2 r3 r4 r5 r6 r7 Zarg_realZarg_imagr9 Zexp_realZexp_imagr; r r r �parse_testfile� s
r> � � c C s� || kr dS d}t | t�r.t |t�r.t|�}n t |t�rJt | t�rJt| �} t | t�r�t |t�r�tj| �rxtj|�rxd}n$tj| �s�tj|�r�nt| |||�}|dk r�d}|j| |�}|dj|�7 }|S dS dS )a� Compare arguments expected and got, as floats, if either
is a float, using a tolerance expressed in multiples of
ulp(expected) or absolutely (if given and greater).
As a convenience, when neither argument is a float, and for
non-finite floats, exact equality is demanded. Also, nan==nan
as far as this function is concerned.
Returns None on success and an error message on failure.
Nz not equalzexpected {!r}, got {!r}z ({}))�
isinstancer �intr r r r, r' )r( r) r* r+ �failure�fail_fmtZfail_msgr r r � result_check� s&
rE c @ s e Zd Zdd� Zdd� ZdS )� MyIndexablec C s
|| _ d S )N)�value)�selfrG r r r �__init__� s zMyIndexable.__init__c C s | j S )N)rG )rH r r r � __index__� s zMyIndexable.__index__N)�__name__�
__module__� __qualname__rI rJ r r r r rF � s rF c @ s� e Zd Zdjdd�Zdd� Zdd� Zd d
� Zd d � Zd
d� Zdd� Z dd� Z
dd� Z dd� Z e
dd� �Zdd� Zdd � Zd d � Zd d � Zd!d"� Zd#d$� Zejd%d&� �Zd'd(� Zd)d*� Zd+d,� Ze
ej e d-�d.d/� ��Z d0d1� Z d2d3� Z d4d5� Z!d6d7� Z"d8d9� Z#e
d:d;� �Z$e
ej%d<d�d=d>� ��Z&d?d@� Z'dAdB� Z(dCdD� Z)dEdF� Z*dGdH� Z+dIdJ� Z,dKdL� Z-dMdN� Z.dOdP� Z/e
ej e0j1dQ�dRkdS�dTdU� ��Z2dVdW� Z3dXdY� Z4dZd[� Z5d\d]� Z6e
d^d_� �Z7e
d`da� �Z8ej9e:db�dcdd� �Z;e
dedf� �Z<e
dgdh� �Z=diS )k� MathTestsr? � c C s, t ||||�}|dk r(| jdj||�� dS )aa Compare arguments expected and got, as floats, if either
is a float, using a tolerance expressed in multiples of
ulp(expected) or absolutely, whichever is greater.
As a convenience, when neither argument is a float, and for
non-finite floats, exact equality is demanded. Also, nan==nan
in this function.
Nz{}: {})rE �failr' )rH �namer) r( r* r+ rC r r r �ftest� s zMathTests.ftestc C s8 | j dtjd� | j dtjd� | jtjdtj � d S )N�pig-DT�! @�egiW�
�@r )rR r rS rT � assertEqualZtau)rH r r r �
testConstants s zMathTests.testConstantsc C s� | j ttj� | jdtjd�tj� | jdtjd�tjd � | jdtjd�d� | j ttjt� | j ttjt� | j ttjdt � | j ttjdt � | j
tj tjt ��� d S ) Nzacos(-1)r zacos(0)r r zacos(1)���rW )
� assertRaises� TypeErrorr ZacosrR rS �
ValueError�INF�NINF�eps�
assertTruer �NAN)rH r r r �testAcos s zMathTests.testAcosc C s� | j ttj� | jdtjd�d� | jdtjd�d� | j ttjd� | j ttjd� | jtjt�t� | j ttjt� | j tj
tjt ��� d S )Nzacosh(1)r r zacosh(2)r g5�qB�?rW ) rX rY r ZacoshrR rZ rU r[ r\ r^ r r_ )rH r r r � testAcosh s zMathTests.testAcoshc C s� | j ttj� | jdtjd�tj d � | jdtjd�d� | jdtjd�tjd � | j ttjt� | j ttjt� | j ttjdt � | j ttjdt � | j
tj tjt ��� d S ) Nzasin(-1)r r zasin(0)r zasin(1)rW rW )
rX rY r ZasinrR rS rZ r[ r\ r] r^ r r_ )rH r r r �testAsin! s zMathTests.testAsinc C s� | j ttj� | jdtjd�d� | jdtjd�d� | jdtjd�d� | jtjt�t� | jtjt�t� | jtj tjt
��� d S ) Nzasinh(0)r zasinh(1)r g'�ya64�?z asinh(-1)rW g'�ya64�) rX rY r ZasinhrR rU r[ r\ r^ r r_ )rH r r r � testAsinh, s zMathTests.testAsinhc C s� | j ttj� | jdtjd
�tj d � | jdtjd�d� | jdtjd�tjd � | jdtjt�tjd � | jd tjt�tj d � | jtj tjt
��� d S ) Nzatan(-1)r r
zatan(0)r zatan(1)z atan(inf)r z
atan(-inf)rW ) rX rY r �atanrR rS r[ r\ r^ r r_ )rH r r r �testAtan5 s zMathTests.testAtanc C s� | j ttj� | jdtjd�d� | jdtjd�d� | jdtjd�d � | j ttjd� | j ttjd
� | j ttjt� | j ttjt� | j tj
tjt ��� d S ) Nzatanh(0)r z
atanh(0.5)g �?g �z��?z atanh(-0.5)r g �g �z��rW ) rX rY r rd rR ZatanhrZ r[ r\ r^ r r_ )rH r r r � testAtanh> s zMathTests.testAtanhc C s | j ttj� | jdtjd%d�tj d � | jdtjd&d�tj d � | jdtjdd�d� | jdtjdd�tjd � | jd tjdd�tjd � | jd
tjd t�tj� | jd tjd d'�tj� | jdtjd d(�tj� | jtjd d �d � | jtjd d
�d � | jtjd t�d � | j tj
tjd t ��� | jdtjd)t�tj � | jdtjd*d+�tj � | jdtjd,d-�tj � | jtjd.d �d/� | jtjd0d
�d1� | jtjd2t�d3� | j tj
tjd4t ��� | jdtjtt�tjd d � | jdtjtd5�tjd � | jdtjtd6�tjd � | jdtjtd �tjd � | jdtjtd
�tjd � | jdtjtt�tjd � | j tj
tjtt ��� | jdtjtt�tj d d � | jdtjtd7�tj d � | jdtjtd8�tj d � | jd tjtd �tj d � | jd tjtd
�tj d � | jd tjtt�tj d � | j tj
tjtt ��� | jd tjd
t�tj� | jd tjd
d9�tjd � | jd!tjd
d �tjd � | jtjd
t�d � | j tj
tjd
t ��� | jd"tjd:t�tj � | jd#tjd;d<�tj d � | jd$tjd=d �tj d � | jtjd>t�d?� | j tj
tjd@t ��� | j tj
tjt t��� | j tj
tjt dA��� | j tj
tjt dB��� | j tj
tjt d ��� | j tj
tjt d
��� | j tj
tjt t��� | j tj
tjt t ��� d S )CNz atan2(-1, 0)r r r z atan2(-1, 1)r
z atan2(0, 1)z atan2(1, 1)z atan2(1, 0)zatan2(0., -inf)g zatan2(0., -2.3)gffffff@zatan2(0., -0.)zatan2(-0., -inf)zatan2(-0., -2.3)zatan2(-0., -0.)zatan2(inf, -inf)� zatan2(inf, -2.3)zatan2(inf, -0.)zatan2(inf, 0.)zatan2(inf, 2.3)zatan2(inf, inf)zatan2(-inf, -inf)zatan2(-inf, -2.3)zatan2(-inf, -0.)zatan2(-inf, 0.)zatan2(-inf, 2.3)zatan2(-inf, inf)zatan2(2.3, -inf)zatan2(2.3, -0.)zatan2(2.3, 0.)zatan2(-2.3, -inf)zatan2(-2.3, -0.)zatan2(-2.3, 0.)rW rW gffffff�g �g �g �gffffff�g �g �g �g �g �g �g �g �g �gffffff�g �gffffff�g �g �gffffff�gffffff�g �gffffff�gffffff�g �gffffff�gffffff�g �) rX rY r Zatan2rR rS r\ rU r[ r^ r r_ )rH r r r � testAtan2I sf " zMathTests.testAtan2c C s | j ttj� | jtttjd��� | jdtjd�d� | jdtjd�d� | jdtjd�d� | jd tjd�d
� | jd tjd�d� | jd tjd�d� G d
d� d�}G dd� d�}| jdtj|� �d� | j ttj|� � |� }dd� |_| j ttj|� | j ttj|d
� d S )Ng �?z ceil(0.5)r z ceil(1.0)g �?z ceil(1.5)g �?r z
ceil(-0.5)r z
ceil(-1.0)z
ceil(-1.5)c @ s e Zd Zdd� ZdS )z$MathTests.testCeil.<locals>.TestCeilc S s dS )N�* r )rH r r r �__ceil__� s z-MathTests.testCeil.<locals>.TestCeil.__ceil__N)rK rL rM rj r r r r �TestCeil� s rk c @ s e Zd ZdS )z&MathTests.testCeil.<locals>.TestNoCeilN)rK rL rM r r r r �
TestNoCeil� s rl zceil(TestCeil())ri c W s | S )Nr )�argsr r r �<lambda>� s z$MathTests.testCeil.<locals>.<lambda>g �g �rW g ��rW ) rX rY r ZceilrU rB �typerR rj )rH rk rl �tr r r �testCeil� s
zMathTests.testCeilc C s� | j tjdd�d� | j tjdd�d� | j tjdd �d
� | j tjdd�d� | j tjdd �d � | jttj� | j tjdd�d� | j tjdd
�d� | j tjtd�t� | j tjtd�t� | j tjtd�t� | j tjtd�t� | j tjdt�d� | j tjdt�d� | j tjtt�t� | j tjtt�t� | j tjtt�t� | j tjtt�t� | jtjtjt d��� | jtjtjt t��� | jtjtjt t��� | jtjtjt t ��� | jtj
tjtt ��� | j t tjdt ��d� d S )Nr ri g �?g rg g @g @g @i���g �g �g �g �g �g �g �g �) rU r �copysignrX rY r[ r\ r^ r r_ r r )rH r r r � testCopysign� s0 zMathTests.testCopysignc C s� | j ttj� | jdtjtj d �dtd�d� | jdtjd�d� | jdtjtjd �dtd�d� | jdtjtj�d � y0| jtjtjt ��� | jtjtjt
��� W n4 t k
r� | j t tjt � | j t tjt
� Y nX | jtjtjt ��� d S )
Nz
cos(-pi/2)r r r )r+ zcos(0)z cos(pi/2)zcos(pi)rW )
rX rY r ZcosrR rS r r^ r r[ r\ rZ r_ )rH r r r �testCos� s $"zMathTests.testCosc C s� | j ttj� | jdtjd�d� | jdtjd�dtjd�d d� | jtjt�t� | jtjt�t� | jtj tjt
��� d S )Nzcosh(0)r r zcosh(2)-2*cosh(1)**2r rW ) rX rY r �coshrR rU r[ r\ r^ r r_ )rH r r r �testCosh� s &zMathTests.testCoshc C sr | j ttj� | jdtjtj�d� | jdtjtjd �d� | jdtjtj d �d � | jd tjd
�d
� d S ) Nz degrees(pi)g �f@z
degrees(pi/2)r g �V@zdegrees(-pi/4)r
g �F@z
degrees(0)r g �F�)rX rY r ZdegreesrR rS )rH r r r � testDegrees� s
zMathTests.testDegreesc C s� | j ttj� | jdtjd�dtj � | jdtjd�d� | jdtjd�tj� | jtjt�t� | jtjt�d� | j tj
tjt ��� | j t tjd� d S ) Nzexp(-1)r zexp(0)r zexp(1)g i@B rW )
rX rY r r: rR rT rU r[ r\ r^ r r_ �
OverflowError)rH r r r �testExp� s zMathTests.testExpc C sN | j ttj� | jdtjd�d� | jdtjd�d� | jdtjd�d� d S )Nzfabs(-1)r zfabs(0)r zfabs(1)rW )rX rY r ZfabsrR )rH r r r �testFabs� s zMathTests.testFabsc C s� | j tjd�d� | j tjd�d� d}xXtdd�D ]J}||9 }| j tj|�|� | j tjt|��|� | j tj|�t|�� q4W | jttjd � | jttjd
� | jttjdd � | jttjd � | jttjtj� d S ) Nr r g i� g �?�
�d g}Ô%�I�TrW g �g}Ô%�I��) rU r � factorialr" r r& rX rZ rS )rH Ztotalr% r r r �
testFactorial� s zMathTests.testFactorialc C s( | j ttjdd � | j ttjd� d S )Nr{ r| g}Ô%�I�T)rX rx r r} )rH r r r �testFactorialHugeInputs� s z!MathTests.testFactorialHugeInputsc C s> | j ttj� | jtttjd��� | jdtjd�d� | jdtjd�d� | jdtjd�d� | jd tjd�d� | jd
tjd�d� | jd tjd �d � | jd
tjd�d� | jdtjd �d � G dd� d�}G dd� d�}| jdtj|� �d� | j ttj|� � |� }dd� |_| j ttj|� | j ttj|d� d S ) Ng �?z
floor(0.5)r z
floor(1.0)g �?r z
floor(1.5)g �?z floor(-0.5)z floor(-1.0)z floor(-1.5)r zfloor(1.23e167)g��Nݯ�bzfloor(-1.23e167)c @ s e Zd Zdd� ZdS )z&MathTests.testFloor.<locals>.TestFloorc S s dS )Nri r )rH r r r � __floor__ s z0MathTests.testFloor.<locals>.TestFloor.__floor__N)rK rL rM r� r r r r � TestFloor s r� c @ s e Zd ZdS )z(MathTests.testFloor.<locals>.TestNoFloorN)rK rL rM r r r r � TestNoFloor s r� zfloor(TestFloor())ri c W s | S )Nr )rm r r r rn s z%MathTests.testFloor.<locals>.<lambda>g �rW g �rW g �����g��Nݯ��g��Nݯ��) rX rY r ZfloorrU rB ro rR r� )rH r� r� rp r r r � testFloor s$
zMathTests.testFloorc C s� | j ttj� | jdtjdd�d� | jdtjdd�d� | jdtjdd�d � | jd
tjdd�d� | jd tjdd�d� | jd tjdd�d� | jtjtjtd ��� | jtjtjd t��� | jtjtjtt��� | j ttjd d� | j ttjt d � | j ttjt
d � | j ttjt d� | j tjd
t �d
� | j tjdt �d� | j tjd
t
�d
� | j tjdt
�d� | j tjdd
�d� | j tjdt
�d� d S )Nz fmod(10, 1)r{ r g z
fmod(10, 0.5)g �?z
fmod(10, 1.5)g �?g �?z fmod(-10, 1)zfmod(-10, 0.5)zfmod(-10, 1.5)g @i����g �i����g �i����g �g �g �g �g �) rX rY r ZfmodrR r^ r r_ rZ r[ r\ rU )rH r r r �testFmod$ s( zMathTests.testFmodc s� � j ttj� � fdd�}|dtjd �d
� |dtjd�d� |dtjd�d� |d tjd
�d� � jtjt�d t� � jtjt�d t� � jtjtjt �d �� d S )Nc sB || \}}\}}t || �tks*||kr>� jd| ||f � d S )Nz%s returned %r, expected %r)r r] rP )rQ �resultr( �mantr: ZemantZeexp)rH r r � testfrexp= s z&MathTests.testFrexp.<locals>.testfrexpz frexp(-1)r � �?zfrexp(0)r zfrexp(1)zfrexp(2)r rW � �)r� r )r r )r� r )r� r )
rX rY r �frexprU r[ r\ r^ r r_ )rH r� r )rH r � testFrexp: s zMathTests.testFrexpz2fsum is not exact on machines with double roundingc s> ddl m} |j�|j� � � �fdd�}g dfdgdfddd.dd d/d0gdfd1d2d5gd7fd8dd:gd<fd>dd@gdBfdDd dFgdHfdd� tdd�D �tjd�fdd� tdd�D �tjd�fdd� td�D �dJg dKfd dd gd fdLdOdQdUgdfd d� tdVd"d#�D �dXg tjd%�fg
}x�t|�D ]|\}\}}ytj |�}W nR t
k
�r^ | j d&|||f � Y n* t k
�r� | j d'|||f � Y nX | j
||� �qW dd(lm}m} m}
x�td�D ]~} d)ddYdZd[d+gd, }d} x8td-�D ],}| d|� �d) | }
| |
7 } |j|
� �q�W |
|� ||�} | j
||�tj |�� �q�W d S )\Nr )�
float_infoc s� d\}}xb| D ]Z}t j|�\}}tt j|���|� }}||krT||| K }|}n ||| K }||7 }qW tttt|���d � � | �}|dkr�d|d > }|d| t||@ o�|d| d @ � }||7 }t j||�S )z�Full precision summation. Compute sum(iterable) without any
intermediate accumulation of error. Based on the 'lsum' function
at http://code.activestate.com/recipes/393090/
r r r rg )r r ) r r� rB �ldexp�max�len�binr �bool)�iterableZtmantZtexpr r� r: �tail�h)�etiny�mant_digr r �msum_ s
" (z MathTests.testFsum.<locals>.msumg g}Ô%�I�Tg �?g0��.�++g�d~�QJg @�5 g �?�6 r| g $@g (@g @g @c S s g | ] }d | �qS )g �?r )�.0r r r r �
<listcomp>� s z&MathTests.testFsum.<locals>.<listcomp>r i� z0x1.df11f45f4e61ap+2c S s g | ]}d| | �qS )g �?g �r )r� r r r r r� � s z-0x1.62a2af1bd3624p-1c S s g | ]}d |d d | �qS )g333333�?r r )r� r% r r r r� � s i� g333333�?g ��7y�ACg��ؗ�Ҝ<g��7y�ACc S s, g | ]$}d | d |d d |d �qS )g @�2 �4 r )r� r r r r r� � s i2 i� r i� z0x1.5555555555555p+970zDtest %d failed: got OverflowError, expected %r for math.fsum(%.100r)zAtest %d failed: got ValueError, expected %r for math.fsum(%.100r))�random�gauss�shuffle� gh�$.5��;g#B� ���;r{ �� g}Ô%�I��g �g�d~�Q�g @Cg �i���g �<g ��g @Cg������?Cg @Ci����g �9g @Cg @Cg @Cg @Ci����g �9g @Cg @Cg @Cg������?Ci���g �<g @Cg������?Cg����-�og����-��g �g��7y�ACi���g �<g�������?g��7y�ACg��7y�A�i���g �<g�������?g�������i���g �g �i����g}Ô%�I��gh�$.5���)�sysr� r� �min_expr" r �fromhex� enumerater Zfsumrx rP rZ rU r� r� r� �append)rH r� r� Z test_valuesr% �valsr( �actualr� r� r� �j�s�vr )r� r� r �testFsumL sT
zMathTests.testFsumc C sz t j}| j|dd�d� | j|dd�d� | j|dd�d� | j|dd�d� | j|dd�d� | j|dd�d� | j|dd�d� | j|dd�d� | j|dd�d� | j|dd�d� | j|d d
�d � d }d
}d}|| }|| }| j|||�|� | j|||�|� | j|| |�|� | j||| �|� | j||| �|� | j|| |�|� | j|| | �|� | j|| | �|� d}|| }|| }| j|||�|� | j|||�|� | j|| |�|� | j||| �|� | j||| �|� | j|| |�|� | j|| | �|� | j|| | �|� | jt|dd� | jt|dd� | j|td�td��d� d S )Nr r r� � � �x �T � l j9W P)�}l "o� `T"-� i� l �8P9�cC�/�[�N
S�9�b�Xl 9@ {�wM�2�75�k l Z��'^%#�>� ; �*Q!/ 61^g ^@g U@rW rW rW i��i����)r �gcdrU rX rY rF )rH r� �cr �y�a�br r r �testGcd� sP
zMathTests.testGcdc C s� | j ttj� | jdtjdd�d� | jdtjdd�d� | jtjtt�t� | jtjtt�t� | jtjtt�t� | jtjtt�t� | j t tjt
t
� | j tj tjdt��� | j tj tjtd ��� d S )
Nz
hypot(0,0)r z
hypot(3,4)rg r
r? g �?g @g �)
rX rY r ZhypotrR rU r_ r[ r\ rx � FLOAT_MAXr^ r )rH r r r � testHypot� s zMathTests.testHypotc C s@ | j ttj� | jdtjdd�d� | jdtjdd�d� | jdtjdd�d� | jdtjdd�d� | j ttjd d
� | j ttjdd
� | jtjd d�d � | jtjdd�d� | jtjtd �t� | jtjtd�t� | j tj
tjt d��� �xFdd d dd gD �]0}| jtjt| �t� | jtjt| �t� | jtjd | �d � | jtjd | �d � | jtjd | �d � | jtjd | �d!� | j tj
tjt | ��� | j ttjd |� | j ttjd"|� | jtjd |�d � | jtjd#|�d$� | jtjt|�t� | jtjt|�t� | j tj
tjt |��� �qW d S )%Nz
ldexp(0,1)r r z
ldexp(1,1)r z ldexp(1,-1)g �?z ldexp(-1,1)g �?i@B g � �� r{ r? � �( rW rW r� g �i���g �i���g �i+���i�� l d( l F�x:^V g �g �g �g �g �g �g �) rX rY r r� rR rx rU r[ r\ r^ r r_ )rH r r r r � testLdexp� s6 zMathTests.testLdexpc C s | j ttj� | jdtjdtj �d� | jdtjd�d� | jdtjtj�d� | jdtjdd�d � | jd
tjd d d �d � | jd
tjd d d�d� | jdtjd d �d� | j ttjd� | j ttjd d � | j ttjt� | jtjt �t � | j
tj tjt ��� d S )Nzlog(1/e)r zlog(1)r zlog(e)z log(32,2)r� r r? zlog(10**40, 10)r{ r� zlog(10**40, 10**20)r� z
log(10**1000)i� gO�+��@g �?rW l F�x:^V g ��)
rX rY r �logrR rT rZ r\ rU r[ r^ r r_ )rH r r r �testLog� s zMathTests.testLogc C sl | j ttj� x6ddd dd gD ] }| jtj|�tjt|��� q"W | j ttjd� | jtjt�t� d S )Nr �Z i, r rW ) rX rY r Zlog1pZassertAlmostEqualr rZ rU r[ )rH r r r r � testLog1p s
zMathTests.testLog1pc C s� | j ttj� | jtjd�d� | jtjd�d� | jtjd�d� | jtjdd �d� | jtjdd �d
� | jtjdd �d � | j ttjd� | j ttjt� | jtjtjt ��� d S )Nr g r g �?r
g @i� g ��@i g �@i� g @�@g �?g ��)
rX rY r �log2rU rZ r\ r^ r r_ )rH r r r �testLog2 s zMathTests.testLog2r{ c C s8 dd� t dd�D �}dd� t dd�D �}| j||� d S )Nc S s g | ]}t jt jd |���qS )g �?)r r� r� )r� r r r r r� , s z+MathTests.testLog2Exact.<locals>.<listcomp>i2 i c S s g | ] }t |��qS r )r )r� r r r r r� - s i���i���)r" rU )rH r� r( r r r �
testLog2Exact' s zMathTests.testLog2Exactc C s� | j ttj� | jdtjd�d � | jdtjd�d� | jdtjd�d� | jdtjdd �d
� | j ttjd
� | j ttjdd � | j ttjt� | jtjt �t � | j
tj tjt ��� d S )Nz
log10(0.1)g�������?r zlog10(1)r z log10(10)r{ zlog10(10**1000)i� g @�@g �?rW g ��)
rX rY r Zlog10rR rZ r\ rU r� r[ r^ r r_ )rH r r r � testLog100 s zMathTests.testLog10c s� � j ttj� � fdd�}|dtjd�d � |dtjd �d� � jtjt�dtf� � jtjt�dtf� tjt�}� jtj |d �� � jtj |d
�� d S )Nc sF || \}}\}}t || �tks.t || �rB� jd| ||f � d S )Nz%s returned %r, expected %r)r r] rP )rQ r� r( Zv1Zv2Ze1Ze2)rH r r �testmodf? s z$MathTests.testModf.<locals>.testmodfz modf(1.5)g �?� �?� �?z
modf(-1.5)g r r )r� r� g ��� �� �)r� r� g �)
rX rY r ZmodfrU r[ r\ r_ r^ r )rH r� Zmodf_nanr )rH r �testModf< s
zMathTests.testModfc C s� | j ttj� | jdtjdd�d� | jdtjdd�d� | jdtjdd�d� | jdtjdd �d� | jtjtd�t� | jtjtd�t� | jtjdt�d � | jtjdt�d � | jtj tjt
d��� | jtj tjdt
��� | jtj tjdt
��� | jtjdt
�d� | jtjd
t�d
� | jtjd
d �d
� | jtjd
d �d
� | jtjd
d
�d
� | jtjd
d
�d � | jtjd
d!�d � | j t tjd
d"� | j t tjd
d#� | j t tjd
d$� | j t tjd
t� | jtj tjd
t
��� | jtjtt�t� | jtjtd �t� | jtjtd �t� | jtjtd
�t� | jtjtd
�d � | jtjtd%�d � | jtjtd&�d
� | jtjtd'�d
� | jtjtd(�d
� | jtjtt�d
� | jtj tjtt
��� | jtjd)t�d
� | jtjd*d �d+� | jtjd,d �d
� | jtjd-d
�d
� | jtjd.d
�d � | jtjd/d0�d � | j t tjd1d2� | j t tjd3d4� | j t tjd5d6� | j t tjd7t� | jtj tjd8t
��� | jtjtt�t� | jtjtd �t� | jtjtd �t� | jtjtd
�t� | jtjtd
�d � | jtjtd9�d � | jtjtd:�d
� | jtjtd;�d
� | jtjtd<�d=� | jtjtt�d
� | jtj tjtt
��� | jtjd>t�d � | jtjd?d �d@� | j t tjdAd � | jtjdBd
�d � | jtjdCd
�d � | jtjdDdE�d � | jtjdFdG�d � | j t tjdHdI� | jtjdJdK�dL� | jtjdMt�d � | jtj tjdNt
��� | jtjd t�d � | jtjd d �d � | jtjd d �d � | jtjd d
�d � | jtjd d
�d � | jtjd dO�d � | jtjd dP�d � | jtjd dQ�d � | jtjd dR�d � | jtjd t�d � | jtjd t
�d � | jtjd d
�d � | jtjdSd
�d � | jtjt
d
�d � | jtjd dT�d � | jtjdUdV�d � | jtjt
dW�d � | j t tjdXd � | j t tjdYdZ� | jtjdt�d
� | jtjdt�d
� | jtjdt�t� | jtjdt�t� | jtjd[t�t� | jtjd\t�t� | jtjd]t�d
� | jtjd^t�d
� | jtjdt�t� | jtjdt�t� | jtjdt�d
� | jtjdt�d
� | jtjd_t�d
� | jtjd`t�d
� | jtjdat�t� | jtjdbt�t� | jdtjdcd �dd� | jdtjded
�d� | jdtjdfd �dg� | jdtjdhd
�d � | jdtjdidj�d � | jdtjdkdl�dm� | jd tjdndo�d � | jd tjdpdq�dr� | j t tjdsdt� | j t tjdud� d S )vNzpow(0,1)r r zpow(1,0)zpow(2,1)r z pow(2,-1)g �?g �?g g @gffffff@g @g .@g������@gffffff�?g�������?g�������?g�������?z (-2.)**3.g @z (-2.)**2.g @z (-2.)**1.z (-2.)**0.z
(-2.)**-0.z
(-2.)**-1.z
(-2.)**-2.g �?z
(-2.)**-3.g �?rW g �g �gffffff�g �g �g �gffffff�g �g �g �g �g �g �g �g �g �g �g �g �gffffff�g �g �g �g �g �g �gffffff�g �g �g �g �g �g �g �g �g �g �g �g �g �gffffff�g �g �g �g �g �g �g �gffffff�g �gffffff�g �gffffff�g �g �g �g .�g�������g��������g�������g�������gffffff��g��������g�������g�������gffffff��g �g �g �g �g �g �g �g �g �g �g �g �g �g �g �g ��g �g �g �) rX rY r �powrR rU r[ r\ r^ r r_ rZ )rH r r r �testPowO s� zMathTests.testPowc C sr | j ttj� | jdtjd�tj� | jdtjd�tjd � | jdtjd �tj d � | jd tjd
�d
� d S ) Nz radians(180)� z radians(90)r� r z radians(-45)�- r
z
radians(0)r i���)rX rY r ZradiansrR rS )rH r r r � testRadians� s
zMathTests.testRadiansc C s� | j ttj� | jdtjd�d� | jdtjtjd �d� | jdtjtj d �d� y0| jtjtjt��� | jtjtjt ��� W n4 t
k
r� | j t
tjt� | j t
tjt � Y nX | jtjtjt ��� d S )Nzsin(0)r z sin(pi/2)r r z
sin(-pi/2)rW ) rX rY r ZsinrR rS r^ r r[ r\ rZ r_ )rH r r r �testSin� s zMathTests.testSinc C s� | j ttj� | jdtjd�d� | jdtjd�d tjd�d d� | jdtjd�tjd� d� | jtjt�t� | jtjt�t� | j tj
tjt ��� d S ) Nzsinh(0)r zsinh(1)**2-cosh(1)**2r r zsinh(1)+sinh(-1)rW rW ) rX rY r ZsinhrR ru rU r[ r\ r^ r r_ )rH r r r �testSinh� s & zMathTests.testSinhc C s� | j ttj� | jdtjd�d� | jdtjd�d� | jdtjd�d� | jtjt�t� | j ttjd� | j ttjt� | j tj
tjt ��� d S ) Nzsqrt(0)r zsqrt(1)r zsqrt(4)r
r rW ) rX rY r �sqrtrR rU r[ rZ r\ r^ r r_ )rH r r r �testSqrt s zMathTests.testSqrtc
C s� | j ttj� | jdtjd�d� | jdtjtjd �d� | jdtjtj d �d� y0| jtjtjt��� | jtjtjt ��� W n, | j t
tjt� | j t
tjt � Y nX | jtjtjt ��� d S )Nztan(0)r z tan(pi/4)r
r z
tan(-pi/4)rW ) rX rY r ZtanrR rS r^ r r[ r\ rZ r_ )rH r r r �testTan s zMathTests.testTanc C s� | j ttj� | jdtjd�d� | jdtjd�tjd� dtd�d� | jdtjt�d� | jdtjt�d � | jtj tjt
��� d S )
Nztanh(0)r ztanh(1)+tanh(-1)r )r+ z tanh(inf)z
tanh(-inf)rW rW ) rX rY r �tanhrR r r[ r\ r^ r r_ )rH r r r �testTanh s zMathTests.testTanh�TANH_PRESERVES_ZERO_SIGNr z,system tanh() function doesn't copy the signc C s8 | j tjd�d� | j tjdtjd��tjdd�� d S )Ng g �?g �g �g �g �)rU r r� rr )rH r r r � testTanhSign! s zMathTests.testTanhSignc C s& | j tjd�d� | j tjd�d� | j ttjd��t� | j ttjd��t� | j tjd�d� | j tjd�d� | j tjd�d� | j tjd�d� | j tjd�d� | j tjd�d� G dd � d t�}G d
d � d t�}| j tj|� �d � | jttj� | jttjdd
� | jttj|� � d S )Nr g �?g����?g� !��?r g-���?Y@r| c @ s e Zd Zdd� ZdS )z'MathTests.test_trunc.<locals>.TestTruncc S s dS )Nr� r )rH r r r � __trunc__7 s z1MathTests.test_trunc.<locals>.TestTrunc.__trunc__N)rK rL rM r� r r r r � TestTrunc6 s r� c @ s e Zd ZdS )z)MathTests.test_trunc.<locals>.TestNoTruncN)rK rL rM r r r r � TestNoTrunc: s r� r� r rW rW g ��rW g�����rW g� !��r g-���?Y�i����)rU r Ztruncro rB �objectrX rY )rH r� r� r r r �
test_trunc* s zMathTests.test_truncc C s� | j tjd�� | j tjd�� | j tjd�� | j tjd�� | jtjtd��� | jtjtd��� | jtjtd��� d S )Ng g �?r r z-infg �g �)r^ r Zisfinite� assertFalser )rH r r r � testIsfiniteC s zMathTests.testIsfinitec C sx | j tjtd��� | j tjtd��� | j tjtd�d �� | jtjtd��� | jtjd�� | jtjd�� d S )Nr z-nanr g g �?)r^ r r r r� )rH r r r � testIsnanL s zMathTests.testIsnanc C s� | j tjtd��� | j tjtd��� | j tjd�� | j tjd�� | jtjtd��� | jtjd�� | jtjd�� d S )Nr z-infg �r g g �?g �)r^ r r r r� )rH r r r � testIsinfT s zMathTests.testIsinfc C s | j tjtj�� d S )N)r^ r r r )rH r r r �test_nan_constant] s zMathTests.test_nan_constantc C sJ | j tjtj�� | jtjd� | jtjtd�� | jtj td�� d S )Ng r z-inf)r^ r r r Z
assertGreaterrU r )rH r r r �test_inf_constanta s zMathTests.test_inf_constantzrequires verbose modec C s� yt jd�}W n | jd� Y nX |dkr8| jd� yt jd�}W n tk
rZ Y n X | jd� yt jd �}W n tk
r� Y n X | jd� d S )
Ni ʚ;z6underflowing exp() should not have raised an exceptionr z)underflowing exp() should have returned 0z.overflowing exp() didn't trigger OverflowErrorg �?z sqrt(-1) didn't raise ValueErrori 6e�g �)r r: rP rx r� rZ )rH r r r r �test_exceptionsm s
zMathTests.test_exceptionsc C st dh}d }t jdkrNtj� d }yttt|jd���}W n tk
rL Y nX d}g }� x�tt �D ]�\}}}} }
} } | dksb| dkr�qb|dkr�qb|d k r�|dk r�||kr�qbt
t |�}
d | ks�d | kr�d
}
n d| kr�d}
y |
|�}W n2 tk
� r� d
}Y n t k
�r d}Y nX d\}}t
|
|||�}|d k�r8qb|j||||�}|j|� qbW |�rp| jddj|� � d S )NZtan0064�darwinr �.z{}: {}({!r}): {}� �rect�polarr{ r? �invalidzdivide-by-zerorZ �overflowrx z Failures in test_testfile:
z
)r� r� )r{ r? )r? r� )r� �platformZmac_ver�tuple�maprB r0 rZ r> � test_file�getattrr rx rE r' r� rP �join)rH Z
SKIP_ON_TIGERZ osx_versionZ version_txtrD �failuresr6 r7 �arZaiZerZeir; �funcr� r* r+ rC �msgr r r �
test_testfile� sL
zMathTests.test_testfilec C s4 d}g }�x
t t�D � ]�\}}}}}tt|�}d|ks>d|krDd}n d|krPd}y ||�} W n. tk
rt d} Y n tk
r� d} Y nX d\}
} |d kr�d
}
n>|d kr�d } n0|d
kr�|dkr�|dk r�d}
n|dk r�d}
nd}
t|| |
| �} | d kr�q|j|||| �}
|j|
� qW |�r0| j ddj
|� � d S )Nz{}: {}({!r}): {}r� zdivide-by-zerorZ r� rx r? � Zgammar� ZlgammagV瞯�<Zerfcg �?r{ g $@r| i� z Failures in test_mtestfile:
z
)r? r� ) r< �math_testcasesr� r rZ rx rE r' r� rP r� )rH rD r� r6 r7 r8 r( r; r� r) r* r+ rC r� r r r �test_mtestfile� sB
zMathTests.test_mtestfileN)r? rO )>rK rL rM rR rV r` ra rb rc re rf rh rq r rs rt rv rw ry rz r~ r Z cpython_onlyr r� r� r� �unittestZskipIf�HAVE_DOUBLE_ROUNDINGr� r� r� r� r� r� r� Zrequires_mac_verr� r� r� r� r� r� r� r� r� r� � sysconfig�get_config_varr� r� r� r� r� r� r� Z
skipUnlessr r� r� r� r r r r rN � sl
= "
Y+
9rN c @ s� e Zd ZejZdd� Zdd� Zdd� Zdd� Zd d
� Z d d � Z
d
d� Z dd� Z dd� Z
dd� Zdd� Zdd� Zdd� Zdd � Zd d � Zd S ) � IsCloseTestsc O s* | j | j||f|�|�d||f d� d S )Nz%s and %s should be close!)r� )r^ �isclose)rH r� r� rm �kwargsr r r �
assertIsClose
s zIsCloseTests.assertIsClosec O s* | j | j||f|�|�d||f d� d S )Nz %s and %s should not be close!)r� )r� r )rH r� r� rm r r r r �assertIsNotClose s z IsCloseTests.assertIsNotClosec O s* x$|D ] \}}| j ||f|�|� qW d S )N)r )rH �examplesrm r r� r� r r r �assertAllClose s zIsCloseTests.assertAllClosec O s* x$|D ] \}}| j ||f|�|� qW d S )N)r )rH r rm r r� r� r r r �assertAllNotClose s z IsCloseTests.assertAllNotClosec C sR | j t�� | jdddd� W d Q R X | j t�� | jddddd� W d Q R X d S )Nr g0��.�++)�rel_tolg _�B)r r+ g0��.�+�g _��)rX rZ r )rH r r r �test_negative_tolerances s z%IsCloseTests.test_negative_tolerancesc C s$ d d
d d ddg}| j |ddd� d S )N� @�N���@i�;����90 � � �@� �NF )r r+ )r
r
)r r )r r )r r � �)r r )r r )r )rH Zidentical_examplesr r r �test_identical$ s zIsCloseTests.test_identicalc C s* d
d
dg}| j |dd� | j|dd� d S )N� �חAr g:�0�yE>g��n�yE>� ѭ��?�fN�ӭ��?)r g��&� .>� �חA)r r �:�0�yE�縦n�yE�)r r )r r )r r
)rH Z eight_decimal_places_examplesr r r �test_eight_decimal_places. s
z&IsCloseTests.test_eight_decimal_placesc C s* dd
d g}| j |dd� | j|dd� d S )
N��&� .>� gu�?j�/� g�������?)r g:�0�yE>)r+ )r r ��&� .�)r r �u�?j�/ʠ)r r )r
r )rH Znear_zero_examplesr r r �test_near_zero6 s
zIsCloseTests.test_near_zeroc C s< | j tt� | j ttdd� | j tt� | j ttdd� d S )Ng )r+ )r r[ r\ )rH r r r �test_identical_infinite@ s z$IsCloseTests.test_identical_infinitec C sR t t ft dfdt ftt ft tfttftdfdtftdfdtfg
}| j|dd� d S )Ng0��.�++g �?g������g�������?)r+ )r_ r[ r\ r
)rH Znot_close_examplesr r r �test_inf_ninf_nanG s
z IsCloseTests.test_inf_ninf_nanc C s4 d
d
dg}| j |dd� dddg}| j|dd� d S )N� �?g333333 @gY�� n�g )r � �?������?�Zb���ti�Tb���ti)r$ r$ �333333 ��333333 �)r) r* �Y�� n���Y�� n��)r+ r, )r$ r% )r& r$ )r' r( )r r
)rH Z zero_tolerance_close_examplesZ!zero_tolerance_not_close_examplesr r r �test_zero_toleranceV s z IsCloseTests.test_zero_tolerancec C s | j ddgdd� d S )N� r{ g�������?)r )r. r{ )r{ r. )r )rH r r r �test_asymmetryb s zIsCloseTests.test_asymmetryc C s( dd g}| j |dd� | j|dd� d S )
N��� ���[��[g:�0�yE>)r g��&� .>)r0 r1 )r2 r3 )r r
)rH Zinteger_examplesr r r �
test_integersf s zIsCloseTests.test_integersc C sh ddl m} |d�|d�f|d�|d�f|d�|d�f|d �|d
�fg}| j|d d � | j|d
d � d S )Nr )�Decimalz
1.00000001z1.0z1.00000001e-20z1.0e-20z1.00000001e-100z1.0e-100z
1.00000001e20z1.0e20g:�0�yE>)r g��&� .>)�decimalr5 r r
)rH r5 Zdecimal_examplesr r r �
test_decimalsn s zIsCloseTests.test_decimalsc C sd ddl m} |dd�d |d�f|d�|d�f|dd�|dd�fg}| j|d
d � | j|d d � d S )Nr )�Fractionr i �i�r{ � � r� g:�0�yE>)r g��&� .>i �i�l � (q�e� l F�x:^V )Z fractionsr8 r r
)rH r8 Zfraction_examplesr r r �test_fractionsy s zIsCloseTests.test_fractionsN)rK rL rM r r r r r r
r r r r! r" r# r- r/ r4 r7 r; r r r r r
s
r c C sN ddl m} tj� }|jtjt�� |jtjt�� |j| d�� t|� d S )Nr )� DocFileSuitez ieee754.txt) Zdoctestr< r� Z TestSuiteZaddTestZ makeSuiterN r r )r< Zsuiter r r � test_main� s r= )r r g��7y�AC)r? r@ )2Z test.supportr r r Ztestr r� r �osr� r r� r r] r r_ r[ r\ r� r� r� r r� r rK �argv�file�__file__�path�dirname�curdirZtest_dirr� r� r� r r r r r& r, r< r> rE r� rF ZTestCaserN r r= r r r r �<module> sX -
. {