Your IP : 216.73.216.213
�
g f�� � @ 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 j( � Z) d# d$ � Z* e d k r�e* � n 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 r< | d } n | 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/python34/lib64/python3.4/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 rD | d | j d � � } n | j � sV q n | j d � \ } } | j � \ } } } | j � } | d } | d d � }
| | t | � t | � |
f Vq WWd QXd 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 n | j d � \ } } | j � \ } } } } | j � } | d | d }
} | d d � } | | t | � t | � t |
� t | � | f Vq WWd QXd 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 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 d# d$ � Z d% d& � Z d' d( � Z e
e j e d) � d* d+ � � � Z d, d- � Z d. d/ � Z d0 d1 � Z d2 d3 � Z e
d4 d5 � � Z e
e! j" d6 d7 � d8 d9 � � � Z# d: d; � Z$ d< d= � Z% d> d? � Z& d@ dA � Z' dB dC � Z( dD dE � Z) dF dG � Z* dH dI � Z+ dJ dK � Z, e
e j e- j. dL � dM k dN � dO dP � � � Z/ dQ dR � Z0 dS dT � Z1 dU dV � Z2 dW dX � Z3 e j4 e5 dY � dZ d[ � � Z6 e
d\ d] � � Z7 e
d^ d_ � � Z8 d` S)a� MathTestsc C s7 t | | � t k r3 | j d | | | f � n d S)Nz%s returned %r, expected %r)r �eps�fail)�self�name�valuer 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�
�@)rC r# rD rE )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# ZacosrC rD �
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�?rG ) rH rI r# ZacoshrC rJ � assertEqualrK rL rM rN rO )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)rG ) rH rI r# ZasinrC rD rJ rK rL rM rN rO )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)rG g'�ya64�) rH rI r# ZasinhrC rQ rK rL rM rN rO )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)rG ) rH rI r# �atanrC rD rK rL rM rN rO )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��rG ) rH rI r# rU rC ZatanhrJ rK rL rM rN rO )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.)rG rG 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 �) rH rI r# Zatan2rC rD rL rQ rK rM rN rO )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__E s z-MathTests.testCeil.<locals>.TestCeil.__ceil__N)�__name__�
__module__� __qualname__r[ r r r r �TestCeilD s r_ c @ s e Z d Z d S)z&MathTests.testCeil.<locals>.TestNoCeilN)r\ r] r^ r r r r �
TestNoCeilG s r` zceil(TestCeil())rZ c W s | S)Nr )�argsr r r �<lambda>M s z$MathTests.testCeil.<locals>.<lambda>g �g �rG g ��rG ) rH rI r# ZceilrQ �int�typerC r[ )r@ r_ r` �tr r r �testCeil7 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 rZ g �?g rX g @g @g @i���g �g �g �g �g �g �g �g �) rQ r# �copysignrH rI rK rL rM rN rO r$ r )r@ r r r � testCopysignQ 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)rG ) rH rI r# ZcosrC rD rM rN rK rL rJ rO )r@ r r r �testCoss 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
rG ) rH rI r# �coshrC rQ rK rL rM rN rO )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�)rH rI r# ZdegreesrC rD )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 rG ) rH rI r# r8 rC rE rQ rK rL rM rN rO )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)rG )rH rI r# ZfabsrC )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 t j � | j t t j t
j d � | j t t j d � d S) Nr r g i� g �?g\�n��TrG g �) rQ r# Z factorialr r. r" rH rJ rD �
OverflowError�sys�maxsize)r@ Ztotalr! r r r �
testFactorial� s
# zMathTests.testFactorialc 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)NrZ r )r@ r r r � __floor__� s z0MathTests.testFloor.<locals>.TestFloor.__floor__N)r\ r] r^ rs r r r r � TestFloor� s rt c @ s e Z d Z d S)z(MathTests.testFloor.<locals>.TestNoFloorN)r\ r] r^ r r r r � TestNoFloor� s ru zfloor(TestFloor())rZ c W s | S)Nr )ra r r r rb � s z%MathTests.testFloor.<locals>.<lambda>g �rG g �rG g �����g��Nݯ��g��Nݯ��) rH rI r# ZfloorrQ rc rd rC rs )r@ rt ru re 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)�
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 �) rH rI r# ZfmodrC rM rN rO rJ rK rL rQ )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 s\ | | \ } } \ } } t | | � t k s; | | k rX � j d | | | f � n d S)Nz%s returned %r, expected %r)r r> r? )rA �resultr �mantr8 ZemantZeexp)r@ r r � testfrexp� s " z&MathTests.testFrexp.<locals>.testfrexpz frexp(-1)r � �?zfrexp(0)r zfrexp(1)zfrexp(2)r
rG � �)r~ r )r r )r} r )r} r
)
rH rI r# �frexprQ rK rL rM rN rO )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} n 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 rX )r r ) r# r rc �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 �d g $@g (@g @g @c S s g | ] } d | � q S)g �?r )�.0r r r r �
<listcomp>) 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� 1 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� ���;rx �� 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���)rp r� r� �min_expr r. �fromhex� enumerater# Zfsumro r? rJ rQ 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 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)rX r � g �?g @g �) rH rI r# ZhypotrC rQ rO rK rL rM rN )r@ r r r � testHypotN 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 � �� rx r� r �( rG rG rv 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 �) rH rI r# r� rC ro rQ rK rL rM rN rO )r@ r r r r � testLdexpY 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� zlog(10**40, 10)rx r� zlog(10**40, 10**20)r z
log(10**1000)i� gO�+��@g �?rG 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 )
rH rI r# �logrC rE rJ rL rQ rK rM rN rO )r@ r r r �testLogy s # zMathTests.testLogc C sE | j t t j � d } | j t j | � t j t | � � � d S)Nr
�Z l )rH rI 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 ��)
rH rI r# �log2rQ rJ rL rM rN rO )r@ r r r �testLog2� s zMathTests.testLog2rx 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 rQ )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)rx zlog10(10**1000)i� g @�@g �?rG 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 )
rH rI r# Zlog10rC rJ rL rQ r� rK rM rN rO )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 | | � r\ � j d | | | f � n d S)Nz%s returned %r, expected %r)r r> r? )rA rz 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 �)
rH rI r# ZmodfrQ rK rL rO rM rN )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 �?rG 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 �) rH rI r# �powrC rQ rK rL rM rN rO rJ )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���)rH rI r# ZradiansrC rD )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)rG ) rH rI r# ZsinrC rD rM rN rK rL rJ rO )r@ r r r �testSinc 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)rG rG ) rH rI r# ZsinhrC rj rQ rK rL rM rN rO )r@ r r r �testSinhp 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
) rH rI r# �sqrtrC rQ rK rJ rL rM rN rO )r@ r r r �testSqrty 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)rG ) rH rI r# ZtanrC rD rM rN rK rL rJ rO )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)rG rG )
rH rI r# �tanhrC rK rL rM rN rO )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 �)rQ r# r� rg )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@r� c @ s e Z d Z d d � Z d S)z'MathTests.test_trunc.<locals>.TestTruncc S s d S)N� r )r@ r r r � __trunc__� s z1MathTests.test_trunc.<locals>.TestTrunc.__trunc__N)r\ r] r^ r� r r r r � TestTrunc� s r� c @ s e Z d Z d S)z)MathTests.test_trunc.<locals>.TestNoTruncN)r\ r] r^ r r r r � TestNoTrunc� s r� r� r
rG rG g ��rG g�����rG g� !��r g-���?Y�i����)rQ r# Ztruncrd rc �objectrH rI )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 �)rM 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 �?)rM r# rN 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 �)rM r# r$ r. r� )r@ r r r � testIsinf� s zMathTests.testIsinfzrequires verbose modec C s� y t j d � } Wn | j d � Yn X| d k rG | j d � n y t j d � } Wn t k
rn 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 r? ro r� rJ )r@ r
r r r �test_exceptions� s
zMathTests.test_exceptionsc C s0 x)t t � D]\ } } } } } } } | d k s
| d k s
| rL q
n | d
k r^ q
n 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# rJ ra r? ro rC ) r@ r4 r5 �arZaiZerZeir9 �funcrz �exc�messager r r �
test_testfile� s" (
zMathTests.test_testfilec C s: d } g } xt t � D]�\ } } } } } t t | � } d | k sU d | k r^ d } n d | k rs d } n 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 n t j | � r�t j | � r�| d k rGt | | d d d
d �}
n3 | d k rht
| | d � }
n t
| | d
� }
|
d k r�q q�q�n t | t � r�t | t � r�| | k r�q q�n | j | | | | | � } |
d k r| d j |
� 7} n | j
| � q W| r6| j d d j | � � n d S)Nz${}:{}({!r}): expected {!r}, got {!r}Zinvalidzdivide-by-zerorJ Zoverflowro Zlgammar& g�+�����<r' Zerfci� r z ({})z Failures in test_mtestfile:
z
)r: �math_testcasesr� r# rJ ro �
isinstancer. rN r) r �strr r� r? �join) r@ Zfail_fmtZfailuresr4 r5 r6 r r9 r� r Zaccuracy_failureZfail_msgr r r �test_mtestfile sL "
zMathTests.test_mtestfileN)9r\ r] r^ rC rF rP rR rS rT rV rW rY rf r rh ri rk rl rm rn rr rw ry 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� Z
skipUnlessr r� r� r� r r r r r= � sb
= " Y ! �
r= c C sS d d l m } t j � } | j t j t � � | j | d � � t | � d S)Nr )� DocFileSuitez ieee754.txt)Zdoctestr� r� Z TestSuiteZaddTestZ makeSuiter= r )r� Zsuiter r r � test_mainK s
r� )+Z test.supportr r r Ztestr r� r# �os�platformrp r r� r> r. rO rK rL r
�yr� r\ �argv�file�__file__�path�dirname�curdirZtest_dirr� r� r� r r r r r" r) r: r< ZTestCaser= r� r r r r �<module> sF ! � � � �