o l hk@s4ddlTddlTddlZddZddZddZdS) )*NcCsdt_tdddks JtdddksJtddd ks!Jtdd d ks+Jtd dd s6Jtd ddsAJtd ddsLJtd d dsWJtd ddsbJtd ddsmJdS)Ng@i@g?gg?ggP? gR[:XL @gXs@gIHb=?g3i?g333333?g[-?g333333ӿgRB3Ϲ ?)mpdpsmpfaer r k/var/www/html/construction_image-detection-poc/venv/lib/python3.10/site-packages/mpmath/tests/test_power.pytest_fractional_powsrcCstddD]A}tdD]:}td|d>d|>}tdd}||}ttt|t||t|ks6Jttt|t||t|ksGJq qdS)zR Test that inexact integer powers are rounded in the right direction. i)r52dN) randomseedrangerandintto_intmpf_powfrom_int round_downround_up)preciababr r r test_pow_integer_directions   "$r#cCsdd}dD]\}}|t_tdtd| d}tdtd| d}| }| }|t_t}||d|dks=J||d|dksGJ||d|dksQJ||d|dks[J||d|dkseJ||d|dksoJ||d|dksyJ||d|dksJ||d |dksJ||d |dksJ||d |dksJ||d |dksJ||d |dksJ||d |dksJ||d |dksJ||d |dksJt}||d|dksJ||d|dksJ||d|dksJ||d|dksJ||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 |dksJt}||d|dksJ||d|dksJ||d|dksJ||d|dksJ||d|dksJ||d|dksJ||d|dksJ||d|dksJ||d |dksJ||d |dksJ||d |dksJ||d |dksJ||d |dksJ||d |dksJ||d |dks(J||d |dks3Jt}||d|dks@J||d|dksKJ||d|dksVJ||d|dksaJ||d|dkslJ||d|dkswJ||d|dksJ||d|dksJ||d |dksJ||d |dksJ||d |dksJ||d |dksJ||d |dksJ||d |dksJ||d |dksJ||d |dksJqd 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)make_mpf mpf_pow_int_mpf_rr)xnrr r r powrCsz'test_pow_epsilon_rounding..powr))r)ii rrirN)rrr rr round_ceiling round_floorr )r*inprecoutprecpos10001pos09999neg10001neg09999r)r r r test_pow_epsilon_rounding#s  r8)mpmath mpmath.libmprrr#r8r r r r s