a d?Z+d@dAZ,dBdCZ-dDdEZ.dFdGZ/dHdIZ0dJdKZ1dLdMZ2dNdOZ3dPdQZ4dRdSZ5dTdUZ6dVdWZ7dXdYZ8dZd[Z9d\d]Z:d^d_Z;d`daZdfdgZ?dhdiZ@djdkZAdldmZBdndoZCdS)p)*NcCs0d}|o|j|j|}|o*|j|j|}|S)NT)realaeimag)abepsresr k/Users/vegardjervell/Documents/master/model/venv/lib/python3.9/site-packages/mpmath/tests/test_functions.pympc_aesr zg3.14159265358979323846264338327950288419716939937510582097494459230781640628620899862803482534211706798zg2.71828182845904523536028747135266249775724709369995957496696762772407663035354759457138217852516642743zi0.0174532925199432957692369076848861271344287188854172545609719144017100911460344944368224156963450948221zh0.577215664901532860606512090082402431042159335939923598805767234884867726777664670936947063291746749516zh0.693147180559945309417232121458176568075500134360255254120680009493393621969694715605863326996418687542zg2.30258509299404568401799145468436420760110148862877297603332790096757260967735248023599720508959829834zh0.915965594177219015054603514932384110774149374281672134266498119621763019776254769479356512926115106249zg2.68545200106530644530971483548179569382038229399446295305115234555721885953715200280114117493184769800zm1.28242712910062263687534256886979172776768892732500119206374002174040630885882646112973649195820237439420646zg1.20205690315959428539973816151144999076498629234049888179227155534183820578631309018645587360933525815zj1.61803398874989484820458683436563811772030917980576286213544862270526046281890244970720720418939113748475zh0.261497212847642783755426838608695859051566648261199206192064213924924510897368209714142631434246651052zh0.660161815846869573927812110014555778432623360284733413319448423335405642304495277143760031413839867912cCsdD]}|t_tttksJtttks.Jtttks>Jt tt ksNJt tt ks^Jt ttksnJtttks~JtttksJtttksJtttksJ|dkrtttksJtttksJqdt_tdksJtdksJtdksJtdksJdS)N)  %Pd2rr )mpdpspimpftpietedegreetdegreeZeulerteulerln2tln2Zln10tln10catalantcatalanZkhinchin tkhinchinZglaisher tglaisherphitphiZmertenstmertensZ twinprime ttwinprimeprecr r r test_constants-s(  r1cCstdD]}tt|||ksJqtddD]T}|t_tdD]@}td|ddd|dd}tt|||ksFJqFq4dt_tdD]`}d D]V}tt||d|ft||fksJtt||d |ft|| fksJqqdS) N Nr,'rrrrr)r3i)rangesqrtrrandomseedrrrandint)ir0Arr r r test_exact_sqrtsCs   $ (rBcCsdD]l}t|}dD]Z}|t_tt|tjtdtjt}tt|tjtdtjt}t||s`Jt ||sJqqt dd}tt t fD]8}tdD]*}t|}t||}t||||ksJqqdt_tdd}tdd}dt_t|d d d ksJt|d d d ksJt|d d d ks(Jt|d d d ks>Jt|d d d ksTJt|d d d ksjJttdtdksJdS)N) rr rr8r r)rrSjirr goK@-rdroundingr nuz7.0503726185518891z2.655253776675949)from_intrrZ mpf_pow_intZmpf_sqrtr0Z round_downZround_upZmpf_ltZmpf_gtr=r>Z round_nearestZ round_ceilingr;Zmpf_randZmpf_mulrr<)r@rrrr0rndr r r test_sqrt_roundingRs4     rTcCsdt_dD] }tt|t|dks Jq tddks皙??r3rr rrCgZd;O?g\( S@rVr?r:y) rrr<rfloatrcmathmathrrxr r r test_float_sqrtps0r_cCsntdddksJtddtdks(Jtddtdks>JtddtdksTJtddtdksjJdS)NrgQ?gQտr rrC)hypotrr r r r test_hypot~s racCstdddD] }tt||||ks Jq tddD]J}|t_td|ddd|dd}tt||||ksrrr?)r@r0rAr r r test_exact_cbrts $rdcCstddksJtdtds&Jtdtdsl_y?:0yEy:0yE>y:0yEy:0yE>?y:0yE?y:0yE>y:0yEy?Ww'&l7ryWw'&l7?gX9v? Fx:^Vg@xDrbgZbti) rrlogrr\r[rldexpmpcinfisnannanrr]r r r test_logs@"& rxcCs:ttdttddD]}|d}tt|t|sBJtt|t|s^Jtt|t|szJtt|t|sJt t|t |sJt t|t |sJqtdt dsJtdt dsJtdt dsJtdt ds6JdS)NrrYrgffffff@y??y ) listr;cosrrr\sintancoshsinhtanhr[)r^tr r r test_trig_hyperb_basicsrcCstdtdksJtdtds*Jtdtds@JtdtdsVJtdtdslJtdtdksJtdtdsJtdtdsJtdtdsJtdtdsJdS)Nrr3Zriih)rzr!rr{r r r r test_degreessrc Csrtdg}t|D]V}tdd}tdd}tdd}tdd}t||}t||}|||fq|S)Nr3r)r=r>r;uniformcomplexappend) Nrr@x1y1x2y2z1z2r r r random_complexess        rcCsdD]@}|t_tdd}|j|jkr>tdddksnJqdt_tdtdD](\}}t|t| ||ds^Jq^t t d d sJd t_t t d d sJdt_dS) N)rrrWrVrrr3r-q=rr) rrrtrrrr=r>rrrr)rrrrr r r test_complex_powerss , "rcCsXdd}tdd}dt_tj}|ddt|D|ddt|Dd t_dS) NcSs|D]\}}t|t|}tt|||s2Jt| t|}tt||| s\Jt|t|}tt|||sJt| t|}tt||| sJqdS)N)rtjr r<)lstrrzr r r test_mpc_sqrts z1test_complex_sqrt_accuracy..test_mpc_sqrtrrrcSs$g|]}tddtddfqS)rr)r=r.0r@r r r z.test_complex_sqrt_accuracy..cSs$g|]}|d|dd|fqS)rU皙?rr rr r r rrr)r=r>rrr;)rrrr r r test_complex_sqrt_accuracys rcCs dt_tdtdsJtddks.JtdtdsDJtddksTJtdt dslJtdd tdsJd D] }|t_d td tsJqdt_td}tttd |sJtttd|sJtttd |sJttd t|sJttdt|s8Jttd t| sTJttt d | srJttt d| sJttt d | sJttd t | sJttdt | sJttd t |sJdS)Nrgffffffrpgd~QJrJzgd~Qrr6)r9Frr5r6rr3rr)rratanrr\rrtru)rZpi2r r r test_atans0rcCspdt_tddtdsJtdddtds:JtdddtdsVJtddt dspJtddt dsJtddtdsJtdddksJttdtdsJtt dt dsJttttsJttt tsJtttt s&Jttdts:JtttdsNJttdtsbJtttdsvJtdtdksJtdt tsJtd tdksJtd tdksJtd t t sJtd t tsJttd tdsJttd tds0Jtt d t dsNJtt d t dslJdS) Nrr3rrr rXrrrr)rratan2rrrurvrwr r r r test_atan2,s6rcCsrttddksJttddks(Jttddksr;rrr\r[rrrr)r@r^rr r r test_areal_inversesIs@         rcCsdt_tdddsJtddds2Jtdd d dsLJtd dd sbJtd ddsxJtdd d dsJtddsJtddsJtdddsJtd ddsJtdddsJtd ddsJdS)Nrgh㈵>ig?绽|=d( r3rprrgh㈵gg|۽rrrqgC6W_G@l Fx:^VgC6W_G)rrrrrrr r r r test_invhyperb_inaccuracyosrcCs:ttdttddD]}ttdttddD]}t||dd}tt|t|slJtt|t|sJtt|t|sJt t|t |sJt t|t |sJt t|t |sJt t|t |sJt t|t |s:Jq:qdS)Nrrrg333333@y{Gz?)ryr;rrjrtrr[rrrzr{r|r~r}r)r^yrr r r test_complex_functions~s rcCsdt_dt_tdD]\}}tt||s2Jt|t|sJJt |t |sbJt |t |szJt |jt |ddsJt |jt |ddsJt d}qtdddD]l}tdddD]Z}d td|d |d|}tt |}||sJtt |}||sJqqt d}d }tdd dD]v}tdd dD]b}d d|td |d|}tt|}|||sJtt|}|||s^Jq^qNdS)Nrrr)Zrel_epsr3rr g?皙?gV瞯JtddtksRJtddtksdJtdddkrt tddt sJdS) Nrg@rgffffff@rrVr9g@) lnrrlog10degreesr!radianspowerrfmodrrr r r r test_aliasess rcCstddksJtdts"Jtttds8Jtt t dsRJtddksbJttdtstJtttdsJttdtsJtttdsJttttsJttdksJtt tsJtttsJtddksJtddksJtddks&Jttdks8Jtt dksLJttts^JtttkspJtdtt ksJtdtdtt dsJdS)Nr rrXrr3r) argrrrrvrrwrusignr<r r r r test_arg_signs,rcCsdt_tddt_dS)Nr6ir)rrrrr r r r test_misc_bugssrc Cstdtdtdtdtdtdtdtdtd td td g ksLJtd d tdtdtdtdtdtdtdtdtdtdg ksJtdddtdtdtdtdtdtdtdtdtdtdg ksJtdd d!td"td#td$td tdtdtdtdtd%g ks6Jtd&dd'tdgksRJtdgksdJtd(d)gksxJtd*d+d,gksJdS)-Nrz0.0z1.02.03.04.05.06.07.08.09.0rCz-5.0z-4.0z-3.0z-2.0z-1.0rr3rUz0.10000000000000001z0.20000000000000001z0.30000000000000004z0.40000000000000002z0.5z0.60000000000000009z0.70000000000000007z0.80000000000000004z0.90000000000000002rrXz17.0z14.0z11.0z-7.0rgr6rgGzgGz @gHz)arangerr r r r test_aranges2     rc Cstdddtdtdtdtdtdtd td gks>Jtdddd d tdtd tdtdtdtdtdgksJtdddtdgksJdS)NrrLrrz3.166666666666667z4.3333333333333339z5.5z6.666666666666667z7.8333333333333339rr)Zendpointrrrrrrr3)Zlinspacerr r r r test_linspace s   rcCsdt_tdddD] }t||||tsJqtddttddsVJtd d}tddd d gD]N}t |d }t|}tjd7_t ||}tjd8_||tstJqtd t_t d dd D]8}d|td}|||}t|}t ||tsJqdt_dS)NrrrrUrrVr rr3g@gj@y?ri-gMbP?r) rrrrcrrrr<rrtpowr;r )rZ one_thirdr1Zr2rPwrrr r r test_float_cbrts&"    rc Csdt_tdtdd}|||}tt|d}||ksBJtddD]"}||}tt|||ksLJqLdt_tdd d D]~}tdtdd fD]\}t||}t|td|}| |sJt|| }t|td |}| |sJqqd D]h}|t_ d D]V}d }tdD]B}|t_ t }d|t_ ||}|t_ t||}||ks$Jq$qqdt_tddD]6}tt t}t|| t|td|sJtt||t|td|sJtdd t t}t||}tjd7_t|td|}tjd8_| |sFJt||tsXJt|| }tjd7_t|td |}tjd8_| |sJt||tsJqdt_tdddksJtdddksJtdddksJttdtksJttdtks$Jttdtks8JttddksLJttddks`JttddkstJtt dt ksJtt ddksJtt dt ksJtttdsJtttdsJtttdsJtttdsJtddtddksJtddddks2JtddddksHJtddddks^JtddddkstJtddddksJtddddksJdS)Nrr3rr7r rr(ir)ZnearestZupZdown)rrXr rCrrrrg?r:rXy@yir)rrr=r>r?nthrootrr;rrrOr0Zrandrr rrurvrwroot) rprrPrrSr0r@rr r r test_root&s          $$  rcCsdD]B}|t_ttdd}|tdddttddsJqd t_td dtdd dttdszJtd ddtdtd sJd t_tdddttdtdsJdt_dS)N)rrz-1e-20rrrr3rrz-1e-3gz-1e-6rz -1e100000000Z 1e25000000r)rrrrrrr<)rrr r r test_issue_136rs,,$(rcCs4dt_td}td}tt|dt|ks0Jttd|tt|t|ksTJtt|dt|ksnJtt|dt|ksJttd|t|ksJttd|tdt|ksJt t|dt |ksJt t|dt |ksJt td| tt |sJt td| tdtt |s@Jt t|d\}}|t |ksdJ|t |ksvJt td|\}}| tt |sJ| tdtt |sJtt|d\}}|t|ksJ|t|ksJttd|\}}|t|ksJ|tdt|ks0JdS)Nrg333333?g?r)rrrrjrtrzr{r}r~cospisinpirrZ cospi_sinpiZcos_sin)r^rcsr r r test_mpcfun_real_imags2$ "(rcCsdt_tdd}t dd}d|}d|}dt_t|dksFJt|dddksZJt|dddksnJt|dddksJt|dddksJt|dksJt|dddksJt|dddksJt|dddksJt|dddksJttttfD]}||| ksJ||dd|ks0J||dd|ksFJ||| ksZJ||dd|kspJ||dd|ksJqt t t t fD]}||| ksJ||| ksJ||dd|ksJ||dd|ksJ||dd|ksJ||dd|ksJqt || ks2Jt || ksFJt |dd|ks\Jt |dd|ksrJt |dd|ksJt |dd|ksJt|dksJt|dksJt|dddksJt|dddksJt|dddksJt|dddksJdS) Nrrrr3rrrNf)rrrrjrzr{rrrrr|r~rrr})rrrrMrr r r test_perturbation_roundingsT rcCsDtddksJtddks Jtddks0Jtddks@JdS)Ng @r ry @@y@@y@@)floorceilr r r r test_integer_partssrcCstddksJtddks Jtddks0Jtddks@JtddksPJtddks`JtddkspJtdd ksJtddksJdS) N3r @@rCy?@r3rry@)fabsreZimZconjr conjugater r r r test_complex_partssrc CsHtddksJtddks Jtddks0Jtddks@JtddksPJtddks`JtddkspJtddksJtddksJtddksJtddksJtddksJtddksJtddksJtd tddsJtd }|jtd tjsJ|jdks.Jtd }|jtd tjsRJ|jdksbJd t _ t d t d }t d t d }t dt d }t dt d }t dt d }t dt d }t dt d }t dt d } dt _ d} t|| t s Jt|| ts$Jt|| tsdksDnJtt tt krddksjnJtdds|JtddsJtddksJtddsJdS) Nrr3rgF?ry?\arng)˿)ZsincZsincpirurr r r r test_sincs""&rcCsdt_ddtddDgdks&Jtdds8Jtd d sJJtd d s\Jttdd dkstJdt_td}|ddksJdt_tttksJtddksJdS)NrcSsg|] }t|qSr ) fibonaccirrPr r r r#rz"test_fibonacci..rr)rCrXrrr3rr3r3rr rCr8rGr"rg3?ryɽ傀C,r6gaڃe@krzv6.24499112864607e+2089876402499787337692720892375554168224592399182109535392875613974104853496745963277658556235103534i4r7rlE! y@r)rrr;Zfibrstrrrurkr r r test_fibonacci!srcCs,dt_ttdddtdddks(JdS)Nrr3r)rrgG^Z9)rrabsrjrr r r r test_call_with_dps0srcCsbdt_tddksJttdks&Jtt dks8JtttsHJttdddks^JdS)Nrrr3rru0)rrrrurvrwrtr r r r test_tanh4s rcCsdt_tddksJtdds(Jtdds:JtdtksJJtdt ks\JtttslJttdts~JttdtsJt t d }tt| sJtt |sJtt td| sJtt td| sJtt td|sJtt dt|s,Jtt dt|sFJtt dt|s`Jtt t d|s|Jtt t d|sJtt t d| sJtt dt | sJtt dt | sJtt dt | sJdS) NrrrVg z?rg zr3rr) rrrrrurvrwrrrrrt)Zjpi2r r r test_atanh<s0 rcCsldt_tddksJtdtdds0Jtttks@JtddsRJtdddshJdS) Nrrr r3rpr _Bgo?)rrexpm1rrjrur r r r test_expm1Ws rcCshdt_tddksJtdtds,Jtttks.)r3rrr rrCr3rrr r3rEr3rGr3r3rrr3r3r3r3r3rCr3r r3rr3cSsg|]}t|dqS)rr rr r r rr)r3r:rr3rr3r r3rr3rCr3r3r3rr3rr3r r3r3r3rEr3r3r3rGr3r3r3r3cSsg|]}t|tqSr r rrr r r rrrr3rrrr rCcSsg|]}t|t qSr rrr r r rriXi@ir/r8rg?3OP@irrng @g#@rgcy@)rrr;rr r<r)rQr r r test_cyclotomics RV r)DZ mpmath.libmpZmpmathr=timer\r[rr rr r"r#r%r&r(r)r*Ztaperyr,r-r.r1rBrTr_rardrlrmrxrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrr rr r r r s   %  & !L , 6