ó ¡¼™\c@ sXdZddlmZmZddlmZdefd„ƒYZd„Zd„Z dS( sñ 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. iÿÿÿÿ(tprint_functiontdivision(t nsimplifytintervalcB søeZdZd„Zed„ƒZed„ƒZd„Zd„Zd„Z d„Z d„Z d „Z d „Z d „Zd „ZeZd „Zd„Zd„Zd„ZeZd„Zd„Zd„ZeZeZd„Zd„Zd„ZRS(s  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 The comparison of two intervals returns a tuple of two 3-valued logic values. The first value determines the comparison as follows: - True: If the comparison is True throughout the intervals. - False: If the comparison is False throughout the intervals. - None: If the comparison is True for some part of the intervals. The second value is determined as follows: - True: If both the intervals in comparison are valid. - False: If at least one of the intervals is False, else - None cO s|jdtƒ|_t|ƒdkr‡t|dtƒr^|dj|dj|_|_q t|dƒ|_t|dƒ|_n„t|ƒdkrÿ|d|dkrÖt|dƒ|_t|dƒ|_q t|dƒ|_t|dƒ|_n t dƒ‚dS(Ntis_validiiis9interval takes a maximum of two float values as arguments( tpoptTrueRtlent isinstanceRtstarttendtfloatt ValueError(tselftargstkwargs((sNlib/python2.7/site-packages/sympy/plotting/intervalmath/interval_arithmetic.pyt__init__Bs$cC s|j|jdS(Ng@(R R (R ((sNlib/python2.7/site-packages/sympy/plotting/intervalmath/interval_arithmetic.pytmidVscC s|j|jS(N(R R (R ((sNlib/python2.7/site-packages/sympy/plotting/intervalmath/interval_arithmetic.pytwidthZscC sd|j|jfS(Nsinterval(%f, %f)(R R (R ((sNlib/python2.7/site-packages/sympy/plotting/intervalmath/interval_arithmetic.pyt__repr__^scC sd|j|jfS(Ns[%f, %f](R R (R ((sNlib/python2.7/site-packages/sympy/plotting/intervalmath/interval_arithmetic.pyt__str__ascC s t|ttfƒr]|j|kr1t|jfS|j|krMt|jfSd|jfSn©t|t ƒr|jtksŠ|jtkr“t}n-|jdks±|jdkrºd}nt}|j|jkrÜt|fS|j|jkrøt|fSd|fSt SdS(N( RtintR R RRR tFalsetNoneRtNotImplemented(R tothertvalid((sNlib/python2.7/site-packages/sympy/plotting/intervalmath/interval_arithmetic.pyt__lt__ds$       cC st|ttfƒr]|j|kr1t|jfS|j|krMt|jfSd|jfSn t|t ƒry|j |ƒSt SdS(N( RRR R RRR RRRRR(R R((sNlib/python2.7/site-packages/sympy/plotting/intervalmath/interval_arithmetic.pyt__gt__|s   cC s2t|ttfƒri|j|kr@|j|kr@t|jfS||krYd|jfSt|jfSnt|t ƒr*|jtks–|jtkrŸt}n-|jdks½|jdkrÆd}nt}|j|jkrú|j|jkrút|fS|j |ƒddk rt|fSd|fSnt SdS(Ni( RRR R R RRRRRRR(R RR((sNlib/python2.7/site-packages/sympy/plotting/intervalmath/interval_arithmetic.pyt__eq__‰s$     $   cC s/t|ttfƒri|j|kr@|j|kr@t|jfS||krYd|jfSt|jfSnt|t ƒr'|jtks–|jtkrŸt}n-|jdks½|jdkrÆd}nt}|j|jkrú|j|jkrút|fS|j |ƒddk rt|fSd|fSt SdS(Ni( RRR R R RRRRRRR(R RR((sNlib/python2.7/site-packages/sympy/plotting/intervalmath/interval_arithmetic.pyt__ne__¢s$     $   cC s t|ttfƒr]|j|kr1t|jfS|j|krMt|jfSd|jfSnt|t ƒr|jtksŠ|jtkr“t}n-|jdks±|jdkrºd}nt}|j|jkrÜt|fS|j|jkrøt|fSd|fSt SdS(N( RRR R RRR RRRR(R RR((sNlib/python2.7/site-packages/sympy/plotting/intervalmath/interval_arithmetic.pyt__le__ºs$       cC s}t|ttfƒr]|j|kr1t|jfS|j|krMt|jfSd|jfSnt|t ƒry|j |ƒSdS(N( RRR R RRR RRRR(R R((sNlib/python2.7/site-packages/sympy/plotting/intervalmath/interval_arithmetic.pyt__ge__Òs  cC s t|ttfƒrl|jr9t|j||j|ƒS|j|}|j|}t||d|jƒSn™t|tƒr|j|j}|j|j}|jrº|jrºt||ƒS|jtksØ|jtkrët||dtƒSt||ddƒSnt SdS(NR( RRR RRR R RRR(R RR R ((sNlib/python2.7/site-packages/sympy/plotting/intervalmath/interval_arithmetic.pyt__add__Ýs    cC söt|ttfƒrE|j|}|j|}t||d|jƒSt|tƒrî|j|j}|j|j}|jr§|jr§t|j|j|j|jƒS|jtksÅ|jtkrØt||dtƒSt||ddƒSnt SdS(NR( RRR R R RRRRR(R RR R ((sNlib/python2.7/site-packages/sympy/plotting/intervalmath/interval_arithmetic.pyt__sub__ôs  !cC sit|ttfƒrE||j}||j}t||d|jƒSt|tƒra|j|ƒStSdS(NR( RRR R R RRR"R(R RR R ((sNlib/python2.7/site-packages/sympy/plotting/intervalmath/interval_arithmetic.pyt__rsub__s   cC s@|jrt|j |j ƒSt|j |j d|jƒSdS(NR(RRR R (R ((sNlib/python2.7/site-packages/sympy/plotting/intervalmath/interval_arithmetic.pyt__neg__s cC sVt|tƒr|jtks-|jtkrMttdƒ tdƒdtƒS|jdksk|jdkr‹ttdƒ tdƒddƒSg}|j|j|jƒ|j|j|jƒ|j|j|jƒ|j|j|jƒt |ƒ}t |ƒ}t||ƒSn=t|t tfƒrNt|j||j|d|jƒSt SdS(NtinfR( RRRRR RtappendR R tmintmaxRR(R RtintersR R ((sNlib/python2.7/site-packages/sympy/plotting/intervalmath/interval_arithmetic.pyt__mul__s     $cC sWt|ttfƒr1|j|ko0|j|kS|j|jkoR|j|jkSdS(N(RRR R R (R R((sNlib/python2.7/site-packages/sympy/plotting/intervalmath/interval_arithmetic.pyt __contains__,scC sRt|ttfƒr.t|ƒ}|j|ƒSt|tƒrJ|j|ƒStSdS(N(RRR Rt__div__R(R R((sNlib/python2.7/site-packages/sympy/plotting/intervalmath/interval_arithmetic.pyt__rdiv__2s    cC sú|js,ttdƒ tdƒd|jƒSt|ttfƒr‹|dkrmttdƒ tdƒdtƒSt|j||j|ƒSnkt|tƒrò|jtks¸|jtkrØttdƒ tdƒdtƒS|jdksö|jdkrttdƒ tdƒddƒSd|krBttdƒ tdƒddƒS|}|jdkrh| }| }ng}|j |j|jƒ|j |j|jƒ|j |j|jƒ|j |j|jƒt |ƒ}t |ƒ}t||ƒSnt SdS(NR%Ri( RRR RRRR R RR&R(R'R(R RtthisR)R R ((sNlib/python2.7/site-packages/sympy/plotting/intervalmath/interval_arithmetic.pyR,;s4 #         cC sµddlm}m}|js#|St|tƒrF||||ƒƒSt|ttfƒr­|dkr~d|jt |ƒƒSt|ƒ|krt ||ƒSt ||ƒSnt SdS(Ni(texptlogi( t lib_intervalR/R0RRRR Rt__pow__tabst_pow_intt _pow_floatR(R RR/R0((sNlib/python2.7/site-packages/sympy/plotting/intervalmath/interval_arithmetic.pyR2gs   cC st|ttfƒrö|js"|S|dkrØ|jdkr]ttdƒ tdƒdtƒSt|jƒ}|j ƒ\}}|ddkr®ttdƒ tdƒdtƒSt |ƒ|j }|}t||ƒSqt||j||j ƒSn t|tƒr|j |ƒSt SdS(NiR%Ri(RR RRRRRRR tas_numer_denomR3R R2R(R Rtpower_rationaltnumtdenomR R ((sNlib/python2.7/site-packages/sympy/plotting/intervalmath/interval_arithmetic.pyt__rpow__ys$    cC st|j|j|jfƒS(N(thashRR R (R ((sNlib/python2.7/site-packages/sympy/plotting/intervalmath/interval_arithmetic.pyt__hash__’s(t__name__t __module__t__doc__RtpropertyRRRRRRRRRR R!t__radd__R"R#R$R*t__rmul__R+R-R,t __truediv__t __rtruediv__R2R:R<(((sNlib/python2.7/site-packages/sympy/plotting/intervalmath/interval_arithmetic.pyR's6             )  cC s”t|ƒ}|jƒ\}}|ddkrŽt|jƒ|}t|jƒ|}|dkr{tdt||ƒƒ}nt||ƒ}|S|ddkr|jdkrÍttdƒ tdƒdtƒS|jdkrötd|j|ddƒSt|j||j|ƒSn||jdkr:t|jƒ| }n |j|}|jdkrmt|jƒ| }n |j|}t||d|j ƒSdS(s1Evaluates an interval raised to a floating point.iiR%RN( RR6R3R R RR(R RRR(tintertpowerR7R8R9R R tret((sNlib/python2.7/site-packages/sympy/plotting/intervalmath/interval_arithmetic.pyR5–s,     cC sžt|ƒ}|d@r1t|j||j|ƒS|jdkr|jdkrd}t|j||j|ƒ}t||ƒSt|j||j|ƒSdS(s0Evaluates an interval raised to an integer poweriiN(RRR R R((RERFR R ((sNlib/python2.7/site-packages/sympy/plotting/intervalmath/interval_arithmetic.pyR4·s   N( R?t __future__RRtsympy.simplify.simplifyRtobjectRR5R4(((sNlib/python2.7/site-packages/sympy/plotting/intervalmath/interval_arithmetic.pyt s ÿp !