a gŗ”hkć@s4ddlTddlTddlZdd„Zdd„Zdd„ZdS) é)Ś*NcCsŽdt_tdƒddksJ‚tdƒddks.J‚tdƒdd ksBJ‚tdƒd d ksVJ‚td ƒd d ”slJ‚td ƒd d”s‚J‚td ƒd d”s˜J‚td ƒd  d”s®J‚td ƒd d”sÄJ‚td ƒd d”sŚJ‚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 śE/var/www/auris/lib/python3.9/site-packages/mpmath/tests/test_power.pyŚtest_fractional_powsr cCs–t d”dD]‚}tdƒD]t}t d|d>d|>”}t dd”}||}ttt|ƒt|ƒ|tƒƒ|kslJ‚ttt|ƒt|ƒ|tƒƒ|ksJ‚qqdS)zR Test that inexact integer powers are rounded in the right direction. iŅ)ré5éČé2ééédN) ŚrandomŚseedŚrangeŚrandintZto_intZmpf_powŚfrom_intŚ round_downŚround_up)ŚprecŚiŚaŚbŚabr r r Śtest_pow_integer_directions   "r cCsīdd„}dD]Ō\}}|t_tdƒtdƒ| d}tdƒtdƒ| d}| }| }|t_t}||d|ƒdkszJ‚||d|ƒdksŽJ‚||d|ƒdks¢J‚||d|ƒdks¶J‚||d|ƒdksŹJ‚||d|ƒdksŽJ‚||d|ƒdksņJ‚||d|ƒdksJ‚||d |ƒdksJ‚||d |ƒdks4J‚||d |ƒdksJJ‚||d |ƒdks`J‚||d |ƒdksvJ‚||d |ƒdksŒJ‚||d |ƒdks¢J‚||d |ƒdksøJ‚t}||d|ƒdksŅJ‚||d|ƒdksčJ‚||d|ƒdksžJ‚||d|ƒdksJ‚||d|ƒdks*J‚||d|ƒdks@J‚||d|ƒdksVJ‚||d|ƒdkslJ‚||d |ƒdks‚J‚||d |ƒdks˜J‚||d |ƒdks®J‚||d |ƒdksÄJ‚||d |ƒdksŚJ‚||d |ƒdksšJ‚||d |ƒdksJ‚||d |ƒdksJ‚t}||d|ƒdks6J‚||d|ƒdksLJ‚||d|ƒdksbJ‚||d|ƒdksxJ‚||d|ƒdksŽJ‚||d|ƒdks¤J‚||d|ƒdksŗJ‚||d|ƒdksŠJ‚||d |ƒdksęJ‚||d |ƒdksüJ‚||d |ƒdksJ‚||d |ƒdks(J‚||d |ƒdks>J‚||d |ƒdksTJ‚||d |ƒdksjJ‚||d |ƒdks€J‚t}||d|ƒdksšJ‚||d|ƒdks°J‚||d|ƒdksĘJ‚||d|ƒdksÜJ‚||d|ƒdksņJ‚||d|ƒdksJ‚||d|ƒdksJ‚||d|ƒdks4J‚||d |ƒdksJJ‚||d |ƒdks`J‚||d |ƒdksvJ‚||d |ƒdksŒJ‚||d |ƒdks¢J‚||d |ƒdksøJ‚||d |ƒdksĪJ‚||d |ƒdks J‚q d 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 r r*r r r r Śs