Rie<TdZddlmZddlmZddlmZGddZdZd Z d S) a Interval Arithmetic for plotting. This module does not implement interval arithmetic accurately and hence cannot be used for purposes other than plotting. If you want to use interval arithmetic, use mpmath's interval arithmetic. The module implements interval arithmetic using numpy and python floating points. The rounding up and down is not handled and hence this is not an accurate implementation of interval arithmetic. The module uses numpy for speed which cannot be achieved with mpmath. ) fuzzy_and) nsimplify)intervalMembershipceZdZdZdddZedZedZdZdZ d Z d Z d Z d Z d ZdZdZeZdZdZdZdZeZdZdZdZdZdZdZdS)intervala: Represents an interval containing floating points as start and end of the interval The is_valid variable tracks whether the interval obtained as the result of the function is in the domain and is continuous. - True: Represents the interval result of a function is continuous and in the domain of the function. - False: The interval argument of the function was not in the domain of the function, hence the is_valid of the result interval is False - None: The function was not continuous over the interval or the function's argument interval is partly in the domain of the function A comparison between an interval and a real number, or a comparison between two intervals may return ``intervalMembership`` of two 3-valued logic values. Tis_validcf||_t|dkrxt|dtr'|dj|djc|_|_dSt |d|_t |d|_dSt|dkr~|d|dkr6t |d|_t |d|_dSt |d|_t |d|_dStd)Nrrz9interval takes a maximum of two float values as arguments)r len isinstancerstartendfloat ValueError)selfr argskwargss Olib/python3.11/site-packages/sympy/plotting/intervalmath/interval_arithmetic.py__init__zinterval.__init__;s  t99>>$q'8,, *'+Aw}d1gk$ DHHH"47^^  a>> YY!^^Awa  "47^^  a>>"47^^  a>>+,, ,c&|j|jzdz S)Ng@rrrs rmidz interval.midOs TX%,,rc |j|jz SN)rrrs rwidthzinterval.widthSsx$*$$rc$d|j|jfzS)Nzinterval(%f, %f)rrs r__repr__zinterval.__repr__Ws!TZ$:::rc$d|j|jfzS)Nz[%f, %f]rrs r__str__zinterval.__str__ZsTZ222rct|ttfrU|j|krt d|jS|j|krt d|jSt d|jSt|trkt|j|jg}|j|jkrt d|S|j|jkrt d|St d|StSNTF rintrrrr rrrNotImplementedrothervalids r__lt__zinterval.__lt__]s ec5\ * * "x%)$ >>>e##)%???)$ >>> x ( ( "t}en=>>Ex%,&&)$666zEI%%)%777%dE22 2! !rcFt|ttfrU|j|krt d|jS|j|krt d|jSt d|jSt|tr||StSr%) rr'rrrr rrr,r(rr*s r__gt__zinterval.__gt__ps ec5\ * * "zE!!)$ >>>E!!)%???)$ >>> x ( ( "<<%% %! !rc0t|ttfrY|j|kr |j|krt d|jS||vrt d|jSt d|jSt|trt|j|jg}|j|jkr |j|jkrt d|S| |dt d|St d|StS)NTFr rr'rrrrr rrr,r(r)s r__eq__zinterval.__eq__}s ec5\ * * @zU""tx5'8'8)$ >>>}})$ >>>)%??? eX & & "t}en=>>EzU[((TX-B-B)$666U##A&2)%777)$666! !rc0t|ttfrY|j|kr |j|krt d|jS||vrt d|jSt d|jSt|trt|j|jg}|j|jkr |j|jkrt d|S| |dt d|St d|StS)NFTrr1r)s r__ne__zinterval.__ne__s ec5\ * * ?zU""tx5'8'8)%???}})$ >>>)$ >>> eX & & "t}en=>>EzU[((TX-B-B)%777;;u%%a(0)$666%dE22 2! !rct|ttfrU|j|krt d|jS|j|krt d|jSt d|jSt|trkt|j|jg}|j|jkrt d|S|j|jkrt d|St d|StSr%r&r)s r__le__zinterval.__le__s ec5\ * * ?x5  )$ >>>zE!!)%???)$ >>> eX & & "t}en=>>Ex5;&&)$666zEI%%)%777%dE22 2! !rc<t|ttfrU|j|krt d|jS|j|krt d|jSt d|jSt|tr||SdSr%) rr'rrrr rrr6r.s r__ge__zinterval.__ge__s ec5\ * * &zU"")$ >>>E!!)%???)$ >>> x ( ( &<<%% % & &rct|ttfrR|jr t |j|z|j|zS|j|z}|j|z}t |||jSt|trK|j|jz}|j|jz}t|j|jg}t |||StSNr ) rr'rr rrrrr(rr*rrr+s r__add__zinterval.__add__s ec5\ * * "} D U 2DHu4DEEE U*h&sT]CCCC x ( ( "J,E(UY&Ct}en=>>EE3777 7! !rc^t|ttfr+|j|z }|j|z }t |||jSt|t rK|j|jz }|j|jz }t|j|jg}t |||StSr:) rr'rrrrr rr(r;s r__sub__zinterval.__sub__s ec5\ * * "J&E(U"CE3??? ? x ( ( "J*E(U[(Ct}en=>>EE3777 7! !rct|ttfr+||jz }||jz }t |||jSt|t r||StSr:) rr'rrrrr r>r()rr*rrs r__rsub__zinterval.__rsub__sq ec5\ * * "DH$E$*$CE3??? ? x ( ( "==&& &! !rc|jrt|j |j St|j |j |jSr:)r rrrrs r__neg__zinterval.__neg__sD = LTXI {33 3TXI {T]KKK Krc&t|tr2|jdus |jdur-ttd tddS|j|j-ttd tddSg}||j|jz||j|jz||j|jz||j|jzt|}t|}t||St|ttfr't|j|z|j|z|jStS)NFinfr ) rrr rappendrrminmaxr'r()rr*intersrrs r__mul__zinterval.__mul__sY eX & & "}%%5)@)@u uU||eLLLL&%.*@u uU||dKKKK dj5;6777 dh4555 dj594555 dh2333F &kks+++ U| , , "DJu,dhunt}UUU U! !rct|ttfr|j|ko |j|kS|j|jko|j|jkSr)rr'rrrr.s r __contains__zinterval.__contains__ sO ec5\ * * G:&<48u+< <:,Fdh1F Frct|ttfr$t|}||St|tr||St Sr)rr'rr __truediv__r(r.s r __rtruediv__zinterval.__rtruediv__sc ec5\ * * "UOOE$$T** * x ( ( "$$T** *! !rcx|js2ttd td|jSt|ttfrS|dkr-ttd tddSt|j|z |j|z St|trv|jdus |jdur-ttd tddS|j|j-ttd tddSd|vr-ttd tddS|}|jdkr| }| }g}||j|jz ||j|jz ||j|jz ||j|jz t|}t|}t||StS)NrDr rF) r rrrr'rrrErGrFr()rr*thisrHrrs rrMzinterval.__truediv__s} QU5\\M5<<$-PPP P ec5\ * *" "zzu uU||eLLLL U 2DHu4DEEE x ( ( "~&&$-5*@*@u uU||eLLLL'4=+@u uU||dKKKK::#U5\\M5<<$OOOO9q== 5D"FE dj5;6777 dh4555 dj594555 dh2333F &kks+++! !rc~ddlm}m}|js|St |t r||||zSt |t tfr^|dkr%d|t|z St||krt||St||StS)Nr)explogr) lib_intervalrRrSr rrrr'__pow__abs_pow_int _pow_floatr()rr*rRrSs rrUzinterval.__pow__As********} K eX & & "3uss4yy()) ) s| , , "qyy4<<E 3333u::&&#D%000%dE222! !rct|ttfr|js|S|dkr|jdkr-t td tddSt |j}|\}}|dzdkr-t td tddSt||jz }|}t ||St ||jz||j zSt|t r| |StS)NrrDFr r ) rrr'r rrrras_numer_denomrVrrUr()rr*power_rationalnumdenomrrs r__rpow__zinterval.__rpow__Ss/ eeS\ * * "= D :>>#U5\\M5<<%PPPP%.tz%:%:N!/!>!>!@!@JCqyA~~'u uU||16 8 8 88"%UTZ!7 7#'s333tz 15$(?CCC x ( ( "==&& &! !rcDt|j|j|jfSr)hashr rrrs r__hash__zinterval.__hash__lsT]DJ9:::rN)__name__ __module__ __qualname____doc__rpropertyrrr!r#r,r/r2r4r6r8r<__radd__r>r@rBrI__rmul__rKrNrMrUr^rarrrr)s"(,,,,,,(--X-%%X%;;;333"""& " " """"("""&"""& & & &""""H " " """"LLL """(HGGG """'"'"'"R"""$"""2;;;;;rrct|}|\}}|dzdkret|j|z}t|j|z}|dkrt dt ||}nt ||}|S|dzdkr}|jdkr-t td tddS|jdkrt d|j|zdSt |j|z|j|zS|jdkrt|j|z }n |j|z}|jdkrt|j|z }n |j|z}t |||jS)z1Evaluates an interval raised to a floating point.r rrDFr N) rrZrVrrrrGrr )interpowerr[r\r]rrrets rrXrXpsu%%N..00JC Qw!||EK  %'%)nne# 1991c%oo..CC5#&&C a 9q==U5\\M5<<%HHH H [1__Auy%/$??? ?EK. 50@AA A ;??%%u,,EEK&E 9q==uy>>5((CC)U"CsU^<<<rrs  D'&&&&&------333333D;D;D;D;D;D;D;D;N ===B B B B B Br