o ¨"`kã@s4ddlTddlTddlZdd„Zdd„Zdd„ZdS) é)Ú*NcCsÞdt_tdƒddks J‚tdƒddksJ‚tdƒdd ks!J‚tdƒd d ks+J‚td ƒd d ¡s6J‚td ƒd d¡sAJ‚td ƒd d¡sLJ‚td ƒd  d¡sWJ‚td ƒd d¡sbJ‚td ƒd d¡smJ‚dS)Néég@ié@gà?égà¿gÀ?gÀgP?é gR[Ú:XL @gX—úí¤Ãs@gIHb=Ô?g·3äÆçi?g333333Ó?g¯¸[-˜ìÿ?g333333Ó¿gRB3Ϲ à?)ÚmpÚdpsÚmpfZae©r r ú9/usr/lib/python3/dist-packages/mpmath/tests/test_power.pyÚtest_fractional_powsr cCs–t d¡dD]A}tdƒD]:}t d|d>d|>¡}t dd¡}||}ttt|ƒt|ƒ|tƒƒ|ks6J‚ttt|ƒt|ƒ|tƒƒ|ksGJ‚q qdS)zR Test that inexact integer powers are rounded in the right direction. iÒ)ré5éÈé2ééédN) ÚrandomZseedÚrangeZrandintZto_intZmpf_powZfrom_intÚ round_downÚround_up)ÚprecÚiÚaÚbÚabr r r Útest_pow_integer_directions   "$úÿrcCsÖdd„}dD]Þ\}}|t_tdƒtdƒ| d}tdƒtdƒ| d}| }| }|t_t}||d|ƒdks=J‚||d|ƒdksGJ‚||d|ƒdksQJ‚||d|ƒdks[J‚||d|ƒdkseJ‚||d|ƒdksoJ‚||d|ƒdksyJ‚||d|ƒdksƒJ‚||d |ƒdksJ‚||d |ƒdks—J‚||d |ƒdks¡J‚||d |ƒdks«J‚||d |ƒdksµJ‚||d |ƒdks¿J‚||d |ƒdksÉJ‚||d |ƒdksÓJ‚t}||d|ƒdksßJ‚||d|ƒdkséJ‚||d|ƒdksóJ‚||d|ƒdksýJ‚||d|ƒdksJ‚||d|ƒdksJ‚||d|ƒdksJ‚||d|ƒdks)J‚||d |ƒdks4J‚||d |ƒdks?J‚||d |ƒdksJJ‚||d |ƒdksUJ‚||d |ƒdks`J‚||d |ƒdkskJ‚||d |ƒdksvJ‚||d |ƒdksJ‚t}||d|ƒdksŽJ‚||d|ƒdks™J‚||d|ƒdks¤J‚||d|ƒdks¯J‚||d|ƒdksºJ‚||d|ƒdksÅJ‚||d|ƒdksÐJ‚||d|ƒdksÛJ‚||d |ƒdksæJ‚||d |ƒdksñJ‚||d |ƒdksüJ‚||d |ƒdksJ‚||d |ƒdksJ‚||d |ƒdksJ‚||d |ƒdks(J‚||d |ƒdks3J‚t}||d|ƒdks@J‚||d|ƒdksKJ‚||d|ƒdksVJ‚||d|ƒdksaJ‚||d|ƒdkslJ‚||d|ƒdkswJ‚||d|ƒdks‚J‚||d|ƒdksJ‚||d |ƒdks˜J‚||d |ƒdks£J‚||d |ƒdks®J‚||d |ƒdks¹J‚||d |ƒdksÄJ‚||d |ƒdksÏJ‚||d |ƒdksÚJ‚||d |ƒdksåJ‚qd 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 rrZ round_ceilingZ round_floorr )r!ZinprecZoutprecZpos10001Zpos09999Zneg10001Zneg09999r r r r Útest_pow_epsilon_rounding#sš  r')ZmpmathZ mpmath.libmprr rr'r r r r Ús