B <Zkã@s4ddlTddlTddlZdd„Zdd„Zdd„ZdS) é)Ú*NcCsÞdt_tdƒddkst‚tdƒddks.t‚tdƒdd ksBt‚tdƒd d ksVt‚td ƒd d ¡slt‚td ƒd d¡s‚t‚td ƒd d¡s˜t‚td ƒd  d¡s®t‚td ƒd d¡sÄt‚td ƒd d¡sÚt‚dS)Néég@ié@gà?égà¿gÀ?gÀgP?é gR[Ú:XL @gX—úí¤Ãs@gIHb=Ô?g·3äÆçi?g333333Ó?g¯¸[-˜ìÿ?g333333Ó¿gRB3Ϲ à?)ÚmpÚdpsÚmpfÚAssertionErrorZae©r r ú6lib/python3.7/site-packages/mpmath/tests/test_power.pyÚtest_fractional_powsrcCsžt d¡xŽdD]†}x€tdƒD]t}t d|d>d|>¡}t dd¡}||}ttt|ƒt|ƒ|tƒƒ|kspt‚ttt|ƒt|ƒ|t ƒƒ|kst‚qWqWdS)zR Test that inexact integer powers are rounded in the right direction. iÒ)ré5éÈé2ééédN) ÚrandomZseedÚrangeZrandintZto_intZmpf_powZfrom_intÚ round_downr Úround_up)ÚprecÚiÚaÚbZabr r r Útest_pow_integer_directions   "rcCsôdd„}xÞdD]Ô\}}|t_tdƒtdƒ| d}tdƒtdƒ| d}| }| }|t_t}||d|ƒdks~t‚||d|ƒdks’t‚||d|ƒdks¦t‚||d|ƒdksºt‚||d|ƒdksÎt‚||d|ƒdksât‚||d|ƒdksöt‚||d|ƒdks t‚||d |ƒdks"t‚||d |ƒdks8t‚||d |ƒdksNt‚||d |ƒdksdt‚||d |ƒdkszt‚||d |ƒdkst‚||d |ƒdks¦t‚||d |ƒdks¼t‚t}||d|ƒdksÖt‚||d|ƒdksìt‚||d|ƒdkst‚||d|ƒdkst‚||d|ƒdks.t‚||d|ƒdksDt‚||d|ƒdksZt‚||d|ƒdkspt‚||d |ƒdks†t‚||d |ƒdksœt‚||d |ƒdks²t‚||d |ƒdksÈt‚||d |ƒdksÞt‚||d |ƒdksôt‚||d |ƒdks t‚||d |ƒdks t‚t}||d|ƒdks:t‚||d|ƒdksPt‚||d|ƒdksft‚||d|ƒdks|t‚||d|ƒdks’t‚||d|ƒdks¨t‚||d|ƒdks¾t‚||d|ƒdksÔt‚||d |ƒdksêt‚||d |ƒdkst‚||d |ƒdkst‚||d |ƒdks,t‚||d |ƒdksBt‚||d |ƒdksXt‚||d |ƒdksnt‚||d |ƒdks„t‚t}||d|ƒdksžt‚||d|ƒdks´t‚||d|ƒdksÊt‚||d|ƒdksàt‚||d|ƒdksöt‚||d|ƒdks t‚||d|ƒdks"t‚||d|ƒdks8t‚||d |ƒdksNt‚||d |ƒdksdt‚||d |ƒdkszt‚||d |ƒdkst‚||d |ƒdks¦t‚||d |ƒdks¼t‚||d |ƒdksÒt‚||d |ƒdkst‚qWd t_d S) aµ Stress test directed rounding for powers with integer exponents. Basically, we look at the following cases: >>> 1.0001 ** -5 # doctest: +SKIP 0.99950014996500702 >>> 0.9999 ** -5 # doctest: +SKIP 1.000500150035007 >>> (-1.0001) ** -5 # doctest: +SKIP -0.99950014996500702 >>> (-0.9999) ** -5 # doctest: +SKIP -1.000500150035007 >>> 1.0001 ** -6 # doctest: +SKIP 0.99940020994401269 >>> 0.9999 ** -6 # doctest: +SKIP 1.0006002100560125 >>> (-1.0001) ** -6 # doctest: +SKIP 0.99940020994401269 >>> (-0.9999) ** -6 # doctest: +SKIP 1.0006002100560125 etc. We run the tests with values a very small epsilon away from 1: small enough that the result is indistinguishable from 1 when rounded to nearest at the output precision. We check that the result is not erroneously rounded to 1 in cases where the rounding should be done strictly away from 1. cSstt|j|tj|ƒƒS)N)Zmake_mpfZ mpf_pow_intZ_mpf_rr)ÚxÚnÚrr r r ÚpowrCsz'test_pow_epsilon_rounding..powr))ré)iˆi¸ rrééÿÿÿÿééûÿÿÿiúÿÿÿrN) rrr rr rZ round_ceilingZ round_floorr )r!ZinprecZoutprecZpos10001Zpos09999Zneg10001Zneg09999r r r r Útest_pow_epsilon_rounding#sš r')ZmpmathZ mpmath.libmprrrr'r r r r Ús