ّRMc@sdZddkZddkZddklZddklZeZdZ dZ dZ xCe e jo5e e e jo"e d9Z e d9Z e d7Z qYWe dZe ddZe Zd ZdedZdedZdZx7de>ejo%ed7Zed jp td qWdMZeeZy8dd klZlZlZlZlZlZlZWndej oXdZdZdZdZ dZ!dZ"dZdZdZdZnXyddkl#Z#Wnej odZ#nXdZ$dZ%dZ&dZ'de'_(dZ)d e)_([ [ [ d!Z*d"Z+d#Z,ei-d$jo d%Z.n d&Z.d'Z/d(Z0d e0_(d)Z1d*Z2de2_(d+Z3d e3_(d,Z4d-e5fd.YZ6d/e6fd0YZ7d1e6fd2YZ8hZ9e:d3Z;d4ei<fd5YZ=d6ei<fd7YZ>e;d8eeZ?e;d9e:eZ@e;d:ed ZAe;d;e:d ZBeAedZCeAe?j o eAZDne ZDd<ZEd=ZFd>ZGd?ZHdZIdZJd@ZKdZLdZMdAZNdBZOddCZPdDZQdEeRfdFYZSdGeRfdHYZTdIei<fdJYZUdKei<fdLYZVdS(Ns This file defines restricted arithmetic: classes and operations to express integer arithmetic, such that before and after translation semantics are consistent r_uint an unsigned integer which has no overflow checking. It is always positive and always truncated to the internal machine word size. intmask mask a possibly long value when running on CPython back to a signed int value ovfcheck check on CPython whether the result of a signed integer operation did overflow ovfcheck_lshift << with oveflow checking catering to 2.3/2.4 differences about << ovfcheck_float_to_int convert to an integer or raise OverflowError r_longlong like r_int but double word size r_ulonglong like r_uint but double word size widen(x) if x is of a type smaller than lltype.Signed or lltype.Unsigned, widen it to lltype.Signed. Useful because the translator doesn't support arithmetic on the smaller types. These are meant to be erased by translation, r_uint in the process should mark unsigned values, ovfcheck should mark where overflow checking is required. iN(t extregistry(t objectmodeliilii@icsLONG_BIT_SHIFT value not found?gZbti(tisinftisnantcopysigntacoshtasinhtatanhtlog1pcCs|tjp |t jS(t NOT_RPYTHON(tINFINITY(tx((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyRFscCs ||jS(R ((tv((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyRJscCs\|djo | }n|djp&|djo!ti|ddjo|Sn| SdS(s(NOT_RPYTHON. Return x with the sign of yggN(tmathtatan2(R ty((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyRNs   3g0>gAg9B.?cCst|otSn|djotdn|tjo*t|o|Sqiti|tSn|djodSn|djo8||}tid|d|ti|dSn|d}t |tid|||S(R g?smath domain errorgg@( RtNANt ValueErrort _2_to_p28RR tlogt_ln2tsqrtR(R tt((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyR[s       . cCst|}t|p t|o|Sn|tjo|Sn|tjoti|t}nq|djo4tid|dti||d|}n0||}t ||dtid|}t ||S(R g@g?( tabsRRt _2_to_m28RR RRRRR(R tabsxtwR((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyRns    4 %cCst|o|Snt|}|djotdn|tjo|Sn|djo.||}dt|||d|}ndt||d|}t||S(R g?smath domain errorg?(RRRRRR(R RR((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyR~s      $cCsddkl}t||idjo|Sn[d|jo djno+d|}ti||d||Sntid|SdS(R i(trfloatg@gg?N(t pypy.rlibRRt DBL_EPSILONR R(R RR((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyRs !(texpm1cCset|djoAti|}|djo|Sn|d|ti|Snti|dS(R gffffff?g?(RR texpR(R tu((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyRs  cCsUt|}|ti|djoti|}nti|}t||S(Ng?(RR tfloortceilR(R Rtr((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt round_aways  cCst|tot|Snt|tio|Snt|t ptt|}|tM}|tjo|dt8}nt|S(Ni( t isinstancetintRtSymbolictfloattAssertionErrortlongt LONG_MASKt LONG_TEST(tn((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pytintmasks   cCsjt|tot|}nt|tpt|tM}|tjo|dt8}nt|S(Ni(R%R&R*R)t LONGLONG_MASKt LONGLONG_TESTt r_longlong(R-((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt longlongmasks  cCs<ddkl}t|i|ot|Sn|SdS(Ni(tlltype(tpypy.rpython.lltypesystemR3t_should_widen_typettypeOfR.(R-R3((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pytwidenssspecialize:argtype(0)cCsddkl}l}||ijotSn||ijotSn|ii|}t |t pt |i t jp|i t jo|iS(Ni(R3trffi(R4R3R8tBooltTruetSignedtFalsetplatformtnumbertype_to_rclasst issubclasstbase_intR)tBITStLONG_BITtSIGNED(ttpR3R8tr_class((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyR5ssspecialize:memocCsut|t p tdt|t p tdt|t p tdt|tjo tdn|S(R s unexpected ovf check on unsigneds$ovfcheck not supported on r_longlongs%ovfcheck not supported on r_ulonglongs&signed integer expression did overflow(R%tr_uintR)R1t r_ulonglongttypeR*t OverflowError(R#((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pytovfchecks  cCs<t|t p tdt|to tdn|S(Ns unexpected ovf check on unsigneds&signed integer expression did overflow(R%RFR)R*RI(R#((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt_local_ovfchecks cCsttt||>S(R (RKR&R*(tatb((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pytovfcheck_lshiftsicCsMt|o tnd|jo djnot|SntdS(Ng gA(RRIR&(R ((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pytovfcheck_float_to_ints   cCsMt|o tnd|jo djnot|SntdS(NggC(RRIR&(R ((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyROs   cCs||jo|tjotSn|Sn|tttfjo|tjotSn|Sn|tttfjo|Sntd|io|it|i|iS(N(tboolR&R*t build_inttNoneRCtmaxRA(t self_typet other_type((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pytcompute_restypes   cCs)|tttfjotSn|iSdS(N(RPR&R*R:RC(R((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt signedtypescCsv|tjotSn|itijotd|itiSn1|itijpttd|itiSdS(N(R&RAtr_intRQRRRCR1R)(R((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pytnormalizedinttypes  cCs#ddkl}t|i|S(Ni(R3(R4R3tmost_neg_value_ofR6(R R3((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pytmost_neg_value_of_same_type#scCsddkl}l}||ijoti dSn|ii|}t|t pt |i o||i d? dSn |dSdS(Ni(R3R8ii( R4R3R8R;tsystmaxintR=R>R?R@R)RCtMASK(RDR3R8RE((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyRZ(s cCsP|o||d@djptd}x |o|d7}|dL}q,W|S(s Calculates the highest set bit in n. This function assumes that n is a power of 2 (and thus only has a single set bit). iii(R)(R-ti((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt highest_bit4s# R@cBseZdZdZdZdZeZdZdZdZ e Z dZ e Z dZ e Zd Zd Zd Zd Zd ZdZdZdZeZdZeZdZeZdZdZdZdZddZ!ddZ"RS(s& fake unsigned integer implementation cCsut|}t|}y|i||f|SWntj onXt||}||i||f<||S(s if one argument is int or long, the other type wins. otherwise, produce the largest class to hold the result. (RHttypemaptKeyErrorRV(tselftothertvalueRTRUtrestype((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt_widenEs  cCs:|tjotdntt|i||SdS(Nsabstract base!(R@t TypeErrortsupert__new__(tklasstval((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyRjTs cCs,t|}t|}|i|||S(N(R*Rg(RcRdR R((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt__add__Zs  cCs,t|}t|}|i|||S(N(R*Rg(RcRdR R((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt__sub__`s  cCs,t|}t|}|i|||S(N(R*Rg(RcRdRR ((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt__rsub__es  cCsNt|}t|ttfp ||Snt|}|i|||S(N(R*R%R&Rg(RcRdR R((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt__mul__js    cCs,t|}t|}|i|||S(N(R*Rg(RcRdR R((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt__div__rs  cCs,t|}t|}|i|||S(N(R*Rg(RcRdRR ((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt__rdiv__ys  cCs,t|}t|}|i|||S(N(R*Rg(RcRdR R((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt__mod__s  cCs,t|}t|}|i|||S(N(R*Rg(RcRdRR ((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt__rmod__s  cCsKt|}t|}t||}|i|d|i|dfS(Nii(R*tdivmodt __class__(RcRdR Rtres((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt __divmod__s  cCs)t|}t|}|i||>S(N(R*Rv(RcR-R R((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt __lshift__s  cCs,t|}t|}|i|||>S(N(R*Rg(RcR-RR ((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt __rlshift__s  cCs,t|}t|}|i|||?S(N(R*Rg(RcR-R R((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt __rshift__s  cCs,t|}t|}|i|||?S(N(R*Rg(RcR-RR ((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt __rrshift__s  cCs,t|}t|}|i|||BS(N(R*Rg(RcRdR R((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt__or__s  cCs,t|}t|}|i|||@S(N(R*Rg(RcRdR R((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt__and__s  cCs,t|}t|}|i|||AS(N(R*Rg(RcRdR R((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt__xor__s  cCst|}|i| S(N(R*Rv(RcR ((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt__neg__s cCst|}|it|S(N(R*RvR(RcR ((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt__abs__s cCs |i|S(N(Rv(Rc((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt__pos__scCst|}|i|S(N(R*Rv(RcR ((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt __invert__s cCs:t|}t|}t|||}|i||S(N(R*tpowRg(RcRdtmR RRw((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt__pow__s  cCs:t|}t|}t|||}|i||S(N(R*RRg(RcRdRRR Rw((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyt__rpow__s  N(#t__name__t __module__t__doc__RgRjRmt__radd__RnRoRpt__rmul__Rqt __floordiv__Rrt __rfloordiv__RsRtRxRyRzR{R|R}t__ror__R~t__rand__Rt__rxor__RRRRRRRR(((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyR@As@                       t signed_intcBs eZeZddZhZRS(icCst|tjot|}n||id?jp||id? djotd||ifn|djo||i@}ntt|i||S(Nis(%s does not fit in signed %d-bit integeri( RHR(R*R^RIRARiRRj(RkRl((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyRjs- (RRR:RCRjRa(((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyRs t unsigned_intcBs eZeZddZhZRS(icCsFt|ttfot|}ntt|i|||i@S(N(R%R(R*RiRRjR^(RkRl((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyRjs(RRR<RCRjRa(((s)/64/home/arigo/u/pypy/rlib/rarithmetic.pyRs c s"t|}|p,yt||fSWq?tj oq?Xn|o t}nt}d|d}|djo$tdddg||fnt||fh|d<|d<|d<|pt||f#s      8                      (   5