ó <Zc@soddlmZddlmZejjZd„Zd„Zd„Zd„Zd„Z d„Z d „Z d S( iÿÿÿÿ(tmp(tlibmpc Cspdt_tjtjƒ}tjdtjƒ¤tjddddƒ}gdd}}}xp|tj||7}|d7}|j|ƒ|j |ƒ\}}||kr²Pn|d krat d ƒ‚qaqaWWdQXtj d tj |ƒƒ}t |tjdd ƒ}||ks t‚tjd „dtjgddddƒ}t ||ƒ}||kslt‚dS(Niitmethodtlevintvarianttuiiièsiteration limit exceededgÍÌÌÌÌÌì?icSs d||S(Ni((tn((s6lib/python2.7/site-packages/mpmath/tests/test_levin.pytst levin_variant(Rtdpstmpftepst extraprectprecRtonetappendt update_psumt RuntimeErrortexptlogtabstpitAssertionErrortnsumtinf( R tLtStsRtvteterrtw((s6lib/python2.7/site-packages/mpmath/tests/test_levin.pyt test_levin_0 s(     *c Csedt_tjtjƒ}tjdtjƒ›tjddddƒ}gd}}xntj|ƒd}|d7}|j|ƒ|j|ƒ\}}||kr©Pn|d krZt d ƒ‚qZqZWWdQXtj d tj |ƒƒ}t |tj dƒƒ}||kst‚tjd „dtjgddddƒ}t ||ƒ}||ksat‚dS(NiiRRRRiy@ièsiteration limit exceededgÍÌÌÌÌÌì?iþÿÿÿcSs|dS(Niy@y@@((R((s6lib/python2.7/site-packages/mpmath/tests/test_levin.pyR3sRy@@yÀÀ(RR R R R R RRtupdateRRRRtzetaRRR( R RtARRRRRR((s6lib/python2.7/site-packages/mpmath/tests/test_levin.pyt test_levin_1"s(      *c s~dt_tjdƒ‰tjtjƒ}tjdtjƒ”tjddddƒ}d}xnd |tj|ƒˆ| }|j|ƒ\}}|d 7}||kr±Pn|d krbt d ƒ‚qbqbWWdQXtj d tj |ƒƒ}tj ‡fd†dtj gƒ}t||ƒ}||ks8t‚tj‡fd†dtj gddddƒ}||kszt‚dS(Nii iRtsidiRttiiÿÿÿÿiièsiteration limit exceededgÍÌÌÌÌÌì?cstj| ƒd|ˆS(Ni(RR(tx(tz(s6lib/python2.7/site-packages/mpmath/tests/test_levin.pyRGscsd|tj|ƒˆ| S(Niÿÿÿÿ(Rtfac(R(R((s6lib/python2.7/site-packages/mpmath/tests/test_levin.pyRLsR(RR R R R R RR)tstepRRRtquadRRRR( R RRRRRtexactRR((R(s6lib/python2.7/site-packages/mpmath/tests/test_levin.pyt test_levin_27s(     $0c sdt_tjdƒ‰tjtjƒ}tjdtjƒÀtjddddƒ}d\}}x”|ˆ |tjd |ƒtj|ƒtjd|ƒd |7}|d 7}|j|ƒ\}}||krÝPn|d krht d ƒ‚qhqhWWdQXtj d tj |ƒƒ}tj ‡fd†dtj gƒdtjdtjƒ}t||ƒ}||ks|t‚tj‡fd†dtj gddddddtjddggtd ƒD] }d ^qɃ} t|| ƒ}||kst‚dS(NiiiRRRR&iiiièsiteration limit exceededgš™™™™™é?cs"tj| |dˆ|dƒS(Nii(RR(R'(R((s6lib/python2.7/site-packages/mpmath/tests/test_levin.pyR^scs@ˆ |tjd|ƒtj|ƒtjd|ƒd|S(Nii(RR)(R(R((s6lib/python2.7/site-packages/mpmath/tests/test_levin.pyRcsRtworkprecitsteps(ii(RR R R R R RR)t step_psumRRRR+RtsqrtRRRRtxrange( R RRRRRR,RR'R((R(s6lib/python2.7/site-packages/mpmath/tests/test_levin.pyt test_levin_3Os*  F   <cc s¢dt_tjtjƒitjdƒd‰tj‡fd†dtjgddƒdˆ}t|tjƒdks€t ‚WdQXtj d tj tj ƒƒ}tjd „dtjgdd ƒ}t|tj d ƒƒ|ksît ‚d‰‡fd†}tj|dtjgdddgt dƒD] }d^q.ƒ}tjd tjdƒdtjdƒdtjdƒˆƒ}t||ƒ|ksžt ‚dS(Nii iöÿÿÿcs |dˆ S(Ni((R(R((s6lib/python2.7/site-packages/mpmath/tests/test_levin.pyRlsiRtlg»½×Ùß|Û=gš™™™™™é?cSsd|d|S(Niÿÿÿÿi((R((s6lib/python2.7/site-packages/mpmath/tests/test_levin.pyRqsR%iyð?csltjdtjdƒ|ƒtjdtjdƒ|ƒˆ|tjdtjdƒ|ƒtj|ƒS(Niiii(RtrfR R)(R(R((s6lib/python2.7/site-packages/mpmath/tests/test_levin.pyRusiRR/ièiiy@ð?(RR R R R RRRteulerRRRR R2thyp2f1(taR tfR'RR,((R(s6lib/python2.7/site-packages/mpmath/tests/test_levin.pyttest_levin_nsumgs 2%$%@?cCsødt_tjƒ}gdd}}}x||d| tj||7}|d7}|j|ƒ|j|ƒ\}}|tjkr‰Pn|dkr,tdƒ‚q,q,Wtjdtj tjƒƒ}t |tj dd ƒ}||ksôt ‚dS( Niiiiÿÿÿÿièsiteration limit exceededgÍÌÌÌÌÌì?ii ( RR t cohen_altRRRR RRRRRR(tACRRRRRR R((s6lib/python2.7/site-packages/mpmath/tests/test_levin.pyttest_cohen_alt_0zs     cCsùdt_g}tjƒ}d}x |jtjdtjd|dƒƒ|jtjdtjd|ƒ ƒ|d7}|j|ƒ\}}|tjkr¥Pn|dkr$tdƒ‚q$q$Wtj |ƒ}t |dƒ}|dksõt ‚dS(Niiiièsiteration limit exceededgƒ­ò‘þð?gê-™—q=( RR R;RtloggammaRR!R RRRR(R#R<RRRR((s6lib/python2.7/site-packages/mpmath/tests/test_levin.pyttest_cohen_alt_1Šs   )&  N( tmpmathRRtbackendR2R R$R-R3R:R=R?(((s6lib/python2.7/site-packages/mpmath/tests/test_levin.pyts