Your IP : 216.73.216.85
��Yf
� � @ s d d l m Z m Z m Z 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 Z d d l Z d Z
e d � Z e d � Z e d � Z d d f \ Z Z e e d d
k Z e d k r� e
j d Z n e Z e j j e � p
e j Z e j j e d � Z e j j e d
� Z d d � Z d d d � Z! d d � Z" d d � Z# d d � Z$ d d d d � Z% d d � Z& d d � Z' Gd! d" � d" e( � Z) Gd# d$ � d$ e j* � Z+ Gd% d& � d& e j* � Z, d' d( � Z- e d k re- � d S))� )� run_unittest�verbose�requires_IEEE_754)�supportNg�h㈵��>�nan�infz-infg ��7y�ACg�):��@� �__main__zmath_testcases.txtzcmath_testcases.txtc C s= t j d t j d | � � d } | d k r9 | d } | S)aV 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 C doubles are represented in IEEE 754 binary64 format.
z<qz<dr � �? l )�structZunpackZpack)�x�n� r �3/opt/alt/python35/lib64/python3.5/test/test_math.py�to_ulps! s
" r � c C s< t | � t | � } t | � | k r, d Sd j | | � S)z�Given non-NaN floats `expected` and `got`,
check that they're equal to within the given number of ulps.
Returns None on success and an error message on failure.Nz*error = {} ulps; permitted error = {} ulps)r �abs�format)�expected�gotZulpsZ
ulps_errorr r r �
ulps_check0 s
r c C s | r d t | | d @� Sd S)zANumber of '1' bits in binary expansion of a nonnnegative integer.� r )�count_set_bits)r r r r r Q s r c C sT | | d ?} | s d S| d k r( | S| | d B} t | | � t | | � Sd S)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 U s r c C sx d } } xY t t | j � � � D]? } | t | | d ?d d B| | ?d d B� 9} | | 9} q# W| | 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_factorialc s
/r" gV瞯�<g
c C sh t j | � r | | k r d S| | } t | | t | � � } t | � | k rX d Sd j | | � S)a Determine whether non-NaN floats a and b are equal to within a
(small) rounding error. The default values for rel_err and
abs_err are chosen to be suitable for platforms where a float is
represented by an IEEE 754 double. They allow an error of between
9 and 19 ulps.Nz error = {}; permitted error = {})�math�isinf�maxr r )r r �rel_err�abs_err�errorZpermitted_errorr r r � acc_checkn s
r) c
c s� t | � �� } x� | D]� } d | k rA | d | j d � � } | j � sP q | j d � \ } } | j � \ } } } | j � } | d } | d d � }
| | t | � t | � |
f Vq WWd QRXd S)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�split�float) �fname�fp�line�lhs�rhs�id�fn�arg�
rhs_pieces�exp�flagsr r r �parse_mtestfile� s
r: c
c s� t | � �� } x� | D]� } | j d � s | j � r; q | j d � \ } } | j � \ } } } } | j � } | d | d }
} | d d � } | | t | � t | � t |
� t | � | f Vq WWd QRXd S)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, r- r. )
r/ r0 r1 r2 r3 r4 r5 Zarg_realZarg_imagr7 Zexp_realZexp_imagr9 r r r �parse_testfile� s
r< c @ s( e Z d Z d d � Z d d � Z d S)� MyIndexablec C s
| | _ d S)N)�value)�selfr> r r r �__init__� s zMyIndexable.__init__c C s | j S)N)r> )r? r r r � __index__� s zMyIndexable.__index__N)�__name__�
__module__� __qualname__r@ rA 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 d
� Z d d � Z d
d � Z d d � Z
d d � Z d d � Z e
d d � � Z d d � Z d d � Z d d � Z d d � Z d d � Z d! d" � Z e j d# d$ � � Z d% d& � Z d' d( � Z d) d* � Z e
e j e d+ � d, d- � � � Z d. d/ � Z d0 d1 � Z d2 d3 � Z! d4 d5 � Z" d6 d7 � Z# e
d8 d9 � � Z$ e
e j% d: d; � d< d= � � � Z& d> d? � Z' d@ dA � Z( dB dC � Z) dD dE � Z* dF dG � Z+ dH dI � Z, dJ dK � Z- dL dM � Z. dN dO � Z/ e
e j e0 j1 dP � dQ k dR � dS dT � � � Z2 dU dV � Z3 dW dX � Z4 dY dZ � Z5 d[ d\ � Z6 e
d] d^ � � Z7 e
d_ d` � � Z8 e j9 e: da � db dc � � Z; e
dd de � � Z< e
df dg � � Z= dh S)i� MathTestsc C s4 t | | � t k r0 | j d | | | f � d S)Nz%s returned %r, expected %r)r �eps�fail)r? �namer> r r r r �ftest� s zMathTests.ftestc C s0 | j d t j d � | j d t j d � d S)N�pigJ�M�! @�egB�
�@)rI r# rJ rK )r? r r r �
testConstants� s zMathTests.testConstantsc C s� | j t t j � | j d t j d � t j � | j d t j d � t j d � | j d t j d � d � | j t t j t � | j t t j t � | j t j
t j t � � � d S)Nzacos(-1)r zacos(0)r r
zacos(1)���) � assertRaises� TypeErrorr# ZacosrI rJ �
ValueError�INF�NINF�
assertTrue�isnan�NAN)r? r r r �testAcos� s # zMathTests.testAcosc C s� | j t t j � | j d t j d � d � | j d t j d � d � | j t t j d � | j t t j d � | j t j t � t � | j t t j t � | j t j
t j t � � � d S)Nzacosh(1)r r zacosh(2)r
g5�qB�?rM ) rN rO r# ZacoshrI rP � assertEqualrQ rR rS rT rU )r? r r r � testAcosh� s zMathTests.testAcoshc C s� | j t t j � | j d t j d � t j d � | j d t j d � d � | j d t j d � t j d � | j t t j t � | j t t j t � | j t j
t j t � � � d S)Nzasin(-1)r r
zasin(0)r zasin(1)rM ) rN rO r# ZasinrI rJ rP rQ rR rS rT rU )r? r r r �testAsin� s $ #zMathTests.testAsinc C s� | j t t j � | j d t j d � d � | j d t j d � d � | j d t j d � d � | j t j t � t � | j t j t � t � | j t j t j t
� � � d S) Nzasinh(0)r zasinh(1)r g'�ya64�?z asinh(-1)rM g'�ya64�) rN rO r# ZasinhrI rW rQ rR rS rT rU )r? r r r � testAsinh� s zMathTests.testAsinhc C s� | j t t j � | j d t j d
� t j d � | j d t j d � d � | j d t j d � t j d � | j d t j t � t j d � | j d t j t � t j d � | j t j t j t
� � � d S) Nzatan(-1)r r zatan(0)r zatan(1)z atan(inf)r
z
atan(-inf)rM ) rN rO r# �atanrI rJ rQ rR rS rT rU )r? r r r �testAtan� s $ ##$zMathTests.testAtanc C s� | j t t j � | j d t j d � d � | j d t j d � d � | j d t j d � d � | j t t j d � | j t t j d
� | j t t j t � | j t t j t � | j t j
t j t � � � d S) Nzatanh(0)r z
atanh(0.5)g �?g �z��?z atanh(-0.5)r g �g �z��rM ) rN rO r# r[ rI ZatanhrP rQ rR rS rT rU )r? r r r � testAtanh� s zMathTests.testAtanhc C s� | j t t j � | j d t j d% d � t j d � | j d t j d& d � t j d � | j d t j d d � d � | j d t j d d � t j d � | j d t j d d � t j d � | j d
t j d t � t j � | j d t j d d' � t j � | j d t j d d( � t j � | j t j d d � d � | j t j d d
� d � | j t j d t � d � | j t j
t j d t � � � | j d t j d) t � t j � | j d t j d* d+ � t j � | j d t j d, d- � t j � | j t j d. d � d/ � | j t j d0 d
� d1 � | j t j d2 t � d3 � | j t j
t j d4 t � � � | j d t j t t � t j d d � | j d t j t d5 � t j d � | j d t j t d6 � t j d � | j d t j t d � t j d � | j d t j t d
� t j d � | j d t j t t � t j d � | j t j
t j t t � � � | j d t j t t � t j d d � | j d t j t d7 � t j d � | j d t j t d8 � t j d � | j d t j t d � t j d � | j d t j t d
� t j d � | j d t j t t � t j d � | j t j
t j t t � � � | j d t j d
t � t j � | j d t j d
d9 � t j d � | j d! t j d
d � t j d � | j t j d
t � d � | j t j
t j d
t � � � | j d" t j d: t � t j � | j d# t j d; d< � t j d � | j d$ t j d= d � t j d � | j t j d> t � d? � | j t j
t j d@ t � � � | j t j
t j t t � � � | j t j
t j t dA � � � | j t j
t j t dB � � � | j t j
t j t d � � � | j t j
t j t d
� � � | j t j
t j t t � � � | j t j
t j t 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.)rM rM 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 �) rN rO r# Zatan2rI rJ rR rW rQ rS rT rU )r? r r r � testAtan2 sf '' &&""" "### "*&&&&&"+'''''""&& "#'' """""""zMathTests.testAtan2c C s� | j t t j � | j t t t j d � � � | j d t j d � d � | j d t j d � d � | j d t j d � d � | j d t j d � d
� | j d t j d � d � | j d t j d � d � Gd
d � d � } Gd d � d � } | j d t j | � � d � | j t t j | � � | � } d d � | _ | j t t j | � | j t t j | 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 Z d Z d d � Z d S)z$MathTests.testCeil.<locals>.TestCeilc S s d S)N�* r )r? r r r �__ceil__M s z-MathTests.testCeil.<locals>.TestCeil.__ceil__N)rB rC rD ra r r r r �TestCeilL s rb c @ s e Z d Z d S)z&MathTests.testCeil.<locals>.TestNoCeilN)rB rC rD r r r r �
TestNoCeilO s rc zceil(TestCeil())r` c W s | S)Nr )�argsr r r �<lambda>U s z$MathTests.testCeil.<locals>.<lambda>g �g �rM g ��rM ) rN rO r# ZceilrW �int�typerI ra )r? rb rc �tr r r �testCeil? s zMathTests.testCeilc C s� | j t j d d � d � | j t j d d � d � | j t j d d � d
� | j t j d d � d � | j t j d d � d � | j t t j � | j t j d d � d � | j t j d d
� d � | j t j t d � t � | j t j t d � t � | j t j t d � t � | j t j t d � t � | j t j d t � d � | j t j d t � d � | j t j t t � t � | j t j t t � t � | j t j t t � t � | j t j t t � t � | j t j t j t d � � � | j t j t j t t � � � | j t j t j t t � � � | j t j t j t t � � � | j t j
t j t t � � � | j t t j d t � � d � d S)Nr r` g �?g r^ g @g @g @i���g �g �g �g �g �g �g �g �) rW r# �copysignrN rO rQ rR rS rT rU r$ r )r? r r r � testCopysignY s0 """""zMathTests.testCopysignc C s; | j t t j � | j d t j t j d � d � | j d t j d � d � | j d t j t j d � d � | j d t j t j � d � yB | j t j t j t � � � | j t j t j t � � � Wn> t
k
r| j t
t j t � | j t
t j t � Yn X| j t j t j t � � � d S) Nz
cos(-pi/2)r
r zcos(0)r z cos(pi/2)zcos(pi)rM ) rN rO r# ZcosrI rJ rS rT rQ rR rP rU )r? r r r �testCos{ s $ # #
zMathTests.testCosc C s� | j t t j � | j d t j d � d � | j d t j d � d t j d � d d � | j t j t � t � | j t j t � t � | j t j t j t
� � � d S)Nzcosh(0)r r zcosh(2)-2*cosh(1)**2r
rM ) rN rO r# �coshrI rW rQ rR rS rT rU )r? r r r �testCosh� s 1zMathTests.testCoshc C s} | j t t j � | j d t j t j � d � | j d t j t j d � d � | j d t j t j d � d � d S)
Nz degrees(pi)g �f@z
degrees(pi/2)r
g �V@zdegrees(-pi/4)r g �F@g �F�)rN rO r# ZdegreesrI rJ )r? r r r � testDegrees� s #zMathTests.testDegreesc C s� | j t t j � | j d t j d � d t j � | j d t j d � d � | j d t j d � t j � | j t j t � t � | j t j t � d � | j t j
t j t � � � d S)Nzexp(-1)r zexp(0)r zexp(1)g rM ) rN rO r# r8 rI rK rW rQ rR rS rT rU )r? r r r �testExp� s # zMathTests.testExpc C sk | j t t j � | j d t j d � d � | j d t j d � d � | j d t j d � d � d S)Nzfabs(-1)r zfabs(0)r zfabs(1)rM )rN rO r# ZfabsrI )r? r r r �testFabs� s zMathTests.testFabsc C s( | j t j d � d � | j t j d � d � d } xx t d d � D]g } | | 9} | j t j | � | � | j t j t | � � | � | j t j | � t | � � qH W| j t t j d � | j t t j d
� | j t t j d � | j t t j d
� | j t t j t j � d S)Nr r g i� g �?�
�d g}Ô%�I�TrM g �l DG��d�r���Y�|g�O�s&YR5I l�� DG��d�r���Y�|g�O�s&YR5I g}Ô%�I��) rW r# � factorialr r. r" rN rP rJ )r? Ztotalr! r r r �
testFactorial� s
#zMathTests.testFactorialc C s0 | j t t j d � | j t t j d � d S)Nrr rs g}Ô%�I�Tl DG��d�r���Y�|g�O�s&YR5I )rN �
OverflowErrorr# rt )r? r r r �testFactorialHugeInputs� s z!MathTests.testFactorialHugeInputsc C s� | j t t j � | j t t t j d � � � | j d t j d � d � | j d t j d � d � | j d t j d � d � | j d t j d � d � | j d
t j d � d � | j d t j d � d � | j d
t j d � d � | j d t j d � d � Gd d � d � } Gd d � d � } | j d t j | � � d � | j t t j | � � | � } d d � | _ | j t t j | � | j t t j | 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 Z d Z d d � Z d S)z&MathTests.testFloor.<locals>.TestFloorc S s d S)Nr` r )r? r r r � __floor__� s z0MathTests.testFloor.<locals>.TestFloor.__floor__N)rB rC rD rx r r r r � TestFloor� s ry c @ s e Z d Z d S)z(MathTests.testFloor.<locals>.TestNoFloorN)rB rC rD r r r r � TestNoFloor� s rz zfloor(TestFloor())r` c W s | S)Nr )rd r r r re � s z%MathTests.testFloor.<locals>.<lambda>g �rM g �rM g �����g��Nݯ��g��Nݯ��) rN rO r# ZfloorrW rf rg rI rx )r? ry rz rh r r r � testFloor� s$ zMathTests.testFloorc C sC | j t t j � | j d t j d d � d � | j d t j d d � d � | j d t j d d � d � | j d
t j d d � d � | j d t j d d � d � | j d t j d d � d � | j t j t j t d � � � | j t j t j d t � � � | j t j t j t t � � � | j t t j d d � | j t t j t d � | j t t j t
d � | j t t j t d � | j t j d
t � d
� | j t j d t � d � | j t j d
t
� d
� | j t j d t
� d � | j t j d d
� d � | j t j d t
� d � d S)Nz fmod(10, 1)rr 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 �) rN rO r# ZfmodrI rS rT rU rP rQ rR rW )r? r r r �testFmod� s( """ zMathTests.testFmodc s� � j t t j � � f d d � } | d t j d � d
� | d t j d � d � | d t j d � d � | d t j d
� d � � j t j t � d t � � j t j t � d t � � j t j t j t � d � � d S)Nc sY | | \ } } \ } } t | | � t k s; | | k rU � j d | | | f � d S)Nz%s returned %r, expected %r)r rF rG )rH �resultr �mantr8 ZemantZeexp)r? r r � testfrexp� s " z&MathTests.testFrexp.<locals>.testfrexpz frexp(-1)r � �?zfrexp(0)r zfrexp(1)zfrexp(2)r
rM � �)r� r )r r )r� r )r� r
)
rN rO r# �frexprW rQ rR rS rT rU )r? r� r )r? r � testFrexp� s zMathTests.testFrexpz2fsum is not exact on machines with double roundingc s� d d l m } | j � | j � � � � f d d � } g d f d g d f d d d. d d d/ d0 g d f d1 d2 d5 g d7 f d8 d d: g d< f d> d d@ g dB f dD d dF g dH f d d � t d d � D� t j d � f d d � t d d � D� t j d � f d d � t d � D� dJ g dK f d d d g d f dL dO dQ dU g d f d d � t dV d" d# � D� dX g t j d% � f g
} x� t | � D]� \ } \ } } y t j | � } WnW t
k
r�| j d& | | | f � Yn, t k
r| j d' | | | f � Yn X| j
| | � q�Wd d( l m } m } m }
x� t d � D]� } d) d dY dZ d[ d+ g d, } d } xE t d- � D]7 } | d | � � d) | }
| |
7} | j |
� q�W|
| � | | � } | j
| | � t j | � � qQWd S)\Nr )�
float_infoc s d \ } } x� | D]y } t j | � \ } } t t j | � � � | � } } | | k rt | | | K} | } n | | | K} | | 7} q Wt t t t | � � � d � � | � } | d k rd | 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 r^ )r r ) r# r� rf �ldexpr% �len�binr �bool)�iterableZtmantZtexpr
r r8 �tail�h)�etiny�mant_digr r �msum s
# - .
z MathTests.testFsum.<locals>.msumg g}Ô%�I�Tg �?g0��.�++g�d~�QJg @�5 g �?�6 rs g $@g (@g @g @c S s g | ] } d | � q S)g �?r )�.0r r r r �
<listcomp>9 s z&MathTests.testFsum.<locals>.<listcomp>r i� z0x1.df11f45f4e61ap+2c S s g | ] } d | | � q S)g �?g �r )r� r r r r r� ; s z-0x1.62a2af1bd3624p-1c S s& g | ] } d | d d | � q S)g333333�?r r )r� r! r r r r� = s i� g333333�?g ��7y�ACg��ؗ�Ҝ<g��7y�ACc S s2 g | ]( } d | d | d d | d � q S)g @�2 �4 r )r� r r r r r� A 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� ���;rr �� 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# Zfsumrv rG rP rW r� r� r� �append)r? r� r� Z test_valuesr! �valsr �actualr� r� r� �j�s�vr )r� r� r �testFsum sT
#
zMathTests.testFsumc C sQ t j } | j | d d � d � | j | d d � d � | j | d d � d � | j | d d � d � | j | d d � d � | j | d d � d � | j | d d � d � | j | d d � d � | j | d d � d � | j | d d � 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 | | | � | � | j t | d d � | j t | d d � | j | t d � t d � � 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@rM rM rM i��i����)r# �gcdrW rN rO r= )r? r� �cr
�y�a�br r r �testGcd^ sP
zMathTests.testGcdc C s | j t t j � | j d t j d d � d � | j d t j d d � d � | j t j t t � t � | j t j t t � t � | j t j t t � t � | j t j t t � t � | j t j
t j d t � � � | j t j
t j t d � � � d S)
Nz
hypot(0,0)r z
hypot(3,4)r^ r � g �?g @g �) rN rO r# ZhypotrI rW rU rQ rR rS rT )r? r r r � testHypot� s "zMathTests.testHypotc C s | j t t j � | j d t j d d � d � | j d t j d d � d � | j d t j d d � d � | j d t j d d � d � | j t t j d d
� | j t t j d d
� | j t j d d � d � | j t j d d � d � | j t j t d � t � | j t j t d � t � | j t j
t j t d � � � x�d d d d g D]�} | j t j t | � t � | j t j t | � t � | j t j d | � d � | j t j d | � d � | j t j d | � d � | j t j d! | � d" � | j t j
t j t | � � � | j t t j d | � | j t t j d# | � | j t j d | � d � | j t j d$ | � d% � | j t j t | � t � | j t j t | � t � | j t j
t j t | � � � qfWd 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 � �� rr r� r �( rM rM r{ g �i���g �i���g �i+���i�� l d( l F�x:^V l �O�z��p�c g �g �g �g �g �g �g �) rN rO r# r� rI rv rW rQ rR rS rT rU )r? r r r r � testLdexp� s6 " # zMathTests.testLdexpc C sh | j t t j � | j d t j d t j � d � | j d t j d � d � | j d t j t j � d � | j d t j d d � d � | j d
t j d d � d � | j d
t j d d � d � | j d t j d � d � | j t t j d � | j t t j d � | j t t j t � | j t j t � t � | j
t j t j t � � � d S)Nzlog(1/e)r zlog(1)r zlog(e)z log(32,2)r� r
r� zlog(10**40, 10)rr r� zlog(10**40, 10**20)r z
log(10**1000)i� gO�+��@g �?rM l �O�z��p�c l �O�z��p�c l F�x:^V l� 3qvN�h�iKa��K�P�" �q�� ��!��a�0mvK0<tI D�[��7� fj �wld: ,"�P�S�^
g�2�$�({P _R�U��vD9�{�)P�Anf�YD8�5 PgH�#�q�EgS�d27�}��W>�Y� �y�{A\< �z�kZ � �p�+�*�eX9� �E��i�/�\;��:�8bG�((ema!� �6�o�}Z s'%COd�
�YP]��[BT�vWL? F�Ch_"!�x"'� �j�j�4�x�o�CD<Y#!*�vYyv�#�e�JU>�+i[�9*-Jc+`�S;?lcy g ��l� 3qvN�h�iKa��K�P�" �q�� ��!��a�0mvK0<tI D�[��7� fj �wld: ,"�P�S�^
g�2�$�({P _R�U��vD9�{�)P�Anf�YD8�5 PgH�#�q�EgS�d27�}��W>�Y� �y�{A\< �z�kZ � �p�+�*�eX9� �E��i�/�\;��:�8bG�((ema!� �6�o�}Z s'%COd�
�YP]��[BT�vWL? F�Ch_"!�x"'� �j�j�4�x�o�CD<Y#!*�vYyv�#�e�JU>�+i[�9*-Jc+`�S;?lcy l"��� 3qvN�h�iKa��K�P�" �q�� ��!��a�0mvK0<tI D�[��7� fj �wld: ,"�P�S�^
g�2�$�({P _R�U��vD9�{�)P�Anf�YD8�5 PgH�#�q�EgS�d27�}��W>�Y� �y�{A\< �z�kZ � �p�+�*�eX9� �E��i�/�\;��:�8bG�((ema!� �6�o�}Z s'%COd�
�YP]��[BT�vWL? F�Ch_"!�x"'� �j�j�4�x�o�CD<Y#!*�vYyv�#�e�JU>�+i[�9*-Jc+`�S;?lcy )
rN rO r# �logrI rK rP rR rW rQ rS rT rU )r? r r r �testLog� s # zMathTests.testLogc C sE | j t t j � d } | j t j | � t j t | � � � d S)Nr
�Z l )rN rO r# Zlog1pZassertAlmostEqualr. )r? r r r r � testLog1p� s zMathTests.testLog1pc C s� | j t t j � | j t j d � d � | j t j d � d � | j t j d � d � | j t j d � d � | j t j d � d
� | j t j d � d � | j t t j d � | j t t j t � | j t j t j t � � � d S)Nr g r
g �?r g @i� g ��@i g �@i� g @�@g �?lE lE l� g ��)
rN rO r# �log2rW rP rR rS rT rU )r? r r r �testLog2� s zMathTests.testLog2rr r� c C sL d d � t d d � D� } d d � t d d � D� } | j | | � d S)Nc S s+ g | ]! } t j t j d | � � � q S)g �?)r# r� r� )r� r r r r r� � s z+MathTests.testLog2Exact.<locals>.<listcomp>i2 i c S s g | ] } t | � � q Sr )r. )r� r r r r r� � s i���i���)r rW )r? r� r r r r �
testLog2Exact� s zMathTests.testLog2Exactc C s | j t t j � | j d t j d � d � | j d t j d � d � | j d t j d � d � | j d t j d
� d
� | j t t j d � | j t t j d � | j t t j t � | j t j t � t � | j
t j t j t � � � d S)Nz
log10(0.1)g�������?r zlog10(1)r z log10(10)rr zlog10(10**1000)i� g @�@g �?rM l� 3qvN�h�iKa��K�P�" �q�� ��!��a�0mvK0<tI D�[��7� fj �wld: ,"�P�S�^
g�2�$�({P _R�U��vD9�{�)P�Anf�YD8�5 PgH�#�q�EgS�d27�}��W>�Y� �y�{A\< �z�kZ � �p�+�*�eX9� �E��i�/�\;��:�8bG�((ema!� �6�o�}Z s'%COd�
�YP]��[BT�vWL? F�Ch_"!�x"'� �j�j�4�x�o�CD<Y#!*�vYyv�#�e�JU>�+i[�9*-Jc+`�S;?lcy g ��l� 3qvN�h�iKa��K�P�" �q�� ��!��a�0mvK0<tI D�[��7� fj �wld: ,"�P�S�^
g�2�$�({P _R�U��vD9�{�)P�Anf�YD8�5 PgH�#�q�EgS�d27�}��W>�Y� �y�{A\< �z�kZ � �p�+�*�eX9� �E��i�/�\;��:�8bG�((ema!� �6�o�}Z s'%COd�
�YP]��[BT�vWL? F�Ch_"!�x"'� �j�j�4�x�o�CD<Y#!*�vYyv�#�e�JU>�+i[�9*-Jc+`�S;?lcy l"��� 3qvN�h�iKa��K�P�" �q�� ��!��a�0mvK0<tI D�[��7� fj �wld: ,"�P�S�^
g�2�$�({P _R�U��vD9�{�)P�Anf�YD8�5 PgH�#�q�EgS�d27�}��W>�Y� �y�{A\< �z�kZ � �p�+�*�eX9� �E��i�/�\;��:�8bG�((ema!� �6�o�}Z s'%COd�
�YP]��[BT�vWL? F�Ch_"!�x"'� �j�j�4�x�o�CD<Y#!*�vYyv�#�e�JU>�+i[�9*-Jc+`�S;?lcy )
rN rO r# Zlog10rI rP rR rW r� rQ rS rT rU )r? r r r � testLog10� s zMathTests.testLog10c s� � j t t j � � f d d � } | d t j d � d � | d t j d � d � � j t j t � d t f � � j t j t � d t f � t j t � } � j t j | d � � � j t j | d
� � d S)Nc s] | | \ } } \ } } t | | � t k s? t | | � rY � j d | | | f � d S)Nz%s returned %r, expected %r)r rF rG )rH r~ r Zv1Zv2Ze1Ze2)r? 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 �)
rN rO r# ZmodfrW rQ rR rU rS rT )r? r� Zmodf_nanr )r? r �testModf� s zMathTests.testModfc C s� | j t t j � | j d t j d d � d � | j d t j d d � d � | j d t j d d � d � | j d t j d d � d � | j t j t d � t � | j t j t d � t � | j t j d t � d � | j t j d t � d � | j t j t j t
d � � � | j t j t j d t
� � � | j t j t j d t
� � � | j t j d t
� d � | j t j d
t � d
� | j t j d
d � d
� | j t j d
d � d
� | j t j d
d
� d
� | j t j d
d
� d � | j t j d
d! � d � | j t t j d
d" � | j t t j d
d# � | j t t j d
d$ � | j t t j d
t � | j t j t j d
t
� � � | j t j t t � t � | j t j t d � t � | j t j t d � t � | j t j t d
� t � | j t j t d
� d � | j t j t d% � d � | j t j t d& � d
� | j t j t d' � d
� | j t j t d( � d
� | j t j t t � d
� | j t j t j t t
� � � | j t j d) t � d
� | j t j d* d � d+ � | j t j d, d � d
� | j t j d- d
� d
� | j t j d. d
� d � | j t j d/ d0 � d � | j t t j d1 d2 � | j t t j d3 d4 � | j t t j d5 d6 � | j t t j d7 t � | j t j t j d8 t
� � � | j t j t t � t � | j t j t d � t � | j t j t d � t � | j t j t d
� t � | j t j t d
� d � | j t j t d9 � d � | j t j t d: � d
� | j t j t d; � d
� | j t j t d< � d= � | j t j t t � d
� | j t j t j t t
� � � | j t j d> t � d � | j t j d? d � d@ � | j t t j dA d � | j t j dB d
� d � | j t j dC d
� d � | j t j dD dE � d � | j t j dF dG � d � | j t t j dH dI � | j t j dJ dK � dL � | j t j dM t � d � | j t j t j dN t
� � � | j t j d t � d � | j t j d d � d � | j t j d d � d � | j t j d d
� d � | j t j d d
� d � | j t j d dO � d � | j t j d dP � d � | j t j d dQ � d � | j t j d dR � d � | j t j d t � d � | j t j d t
� d � | j t j d d
� d � | j t j dS d
� d � | j t j t
d
� d � | j t j d dT � d � | j t j dU dV � d � | j t j t
dW � d � | j t t j dX d � | j t t j dY dZ � | j t j d t � d
� | j t j d t � d
� | j t j d t � t � | j t j d t � t � | j t j d[ t � t � | j t j d\ t � t � | j t j d] t � d
� | j t j d^ t � d
� | j t j d t � t � | j t j d t � t � | j t j d t � d
� | j t j d t � d
� | j t j d_ t � d
� | j t j d` t � d
� | j t j da t � t � | j t j db t � t � | j d t j dc d � dd � | j d t j de d
� d � | j d t j df d � dg � | j d t j dh d
� d � | j d t j di dj � d � | j d t j dk dl � dm � | j d t j dn do � d � | j d t j dp dq � dr � | j t t j ds dt � | j t t j du d � 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 �?rM 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 �) rN rO r# �powrI rW rQ rR rS rT rU rP )r? r r r �testPow s� """ " " " " " zMathTests.testPowc C s} | j t t j � | j d t j d � t j � | j d t j d � t j d � | j d t j d � t j d � d S)
Nz radians(180)� z radians(90)r� r
z radians(-45)�- r i���)rN rO r# ZradiansrI rJ )r? r r r � testRadians� s #zMathTests.testRadiansc C s | j t t j � | j d t j d � d � | j d t j t j d � d � | j d t j t j d � d � yB | j t j t j t � � � | j t j t j t � � � Wn> t
k
r� | j t
t j t � | j t
t j t � Yn X| j t j t j t � � � d S)Nzsin(0)r z sin(pi/2)r
r z
sin(-pi/2)rM ) rN rO r# ZsinrI rJ rS rT rQ rR rP rU )r? r r r �testSin� s #$ #
zMathTests.testSinc C s� | j t t j � | j d t j d � d � | j d t j d � d t j d � d d � | j d t j d � t j d � d � | j t j t � t � | j t j t � t � | j t j
t j t � � � d S) Nzsinh(0)r zsinh(1)**2-cosh(1)**2r r
zsinh(1)+sinh(-1)rM rM ) rN rO r# ZsinhrI rm rW rQ rR rS rT rU )r? r r r �testSinh� s 1)zMathTests.testSinhc C s� | j t t j � | j d t j d � d � | j d t j d � d � | j d t j d � d � | j t j t � t � | j t t j t � | j t j
t j t � � � d S)Nzsqrt(0)r zsqrt(1)r zsqrt(4)r r
) rN rO r# �sqrtrI rW rQ rP rR rS rT rU )r? r r r �testSqrt� s zMathTests.testSqrtc
C s | j t t j � | j d t j d � d � | j d t j t j d � d � | j d t j t j d � d � yB | j t j t j t � � � | j t j t j t � � � Wn4 | j t
t j t � | j t
t j t � Yn X| j t j t j t � � � d S)Nztan(0)r z tan(pi/4)r r z
tan(-pi/4)rM ) rN rO r# ZtanrI rJ rS rT rQ rR rP rU )r? r r r �testTan� s #$ #zMathTests.testTanc C s� | j t t j � | j d t j d � d � | j d t j d � t j d � d � | j d t j t � d � | j d t j t � d � | j t j t j t � � � d S) Nztanh(0)r ztanh(1)+tanh(-1)r z tanh(inf)z
tanh(-inf)rM rM )
rN rO r# �tanhrI rQ rR rS rT rU )r? r r r �testTanh� s ) zMathTests.testTanh�TANH_PRESERVES_ZERO_SIGNr z,system tanh() function doesn't copy the signc C sN | j t j d � d � | j t j d t j d � � t j d d � � d S)Ng g �?g �g �g �g �)rW r# r� rj )r? r r r � testTanhSign� s zMathTests.testTanhSignc C s� | j t j d � d � | j t j d � d � | j t t j d � � t � | j t t j d � � t � | j t j d � d � | j t j d � d � | j t j d � d � | j t j d � d � | j t j d � d � | j t j d � d � Gd d � d t � } Gd
d � d t � } | j t j | � � d � | j t t j � | j t t j d d
� | j t t j | � � d S)Nr g �?g����?g� !��?r g-���?Y@rs c @ s e Z d Z d d � Z d S)z'MathTests.test_trunc.<locals>.TestTruncc S s d S)Nr� r )r? r r r � __trunc__� s z1MathTests.test_trunc.<locals>.TestTrunc.__trunc__N)rB rC rD r� r r r r � TestTrunc� s r� c @ s e Z d Z d S)z)MathTests.test_trunc.<locals>.TestNoTruncN)rB rC rD r r r r � TestNoTrunc� s r� r� r
rM rM g ��rM g�����rM g� !��r g-���?Y�i����)rW r# Ztruncrg rf �objectrN rO )r? r� r� r r r �
test_trunc� s zMathTests.test_truncc C s� | j t j d � � | j t j d � � | j t j d � � | j t j d � � | j t j t d � � � | j t j t d � � � | j t j t d � � � d S)Ng g �?r r z-infg �g �)rS r# Zisfinite� assertFalser. )r? r r r � testIsfinite� s zMathTests.testIsfinitec C s� | j t j t d � � � | j t j t d � d � � | j t j t d � � � | j t j d � � | j t j d � � d S)Nr r g g �?)rS r# rT r. r� )r? r r r � testIsnan� s
zMathTests.testIsnanc C s� | j t j t d � � � | j t j t d � � � | j t j d � � | j t j d � � | j t j t d � � � | j t j d � � | j t j d � � d S)Nr z-infg �r g g �?g �)rS r# r$ r. r� )r? r r r � testIsinf s zMathTests.testIsinfc C s | j t j t j � � d S)N)rS r# rT r )r? r r r �test_nan_constant
s zMathTests.test_nan_constantc C sc | j t j t j � � | j t j d � | j t j t d � � | j t j t d � � d S)Ng r z-inf)rS r# r$ r Z
assertGreaterrW r. )r? r r r �test_inf_constant s zMathTests.test_inf_constantzrequires verbose modec C s� y t j d � } Wn | j d � Yn X| d k rD | j d � y t j d � } Wn t k
rk Yn X| j d � y t j d � } Wn t k
r� Yn X| j d � 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# r8 rG rv r� rP )r? r
r r r �test_exceptions s
zMathTests.test_exceptionsc C s* x#t t � D]\ } } } } } } } | d k s
| d k s
| rI q
| d
k rX q
t t | � } y | | � } Wn� t k
r� }
z3 d d |
j d | | | f } | j | � WYd d }
~
Xn6 t k
rd d | | | f } | j | � Yn X| j d | | | f | | � q
Wd S) Ng �rect�polarz"Unexpected ValueError: %s
zin test %s:%s(%r)
r z Unexpected OverflowError in ztest %s:%s(%r)
z %s:%s(%r))r� r� ) r< � test_file�getattrr# rP rd rG rv rI ) r? r4 r5 �arZaiZerZeir9 �funcr~ �exc�messager r r �
test_testfile= s" (
zMathTests.test_testfilec C s d } g } x�t t � D]�\ } } } } } t t | � } d | k sU d | k r^ d } n d | k rp d } y | | � } Wn/ t k
r� d } Yn t k
r� d } Yn Xd }
t | t � r�t | t � r�t j | � r� t j | � r� q t j | � r�t j | � r�| d k rAt | | d d d
d �}
n3 | d k rbt
| | d � }
n t
| | d
� }
|
d k r�q t | t � r�t | t � r�| | k r�q | j | | | | | � } |
d k r�| d j |
� 7} | j
| � q W| r| j d d j | � � d S)Nz${}:{}({!r}): expected {!r}, got {!r}Zinvalidzdivide-by-zerorP Zoverflowrv Zlgammar& g�+�����<r' Zerfci� r z ({})z Failures in test_mtestfile:
z
)r: �math_testcasesr� r# rP rv �
isinstancer. rT r) r �strr r� rG �join) r? Zfail_fmtZfailuresr4 r5 r6 r r9 r� r Zaccuracy_failureZfail_msgr r r �test_mtestfileT sL "
zMathTests.test_mtestfileN)>rB rC rD rI rL rV rX rY rZ r\ r] r_ ri r rk rl rn ro rp rq ru r Z cpython_onlyrw 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 rE � sj
= " Y + ! �
rE c @ s� e Z d Z e j Z d d � Z d d � Z d d � Z d d � Z d d
� Z d d � Z
d
d � Z d d � Z d d � Z
d d � Z d d � Z d d � Z d d � Z d d � Z d d � Z d S) � IsCloseTestsc O s3 | j | j | | | | � d d | | f �d S)N�msgz%s and %s should be close!)rS �isclose)r? r� r� rd �kwargsr r r �
assertIsClose� s zIsCloseTests.assertIsClosec O s3 | j | j | | | | � d d | | f �d S)Nr� z %s and %s should not be close!)r� r� )r? r� r� rd r� r r r �assertIsNotClose� s z IsCloseTests.assertIsNotClosec O s1 x* | D]" \ } } | j | | | | � q Wd S)N)r� )r? �examplesrd r� r� r� r r r �assertAllClose� s zIsCloseTests.assertAllClosec O s1 x* | D]" \ } } | j | | | | � q Wd S)N)r� )r? r� rd r� r� r� r r r �assertAllNotClose� s z IsCloseTests.assertAllNotClosec C sd | j t � � | j d d d d �Wd QRX| j t � �! | j d d d d d d �Wd QRXd S)Nr �rel_tolg0��.�++�abs_tolg _�Bg0��.�+�g _��)rN rP r� )r? r r r �test_negative_tolerances� s z%IsCloseTests.test_negative_tolerancesc C s5 d
d d d
d d g } | j | d d d d �d S)N� @�N���@i�;����90 � � �@� �NF r� r� )r� r� )r� r� )r r )r r � �)r r )r r )r� )r? Zidentical_examplesr r r �test_identical� s zIsCloseTests.test_identicalc C s9 d
d
d g } | j | d d �| j | d d �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� )r? Z eight_decimal_places_examplesr r r �test_eight_decimal_places� s
z&IsCloseTests.test_eight_decimal_placesc C s9 d d
d g } | j | d d �| j | d d �d S)
N��&� .>� gu�?j�/� r� g�������?r� g:�0�yE>)r r ��&� .�)r r �u�?j�/ʠ)r r )r� r� )r? Znear_zero_examplesr r r �test_near_zero� s
zIsCloseTests.test_near_zeroc C sP | j t t � | j t t d d �| j t t � | j t t d d �d S)Nr� g )r� rQ rR )r? r r r �test_identical_infinite� s z$IsCloseTests.test_identical_infinitec C sw t t f t d f d t f t t f t t f t t f t d f d t f t d f d t f g
} | j | d d �d S)Ng0��.�++g �?g������r� g�������?)rU rQ rR r� )r? Znot_close_examplesr r r �test_inf_ninf_nan� s z IsCloseTests.test_inf_ninf_nanc C sH d
d
d g } | j | d d �d d d g } | j | d d �d S)N� �?g333333 @gY�� n�r� g � �?������?�Zb���ti�Tb���ti)r r �333333 ��333333 �)r r �Y�� n���Y�� n��)r r )r r )r r )r r )r� r� )r? Z zero_tolerance_close_examplesZ!zero_tolerance_not_close_examplesr r r �test_zero_tolerance� s z IsCloseTests.test_zero_tolerancec C s | j d d g d d �d S)N� rr r� g�������?)r rr )rr r )r� )r? r r r �test_asymmetry� s zIsCloseTests.test_asymmetryc C s6 d d g } | j | d d �| j | d d �d S)
N��� ���[��[r� g:�0�yE>g��&� .>)r! r" )r# r$ )r� r� )r? Zinteger_examplesr r r �
test_integers� s zIsCloseTests.test_integersc C s d d l m } | d � | d � f | d � | d � f | d � | d � f g } | j | d d
�| j | d d �d S) Nr )�Decimalz
1.00000001z1.0z1.00000001e-20z1.0e-20z1.00000001e-100z1.0e-100r� g:�0�yE>g��&� .>)�decimalr& r� r� )r? r& Zdecimal_examplesr r r �
test_decimals� s zIsCloseTests.test_decimalsc C s\ d d l m } | d d � d | d � f g } | j | d d �| j | d d �d S)Nr )�Fractionr i �r� g:�0�yE>g��&� .>)Z fractionsr) r� r� )r? r) Zfraction_examplesr r r �test_fractions� s "zIsCloseTests.test_fractionsN)rB rC rD r# r� r� r� r� r� r� r r
r r r r r r% r( r* r r r r r� � s
r� c C si d d l m } t j � } | j t j t � � | j t j t � � | j | d � � t | � d S)Nr )� DocFileSuitez ieee754.txt) Zdoctestr+ r� Z TestSuiteZaddTestZ makeSuiterE r� r )r+ Zsuiter r r � test_main s r, ).Z test.supportr r r Ztestr r� r# �os�platformr� r r� rF r. rU rQ rR r
r� r� rB �argv�file�__file__�path�dirname�curdirZtest_dirr� r� r� r r r r r" r) r: r< r� r= ZTestCaserE r� r, r r r r �<module> sJ ! � � � �x