3 \Tƒã$@sJdZddddddddd d ddd ddddddddddddddddddd d!d"d#d$g$ZeZd%Zd&Zd'Zd(d)lZd(d)lZ d(d)l Z yd(d*lmZ e dd+ƒZWnek r®d,d-„ZYnXdZdZdZdZdZdZdZdZd.Ze jd•krìd2Zd2Zd–Znd3Zd3Zd—Zeed1ZGd4d„deƒZGd5d„deƒZ Gd6d „d eƒZ!Gd7d„de!ƒZ"Gd8d „d ee#ƒZ$Gd9d„de!ƒZ%Gd:d„de!e#ƒZ&Gd;d„deƒZ'Gdd „d eƒZ*Gd?d„de'e)ƒZ+Gd@d„de'e)e*ƒZ,GdAd„dee-ƒZ.e e$e'e+e)e,e!e*e.g Z/e"e!e%e!e&e!e(e!iZ0eeeeeeeefZ1yd(d)l2Z2Wn.ek rhGdBdC„dCe3ƒZ4e4ƒZ2[4YnXy e2j5Wn>e6k r²e7e2j8ƒdDƒrže2j8ƒ`9dEd„Z:dFd„Z;Yn6Xe2j5ƒZ5e7e5dDƒrÌe5`9e5fdGd„Z;e5fdHd„Z:[2[5d˜dId„Ze j?j@e=ƒGdNdO„dOe3ƒZAGdPd„de3ƒZBGdQdR„dRe3ƒZCdšdSdT„ZDeEjFZGdUdV„ZHdWdX„ZIdYdZ„ZJd[d\„ZKd›d^d_„ZLd`da„ZMdbdc„ZNGddde„dee3ƒZOeOƒjPZQdœdfdg„ZRdhdi„ZSdjdk„ZTdldmdndodpdqdrdsdtduœ fdvdw„ZUddxdy„ZVdždzd{„ZWeBd|ee$e+e!ggd}dŸd1d(d~ZXeBdee$e+e!e e,ggd€ZYeBdeggd€ZZd(d)l[Z[e[j\de[j]e[j^Bƒj_Z`e[j\d‚ƒj_Zae[j\dƒƒj_Zbe[j\d„e[j]e[jcBƒZd[[yd(d)leZfWnek r®YnXd d…d†„Zgd‡dˆ„Zhd‰dŠ„Zid¡d‹dŒ„ZjddŽ„Zkdd„Zle=d‘ƒZme=d’ƒZne=d“ƒZoe=d(ƒZpe=d1ƒZqe=d¢ƒZremenfZse jtjuZve jtjwZxe jtjyZze{dsevd/evƒZ|[ d)S)£a¿ This is an implementation of decimal floating point arithmetic based on the General Decimal Arithmetic Specification: http://speleotrove.com/decimal/decarith.html and IEEE standard 854-1987: http://en.wikipedia.org/wiki/IEEE_854-1987 Decimal floating point has finite precision with arbitrarily large bounds. The purpose of this module is to support arithmetic using familiar "schoolhouse" rules and to avoid some of the tricky representation issues associated with binary floating point. The package is especially useful for financial applications or for contexts where users have expectations that are at odds with binary floating point (for instance, in binary floating point, 1.00 % 0.1 gives 0.09999999999999995 instead of 0.0; Decimal('1.00') % Decimal('0.1') returns the expected Decimal('0.00')). Here are some examples of using the decimal module: >>> from decimal import * >>> setcontext(ExtendedContext) >>> Decimal(0) Decimal('0') >>> Decimal('1') Decimal('1') >>> Decimal('-.0123') Decimal('-0.0123') >>> Decimal(123456) Decimal('123456') >>> Decimal('123.45e12345678') Decimal('1.2345E+12345680') >>> Decimal('1.33') + Decimal('1.27') Decimal('2.60') >>> Decimal('12.34') + Decimal('3.87') - Decimal('18.41') Decimal('-2.20') >>> dig = Decimal(1) >>> print(dig / Decimal(3)) 0.333333333 >>> getcontext().prec = 18 >>> print(dig / Decimal(3)) 0.333333333333333333 >>> print(dig.sqrt()) 1 >>> print(Decimal(3).sqrt()) 1.73205080756887729 >>> print(Decimal(3) ** 123) 4.85192780976896427E+58 >>> inf = Decimal(1) / Decimal(0) >>> print(inf) Infinity >>> neginf = Decimal(-1) / Decimal(0) >>> print(neginf) -Infinity >>> print(neginf + inf) NaN >>> print(neginf * inf) -Infinity >>> print(dig / 0) Infinity >>> getcontext().traps[DivisionByZero] = 1 >>> print(dig / 0) Traceback (most recent call last): ... ... ... decimal.DivisionByZero: x / 0 >>> c = Context() >>> c.traps[InvalidOperation] = 0 >>> print(c.flags[InvalidOperation]) 0 >>> c.divide(Decimal(0), Decimal(0)) Decimal('NaN') >>> c.traps[InvalidOperation] = 1 >>> print(c.flags[InvalidOperation]) 1 >>> c.flags[InvalidOperation] = 0 >>> print(c.flags[InvalidOperation]) 0 >>> print(c.divide(Decimal(0), Decimal(0))) Traceback (most recent call last): ... ... ... decimal.InvalidOperation: 0 / 0 >>> print(c.flags[InvalidOperation]) 1 >>> c.flags[InvalidOperation] = 0 >>> c.traps[InvalidOperation] = 0 >>> print(c.divide(Decimal(0), Decimal(0))) NaN >>> print(c.flags[InvalidOperation]) 1 >>> ÚDecimalÚContextÚDecimalTupleÚDefaultContextÚBasicContextÚExtendedContextÚDecimalExceptionÚClampedÚInvalidOperationÚDivisionByZeroÚInexactÚRoundedÚ SubnormalÚOverflowÚ UnderflowÚFloatOperationÚDivisionImpossibleÚInvalidContextÚConversionSyntaxÚDivisionUndefinedÚ ROUND_DOWNÚ ROUND_HALF_UPÚROUND_HALF_EVENÚ ROUND_CEILINGÚROUND_FLOORÚROUND_UPÚROUND_HALF_DOWNÚ ROUND_05UPÚ setcontextÚ getcontextÚlocalcontextÚMAX_PRECÚMAX_EMAXÚMIN_EMINÚ MIN_ETINYÚHAVE_THREADSÚdecimalz1.70z2.4.2éN)Ú namedtuplezsign digits exponentcGs|S)N©)Úargsr(r(ú"/usr/lib64/python3.6/_pydecimal.pyÚ¡sr+Téé?élÿÇNÎZoi@üTc@seZdZdZdd„ZdS)ra1Base exception class. Used exceptions derive from this. If an exception derives from another exception besides this (such as Underflow (Inexact, Rounded, Subnormal) that indicates that it is only called if the others are present. This isn't actually used for anything, though. handle -- Called when context._raise_error is called and the trap_enabler is not set. First argument is self, second is the context. More arguments can be given, those being after the explanation in _raise_error (For example, context._raise_error(NewError, '(-x)!', self._sign) would call NewError().handle(context, self._sign).) To define a new exception, it should be sufficient to have it derive from DecimalException. cGsdS)Nr()ÚselfÚcontextr)r(r(r*ÚhandleÏszDecimalException.handleN)Ú__name__Ú __module__Ú__qualname__Ú__doc__r1r(r(r(r*r¼sc@seZdZdZdS)ra)Exponent of a 0 changed to fit bounds. This occurs and signals clamped if the exponent of a result has been altered in order to fit the constraints of a specific concrete representation. This may occur when the exponent of a zero result would be outside the bounds of a representation, or when a large normal number would have an encoded exponent that cannot be represented. In this latter case, the exponent is reduced to fit and the corresponding number of zero digits are appended to the coefficient ("fold-down"). N)r2r3r4r5r(r(r(r*rÓs c@seZdZdZdd„ZdS)r a0An invalid operation was performed. Various bad things cause this: Something creates a signaling NaN -INF + INF 0 * (+-)INF (+-)INF / (+-)INF x % 0 (+-)INF % x x._rescale( non-integer ) sqrt(-x) , x > 0 0 ** 0 x ** (non-integer) x ** (+-)INF An operand is invalid The result of the operation after these is a quiet positive NaN, except when the cause is a signaling NaN, in which case the result is also a quiet NaN, but with the original sign, and an optional diagnostic information. cGs,|r(t|dj|djddƒ}|j|ƒStS)Nr&ÚnT)Ú_dec_from_tripleÚ_signÚ_intÚ_fix_nanÚ_NaN)r/r0r)Úansr(r(r*r1ös zInvalidOperation.handleN)r2r3r4r5r1r(r(r(r*r ßsc@seZdZdZdd„ZdS)rzÜTrying to convert badly formed string. This occurs and signals invalid-operation if a string is being converted to a number and it does not conform to the numeric string syntax. The result is [0,qNaN]. cGstS)N)r;)r/r0r)r(r(r*r1szConversionSyntax.handleN)r2r3r4r5r1r(r(r(r*rüsc@seZdZdZdd„ZdS)r a²Division by 0. This occurs and signals division-by-zero if division of a finite number by zero was attempted (during a divide-integer or divide operation, or a power operation with negative right-hand operand), and the dividend was not zero. The result of the operation is [sign,inf], where sign is the exclusive or of the signs of the operands for divide, or is 1 for an odd power of -0, for power. cGst|S)N)Ú_SignedInfinity)r/r0Úsignr)r(r(r*r1szDivisionByZero.handleN)r2r3r4r5r1r(r(r(r*r sc@seZdZdZdd„ZdS)rzóCannot perform the division adequately. This occurs and signals invalid-operation if the integer result of a divide-integer or remainder operation had too many digits (would be longer than precision). The result is [0,qNaN]. cGstS)N)r;)r/r0r)r(r(r*r1szDivisionImpossible.handleN)r2r3r4r5r1r(r(r(r*rsc@seZdZdZdd„ZdS)rzîUndefined result of division. This occurs and signals invalid-operation if division by zero was attempted (during a divide-integer, divide, or remainder operation), and the dividend is also zero. The result is [0,qNaN]. cGstS)N)r;)r/r0r)r(r(r*r1)szDivisionUndefined.handleN)r2r3r4r5r1r(r(r(r*r!sc@seZdZdZdS)raHad to round, losing information. This occurs and signals inexact whenever the result of an operation is not exact (that is, it needed to be rounded and any discarded digits were non-zero), or if an overflow or underflow condition occurs. The result in all cases is unchanged. The inexact signal may be tested (or trapped) to determine if a given operation (or sequence of operations) was inexact. N)r2r3r4r5r(r(r(r*r,s c@seZdZdZdd„ZdS)raìInvalid context. Unknown rounding, for example. This occurs and signals invalid-operation if an invalid context was detected during an operation. This can occur if contexts are not checked on creation and either the precision exceeds the capability of the underlying concrete representation or an unknown or unsupported rounding was specified. These aspects of the context need only be checked when the values are required to be used. The result is [0,qNaN]. cGstS)N)r;)r/r0r)r(r(r*r1CszInvalidContext.handleN)r2r3r4r5r1r(r(r(r*r8s c@seZdZdZdS)raÙNumber got rounded (not necessarily changed during rounding). This occurs and signals rounded whenever the result of an operation is rounded (that is, some zero or non-zero digits were discarded from the coefficient), or if an overflow or underflow condition occurs. The result in all cases is unchanged. The rounded signal may be tested (or trapped) to determine if a given operation (or sequence of operations) caused a loss of precision. N)r2r3r4r5r(r(r(r*rFs c@seZdZdZdS)r a˜Exponent < Emin before rounding. This occurs and signals subnormal whenever the result of a conversion or operation is subnormal (that is, its adjusted exponent is less than Emin, before any rounding). The result in all cases is unchanged. The subnormal signal may be tested (or trapped) to determine if a given or operation (or sequence of operations) yielded a subnormal result. N)r2r3r4r5r(r(r(r*r Rs c@seZdZdZdd„ZdS)raNumerical overflow. This occurs and signals overflow if the adjusted exponent of a result (from a conversion or from an operation that is not an attempt to divide by zero), after rounding, would be greater than the largest value that can be handled by the implementation (the value Emax). The result depends on the rounding mode: For round-half-up and round-half-even (and for round-half-down and round-up, if implemented), the result of the operation is [sign,inf], where sign is the sign of the intermediate result. For round-down, the result is the largest finite number that can be represented in the current precision, with the sign of the intermediate result. For round-ceiling, the result is the same as for round-down if the sign of the intermediate result is 1, or is [0,inf] otherwise. For round-floor, the result is the same as for round-down if the sign of the intermediate result is 0, or is [1,inf] otherwise. In all cases, Inexact and Rounded will also be raised. cGsŽ|jttttfkrt|S|dkrR|jtkr4t|St|d|j|j |jdƒS|dkrŠ|jt krlt|St|d|j|j |jdƒSdS)Nr&Ú9r.)Úroundingrrrrr=rr7ÚprecÚEmaxr)r/r0r>r)r(r(r*r1ss zOverflow.handleN)r2r3r4r5r1r(r(r(r*r]sc@seZdZdZdS)raxNumerical underflow with result rounded to 0. This occurs and signals underflow if a result is inexact and the adjusted exponent of the result would be smaller (more negative) than the smallest value that can be handled by the implementation (the value Emin). That is, the result is both inexact and subnormal. The result after an underflow will be a subnormal number rounded, if necessary, so that its exponent is not less than Etiny. This may result in 0 with the sign of the intermediate result and an exponent of Etiny. In all cases, Inexact, Rounded, and Subnormal will also be raised. N)r2r3r4r5r(r(r(r*rƒs c@seZdZdZdS)raœEnable stricter semantics for mixing floats and Decimals. If the signal is not trapped (default), mixing floats and Decimals is permitted in the Decimal() constructor, context.create_decimal() and all comparison operators. Both conversion and comparisons are exact. Any occurrence of a mixed operation is silently recorded by setting FloatOperation in the context flags. Explicit conversions with Decimal.from_float() or context.create_decimal_from_float() do not set the flag. Otherwise (the signal is trapped), only equality comparisons and explicit conversions are silent. All other mixed operations raise FloatOperation. N)r2r3r4r5r(r(r(r*r’s c@seZdZefdd„ZdS)Ú MockThreadingcCs |jtS)N)ÚmodulesÚ __xname__)r/Úsysr(r(r*Úlocal¼szMockThreading.localN)r2r3r4rFrGr(r(r(r*rC»srCÚ__decimal_context__cCs,|tttfkr|jƒ}|jƒ|tjƒ_dS)z%Set this thread's context to context.N)rrrÚcopyÚclear_flagsÚ threadingÚcurrent_threadrH)r0r(r(r*rËscCs4y tjƒjStk r.tƒ}|tjƒ_|SXdS)zÍReturns this thread's context. If this thread does not yet have a context, returns a new context and sets this thread's context. New contexts are copies of DefaultContext. N)rKrLrHÚAttributeErrorr)r0r(r(r*rÒs cCs,y|jStk r&tƒ}||_|SXdS)zÍReturns this thread's context. If this thread does not yet have a context, returns a new context and sets this thread's context. New contexts are copies of DefaultContext. N)rHrMr)Ú_localr0r(r(r*ræscCs(|tttfkr|jƒ}|jƒ||_dS)z%Set this thread's context to context.N)rrrrIrJrH)r0rNr(r(r*rôscCs|dkrtƒ}t|ƒS)abReturn a context manager for a copy of the supplied context Uses a copy of the current context if no context is specified The returned context manager creates a local decimal context in a with statement: def sin(x): with localcontext() as ctx: ctx.prec += 2 # Rest of sin calculation algorithm # uses a precision 2 greater than normal return +s # Convert result to normal precision def sin(x): with localcontext(ExtendedContext): # Rest of sin calculation algorithm # uses the Extended Context from the # General Decimal Arithmetic Specification return +s # Convert result to normal context >>> setcontext(DefaultContext) >>> print(getcontext().prec) 28 >>> with localcontext(): ... ctx = getcontext() ... ctx.prec += 2 ... print(ctx.prec) ... 30 >>> with localcontext(ExtendedContext): ... print(getcontext().prec) ... 9 >>> print(getcontext().prec) 28 N)rÚ_ContextManager)Zctxr(r(r*rýs$c @sàeZdZdZdõZdödd „Zed d„ƒZdd „Zdd„Z d÷dd„Z dd„Zdd„Zdd„Z dødd„Zdùdd„Zdúdd„Zdûdd„Zdüd d!„Zdýd"d#„Zd$d%„Zd&d'„Zd(d)„Zd*d+„Zdþd-d.„Zdÿd/d0„Zdd1d2„Zdd3d4„Zdd6d7„Zdd8d9„ZeZdd:d;„Zddd?„Z!e!Z"dd@dA„Z#dBdC„Z$ddDdE„Z%d dFdG„Z&d dHdI„Z'ddJdK„Z(ddLdM„Z)d dNdO„Z*ddPdQ„Z+ddRdS„Z,dTdU„Z-dVdW„Z.e.Z/dXdY„Z0e1e0ƒZ0dZd[„Z2e1e2ƒZ2d\d]„Z3d^d_„Z4d`da„Z5dbdc„Z6ddde„Z7dfdg„Z8dhdi„Z9djdk„Z:dldm„Z;dndo„Ze?e7e8e9e:e;edtZ@ddudv„ZAdwdx„ZBdydz„ZCdd{d|„ZDdd}d~„ZEdd€„ZFddd‚„ZGddƒd„„ZHdd…d†„ZIdd‡dˆ„ZJdd‰dŠ„ZKd‹dŒ„ZLddŽ„ZMddd„ZNdd‘d’„ZOeOZPdd“d”„ZQdd•d–„ZRdd—d˜„ZSd™dš„ZTd›dœ„ZUddž„ZVdŸd „ZWdd¡d¢„ZXdd£d¤„ZYdd¥d¦„ZZd§d¨„Z[d©dª„Z\d d«d¬„Z]d!dd®„Z^d¯d°„Z_d±d²„Z`d³d´„Zadµd¶„Zbd"d·d¸„Zcd¹dº„Zdd»d¼„Zed½d¾„Zfd#d¿dÀ„ZgdÁdÂ„ZhdÃdÄ„Zid$dÅdÆ„ZjdÇdÈ„Zkd%dÉdÊ„Zld&dËdÌ„ZmdÍdÎ„ZndÏdÐ„Zod'dÑdÒ„Zpd(dÓdÔ„Zqd)dÕdÖ„Zrd*d×dØ„Zsd+dÙdÚ„Ztd,dÛdÜ„Zud-dÝdÞ„Zvd.dßdà„Zwd/dádâ„Zxd0dãdä„Zydådæ„Zzd1dçdè„Z{d2dédê„Z|d3dëdì„Z}dídî„Z~dïdð„Zdñdò„Z€d4dódô„ZdS(5rz,Floating point class for decimal arithmetic.Ú_expr9r8Ú_is_specialÚ0NcCsŒtj|ƒ}t|tƒr$t|jƒjddƒƒ}|dkrP|dkr@tƒ}|jt d|ƒS|j dƒdkrfd|_nd|_|j d ƒ}|dk rÆ|j d ƒpŠd}t|j dƒpšdƒ}tt||ƒƒ|_ |t|ƒ|_d |_nZ|j dƒ}|dk rtt|pädƒƒjdƒ|_ |j dƒrd|_nd|_nd|_ d|_d|_|St|tƒrf|dkrBd|_nd|_d|_tt|ƒƒ|_ d |_|St|tƒr–|j|_|j|_|j |_ |j|_|St|tƒrÌ|j|_t|jƒ|_ t|jƒ|_d |_|St|ttfƒr&t|ƒdkròtdƒ‚t|dtƒo|ddkstdƒ‚|d|_|ddkrHd|_ |d|_d|_nÚg} x^|dD]R} t| tƒržd| kozdknrž| s’| dkr¦| j| ƒntdƒ‚qVW|ddkrÞdjtt| ƒƒ|_ |d|_d|_nDt|dtƒrdjtt| pdgƒƒ|_ |d|_d |_ntdƒ‚|St|tƒr||dkrBtƒ}|jtdƒtj|ƒ}|j|_|j|_|j |_ |j|_|St d|ƒ‚dS)aêCreate a decimal point instance. >>> Decimal('3.14') # string input Decimal('3.14') >>> Decimal((0, (3, 1, 4), -2)) # tuple (sign, digit_tuple, exponent) Decimal('3.14') >>> Decimal(314) # int Decimal('314') >>> Decimal(Decimal(314)) # another decimal instance Decimal('314') >>> Decimal(' 3.14 \n') # leading and trailing whitespace okay Decimal('3.14') Ú_ÚNzInvalid literal for Decimal: %rr>ú-r.r&ÚintZfracÚexprRFÚdiagÚsignalÚNr6ÚFTéztInvalid tuple size in creation of Decimal from list or tuple. The list or tuple should have exactly three elements.z|Invalid sign. The first value in the tuple should be an integer; either 0 for a positive number or 1 for a negative number.r,é zTThe second value in the tuple must be composed of integers in the range 0 through 9.zUThe third value in the tuple must be an integer, or one of the strings 'F', 'n', 'N'.z;strict semantics for mixing floats and Decimals are enabledzCannot convert %r to Decimal)r&r.)r6rZ)!ÚobjectÚ__new__Ú isinstanceÚstrÚ_parserÚstripÚreplacerÚ_raise_errorrÚgroupr8rVr9ÚlenrPrQÚlstripÚabsrÚ_WorkRepr>rWÚlistÚtupleÚ ValueErrorÚappendÚjoinÚmapÚfloatrÚ from_floatÚ TypeError)ÚclsÚvaluer0r/ÚmÚintpartÚfracpartrWrXÚdigitsÚdigitr(r(r*r_4s¬ ( zDecimal.__new__cCs¬t|tƒr||ƒSt|tƒs$tdƒ‚tj|ƒs8tj|ƒrD|t|ƒƒStjd|ƒdkrZd}nd}t |ƒj ƒ\}}|jƒd}t|t |d|ƒ|ƒ}|tkr |S||ƒSdS)a.Converts a float to a decimal number, exactly. Note that Decimal.from_float(0.1) is not the same as Decimal('0.1'). Since 0.1 is not exactly representable in binary floating point, the value is stored as the nearest representable value which is 0x1.999999999999ap-4. The exact equivalent of the value in decimal is 0.1000000000000000055511151231257827021181583404541015625. >>> Decimal.from_float(0.1) Decimal('0.1000000000000000055511151231257827021181583404541015625') >>> Decimal.from_float(float('nan')) Decimal('NaN') >>> Decimal.from_float(float('inf')) Decimal('Infinity') >>> Decimal.from_float(-float('inf')) Decimal('-Infinity') >>> Decimal.from_float(-0.0) Decimal('-0') zargument must be int or float.gð?r&r.éN)r`rVrqrsÚ_mathZisinfZisnanÚreprZcopysignriÚas_integer_ratioÚ bit_lengthr7rar)rtÚfr>r6ÚdÚkÚresultr(r(r*rrÊs zDecimal.from_floatcCs(|jr$|j}|dkrdS|dkr$dSdS)zrReturns whether the number is not actually one. 0 if a number 1 if NaN 2 if sNaN r6r.rZr,r&)rQrP)r/rWr(r(r*Ú_isnanòszDecimal._isnancCs|jdkr|jrdSdSdS)zyReturns whether the number is infinite 0 if finite or not a number 1 if +INF -1 if -INF r[r.r&éÿÿÿÿ)rPr8)r/r(r(r*Ú_isinfinitys zDecimal._isinfinitycCs||jƒ}|dkrd}n|jƒ}|s&|rx|dkr4tƒ}|dkrJ|jtd|ƒS|dkr`|jtd|ƒS|rn|j|ƒS|j|ƒSdS)z½Returns whether the number is not actually one. if self, other are sNaN, signal if self, other are NaN return nan return 0 Done before operations. NFr,ÚsNaNr&)r„rrer r:)r/Úotherr0Úself_is_nanÚother_is_nanr(r(r*Ú_check_nanss" zDecimal._check_nanscCsv|dkrtƒ}|js|jrr|jƒr0|jtd|ƒS|jƒrF|jtd|ƒS|jƒr\|jtd|ƒS|jƒrr|jtd|ƒSdS)aCVersion of _check_nans used for the signaling comparisons compare_signal, __le__, __lt__, __ge__, __gt__. Signal InvalidOperation if either self or other is a (quiet or signaling) NaN. Signaling NaNs take precedence over quiet NaNs. Return 0 if neither operand is a NaN. Nzcomparison involving sNaNzcomparison involving NaNr&)rrQÚis_snanrer Úis_qnan)r/rˆr0r(r(r*Ú_compare_check_nans.s(zDecimal._compare_check_nanscCs|jp|jdkS)zuReturn True if self is nonzero; otherwise return False. NaNs and infinities are considered nonzero. rR)rQr9)r/r(r(r*Ú__bool__OszDecimal.__bool__cCs|js|jr8|jƒ}|jƒ}||kr(dS||kr4dSdS|sP|sDdSd|jS|s^d|jS|j|jkrndS|j|jkr~dS|jƒ}|jƒ}||krð|jd|j|j}|jd|j|j}||krÐdS||kräd |jSd |jSn ||krd|jSd|jSdS) z¸Compare the two non-NaN decimal instances self and other. Returns -1 if self < other, 0 if self == other and 1 if self > other. This routine is for internal use only.r&r.rRNr…r…r…r…r…r…r…r…)rQr†r8Úadjustedr9rP)r/rˆZself_infZ other_infÚ self_adjustedZother_adjustedÚself_paddedZother_paddedr(r(r*Ú_cmpVs> zDecimal._cmpcCs<t||dd\}}|tkr|S|j||ƒr.dS|j|ƒdkS)NT)Úequality_opFr&)Ú_convert_for_comparisonÚNotImplementedr‹r“)r/rˆr0r(r(r*Ú__eq__–szDecimal.__eq__cCs<t||ƒ\}}|tkr|S|j||ƒ}|r.dS|j|ƒdkS)NFr&)r•r–rŽr“)r/rˆr0r<r(r(r*Ú__lt__žszDecimal.__lt__cCs<t||ƒ\}}|tkr|S|j||ƒ}|r.dS|j|ƒdkS)NFr&)r•r–rŽr“)r/rˆr0r<r(r(r*Ú__le__§szDecimal.__le__cCs<t||ƒ\}}|tkr|S|j||ƒ}|r.dS|j|ƒdkS)NFr&)r•r–rŽr“)r/rˆr0r<r(r(r*Ú__gt__°szDecimal.__gt__cCs<t||ƒ\}}|tkr|S|j||ƒ}|r.dS|j|ƒdkS)NFr&)r•r–rŽr“)r/rˆr0r<r(r(r*Ú__ge__¹szDecimal.__ge__cCs>t|dd}|js|r0|jr0|j||ƒ}|r0|St|j|ƒƒS)zàCompare self to other. Return a decimal value: a or b is a NaN ==> Decimal('NaN') a < b ==> Decimal('-1') a == b ==> Decimal('0') a > b ==> Decimal('1') T)Úraiseit)Ú_convert_otherrQr‹rr“)r/rˆr0r<r(r(r*ÚcompareÂszDecimal.comparecCs’|jr4|jƒrtdƒ‚n|jƒr$tS|jr0tStS|jdkrNtd|jt ƒ}ntt |jt ƒ}t|jƒ|t }|dkr||n|}|dkrŽdS|S)zx.__hash__() <==> hash(x)z"Cannot hash a signaling NaN value.r&é r.r,r…éþÿÿÿ) rQrŒrsÚis_nanÚ_PyHASH_NANr8Ú_PyHASH_INFrPÚpowÚ_PyHASH_MODULUSÚ _PyHASH_10INVrVr9)r/Zexp_hashZhash_r<r(r(r*Ú__hash__Ôs zDecimal.__hash__cCst|jttt|jƒƒ|jƒS)zeRepresents the number as a triple tuple. To show the internals exactly as they are. )rr8rlrprVr9rP)r/r(r(r*Úas_tupleîszDecimal.as_tuplecCsØ|jr |jƒrtdƒ‚ntdƒ‚|s(dSt|jƒ}|jdkrR|d|jd}}nr|j}x(|dkr‚|ddkr‚|d}|d8}q\W|j}t||@jƒd|ƒ}|r¸||L}||8}d||>}|j rÐ|}||fS)aExpress a finite Decimal instance in the form n / d. Returns a pair (n, d) of integers. When called on an infinity or NaN, raises OverflowError or ValueError respectively. >>> Decimal('3.14').as_integer_ratio() (157, 50) >>> Decimal('-123e5').as_integer_ratio() (-12300000, 1) >>> Decimal('0.00').as_integer_ratio() (0, 1) z#cannot convert NaN to integer ratioz(cannot convert Infinity to integer ratior&r.rŸr{)r&r.) rQr¡rmÚ OverflowErrorrVr9rPÚminrr8)r/r6rZd5Zd2Zshift2r(r(r*r~õs, zDecimal.as_integer_ratiocCsdt|ƒS)z0Represents the number as an instance of Decimal.z Decimal('%s'))ra)r/r(r(r*Ú__repr__'szDecimal.__repr__Fc Csdddg|j}|jrL|jdkr&|dS|jdkr>|d|jS|d|jS|jt|jƒ}|jdkrt|dkrt|}n6|s~d }n,|jdkrš|d dd }n|d dd }|dkrÌd}d d||j}nN|t|jƒkrú|jd|t|jƒ}d}n |jd|…}d |j|d…}||kr*d}n*|dkr:tƒ}ddg|jd||}||||S)z–Return string representation of the number in scientific notation. Captures all of the information in the underlying representation. rTrUr[ZInfinityr6ÚNaNr‡r&ér.rRr\Ú.NÚeÚEz%+diúÿÿÿ)r8rQrPr9rgrÚcapitals) r/Úengr0r>Ú leftdigitsÚdotplacerwrxrWr(r(r*Ú__str__,s: zDecimal.__str__cCs|jd|dS)a,Convert to a string, using engineering notation if an exponent is needed. Engineering notation has an exponent which is a multiple of 3. This can leave up to 3 digits to the left of the decimal place and may require the addition of either one or two trailing zeros. T)r²r0)rµ)r/r0r(r(r*Ú to_eng_string`szDecimal.to_eng_stringcCsT|jr|j|d}|r|S|dkr(tƒ}|rB|jtkrB|jƒ}n|jƒ}|j|ƒS)zRReturns a copy with the sign switched. Rounds, if it has reason. )r0N)rQr‹rr@rÚcopy_absÚcopy_negateÚ_fix)r/r0r<r(r(r*Ú__neg__is zDecimal.__neg__cCsT|jr|j|d}|r|S|dkr(tƒ}|rB|jtkrB|jƒ}nt|ƒ}|j|ƒS)zhReturns a copy, unless it is a sNaN. Rounds the number (if more than precision digits) )r0N)rQr‹rr@rr·rr¹)r/r0r<r(r(r*Ú__pos__s zDecimal.__pos__TcCsJ|s|jƒS|jr&|j|d}|r&|S|jr:|j|d}n|j|d}|S)zÉReturns the absolute value of self. If the keyword argument 'round' is false, do not round. The expression self.__abs__(round=False) is equivalent to self.copy_abs(). )r0)r·rQr‹r8rºr»)r/Úroundr0r<r(r(r*Ú__abs__”szDecimal.__abs__c Cslt|ƒ}|tkr|S|dkr"tƒ}|js.|jr‚|j||ƒ}|rB|S|jƒrr|j|jkrj|jƒrj|jtdƒSt |ƒS|jƒr‚t |ƒSt |j|jƒ}d}|jt kr®|j|jkr®d}|rê|rêt |j|jƒ}|rÐd}t|d|ƒ}|j|ƒ}|S|s"t||j|jdƒ}|j||jƒ}|j|ƒ}|S|sZt||j|jdƒ}|j||jƒ}|j|ƒ}|St|ƒ}t|ƒ}t|||jƒ\}}tƒ} |j|jkrþ|j|jkr¸t|d|ƒ}|j|ƒ}|S|j|jkrÐ||}}|jdkröd| _|j|j|_|_nd| _n&|jdkrd| _d\|_|_nd| _|jdkr@|j|j| _n|j|j| _|j| _t | ƒ}|j|ƒ}|S)zbReturns self + other. -INF + INF (or the reverse) cause InvalidOperation errors. Nz -INF + INFr&r.rR)r&r&)rr–rrQr‹r†r8rer rrªrPr@rr7r¹ÚmaxrAÚ_rescalerjÚ _normalizer>rVrW) r/rˆr0r<rWZnegativezeror>Úop1Úop2rƒr(r(r*Ú__add__ªs| zDecimal.__add__cCsHt|ƒ}|tkr|S|js |jr6|j||d}|r6|S|j|jƒ|dS)zReturn self - other)r0)rr–rQr‹rÃr¸)r/rˆr0r<r(r(r*Ú__sub__szDecimal.__sub__cCs"t|ƒ}|tkr|S|j||dS)zReturn other - self)r0)rr–rÄ)r/rˆr0r(r(r*Ú__rsub__szDecimal.__rsub__cCsDt|ƒ}|tkr|S|dkr"tƒ}|j|jA}|js:|jrŽ|j||ƒ}|rN|S|jƒrn|sf|jtdƒSt |S|jƒrŽ|s†|jtdƒSt |S|j |j }|s¦|rÀt|d|ƒ}|j|ƒ}|S|j dkræt||j |ƒ}|j|ƒ}|S|j dkrt||j |ƒ}|j|ƒ}|St|ƒ}t|ƒ}t|t|j|jƒ|ƒ}|j|ƒ}|S)z\Return self * other. (+-) INF * 0 (or its reverse) raise InvalidOperation. Nz(+-)INF * 0z0 * (+-)INFrRÚ1)rr–rr8rQr‹r†rer r=rPr7r¹r9rjrarV)r/rˆr0Z resultsignr<Z resultexprÁrÂr(r(r*Ú__mul__sH zDecimal.__mul__cCsÊt|ƒ}|tkrtS|dkr"tƒ}|j|jA}|js:|jrž|j||ƒ}|rN|S|jƒrj|jƒrj|jtdƒS|jƒrzt |S|jƒrž|jt dƒt|d|jƒƒS|sÀ|s²|jt dƒS|jtd|ƒS|sÖ|j|j}d}nÚt|jƒt|jƒ|jd}|j|j|}t|ƒ}t|ƒ} |dkr:t|jd || jƒ\}} nt|j| jd |ƒ\}} | rt|d dkr°|d7}n<|j|j}x.||kr®|d dkr®|d }|d7}q‚Wt|t|ƒ|ƒ}|j|ƒS)zReturn self / other.Nz(+-)INF/(+-)INFzDivision by infinityrRz0 / 0zx / 0r&r.rŸr{)rr–rr8rQr‹r†rer r=rr7ÚEtinyrr rPrgr9rArjÚdivmodrVrar¹)r/rˆr0r>r<rWÚcoeffÚshiftrÁrÂÚ remainderÚ ideal_expr(r(r*Ú__truediv__QsP zDecimal.__truediv__cCs |j|jA}|jƒr|j}nt|j|jƒ}|jƒ|jƒ}|sP|jƒsP|dkrjt|ddƒ|j||jƒfS||jkrt |ƒ}t |ƒ}|j |j kr®|jd|j |j 9_n|jd|j |j 9_t|j|jƒ\}} |d|jkrt|t |ƒdƒt|jt | ƒ|ƒfS|jtdƒ} | | fS)z½Return (self // other, self % other), to context.prec precision. Assumes that neither self nor other is a NaN, that self is not infinite and that other is nonzero. r,rRr&rŸz%quotient too large in //, % or divmodr )r8r†rPrªrr7r¿r@rArjrWrVrÉrarer)r/rˆr0r>rÍÚexpdiffrÁrÂÚqÚrr<r(r(r*Ú_divideŒs* zDecimal._dividecCs"t|ƒ}|tkr|S|j||dS)z)Swaps self/other and returns __truediv__.)r0)rr–rÎ)r/rˆr0r(r(r*Ú__rtruediv__szDecimal.__rtruediv__cCsÖt|ƒ}|tkr|S|dkr"tƒ}|j||ƒ}|r:||fS|j|jA}|jƒr~|jƒrj|jtdƒ}||fSt||jtdƒfS|s´|sš|jt dƒ}||fS|jt d|ƒ|jtdƒfS|j||ƒ\}}|j|ƒ}||fS)z6 Return (self // other, self % other) Nzdivmod(INF, INF)zINF % xzdivmod(0, 0)zx // 0zx % 0) rr–rr‹r8r†rer r=rr rÒr¹)r/rˆr0r<r>ZquotientrÌr(r(r*Ú __divmod__´s0 zDecimal.__divmod__cCs"t|ƒ}|tkr|S|j||dS)z(Swaps self/other and returns __divmod__.)r0)rr–rÔ)r/rˆr0r(r(r*Ú__rdivmod__ØszDecimal.__rdivmod__cCsˆt|ƒ}|tkr|S|dkr"tƒ}|j||ƒ}|r6|S|jƒrJ|jtdƒS|sj|r^|jtdƒS|jtdƒS|j||ƒd}|j |ƒ}|S)z self % other NzINF % xzx % 0z0 % 0r.) rr–rr‹r†rer rrÒr¹)r/rˆr0r<rÌr(r(r*Ú__mod__ßs" zDecimal.__mod__cCs"t|ƒ}|tkr|S|j||dS)z%Swaps self/other and returns __mod__.)r0)rr–rÖ)r/rˆr0r(r(r*Ú__rmod__úszDecimal.__rmod__cCsÐ|dkrtƒ}t|dd}|j||ƒ}|r.|S|jƒrB|jtdƒS|sb|rV|jtdƒS|jtdƒS|jƒr|t|ƒ}|j|ƒSt |j |j ƒ}|s¦t|jd|ƒ}|j|ƒS|j ƒ|j ƒ}||jdkrÎ|jtƒS|dkrî|j||jƒ}|j|ƒSt|ƒ}t|ƒ}|j|jkr(|jd |j|j9_n|jd |j|j9_t|j|jƒ\}} d | |d@|jkr~| |j8} |d7}|d |jkr˜|jtƒS|j} | dkr¶d| } | } t| t| ƒ|ƒ}|j|ƒS) zI Remainder nearest to 0- abs(remainder-near) <= other/2 NT)rœzremainder_near(infinity, x)zremainder_near(x, 0)zremainder_near(0, 0)rRr.r,rŸr&r )rrr‹r†rer rrr¹rªrPr7r8rrArr¿r@rjrWrVrÉra)r/rˆr0r<Úideal_exponentrÏrÁrÂrÐrÑr>r(r(r*Úremainder_nearsZ zDecimal.remainder_nearcCsœt|ƒ}|tkr|S|dkr"tƒ}|j||ƒ}|r6|S|jƒrb|jƒrR|jtdƒSt|j|jAS|sŒ|r€|jt d|j|jAƒS|jt dƒS|j||ƒdS)z self // otherNz INF // INFzx // 0z0 // 0r&)rr–rr‹r†rer r=r8r rrÒ)r/rˆr0r<r(r(r*Ú__floordiv__Ls$zDecimal.__floordiv__cCs"t|ƒ}|tkr|S|j||dS)z*Swaps self/other and returns __floordiv__.)r0)rr–rÚ)r/rˆr0r(r(r*Ú __rfloordiv__hszDecimal.__rfloordiv__cCs8|jƒr(|jƒrtdƒ‚|jr"dnd}nt|ƒ}t|ƒS)zFloat representation.z%Cannot convert signaling NaN to floatz-nanÚnan)r„rŒrmr8rarq)r/Úsr(r(r*Ú __float__oszDecimal.__float__cCst|jr(|jƒrtdƒ‚n|jƒr(tdƒ‚d|j}|jdkrT|t|jƒd|jS|t|jd|j…pjdƒSdS) z1Converts self to an int, truncating if necessary.zCannot convert NaN to integerz"Cannot convert infinity to integerr.r&rŸNrRr…) rQr„rmr†r©r8rPrVr9)r/rÝr(r(r*Ú__int__ys zDecimal.__int__cCs|S)Nr()r/r(r(r*ÚrealˆszDecimal.realcCstdƒS)Nr&)r)r/r(r(r*ÚimagŒszDecimal.imagcCs|S)Nr()r/r(r(r*Ú conjugateszDecimal.conjugatecCstt|ƒƒS)N)Úcomplexrq)r/r(r(r*Ú__complex__“szDecimal.__complex__cCsR|j}|j|j}t|ƒ|krJ|t|ƒ|d…jdƒ}t|j||jdƒSt|ƒS)z2Decapitate the payload of a NaN to fit the contextNrRT) r9rAÚclamprgrhr7r8rPr)r/r0ZpayloadZmax_payload_lenr(r(r*r:–szDecimal._fix_nancCsX|jr |jƒr|j|ƒSt|ƒS|jƒ}|jƒ}|s€|j|g|j}tt |j |ƒ|ƒ}||j krx|jtƒt |jd|ƒSt|ƒSt|jƒ|j |j}||krÆ|jtd|jƒ}|jtƒ|jtƒ|S||k}|rÖ|}|j |krüt|jƒ|j |} | dkrt |jd|dƒ}d} |j|j} | || ƒ}|jd| …p>d}|dkr~tt|ƒdƒ}t|ƒ|jkr~|dd…}|d7}||krš|jtd|jƒ}nt |j||ƒ}|r¾|r¾|jtƒ|rÎ|jtƒ|rÞ|jtƒ|jtƒ|sø|jtƒ|S|r|jtƒ|jdkrP|j |krP|jtƒ|jd|j |} t |j| |ƒSt|ƒS)zÜRound if it is necessary to keep self within prec precision. Rounds and fixes the exponent. Does not raise on a sNaN. Arguments: self - Decimal instance context - context used. rRz above Emaxr&rÆr.Nr…)rQr„r:rrÈÚEtoprBrårªr¾rPrerr7r8rgr9rArrrÚ_pick_rounding_functionr@rarVrr )r/r0rÈræÚexp_maxZnew_expZexp_minr<Zself_is_subnormalryZrounding_methodÚchangedrÊr’r(r(r*r¹¢sn zDecimal._fixcCst|j|ƒrdSdSdS)z(Also known as round-towards-0, truncate.r&r.Nr…)Ú _all_zerosr9)r/rAr(r(r*Ú_round_downszDecimal._round_downcCs|j|ƒS)zRounds away from 0.)rë)r/rAr(r(r*Ú _round_upszDecimal._round_upcCs*|j|dkrdSt|j|ƒr"dSdSdS)zRounds 5 up (away from 0)Z56789r.r&Nr…)r9rê)r/rAr(r(r*Ú_round_half_ups zDecimal._round_half_upcCst|j|ƒrdS|j|ƒSdS)zRound 5 downr.Nr…)Ú_exact_halfr9rí)r/rAr(r(r*Ú_round_half_downszDecimal._round_half_downcCs8t|j|ƒr*|dks&|j|ddkr*dS|j|ƒSdS)z!Round 5 to even, rest to nearest.r&r.Ú02468Nr…)rîr9rí)r/rAr(r(r*Ú_round_half_even#szDecimal._round_half_evencCs |jr|j|ƒS|j|ƒSdS)z(Rounds up (not away from 0 if negative.)N)r8rë)r/rAr(r(r*Ú_round_ceiling+s zDecimal._round_ceilingcCs |js|j|ƒS|j|ƒSdS)z'Rounds down (not towards 0 if negative)N)r8rë)r/rAr(r(r*Ú_round_floor2s zDecimal._round_floorcCs0|r |j|ddkr |j|ƒS|j|ƒSdS)z)Round down unless digit prec-1 is 0 or 5.r.Z05N)r9rë)r/rAr(r(r*Ú_round_05up9s zDecimal._round_05up)rrrrrrrrcCsb|dk r2t|tƒstdƒ‚tdd|ƒ}|j|ƒS|jrR|jƒrJtdƒ‚ntdƒ‚t|j dt ƒƒS)aÊRound self to the nearest integer, or to a given precision. If only one argument is supplied, round a finite Decimal instance self to the nearest integer. If self is infinite or a NaN then a Python exception is raised. If self is finite and lies exactly halfway between two integers then it is rounded to the integer with even last digit. >>> round(Decimal('123.456')) 123 >>> round(Decimal('-456.789')) -457 >>> round(Decimal('-3.0')) -3 >>> round(Decimal('2.5')) 2 >>> round(Decimal('3.5')) 4 >>> round(Decimal('Inf')) Traceback (most recent call last): ... OverflowError: cannot round an infinity >>> round(Decimal('NaN')) Traceback (most recent call last): ... ValueError: cannot round a NaN If a second argument n is supplied, self is rounded to n decimal places using the rounding mode for the current context. For an integer n, round(self, -n) is exactly equivalent to self.quantize(Decimal('1En')). >>> round(Decimal('123.456'), 0) Decimal('123') >>> round(Decimal('123.456'), 2) Decimal('123.46') >>> round(Decimal('123.456'), -2) Decimal('1E+2') >>> round(Decimal('-Infinity'), 37) Decimal('NaN') >>> round(Decimal('sNaN123'), 0) Decimal('NaN123') Nz+Second argument to round should be integralr&rÆzcannot round a NaNzcannot round an infinity)r`rVrsr7ÚquantizerQr¡rmr©r¿r)r/r6rWr(r(r*Ú __round__Ks/ zDecimal.__round__cCs0|jr |jƒrtdƒ‚ntdƒ‚t|jdtƒƒS)zãReturn the floor of self, as an integer. For a finite Decimal instance self, return the greatest integer n such that n <= self. If self is infinite or a NaN then a Python exception is raised. zcannot round a NaNzcannot round an infinityr&)rQr¡rmr©rVr¿r)r/r(r(r*Ú __floor__‰s zDecimal.__floor__cCs0|jr |jƒrtdƒ‚ntdƒ‚t|jdtƒƒS)zâReturn the ceiling of self, as an integer. For a finite Decimal instance self, return the least integer n such that n >= self. If self is infinite or a NaN then a Python exception is raised. zcannot round a NaNzcannot round an infinityr&)rQr¡rmr©rVr¿r)r/r(r(r*Ú__ceil__˜s zDecimal.__ceil__cCst|dd}t|dd}|js$|jrÚ|dkr2tƒ}|jdkrJ|jtd|ƒS|jdkrb|jtd|ƒS|jdkrr|}nf|jdkr‚|}nV|jdkr®|sœ|jtdƒSt|j|jA}n*|jdkrØ|sÈ|jtd ƒSt|j|jA}n0t|j|jAt t |jƒt |jƒƒ|j|jƒ}|j||ƒS) a:Fused multiply-add. Returns self*other+third with no rounding of the intermediate product self*other. self and other are multiplied together, with no rounding of the result. The third operand is then added to the result, and a single final rounding is performed. T)rœNrZr‡r6r[zINF * 0 in fmaz0 * INF in fma) rrQrrPrer r=r8r7rarVr9rÃ)r/rˆZthirdr0Úproductr(r(r*Úfma§s6 zDecimal.fmacCsÚt|ƒ}|tkr|St|ƒ}|tkr(|S|dkr6tƒ}|jƒ}|jƒ}|jƒ}|sZ|sZ|rÂ|dkrp|jtd|ƒS|dkr†|jtd|ƒS|dkrœ|jtd|ƒS|rª|j|ƒS|r¸|j|ƒS|j|ƒS|jƒoØ|jƒoØ|jƒsæ|jtdƒS|dkrú|jtdƒS|s|jtdƒS|jƒ|j kr(|jtdƒS|rD|rD|jtd ƒS|j ƒrTd}n|j}tt |ƒƒ}t|jƒƒ}t|jƒƒ} |j |td |j|ƒ|}x t| jƒD]} t|d |ƒ}q¦Wt|| j |ƒ}t|t|ƒdƒS)z!Three argument version of __pow__Nr,r‡z@pow() 3rd argument not allowed unless all arguments are integersr&zApow() 2nd argument cannot be negative when 3rd argument specifiedzpow() 3rd argument cannot be 0zSinsufficient precision: pow() 3rd argument must not have more than precision digitszXat least one of pow() 1st argument and 2nd argument must be nonzero; 0**0 is not definedrŸ)rr–rr„rer r:Ú _isintegerrrAÚ_isevenr8rirVrjÚto_integral_valuer¤rWÚranger7ra)r/rˆÚmodulor0r‰rŠZ modulo_is_nanr>ÚbaseÚexponentÚir(r(r*Ú _power_moduloÓsl zDecimal._power_modulocCst|ƒ}|j|j}}x |ddkr6|d}|d7}qWt|ƒ}|j|j}}x |ddkrn|d}|d7}qPW|dkr||9}x |ddkr¢|d}|d7}q„W|dkr°dS|d|} |jdkrÌ| } |jƒoÜ|jdkr|jt|ƒ} t| | |dƒ}nd}tddd|| |ƒS|jdkrÄ|d}|dkrÊ||@|krPdSt |ƒd} |dd}|t t|ƒƒkr~dSt| ||ƒ} t|||ƒ}| dks®|dkr²dS| |krÀdSd | }nÎ|d kr”t |ƒdd} t d | |ƒ\}}|rdSx$|d dkr$|d }| d8} qW|dd}|t t|ƒƒkrHdSt| ||ƒ} t|||ƒ}| dksx|dkr|dS| |krŠdSd| }ndS|d|krªdS| |}tdt|ƒ|ƒS|dkrâ|d|d}}nä|dkrt tt||ƒƒƒ|krdSt |ƒ}|dkr>t tt|ƒ|ƒƒ|kr>dS|d|}}x:|d|dkoldknrˆ|d}|d}qPWx:|d |d ko¨dknrÄ|d }|d }qŒW|dkrt|dkrè||krèdSt ||ƒ\}}|dkrdSdt |ƒ|>}x>t |||dƒ\}}||kr>Pn||d||}qW||koh|dkspdS|}|dkr˜||dt|ƒkr˜dS||}||9}|d|krºdSt|ƒ}|jƒrþ|jdkrþ|jt|ƒ} t|| |t |ƒƒ}nd}td|d|||ƒS)ahAttempt to compute self**other exactly. Given Decimals self and other and an integer p, attempt to compute an exact result for the power self**other, with p digits of precision. Return None if self**other is not exactly representable in p digits. Assumes that elimination of special cases has already been performed: self and other must both be nonspecial; self must be positive and not numerically equal to 1; other must be nonzero. For efficiency, other._exp should not be too large, so that 10**abs(other._exp) is a feasible calculation.rŸr&r.NrÆrRr,éréé]éAr{ér\éd)r,rrr)rjrVrWr>rûr8rPrªr7Ú_nbitsrgraÚ_decimal_lshift_exactrÉriÚ _log10_lb)r/rˆÚpÚxÚxcÚxeÚyÚycÚyerrØZzerosZ last_digitr¯ZemaxrÌrvr6Zxc_bitsÚremÚarÐrÑZstr_xcr(r(r*Ú_power_exact(sÒ: &&&& zDecimal._power_exactcCs@|dk r|j|||ƒSt|ƒ}|tkr*|S|dkr8tƒ}|j||ƒ}|rL|S|sd|s`|jtdƒStSd}|jdkr |j ƒrˆ|j ƒs˜d}n|r˜|jtdƒS|jƒ}|sÂ|jdkrºt|ddƒSt |S|jƒrè|jdkrÜt |St|ddƒS|tkrŽ|j ƒrZ|jdkrd}n||jkr"|j}nt|ƒ}|j|}|d|jkrxd|j}|jtƒn|jtƒ|jtƒd|j}t|dd||ƒS|jƒ}|jƒrÈ|jdk|dkkrÀt|ddƒSt |Sd}d} |jƒ|jƒ} |dk|jdkkr| tt|jƒƒkrHt|d|jdƒ}n,|jƒ}| tt|ƒƒkrHt|d|dƒ}|dkrŒ|j||jdƒ}|dk rŒ|dkrˆtd|j|jƒ}d } |dkr:|j}t|ƒ} | j| j}}t|ƒ}|j|j}}|jdkrÚ|}d }xJt||||||ƒ\}}|ddtt|ƒƒ|drP|d 7}qàWt|t|ƒ|ƒ}| r2|j ƒr2t|jƒ|jkr|jdt|jƒ}t|j|jd||j|ƒ}|j ƒ}|j!ƒxt"D]}d|j#|<q¦W|j$|ƒ}|jtƒ|j%t&rä|jt'ƒ|j%t(r|jt(d |jƒx:t't&ttt)fD]}|j%|r|j|ƒqWn |j$|ƒ}|S)aHReturn self ** other [ % modulo]. With two arguments, compute self**other. With three arguments, compute (self**other) % modulo. For the three argument form, the following restrictions on the arguments hold: - all three arguments must be integral - other must be nonnegative - either self or other (or both) must be nonzero - modulo must be nonzero and must have at most p digits, where p is the context precision. If any of these restrictions is violated the InvalidOperation flag is raised. The result of pow(self, other, modulo) is identical to the result that would be obtained by computing (self**other) % modulo with unbounded precision, but is computed more efficiently. It is always exact. Nz0 ** 0r&r.z+x ** y with x negative and y not an integerrRrÆFTr\r{rŸz above Emax)*rrr–rr‹rer Ú_Oner8rûrür¸r7r=r†rArVrPrrrÚ_log10_exp_boundrgrarBrÈrr9rjrWr>Ú_dpowerrIrJÚ_signalsÚtrapsr¹Úflagsr rrr)r/rˆrÿr0r<Zresult_signZ multiplierrWZself_adjÚexactZboundrÈr rrrrrrÚextrarÊrÏZ newcontextZ exceptionr(r(r*Ú__pow__ sÊ " zDecimal.__pow__cCs"t|ƒ}|tkr|S|j||dS)z%Swaps self/other and returns __pow__.)r0)rr–r)r/rˆr0r(r(r*Ú__rpow__í szDecimal.__rpow__cCs¼|dkrtƒ}|jr(|j|d}|r(|S|j|ƒ}|jƒr>|S|sPt|jddƒS|j|jƒg|j }t |jƒ}|j}x.|j|ddkr¢||kr¢|d7}|d8}qvWt|j|jd|…|ƒS)z?Normalize- strip trailing 0s, change anything equal to 0 to 0e0N)r0rRr&r.) rrQr‹r¹r†r7r8rBrærårgr9rP)r/r0r<ÚduprèÚendrWr(r(r*Ú normalizeô s$ zDecimal.normalizecCs¨t|dd}|dkrtƒ}|dkr(|j}|js4|jr||j||ƒ}|rH|S|jƒsX|jƒr||jƒrp|jƒrpt|ƒS|jtdƒS|j ƒ|j ko–|jkns¨|jtdƒS|sÆt|j d|j ƒ}|j|ƒS|jƒ}||jkrä|jtdƒS||j d|jkr|jtd ƒS|j|j |ƒ}|jƒ|jkr0|jtdƒSt|jƒ|jkrN|jtd ƒS|rn|jƒ|jkrn|jtƒ|j |j krš||kr|jtƒ|jtƒ|j|ƒ}|S) z‡Quantize self so its exponent is the same as that of exp. Similar to self._rescale(exp._exp) but with error checking. T)rœNzquantize with one INFz)target exponent out of bounds in quantizerRz9exponent of quantize result too large for current contextr.z7quantize result has too many digits for current context)rrr@rQr‹r†rrer rÈrPrBr7r8r¹rrAr¿rgr9ÚEminr rr)r/rWr@r0r<r‘r(r(r*rõ sT zDecimal.quantizecCsDt|dd}|js|jr8|jƒr(|jƒp6|jƒo6|jƒS|j|jkS)a=Return True if self and other have the same exponent; otherwise return False. If either operand is a special value, the following rules are used: * return True if both operands are infinities * return True if both operands are NaNs * otherwise, return False. T)rœ)rrQr¡Úis_infiniterP)r/rˆr0r(r(r*Úsame_quantumJ s zDecimal.same_quantumcCsÆ|jrt|ƒS|s t|jd|ƒS|j|krHt|j|jd|j||ƒSt|jƒ|j|}|dkrzt|jd|dƒ}d}|j|}|||ƒ}|jd|…pžd}|dkr¸tt |ƒdƒ}t|j||ƒS)asRescale self so that the exponent is exp, either by padding with zeros or by truncating digits, using the given rounding mode. Specials are returned without change. This operation is quiet: it raises no flags, and uses no information from the context. exp = exp to scale to (an integer) rounding = rounding mode rRr&rÆr.N) rQrr7r8rPr9rgrçrarV)r/rWr@ryZ this_functionrérÊr(r(r*r¿Y s" zDecimal._rescalecCsh|dkrtdƒ‚|js|r$t|ƒS|j|jƒd||ƒ}|jƒ|jƒkrd|j|jƒd||ƒ}|S)a"Round a nonzero, nonspecial Decimal to a fixed number of significant figures, using the given rounding mode. Infinities, NaNs and zeros are returned unaltered. This operation is quiet: it raises no flags, and uses no information from the context. r&z'argument should be at least 1 in _roundr.)rmrQrr¿r)r/Úplacesr@r<r(r(r*Ú_round{ s zDecimal._roundcCsŽ|jr"|j|d}|r|St|ƒS|jdkr4t|ƒS|sFt|jddƒS|dkrTtƒ}|dkrb|j}|jd|ƒ}||kr€|j t ƒ|j tƒ|S)aVRounds to a nearby integer. If no rounding mode is specified, take the rounding mode from the context. This method raises the Rounded and Inexact flags when appropriate. See also: to_integral_value, which does exactly the same as this method except that it doesn't raise Inexact or Rounded. )r0r&rRN)rQr‹rrPr7r8rr@r¿rerr)r/r@r0r<r(r(r*Úto_integral_exact’ s$ zDecimal.to_integral_exactcCs`|dkrtƒ}|dkr|j}|jr>|j|d}|r6|St|ƒS|jdkrPt|ƒS|jd|ƒSdS)z@Rounds to the nearest integer, without raising inexact, rounded.N)r0r&)rr@rQr‹rrPr¿)r/r@r0r<r(r(r*rý¯ s zDecimal.to_integral_valuecCsà|dkrtƒ}|jrB|j|d}|r(|S|jƒrB|jdkrBt|ƒS|sdt|jd|jdƒ}|j|ƒS|jdkrz|j t dƒS|jd}t|ƒ}|j d?}|j d@r¾|jd}t|jƒd?d}n|j}t|jƒdd?}||}|dkrø|d |9}d } nt|d |ƒ\}} | } ||8}d|}x(||}||kr:Pn||d?}q$W| o\|||k} | r”|dkr||d|}n|d|9}||7}n|ddkrª|d7}tdt|ƒ|ƒ}|jƒ}|jtƒ} |j|ƒ}| |_|S)zReturn the square root of self.N)r0r&rRr,r.zsqrt(-x), x > 0rŸr Tr{)rrQr‹r†r8rr7rPr¹rer rArjrWrVrgr9rÉraÚ _shallow_copyÚ _set_roundingrr@)r/r0r<rAÚopr¯ÚcÚlrËrrÌr6rÐr@r(r(r*ÚsqrtÂ s` zDecimal.sqrtcCs¶t|dd}|dkrtƒ}|js&|jr~|jƒ}|jƒ}|s>|r~|dkrX|dkrX|j|ƒS|dkrr|dkrr|j|ƒS|j||ƒS|j|ƒ}|dkrš|j|ƒ}|dkr¨|}n|}|j|ƒS)z Returns the larger value. Like max(self, other) except if one is not a number, returns NaN (and signals if one is sNaN). Also rounds. T)rœNr.r&r…)rrrQr„r¹r‹r“Ú compare_total)r/rˆr0ÚsnÚonr-r<r(r(r*r¾%s& zDecimal.maxcCs¶t|dd}|dkrtƒ}|js&|jr~|jƒ}|jƒ}|s>|r~|dkrX|dkrX|j|ƒS|dkrr|dkrr|j|ƒS|j||ƒS|j|ƒ}|dkrš|j|ƒ}|dkr¨|}n|}|j|ƒS)z¡Returns the smaller value. Like min(self, other) except if one is not a number, returns NaN (and signals if one is sNaN). Also rounds. T)rœNr.r&r…)rrrQr„r¹r‹r“r0)r/rˆr0r1r2r-r<r(r(r*rªOs& zDecimal.mincCs8|jr dS|jdkrdS|j|jd…}|dt|ƒkS)z"Returns whether self is an integerFr&TNrR)rQrPr9rg)r/Úrestr(r(r*rûqs zDecimal._isintegercCs(|s|jdkrdS|jd|jdkS)z:Returns True if self is even. Assumes self is an integer.r&Tr.rðr…)rPr9)r/r(r(r*rüzszDecimal._isevencCs.y|jt|jƒdStk r(dSXdS)z$Return the adjusted exponent of selfr.r&N)rPrgr9rs)r/r(r(r*r€szDecimal.adjustedcCs|S)z«Returns the same Decimal object. As we do not have different encodings for the same number, the received object already is in its canonical form. r()r/r(r(r*Ú canonicalˆszDecimal.canonicalcCs.t|dd}|j||ƒ}|r |S|j||dS)z¶Compares self to the other operand numerically. It's pretty much like compare(), but all NaNs signal, with signaling NaNs taking precedence over quiet NaNs. T)rœ)r0)rrŽrž)r/rˆr0r<r(r(r*Úcompare_signals zDecimal.compare_signalcCsft|dd}|jr|jrtS|jr0|jr0tS|j}|jƒ}|jƒ}|sP|r||kr¤t|jƒ|jf}t|jƒ|jf}||krŒ|rˆtStS||kr |rœtStStS|rÚ|dkr´tS|dkrÀtS|dkrÌtS|dkrØtSn4|dkrætS|dkròtS|dkrtS|dkrtS||krtS||kr*tS|j|jkrF|rBtStS|j|jkrb|r^tStStS)zõCompares self to other using the abstract representations. This is not like the standard compare, which use their numerical value. Note that a total ordering is defined for all possible abstract representations. T)rœr.r,) rr8Ú_NegativeOnerr„rgr9Ú_ZerorP)r/rˆr0r>Zself_nanZ other_nanZself_keyZ other_keyr(r(r*r0œsf zDecimal.compare_totalcCs&t|dd}|jƒ}|jƒ}|j|ƒS)z–Compares self to other using abstract repr., ignoring sign. Like compare_total, but with operand's sign ignored and assumed to be 0. T)rœ)rr·r0)r/rˆr0rÝÚor(r(r*Úcompare_total_magåszDecimal.compare_total_magcCstd|j|j|jƒS)z'Returns a copy with the sign set to 0. r&)r7r9rPrQ)r/r(r(r*r·ðszDecimal.copy_abscCs2|jrtd|j|j|jƒStd|j|j|jƒSdS)z&Returns a copy with the sign inverted.r&r.N)r8r7r9rPrQ)r/r(r(r*r¸ôszDecimal.copy_negatecCs"t|dd}t|j|j|j|jƒS)z$Returns self with the sign of other.T)rœ)rr7r8r9rPrQ)r/rˆr0r(r(r*Ú copy_signûs zDecimal.copy_signcCsì|dkrtƒ}|j|d}|r"|S|jƒdkr2tS|s:tS|jƒdkrNt|ƒS|j}|jƒ}|jdkr–|t t |jddƒƒkr–tdd|jdƒ}n0|jdkrÔ|t t |j ƒddƒƒkrÔtdd|j ƒdƒ}nò|jdkoæ||kr tddd|dd|ƒ}n¼|jdkrB||dkrBtdd|d|dƒ}n„t|ƒ}|j|j}}|jdkrj|}d}xFt||||ƒ\} } | d d t t | ƒƒ|dr¨P|d7}qpWtdt | ƒ| ƒ}|jƒ}|jtƒ}|j|ƒ}||_|S)zReturns e ** self.N)r0r.r&r\rÆrRr?r{rŸr…)rr‹r†r7rrrArr8rgrarBr7rÈrjrVrWr>Ú_dexpr*r+rr¹r@)r/r0r<r Úadjr,r-r¯rrÊrWr@r(r(r*rWsJ$( " zDecimal.expcCsdS)zÃReturn True if self is canonical; otherwise return False. Currently, the encoding of a Decimal instance is always canonical, so this method returns True for any Decimal. Tr()r/r(r(r*Úis_canonicalLszDecimal.is_canonicalcCs|jS)zReturn True if self is finite; otherwise return False. A Decimal instance is considered finite if it is neither infinite nor a NaN. )rQ)r/r(r(r*Ú is_finiteTszDecimal.is_finitecCs |jdkS)z8Return True if self is infinite; otherwise return False.r[)rP)r/r(r(r*r%\szDecimal.is_infinitecCs |jdkS)z>Return True if self is a qNaN or sNaN; otherwise return False.r6rZ)r6rZ)rP)r/r(r(r*r¡`szDecimal.is_nancCs,|js|rdS|dkrtƒ}|j|jƒkS)z?Return True if self is a normal number; otherwise return False.FN)rQrr$r)r/r0r(r(r*Ú is_normalds zDecimal.is_normalcCs |jdkS)z;Return True if self is a quiet NaN; otherwise return False.r6)rP)r/r(r(r*rlszDecimal.is_qnancCs |jdkS)z8Return True if self is negative; otherwise return False.r.)r8)r/r(r(r*Ú is_signedpszDecimal.is_signedcCs |jdkS)z?Return True if self is a signaling NaN; otherwise return False.rZ)rP)r/r(r(r*rŒtszDecimal.is_snancCs,|js|rdS|dkrtƒ}|jƒ|jkS)z9Return True if self is subnormal; otherwise return False.FN)rQrrr$)r/r0r(r(r*Úis_subnormalxs zDecimal.is_subnormalcCs|jo|jdkS)z6Return True if self is a zero; otherwise return False.rR)rQr9)r/r(r(r*Úis_zero€szDecimal.is_zerocCsÆ|jt|jƒd}|dkr4tt|ddƒƒdS|dkrXttd|ddƒƒdSt|ƒ}|j|j}}|dkr¨t|d|ƒ}t|ƒ}t|ƒt|ƒ||kS|ttd||ƒƒdS)zÌCompute a lower bound for the adjusted exponent of self.ln(). In other words, compute r such that self.ln() >= 10**r. Assumes that self is finite and positive and that self != 1. r.érŸr,r&r r…)rPrgr9rarjrVrW)r/r<r,r-r¯ÚnumÚdenr(r(r*Ú _ln_exp_bound„szDecimal._ln_exp_boundc Cs|dkrtƒ}|j|d}|r"|S|s*tS|jƒdkr:tS|tkrFtS|jdkr\|jt dƒSt |ƒ}|j|j}}|j }||jƒd}x>t|||ƒ}|ddttt|ƒƒƒ|dr¼P|d7}qŠWtt|d kƒtt|ƒƒ|ƒ}|jƒ}|jtƒ} |j|ƒ}| |_|S) z/Returns the natural (base e) logarithm of self.N)r0r.zln of a negative valuer,r{rŸr\r&)rr‹Ú_NegativeInfinityr†Ú _Infinityrr7r8rer rjrVrWrArFÚ_dlogrgrarir7r*r+rr¹r@) r/r0r<r,r-r¯r r'rÊr@r(r(r*Úlns: $ z Decimal.lncCsÊ|jt|jƒd}|dkr,tt|ƒƒdS|dkrHttd|ƒƒdSt|ƒ}|j|j}}|dkr t|d|ƒ}td|ƒ}t|ƒt|ƒ||kdStd||ƒ}t|ƒ||dkdS) zÎCompute a lower bound for the adjusted exponent of self.log10(). In other words, find r such that self.log10() >= 10**r. Assumes that self is finite and positive and that self != 1. r.r,r&rŸéçZ231r r…)rPrgr9rarjrVrW)r/r<r,r-r¯rDrEr(r(r*rÏszDecimal._log10_exp_boundc CsH|dkrtƒ}|j|d}|r"|S|s*tS|jƒdkr:tS|jdkrP|jtdƒS|jddkr˜|jdd…dt |jƒdkr˜t |jt |jƒdƒ}nŠt|ƒ}|j |j}}|j}||jƒd}x>t|||ƒ}|d d t tt|ƒƒƒ|drøP|d7}qÆWtt |dkƒtt|ƒƒ|ƒ}|jƒ}|jtƒ} |j|ƒ}| |_|S)z&Returns the base 10 logarithm of self.N)r0r.zlog10 of a negative valuer&rÆrRr,r{rŸr\)rr‹rGr†rHr8rer r9rgrrPrjrVrWrArÚ_dlog10rarir7r*r+rr¹r@) r/r0r<r,r-r¯r r'rÊr@r(r(r*Úlog10ís: .$ z Decimal.log10cCsV|j|d}|r|S|dkr"tƒ}|jƒr.tS|s@|jtddƒSt|jƒƒ}|j|ƒS)aM Returns the exponent of the magnitude of self's MSD. The result is the integer which is the exponent of the magnitude of the most significant digit of self (as though it were truncated to a single digit while maintaining the value of that digit and without limiting the resulting exponent). )r0Nzlogb(0)r.) r‹rr†rHrer rrr¹)r/r0r<r(r(r*Úlogb s zDecimal.logbcCs8|jdks|jdkrdSx|jD]}|dkr dSq WdS)z×Return True if self is a logical operand. For being logical, it must be a finite number with a sign of 0, an exponent of 0, and a coefficient whose digits must all be either 0 or 1. r&FZ01T)r8rPr9)r/Údigr(r(r*Ú _islogical> szDecimal._islogicalcCs€|jt|ƒ}|dkr$d||}n|dkr<||jd…}|jt|ƒ}|dkr`d||}n|dkrx||jd…}||fS)Nr&rR)rArg)r/r0ÚopaÚopbZdifr(r(r*Ú _fill_logicalL szDecimal._fill_logicalcCs~|dkrtƒ}t|dd}|jƒs.|jƒr8|jtƒS|j||j|jƒ\}}djdd„t||ƒDƒƒ}t d|j dƒpxddƒS) z;Applies an 'and' operation between self and other's digits.NT)rœrTcSs$g|]\}}tt|ƒt|ƒ@ƒ‘qSr()rarV)Ú.0rÚbr(r(r*ú g sz'Decimal.logical_and..r&rR)rrrPrer rSr9roÚzipr7rh)r/rˆr0rQrRrƒr(r(r*Úlogical_andY s zDecimal.logical_andcCs(|dkrtƒ}|jtdd|jdƒ|ƒS)zInvert all its digits.Nr&rÆ)rÚlogical_xorr7rA)r/r0r(r(r*Úlogical_invertj szDecimal.logical_invertcCs~|dkrtƒ}t|dd}|jƒs.|jƒr8|jtƒS|j||j|jƒ\}}djdd„t||ƒDƒƒ}t d|j dƒpxddƒS) z:Applies an 'or' operation between self and other's digits.NT)rœrTcSs$g|]\}}tt|ƒt|ƒBƒ‘qSr()rarV)rTrrUr(r(r*rV sz&Decimal.logical_or..r&rR)rrrPrer rSr9rorWr7rh)r/rˆr0rQrRrƒr(r(r*Ú logical_orq s zDecimal.logical_orcCs~|dkrtƒ}t|dd}|jƒs.|jƒr8|jtƒS|j||j|jƒ\}}djdd„t||ƒDƒƒ}t d|j dƒpxddƒS) z;Applies an 'xor' operation between self and other's digits.NT)rœrTcSs$g|]\}}tt|ƒt|ƒAƒ‘qSr()rarV)rTrrUr(r(r*rV sz'Decimal.logical_xor..r&rR)rrrPrer rSr9rorWr7rh)r/rˆr0rQrRrƒr(r(r*rY‚ s zDecimal.logical_xorcCs¾t|dd}|dkrtƒ}|js&|jr~|jƒ}|jƒ}|s>|r~|dkrX|dkrX|j|ƒS|dkrr|dkrr|j|ƒS|j||ƒS|jƒj|jƒƒ}|dkr¢|j|ƒ}|dkr°|}n|}|j|ƒS)z8Compares the values numerically with their sign ignored.T)rœNr.r&r…) rrrQr„r¹r‹r·r“r0)r/rˆr0r1r2r-r<r(r(r*Úmax_mag“ s& zDecimal.max_magcCs¾t|dd}|dkrtƒ}|js&|jr~|jƒ}|jƒ}|s>|r~|dkrX|dkrX|j|ƒS|dkrr|dkrr|j|ƒS|j||ƒS|jƒj|jƒƒ}|dkr¢|j|ƒ}|dkr°|}n|}|j|ƒS)z8Compares the values numerically with their sign ignored.T)rœNr.r&r…) rrrQr„r¹r‹r·r“r0)r/rˆr0r1r2r-r<r(r(r*Úmin_mag± s& zDecimal.min_magcCs |dkrtƒ}|j|d}|r"|S|jƒdkr2tS|jƒdkrTtdd|j|jƒƒS|jƒ}|jt ƒ|j ƒ|j|ƒ}||kr„|S|jtdd|j ƒdƒ|ƒS)z=Returns the largest representable number smaller than itself.N)r0r.r&r?rÆr…)rr‹r†rGr7rArærIr+rÚ_ignore_all_flagsr¹rÄrÈ)r/r0r<Únew_selfr(r(r*Ú next_minusÏ s" zDecimal.next_minuscCs |dkrtƒ}|j|d}|r"|S|jƒdkr2tS|jƒdkrTtdd|j|jƒƒS|jƒ}|jt ƒ|j ƒ|j|ƒ}||kr„|S|jtdd|j ƒdƒ|ƒS)z=Returns the smallest representable number larger than itself.N)r0r.r?r&rÆr…)rr‹r†rHr7rArærIr+rr^r¹rÃrÈ)r/r0r<r_r(r(r*Ú next_plusæ s" zDecimal.next_pluscCsÞt|dd}|dkrtƒ}|j||ƒ}|r.|S|j|ƒ}|dkrJ|j|ƒS|dkr^|j|ƒ}n |j|ƒ}|jƒr–|jt d|j ƒ|jtƒ|jtƒnD|j ƒ|jkrÚ|jtƒ|jtƒ|jtƒ|jtƒ|sÚ|jtƒ|S)a‹Returns the number closest to self, in the direction towards other. The result is the closest representable number to self (excluding self) that is in the direction towards other, unless both have the same value. If the two operands are numerically equal, then the result is a copy of self with the sign set to be the same as the sign of other. T)rœNr&r.z Infinite result from next_towardr…)rrr‹r“r:rar`r†rerr8rrrr$rr r)r/rˆr0r<Z comparisonr(r(r*Únext_towardý s4 zDecimal.next_towardcCsˆ|jƒrdS|jƒrdS|jƒ}|dkr,dS|dkr8dS|jƒrN|jrJdSdS|dkr\tƒ}|j|d rv|jrrd SdS|jr€dSd SdS)aReturns an indication of the class of self. The class is one of the following strings: sNaN NaN -Infinity -Normal -Subnormal -Zero +Zero +Subnormal +Normal +Infinity r‡r¬r.z +Infinityz -Infinityz-Zeroz+ZeroN)r0z -Subnormalz +Subnormalz-Normalz+Normalr…)rŒrr†rBr8rrA)r/r0Úinfr(r(r*Únumber_class+s,zDecimal.number_classcCstdƒS)z'Just returns 10, as this is Decimal, :)rŸ)r)r/r(r(r*ÚradixUsz Decimal.radixcCsø|dkrtƒ}t|dd}|j||ƒ}|r.|S|jdkrB|jtƒS|jt|ƒko^|jknsn|jtƒS|jƒr~t |ƒSt|ƒ}|j }|jt|ƒ}|dkr°d||}n|dkrÆ||d…}||d…|d|…}t|j |jdƒpðd|jƒS)z5Returns a rotated copy of self, value-of-other times.NT)rœr&rR)rrr‹rPrer rArVr†rr9rgr7r8rh)r/rˆr0r<ÚtorotÚrotdigÚtopadZrotatedr(r(r*ÚrotateYs, " zDecimal.rotatecCsÀ|dkrtƒ}t|dd}|j||ƒ}|r.|S|jdkrB|jtƒSd|j|j}d|j|j}|t|ƒkox|knsˆ|jtƒS|j ƒr˜t |ƒSt|j|j |jt|ƒƒ}|j|ƒ}|S)z>Returns self operand after adding the second value to its exp.NT)rœr&r,r )rrr‹rPrer rBrArVr†rr7r8r9r¹)r/rˆr0r<ZliminfZlimsuprr(r(r*Úscalebzs" zDecimal.scalebcCs|dkrtƒ}t|dd}|j||ƒ}|r.|S|jdkrB|jtƒS|jt|ƒko^|jknsn|jtƒS|jƒr~t |ƒSt|ƒ}|j }|jt|ƒ}|dkr°d||}n|dkrÆ||d…}|dkrÜ|d|…}n|d|}||jd…}t|j |jdƒpd|jƒS)z5Returns a shifted copy of self, value-of-other times.NT)rœr&rR)rrr‹rPrer rArVr†rr9rgr7r8rh)r/rˆr0r<rfrgrhZshiftedr(r(r*rË“s2 " z Decimal.shiftcCs|jt|ƒffS)N)Ú __class__ra)r/r(r(r*Ú __reduce__ºszDecimal.__reduce__cCst|ƒtkr|S|jt|ƒƒS)N)Útyperrkra)r/r(r(r*Ú__copy__½szDecimal.__copy__cCst|ƒtkr|S|jt|ƒƒS)N)rmrrkra)r/Úmemor(r(r*Ú__deepcopy__ÂszDecimal.__deepcopy__cCsN|dkrtƒ}t||d}|jrXt|j|ƒ}t|jƒƒ}|ddkrL|d7}t|||ƒS|ddkrvddg|j|d<|ddkr˜t |j|j |jdƒ}|j}|d}|dk r|dd krÎ|j |d |ƒ}nF|ddkrê|j||ƒ}n*|ddkrt|j ƒ|kr|j ||ƒ}|rB|jd krB|ddkrB|jd |ƒ}|jt|j ƒ} |dd kr‚|r||dk r|d |} nd } nB|ddkr–| } n.|ddkrÄ|jd krÀ| dkrÀ| } nd } | d kräd}d| |j }nP| t|j ƒkr|j d| t|j ƒ}d}n"|j d| …p$d}|j | d…}| | } t|j||| |ƒS)a|Format a Decimal instance according to the given specifier. The specifier should be a standard format specifier, with the form described in PEP 3101. Formatting types 'e', 'E', 'f', 'F', 'g', 'G', 'n' and '%' are supported. If the formatting type is omitted it defaults to 'g' or 'G', depending on the value of context.capitals. N)Ú_localeconvrmú%ÚgÚGr,Ú precisionÚeEr.zfF%ZgGr&rrRrTiúÿÿÿ)rÚ_parse_format_specifierrQÚ_format_signr8rar·Ú _format_alignr±r7r9rPr@r(r¿rgÚ_format_number)r/Z specifierr0rqÚspecr>Úbodyr@rur³r´rwrxrWr(r(r*Ú __format__ÉsZ " zDecimal.__format__)rPr9r8rQ)rRN)NN)N)N)N)N)N)N)FN)N)N)N)TN)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)NN)N)N)NN)N)NN)NN)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)N)NN)‚r2r3r4r5Ú __slots__r_Úclassmethodrrr„r†r‹rŽrr“r—r˜r™ršr›ržr§r¨r~r«rµr¶rºr»r½rÃÚ__radd__rÄrÅrÇÚ__rmul__rÎrÒrÓrÔrÕrÖr×rÙrÚrÛrÞrßÚ __trunc__ràÚpropertyrárârär:r¹rërìrírïrñròrórôÚdictrçrör÷rørúrrrr r#rõr&r¿r(r)rýÚto_integralr/r¾rªrûrürr4r5r0r9r·r¸r:rWr=r>r%r¡r?rr@rŒrArBrFrJrrMrNrPrSrXrZr[rYr\r]r`rarbrdrerirjrËrlrnrpr}r(r(r(r*r+s ( !@ 2 4 V7;!$K f >,UnY="c*" IK23 .*!'FcCs&tjtƒ}||_||_||_||_|S)z½Create a decimal instance directly, without any validation, normalization (e.g. removal of leading zeros) or argument conversion. This function is for *internal use only*. )r^r_rr8r9rPrQ)r>ZcoefficientrZspecialr/r(r(r*r7s r7c@s(eZdZdZdd„Zdd„Zdd„ZdS) rOzContext manager class to support localcontext(). Sets a copy of the supplied context in __enter__() and restores the previous decimal context in __exit__() cCs|jƒ|_dS)N)rIÚnew_context)r/r†r(r(r*Ú__init__9sz_ContextManager.__init__cCstƒ|_t|jƒ|jS)N)rÚ saved_contextrr†)r/r(r(r*Ú __enter__;s z_ContextManager.__enter__cCst|jƒdS)N)rrˆ)r/ÚtÚvÚtbr(r(r*Ú__exit__?sz_ContextManager.__exit__N)r2r3r4r5r‡r‰rr(r(r(r*rO3srOc @s¬eZdZdZd¦dd„Zdd„Zdd„Zd d „Zdd„Zd d„Z dd„Z dd„Zdd„Zdd„Z dd„ZeZd§dd„Zdd„Zdd„Zdd „ZdZd!d"„Zd#d$„Zd%d&„Zd¨d(d)„Zd*d+„Zd,d-„Zd.d/„Zd0d1„Zd2d3„Zd4d5„Zd6d7„Zd8d9„Z d:d;„Z!dd?„Z#d@dA„Z$dBdC„Z%dDdE„Z&dFdG„Z'dHdI„Z(dJdK„Z)dLdM„Z*dNdO„Z+dPdQ„Z,dRdS„Z-dTdU„Z.dVdW„Z/dXdY„Z0dZd[„Z1d\d]„Z2d^d_„Z3d`da„Z4dbdc„Z5ddde„Z6dfdg„Z7dhdi„Z8djdk„Z9dldm„Z:dndo„Z;dpdq„Zdvdw„Z?dxdy„Z@dzd{„ZAd|d}„ZBd~d„ZCd€d„ZDd‚dƒ„ZEd„d…„ZFd†d‡„ZGd©dˆd‰„ZHdŠd‹„ZIdŒd„ZJdŽd„ZKdd‘„ZLd’d“„ZMd”d•„ZNd–d—„ZOd˜d™„ZPdšd›„ZQdœd„ZRdždŸ„ZSd d¡„ZTd¢d£„ZUd¤d¥„ZVeVZWdS)ªraßContains the context for a Decimal instance. Contains: prec - precision (for use in rounding, division, square roots..) rounding - rounding type (how you round) traps - If traps[exception] = 1, then the exception is raised when it is caused. Otherwise, a value is substituted in. flags - When an exception is caused, flags[exception] is set. (Whether or not the trap_enabler is set) Should be reset by user of Decimal instance. Emin - Minimum exponent Emax - Maximum exponent capitals - If 1, 1*10^1 is printed as 1E+1. If 0, printed as 1e1 clamp - If 1, change exponents if too high (Default 0) Nc s>yt} Wntk rYnX|dk r*|n| j|_|dk r>|n| j|_|dk rR|n| j|_|dk rf|n| j|_|dk rz|n| j|_|dk rŽ|n| j|_| dkr¦g|_n| |_ˆdkrÂ| j j ƒ|_ n.tˆtƒsêt‡fdd„t ˆDƒƒ|_ nˆ|_ ˆdkr tjt dƒ|_n0tˆtƒs4t‡fdd„t ˆDƒƒ|_nˆ|_dS)Nc3s|]}|t|ˆkƒfVqdS)N)rV)rTrÝ)rr(r*ú nsz#Context.__init__..r&c3s|]}|t|ˆkƒfVqdS)N)rV)rTrÝ)rr(r*rŽus)rÚ NameErrorrAr@r$rBr±råÚ_ignored_flagsrrIr`r„rÚfromkeysr)r/rAr@r$rBr±rårrrZdcr()rrr*r‡Us. zContext.__init__cCs”t|tƒstd|ƒ‚|dkr<||kr†td||||fƒ‚nJ|dkrb||kr†td||||fƒ‚n$||ksr||kr†td||||fƒ‚tj|||ƒS)Nz%s must be an integerz-infz%s must be in [%s, %d]. got: %srcz%s must be in [%d, %s]. got: %sz%s must be in [%d, %d]. got %s)r`rVrsrmr^Ú__setattr__)r/ÚnameruZvminZvmaxr(r(r*Ú_set_integer_checkys zContext._set_integer_checkcCsht|tƒstd|ƒ‚x |D]}|tkrtd|ƒ‚qWx tD]}||kr>td|ƒ‚q>Wtj|||ƒS)Nz%s must be a signal dictz%s is not a valid signal dict)r`r„rsrÚKeyErrorr^r’)r/r“rÚkeyr(r(r*Ú_set_signal_dict‡s zContext._set_signal_dictcCsä|dkr|j||ddƒS|dkr0|j||ddƒS|dkrH|j||ddƒS|dkr`|j||ddƒS|d krx|j||ddƒS|d kr¢|tkr”td|ƒ‚tj|||ƒS|dks²|d kr¾|j||ƒS|dkrÔtj|||ƒStd|ƒ‚dS)NrAr.rcr$z-infr&rBr±rår@z%s: invalid rounding moderrrz.'decimal.Context' object has no attribute '%s')r”Ú_rounding_modesrsr^r’r—rM)r/r“rur(r(r*r’’s(zContext.__setattr__cCstd|ƒ‚dS)Nz%s cannot be deleted)rM)r/r“r(r(r*Ú__delattr__«szContext.__delattr__c CsNdd„|jjƒDƒ}dd„|jjƒDƒ}|j|j|j|j|j|j|j ||ffS)NcSsg|]\}}|r|‘qSr(r()rTÚsigr‹r(r(r*rV°sz&Context.__reduce__..cSsg|]\}}|r|‘qSr(r()rTršr‹r(r(r*rV±s) rÚitemsrrkrAr@r$rBr±rå)r/rrr(r(r*rl¯s zContext.__reduce__cCs|g}|jdt|ƒƒdd„|jjƒDƒ}|jddj|ƒdƒdd„|jjƒDƒ}|jddj|ƒdƒdj|ƒd S) zShow the current context.zrContext(prec=%(prec)d, rounding=%(rounding)s, Emin=%(Emin)d, Emax=%(Emax)d, capitals=%(capitals)d, clamp=%(clamp)dcSsg|]\}}|r|j‘qSr()r2)rTr€r‹r(r(r*rV½sz$Context.__repr__..zflags=[z, ú]cSsg|]\}}|r|j‘qSr()r2)rTrŠr‹r(r(r*rV¿sztraps=[ú))rnÚvarsrr›ror)r/rÝÚnamesr(r(r*r«¶szContext.__repr__cCsx|jD]}d|j|<qWdS)zReset all flags to zeror&N)r)r/Úflagr(r(r*rJÃszContext.clear_flagscCsx|jD]}d|j|<qWdS)zReset all traps to zeror&N)r)r/r r(r(r*Úclear_trapsÈszContext.clear_trapsc Cs.t|j|j|j|j|j|j|j|j|j ƒ }|S)z!Returns a shallow copy from self.) rrAr@r$rBr±rårrr)r/Úncr(r(r*r*ÍszContext._shallow_copyc Cs6t|j|j|j|j|j|j|jjƒ|j jƒ|j ƒ }|S)zReturns a deep copy from self.)rrAr@r$rBr±rårrIrr)r/r¢r(r(r*rIÔs zContext.copycGsZtj||ƒ}||jkr(|ƒj|f|žŽSd|j|<|j|sN|ƒj|f|žŽS||ƒ‚dS)a#Handles an error If the flag is in _ignored_flags, returns the default response. Otherwise, it sets the flag, then, if the corresponding trap_enabler is set, it reraises the exception. Otherwise, it returns the default value after setting the flag. r.N)Ú_condition_mapÚgetrr1rr)r/Z conditionZexplanationr)Úerrorr(r(r*reÝs zContext._raise_errorcCs |jtŽS)z$Ignore all flags, if they are raised)Ú _ignore_flagsr)r/r(r(r*r^ószContext._ignore_all_flagscGs|jt|ƒ|_t|ƒS)z$Ignore the flags, if they are raised)rrk)r/rr(r(r*r¦÷szContext._ignore_flagscGs<|rt|dttfƒr|d}x|D]}|jj|ƒq$WdS)z+Stop ignoring the flags, if they are raisedr&N)r`rlrkrÚremove)r/rr r(r(r*Ú _regard_flagsþs zContext._regard_flagscCst|j|jdƒS)z!Returns Etiny (= Emin - prec + 1)r.)rVr$rA)r/r(r(r*rÈsz Context.EtinycCst|j|jdƒS)z,Returns maximum exponent (= Emax - prec + 1)r.)rVrBrA)r/r(r(r*ræszContext.EtopcCs|j}||_|S)aÓSets the rounding type. Sets the rounding type, and returns the current (previous) rounding type. Often used like: context = context.copy() # so you don't change the calling context # if an error occurs in the middle. rounding = context._set_rounding(ROUND_UP) val = self.__sub__(other, context=context) context._set_rounding(rounding) This will make it round up for that operation. )r@)r/rmr@r(r(r*r+szContext._set_roundingrRcCsjt|tƒr*||jƒksd|kr*|jtdƒSt||d}|jƒr`t|jƒ|j |j kr`|jtdƒS|j|ƒS)z›Creates a new Decimal instance but using self as context. This method implements the to-number operation of the IBM Decimal specification.rSzAtrailing or leading whitespace and underscores are not permitted.)r0zdiagnostic info too long in NaN)r`rarcrerrr„rgr9rArår¹)r/rDrr(r(r*Úcreate_decimal#szContext.create_decimalcCstj|ƒ}|j|ƒS)aÏCreates a new Decimal instance from a float but rounding using self as the context. >>> context = Context(prec=5, rounding=ROUND_DOWN) >>> context.create_decimal_from_float(3.1415926535897932) Decimal('3.1415') >>> context = Context(prec=5, traps=[Inexact]) >>> context.create_decimal_from_float(3.1415926535897932) Traceback (most recent call last): ... decimal.Inexact: None )rrrr¹)r/r€rr(r(r*Úcreate_decimal_from_float4s z!Context.create_decimal_from_floatcCst|dd}|j|dS)a[Returns the absolute value of the operand. If the operand is negative, the result is the same as using the minus operation on the operand. Otherwise, the result is the same as using the plus operation on the operand. >>> ExtendedContext.abs(Decimal('2.1')) Decimal('2.1') >>> ExtendedContext.abs(Decimal('-100')) Decimal('100') >>> ExtendedContext.abs(Decimal('101.5')) Decimal('101.5') >>> ExtendedContext.abs(Decimal('-101.5')) Decimal('101.5') >>> ExtendedContext.abs(-1) Decimal('1') T)rœ)r0)rr½)r/rr(r(r*riFszContext.abscCs8t|dd}|j||d}|tkr0td|ƒ‚n|SdS)a«Return the sum of the two operands. >>> ExtendedContext.add(Decimal('12'), Decimal('7.00')) Decimal('19.00') >>> ExtendedContext.add(Decimal('1E+2'), Decimal('1.01E+4')) Decimal('1.02E+4') >>> ExtendedContext.add(1, Decimal(2)) Decimal('3') >>> ExtendedContext.add(Decimal(8), 5) Decimal('13') >>> ExtendedContext.add(5, 5) Decimal('10') T)rœ)r0zUnable to convert %s to DecimalN)rrÃr–rs)r/rrUrÑr(r(r*Úadd[s zContext.addcCst|j|ƒƒS)N)rar¹)r/rr(r(r*Ú_applypszContext._applycCst|tƒstdƒ‚|jƒS)zûReturns the same Decimal object. As we do not have different encodings for the same number, the received object already is in its canonical form. >>> ExtendedContext.canonical(Decimal('2.50')) Decimal('2.50') z,canonical requires a Decimal as an argument.)r`rrsr4)r/rr(r(r*r4ss zContext.canonicalcCst|dd}|j||dS)a…Compares values numerically. If the signs of the operands differ, a value representing each operand ('-1' if the operand is less than zero, '0' if the operand is zero or negative zero, or '1' if the operand is greater than zero) is used in place of that operand for the comparison instead of the actual operand. The comparison is then effected by subtracting the second operand from the first and then returning a value according to the result of the subtraction: '-1' if the result is less than zero, '0' if the result is zero or negative zero, or '1' if the result is greater than zero. >>> ExtendedContext.compare(Decimal('2.1'), Decimal('3')) Decimal('-1') >>> ExtendedContext.compare(Decimal('2.1'), Decimal('2.1')) Decimal('0') >>> ExtendedContext.compare(Decimal('2.1'), Decimal('2.10')) Decimal('0') >>> ExtendedContext.compare(Decimal('3'), Decimal('2.1')) Decimal('1') >>> ExtendedContext.compare(Decimal('2.1'), Decimal('-3')) Decimal('1') >>> ExtendedContext.compare(Decimal('-3'), Decimal('2.1')) Decimal('-1') >>> ExtendedContext.compare(1, 2) Decimal('-1') >>> ExtendedContext.compare(Decimal(1), 2) Decimal('-1') >>> ExtendedContext.compare(1, Decimal(2)) Decimal('-1') T)rœ)r0)rrž)r/rrUr(r(r*rž€s!zContext.comparecCst|dd}|j||dS)aCompares the values of the two operands numerically. It's pretty much like compare(), but all NaNs signal, with signaling NaNs taking precedence over quiet NaNs. >>> c = ExtendedContext >>> c.compare_signal(Decimal('2.1'), Decimal('3')) Decimal('-1') >>> c.compare_signal(Decimal('2.1'), Decimal('2.1')) Decimal('0') >>> c.flags[InvalidOperation] = 0 >>> print(c.flags[InvalidOperation]) 0 >>> c.compare_signal(Decimal('NaN'), Decimal('2.1')) Decimal('NaN') >>> print(c.flags[InvalidOperation]) 1 >>> c.flags[InvalidOperation] = 0 >>> print(c.flags[InvalidOperation]) 0 >>> c.compare_signal(Decimal('sNaN'), Decimal('2.1')) Decimal('NaN') >>> print(c.flags[InvalidOperation]) 1 >>> c.compare_signal(-1, 2) Decimal('-1') >>> c.compare_signal(Decimal(-1), 2) Decimal('-1') >>> c.compare_signal(-1, Decimal(2)) Decimal('-1') T)rœ)r0)rr5)r/rrUr(r(r*r5¤s zContext.compare_signalcCst|dd}|j|ƒS)a+Compares two operands using their abstract representation. This is not like the standard compare, which use their numerical value. Note that a total ordering is defined for all possible abstract representations. >>> ExtendedContext.compare_total(Decimal('12.73'), Decimal('127.9')) Decimal('-1') >>> ExtendedContext.compare_total(Decimal('-127'), Decimal('12')) Decimal('-1') >>> ExtendedContext.compare_total(Decimal('12.30'), Decimal('12.3')) Decimal('-1') >>> ExtendedContext.compare_total(Decimal('12.30'), Decimal('12.30')) Decimal('0') >>> ExtendedContext.compare_total(Decimal('12.3'), Decimal('12.300')) Decimal('1') >>> ExtendedContext.compare_total(Decimal('12.3'), Decimal('NaN')) Decimal('-1') >>> ExtendedContext.compare_total(1, 2) Decimal('-1') >>> ExtendedContext.compare_total(Decimal(1), 2) Decimal('-1') >>> ExtendedContext.compare_total(1, Decimal(2)) Decimal('-1') T)rœ)rr0)r/rrUr(r(r*r0ÇszContext.compare_totalcCst|dd}|j|ƒS)z£Compares two operands using their abstract representation ignoring sign. Like compare_total, but with operand's sign ignored and assumed to be 0. T)rœ)rr9)r/rrUr(r(r*r9äszContext.compare_total_magcCst|dd}|jƒS)aReturns a copy of the operand with the sign set to 0. >>> ExtendedContext.copy_abs(Decimal('2.1')) Decimal('2.1') >>> ExtendedContext.copy_abs(Decimal('-100')) Decimal('100') >>> ExtendedContext.copy_abs(-1) Decimal('1') T)rœ)rr·)r/rr(r(r*r·ìs zContext.copy_abscCst|dd}t|ƒS)aReturns a copy of the decimal object. >>> ExtendedContext.copy_decimal(Decimal('2.1')) Decimal('2.1') >>> ExtendedContext.copy_decimal(Decimal('-1.00')) Decimal('-1.00') >>> ExtendedContext.copy_decimal(1) Decimal('1') T)rœ)rr)r/rr(r(r*Úcopy_decimalùs zContext.copy_decimalcCst|dd}|jƒS)a(Returns a copy of the operand with the sign inverted. >>> ExtendedContext.copy_negate(Decimal('101.5')) Decimal('-101.5') >>> ExtendedContext.copy_negate(Decimal('-101.5')) Decimal('101.5') >>> ExtendedContext.copy_negate(1) Decimal('-1') T)rœ)rr¸)r/rr(r(r*r¸s zContext.copy_negatecCst|dd}|j|ƒS)aCopies the second operand's sign to the first one. In detail, it returns a copy of the first operand with the sign equal to the sign of the second operand. >>> ExtendedContext.copy_sign(Decimal( '1.50'), Decimal('7.33')) Decimal('1.50') >>> ExtendedContext.copy_sign(Decimal('-1.50'), Decimal('7.33')) Decimal('1.50') >>> ExtendedContext.copy_sign(Decimal( '1.50'), Decimal('-7.33')) Decimal('-1.50') >>> ExtendedContext.copy_sign(Decimal('-1.50'), Decimal('-7.33')) Decimal('-1.50') >>> ExtendedContext.copy_sign(1, -2) Decimal('-1') >>> ExtendedContext.copy_sign(Decimal(1), -2) Decimal('-1') >>> ExtendedContext.copy_sign(1, Decimal(-2)) Decimal('-1') T)rœ)rr:)r/rrUr(r(r*r:szContext.copy_signcCs8t|dd}|j||d}|tkr0td|ƒ‚n|SdS)aˆDecimal division in a specified context. >>> ExtendedContext.divide(Decimal('1'), Decimal('3')) Decimal('0.333333333') >>> ExtendedContext.divide(Decimal('2'), Decimal('3')) Decimal('0.666666667') >>> ExtendedContext.divide(Decimal('5'), Decimal('2')) Decimal('2.5') >>> ExtendedContext.divide(Decimal('1'), Decimal('10')) Decimal('0.1') >>> ExtendedContext.divide(Decimal('12'), Decimal('12')) Decimal('1') >>> ExtendedContext.divide(Decimal('8.00'), Decimal('2')) Decimal('4.00') >>> ExtendedContext.divide(Decimal('2.400'), Decimal('2.0')) Decimal('1.20') >>> ExtendedContext.divide(Decimal('1000'), Decimal('100')) Decimal('10') >>> ExtendedContext.divide(Decimal('1000'), Decimal('1')) Decimal('1000') >>> ExtendedContext.divide(Decimal('2.40E+6'), Decimal('2')) Decimal('1.20E+6') >>> ExtendedContext.divide(5, 5) Decimal('1') >>> ExtendedContext.divide(Decimal(5), 5) Decimal('1') >>> ExtendedContext.divide(5, Decimal(5)) Decimal('1') T)rœ)r0zUnable to convert %s to DecimalN)rrÎr–rs)r/rrUrÑr(r(r*Údivide+s zContext.dividecCs8t|dd}|j||d}|tkr0td|ƒ‚n|SdS)a/Divides two numbers and returns the integer part of the result. >>> ExtendedContext.divide_int(Decimal('2'), Decimal('3')) Decimal('0') >>> ExtendedContext.divide_int(Decimal('10'), Decimal('3')) Decimal('3') >>> ExtendedContext.divide_int(Decimal('1'), Decimal('0.3')) Decimal('3') >>> ExtendedContext.divide_int(10, 3) Decimal('3') >>> ExtendedContext.divide_int(Decimal(10), 3) Decimal('3') >>> ExtendedContext.divide_int(10, Decimal(3)) Decimal('3') T)rœ)r0zUnable to convert %s to DecimalN)rrÚr–rs)r/rrUrÑr(r(r*Ú divide_intPs zContext.divide_intcCs8t|dd}|j||d}|tkr0td|ƒ‚n|SdS)aÝReturn (a // b, a % b). >>> ExtendedContext.divmod(Decimal(8), Decimal(3)) (Decimal('2'), Decimal('2')) >>> ExtendedContext.divmod(Decimal(8), Decimal(4)) (Decimal('2'), Decimal('0')) >>> ExtendedContext.divmod(8, 4) (Decimal('2'), Decimal('0')) >>> ExtendedContext.divmod(Decimal(8), 4) (Decimal('2'), Decimal('0')) >>> ExtendedContext.divmod(8, Decimal(4)) (Decimal('2'), Decimal('0')) T)rœ)r0zUnable to convert %s to DecimalN)rrÔr–rs)r/rrUrÑr(r(r*rÉgs zContext.divmodcCst|dd}|j|dS)a#Returns e ** a. >>> c = ExtendedContext.copy() >>> c.Emin = -999 >>> c.Emax = 999 >>> c.exp(Decimal('-Infinity')) Decimal('0') >>> c.exp(Decimal('-1')) Decimal('0.367879441') >>> c.exp(Decimal('0')) Decimal('1') >>> c.exp(Decimal('1')) Decimal('2.71828183') >>> c.exp(Decimal('0.693147181')) Decimal('2.00000000') >>> c.exp(Decimal('+Infinity')) Decimal('Infinity') >>> c.exp(10) Decimal('22026.4658') T)rœ)r0)rrW)r/rr(r(r*rW|szContext.expcCst|dd}|j|||dS)aReturns a multiplied by b, plus c. The first two operands are multiplied together, using multiply, the third operand is then added to the result of that multiplication, using add, all with only one final rounding. >>> ExtendedContext.fma(Decimal('3'), Decimal('5'), Decimal('7')) Decimal('22') >>> ExtendedContext.fma(Decimal('3'), Decimal('-5'), Decimal('7')) Decimal('-8') >>> ExtendedContext.fma(Decimal('888565290'), Decimal('1557.96930'), Decimal('-86087.7578')) Decimal('1.38435736E+12') >>> ExtendedContext.fma(1, 3, 4) Decimal('7') >>> ExtendedContext.fma(1, Decimal(3), 4) Decimal('7') >>> ExtendedContext.fma(1, 3, Decimal(4)) Decimal('7') T)rœ)r0)rrú)r/rrUr-r(r(r*rú”szContext.fmacCst|tƒstdƒ‚|jƒS)aReturn True if the operand is canonical; otherwise return False. Currently, the encoding of a Decimal instance is always canonical, so this method returns True for any Decimal. >>> ExtendedContext.is_canonical(Decimal('2.50')) True z/is_canonical requires a Decimal as an argument.)r`rrsr=)r/rr(r(r*r=«s zContext.is_canonicalcCst|dd}|jƒS)a,Return True if the operand is finite; otherwise return False. A Decimal instance is considered finite if it is neither infinite nor a NaN. >>> ExtendedContext.is_finite(Decimal('2.50')) True >>> ExtendedContext.is_finite(Decimal('-0.3')) True >>> ExtendedContext.is_finite(Decimal('0')) True >>> ExtendedContext.is_finite(Decimal('Inf')) False >>> ExtendedContext.is_finite(Decimal('NaN')) False >>> ExtendedContext.is_finite(1) True T)rœ)rr>)r/rr(r(r*r>¸szContext.is_finitecCst|dd}|jƒS)aUReturn True if the operand is infinite; otherwise return False. >>> ExtendedContext.is_infinite(Decimal('2.50')) False >>> ExtendedContext.is_infinite(Decimal('-Inf')) True >>> ExtendedContext.is_infinite(Decimal('NaN')) False >>> ExtendedContext.is_infinite(1) False T)rœ)rr%)r/rr(r(r*r%ÎszContext.is_infinitecCst|dd}|jƒS)aOReturn True if the operand is a qNaN or sNaN; otherwise return False. >>> ExtendedContext.is_nan(Decimal('2.50')) False >>> ExtendedContext.is_nan(Decimal('NaN')) True >>> ExtendedContext.is_nan(Decimal('-sNaN')) True >>> ExtendedContext.is_nan(1) False T)rœ)rr¡)r/rr(r(r*r¡Ýs zContext.is_nancCst|dd}|j|dS)aïReturn True if the operand is a normal number; otherwise return False. >>> c = ExtendedContext.copy() >>> c.Emin = -999 >>> c.Emax = 999 >>> c.is_normal(Decimal('2.50')) True >>> c.is_normal(Decimal('0.1E-999')) False >>> c.is_normal(Decimal('0.00')) False >>> c.is_normal(Decimal('-Inf')) False >>> c.is_normal(Decimal('NaN')) False >>> c.is_normal(1) True T)rœ)r0)rr?)r/rr(r(r*r?íszContext.is_normalcCst|dd}|jƒS)aHReturn True if the operand is a quiet NaN; otherwise return False. >>> ExtendedContext.is_qnan(Decimal('2.50')) False >>> ExtendedContext.is_qnan(Decimal('NaN')) True >>> ExtendedContext.is_qnan(Decimal('sNaN')) False >>> ExtendedContext.is_qnan(1) False T)rœ)rr)r/rr(r(r*rszContext.is_qnancCst|dd}|jƒS)aReturn True if the operand is negative; otherwise return False. >>> ExtendedContext.is_signed(Decimal('2.50')) False >>> ExtendedContext.is_signed(Decimal('-12')) True >>> ExtendedContext.is_signed(Decimal('-0')) True >>> ExtendedContext.is_signed(8) False >>> ExtendedContext.is_signed(-8) True T)rœ)rr@)r/rr(r(r*r@szContext.is_signedcCst|dd}|jƒS)aTReturn True if the operand is a signaling NaN; otherwise return False. >>> ExtendedContext.is_snan(Decimal('2.50')) False >>> ExtendedContext.is_snan(Decimal('NaN')) False >>> ExtendedContext.is_snan(Decimal('sNaN')) True >>> ExtendedContext.is_snan(1) False T)rœ)rrŒ)r/rr(r(r*rŒ$s zContext.is_snancCst|dd}|j|dS)aôReturn True if the operand is subnormal; otherwise return False. >>> c = ExtendedContext.copy() >>> c.Emin = -999 >>> c.Emax = 999 >>> c.is_subnormal(Decimal('2.50')) False >>> c.is_subnormal(Decimal('0.1E-999')) True >>> c.is_subnormal(Decimal('0.00')) False >>> c.is_subnormal(Decimal('-Inf')) False >>> c.is_subnormal(Decimal('NaN')) False >>> c.is_subnormal(1) False T)rœ)r0)rrA)r/rr(r(r*rA4szContext.is_subnormalcCst|dd}|jƒS)auReturn True if the operand is a zero; otherwise return False. >>> ExtendedContext.is_zero(Decimal('0')) True >>> ExtendedContext.is_zero(Decimal('2.50')) False >>> ExtendedContext.is_zero(Decimal('-0E+2')) True >>> ExtendedContext.is_zero(1) False >>> ExtendedContext.is_zero(0) True T)rœ)rrB)r/rr(r(r*rBJszContext.is_zerocCst|dd}|j|dS)aþReturns the natural (base e) logarithm of the operand. >>> c = ExtendedContext.copy() >>> c.Emin = -999 >>> c.Emax = 999 >>> c.ln(Decimal('0')) Decimal('-Infinity') >>> c.ln(Decimal('1.000')) Decimal('0') >>> c.ln(Decimal('2.71828183')) Decimal('1.00000000') >>> c.ln(Decimal('10')) Decimal('2.30258509') >>> c.ln(Decimal('+Infinity')) Decimal('Infinity') >>> c.ln(1) Decimal('0') T)rœ)r0)rrJ)r/rr(r(r*rJ[sz Context.lncCst|dd}|j|dS)a§Returns the base 10 logarithm of the operand. >>> c = ExtendedContext.copy() >>> c.Emin = -999 >>> c.Emax = 999 >>> c.log10(Decimal('0')) Decimal('-Infinity') >>> c.log10(Decimal('0.001')) Decimal('-3') >>> c.log10(Decimal('1.000')) Decimal('0') >>> c.log10(Decimal('2')) Decimal('0.301029996') >>> c.log10(Decimal('10')) Decimal('1') >>> c.log10(Decimal('70')) Decimal('1.84509804') >>> c.log10(Decimal('+Infinity')) Decimal('Infinity') >>> c.log10(0) Decimal('-Infinity') >>> c.log10(1) Decimal('0') T)rœ)r0)rrM)r/rr(r(r*rMqsz Context.log10cCst|dd}|j|dS)a4 Returns the exponent of the magnitude of the operand's MSD. The result is the integer which is the exponent of the magnitude of the most significant digit of the operand (as though the operand were truncated to a single digit while maintaining the value of that digit and without limiting the resulting exponent). >>> ExtendedContext.logb(Decimal('250')) Decimal('2') >>> ExtendedContext.logb(Decimal('2.50')) Decimal('0') >>> ExtendedContext.logb(Decimal('0.03')) Decimal('-2') >>> ExtendedContext.logb(Decimal('0')) Decimal('-Infinity') >>> ExtendedContext.logb(1) Decimal('0') >>> ExtendedContext.logb(10) Decimal('1') >>> ExtendedContext.logb(100) Decimal('2') T)rœ)r0)rrN)r/rr(r(r*rNszContext.logbcCst|dd}|j||dS)a”Applies the logical operation 'and' between each operand's digits. The operands must be both logical numbers. >>> ExtendedContext.logical_and(Decimal('0'), Decimal('0')) Decimal('0') >>> ExtendedContext.logical_and(Decimal('0'), Decimal('1')) Decimal('0') >>> ExtendedContext.logical_and(Decimal('1'), Decimal('0')) Decimal('0') >>> ExtendedContext.logical_and(Decimal('1'), Decimal('1')) Decimal('1') >>> ExtendedContext.logical_and(Decimal('1100'), Decimal('1010')) Decimal('1000') >>> ExtendedContext.logical_and(Decimal('1111'), Decimal('10')) Decimal('10') >>> ExtendedContext.logical_and(110, 1101) Decimal('100') >>> ExtendedContext.logical_and(Decimal(110), 1101) Decimal('100') >>> ExtendedContext.logical_and(110, Decimal(1101)) Decimal('100') T)rœ)r0)rrX)r/rrUr(r(r*rX§szContext.logical_andcCst|dd}|j|dS)aInvert all the digits in the operand. The operand must be a logical number. >>> ExtendedContext.logical_invert(Decimal('0')) Decimal('111111111') >>> ExtendedContext.logical_invert(Decimal('1')) Decimal('111111110') >>> ExtendedContext.logical_invert(Decimal('111111111')) Decimal('0') >>> ExtendedContext.logical_invert(Decimal('101010101')) Decimal('10101010') >>> ExtendedContext.logical_invert(1101) Decimal('111110010') T)rœ)r0)rrZ)r/rr(r(r*rZÂszContext.logical_invertcCst|dd}|j||dS)aApplies the logical operation 'or' between each operand's digits. The operands must be both logical numbers. >>> ExtendedContext.logical_or(Decimal('0'), Decimal('0')) Decimal('0') >>> ExtendedContext.logical_or(Decimal('0'), Decimal('1')) Decimal('1') >>> ExtendedContext.logical_or(Decimal('1'), Decimal('0')) Decimal('1') >>> ExtendedContext.logical_or(Decimal('1'), Decimal('1')) Decimal('1') >>> ExtendedContext.logical_or(Decimal('1100'), Decimal('1010')) Decimal('1110') >>> ExtendedContext.logical_or(Decimal('1110'), Decimal('10')) Decimal('1110') >>> ExtendedContext.logical_or(110, 1101) Decimal('1111') >>> ExtendedContext.logical_or(Decimal(110), 1101) Decimal('1111') >>> ExtendedContext.logical_or(110, Decimal(1101)) Decimal('1111') T)rœ)r0)rr[)r/rrUr(r(r*r[ÕszContext.logical_orcCst|dd}|j||dS)a˜Applies the logical operation 'xor' between each operand's digits. The operands must be both logical numbers. >>> ExtendedContext.logical_xor(Decimal('0'), Decimal('0')) Decimal('0') >>> ExtendedContext.logical_xor(Decimal('0'), Decimal('1')) Decimal('1') >>> ExtendedContext.logical_xor(Decimal('1'), Decimal('0')) Decimal('1') >>> ExtendedContext.logical_xor(Decimal('1'), Decimal('1')) Decimal('0') >>> ExtendedContext.logical_xor(Decimal('1100'), Decimal('1010')) Decimal('110') >>> ExtendedContext.logical_xor(Decimal('1111'), Decimal('10')) Decimal('1101') >>> ExtendedContext.logical_xor(110, 1101) Decimal('1011') >>> ExtendedContext.logical_xor(Decimal(110), 1101) Decimal('1011') >>> ExtendedContext.logical_xor(110, Decimal(1101)) Decimal('1011') T)rœ)r0)rrY)r/rrUr(r(r*rYðszContext.logical_xorcCst|dd}|j||dS)a³max compares two values numerically and returns the maximum. If either operand is a NaN then the general rules apply. Otherwise, the operands are compared as though by the compare operation. If they are numerically equal then the left-hand operand is chosen as the result. Otherwise the maximum (closer to positive infinity) of the two operands is chosen as the result. >>> ExtendedContext.max(Decimal('3'), Decimal('2')) Decimal('3') >>> ExtendedContext.max(Decimal('-10'), Decimal('3')) Decimal('3') >>> ExtendedContext.max(Decimal('1.0'), Decimal('1')) Decimal('1') >>> ExtendedContext.max(Decimal('7'), Decimal('NaN')) Decimal('7') >>> ExtendedContext.max(1, 2) Decimal('2') >>> ExtendedContext.max(Decimal(1), 2) Decimal('2') >>> ExtendedContext.max(1, Decimal(2)) Decimal('2') T)rœ)r0)rr¾)r/rrUr(r(r*r¾szContext.maxcCst|dd}|j||dS)aÇCompares the values numerically with their sign ignored. >>> ExtendedContext.max_mag(Decimal('7'), Decimal('NaN')) Decimal('7') >>> ExtendedContext.max_mag(Decimal('7'), Decimal('-10')) Decimal('-10') >>> ExtendedContext.max_mag(1, -2) Decimal('-2') >>> ExtendedContext.max_mag(Decimal(1), -2) Decimal('-2') >>> ExtendedContext.max_mag(1, Decimal(-2)) Decimal('-2') T)rœ)r0)rr\)r/rrUr(r(r*r\&szContext.max_magcCst|dd}|j||dS)a¸min compares two values numerically and returns the minimum. If either operand is a NaN then the general rules apply. Otherwise, the operands are compared as though by the compare operation. If they are numerically equal then the left-hand operand is chosen as the result. Otherwise the minimum (closer to negative infinity) of the two operands is chosen as the result. >>> ExtendedContext.min(Decimal('3'), Decimal('2')) Decimal('2') >>> ExtendedContext.min(Decimal('-10'), Decimal('3')) Decimal('-10') >>> ExtendedContext.min(Decimal('1.0'), Decimal('1')) Decimal('1.0') >>> ExtendedContext.min(Decimal('7'), Decimal('NaN')) Decimal('7') >>> ExtendedContext.min(1, 2) Decimal('1') >>> ExtendedContext.min(Decimal(1), 2) Decimal('1') >>> ExtendedContext.min(1, Decimal(29)) Decimal('1') T)rœ)r0)rrª)r/rrUr(r(r*rª7szContext.mincCst|dd}|j||dS)aÄCompares the values numerically with their sign ignored. >>> ExtendedContext.min_mag(Decimal('3'), Decimal('-2')) Decimal('-2') >>> ExtendedContext.min_mag(Decimal('-3'), Decimal('NaN')) Decimal('-3') >>> ExtendedContext.min_mag(1, -2) Decimal('1') >>> ExtendedContext.min_mag(Decimal(1), -2) Decimal('1') >>> ExtendedContext.min_mag(1, Decimal(-2)) Decimal('1') T)rœ)r0)rr])r/rrUr(r(r*r]RszContext.min_magcCst|dd}|j|dS)aÎMinus corresponds to unary prefix minus in Python. The operation is evaluated using the same rules as subtract; the operation minus(a) is calculated as subtract('0', a) where the '0' has the same exponent as the operand. >>> ExtendedContext.minus(Decimal('1.3')) Decimal('-1.3') >>> ExtendedContext.minus(Decimal('-1.3')) Decimal('1.3') >>> ExtendedContext.minus(1) Decimal('-1') T)rœ)r0)rrº)r/rr(r(r*Úminuscsz Context.minuscCs8t|dd}|j||d}|tkr0td|ƒ‚n|SdS)aàmultiply multiplies two operands. If either operand is a special value then the general rules apply. Otherwise, the operands are multiplied together ('long multiplication'), resulting in a number which may be as long as the sum of the lengths of the two operands. >>> ExtendedContext.multiply(Decimal('1.20'), Decimal('3')) Decimal('3.60') >>> ExtendedContext.multiply(Decimal('7'), Decimal('3')) Decimal('21') >>> ExtendedContext.multiply(Decimal('0.9'), Decimal('0.8')) Decimal('0.72') >>> ExtendedContext.multiply(Decimal('0.9'), Decimal('-0')) Decimal('-0.0') >>> ExtendedContext.multiply(Decimal('654321'), Decimal('654321')) Decimal('4.28135971E+11') >>> ExtendedContext.multiply(7, 7) Decimal('49') >>> ExtendedContext.multiply(Decimal(7), 7) Decimal('49') >>> ExtendedContext.multiply(7, Decimal(7)) Decimal('49') T)rœ)r0zUnable to convert %s to DecimalN)rrÇr–rs)r/rrUrÑr(r(r*Úmultiplyts zContext.multiplycCst|dd}|j|dS)a"Returns the largest representable number smaller than a. >>> c = ExtendedContext.copy() >>> c.Emin = -999 >>> c.Emax = 999 >>> ExtendedContext.next_minus(Decimal('1')) Decimal('0.999999999') >>> c.next_minus(Decimal('1E-1007')) Decimal('0E-1007') >>> ExtendedContext.next_minus(Decimal('-1.00000003')) Decimal('-1.00000004') >>> c.next_minus(Decimal('Infinity')) Decimal('9.99999999E+999') >>> c.next_minus(1) Decimal('0.999999999') T)rœ)r0)rr`)r/rr(r(r*r`”szContext.next_minuscCst|dd}|j|dS)aReturns the smallest representable number larger than a. >>> c = ExtendedContext.copy() >>> c.Emin = -999 >>> c.Emax = 999 >>> ExtendedContext.next_plus(Decimal('1')) Decimal('1.00000001') >>> c.next_plus(Decimal('-1E-1007')) Decimal('-0E-1007') >>> ExtendedContext.next_plus(Decimal('-1.00000003')) Decimal('-1.00000002') >>> c.next_plus(Decimal('-Infinity')) Decimal('-9.99999999E+999') >>> c.next_plus(1) Decimal('1.00000001') T)rœ)r0)rra)r/rr(r(r*ra¨szContext.next_pluscCst|dd}|j||dS)a´Returns the number closest to a, in direction towards b. The result is the closest representable number from the first operand (but not the first operand) that is in the direction towards the second operand, unless the operands have the same value. >>> c = ExtendedContext.copy() >>> c.Emin = -999 >>> c.Emax = 999 >>> c.next_toward(Decimal('1'), Decimal('2')) Decimal('1.00000001') >>> c.next_toward(Decimal('-1E-1007'), Decimal('1')) Decimal('-0E-1007') >>> c.next_toward(Decimal('-1.00000003'), Decimal('0')) Decimal('-1.00000002') >>> c.next_toward(Decimal('1'), Decimal('0')) Decimal('0.999999999') >>> c.next_toward(Decimal('1E-1007'), Decimal('-100')) Decimal('0E-1007') >>> c.next_toward(Decimal('-1.00000003'), Decimal('-10')) Decimal('-1.00000004') >>> c.next_toward(Decimal('0.00'), Decimal('-0.0000')) Decimal('-0.00') >>> c.next_toward(0, 1) Decimal('1E-1007') >>> c.next_toward(Decimal(0), 1) Decimal('1E-1007') >>> c.next_toward(0, Decimal(1)) Decimal('1E-1007') T)rœ)r0)rrb)r/rrUr(r(r*rb¼s zContext.next_towardcCst|dd}|j|dS)a³normalize reduces an operand to its simplest form. Essentially a plus operation with all trailing zeros removed from the result. >>> ExtendedContext.normalize(Decimal('2.1')) Decimal('2.1') >>> ExtendedContext.normalize(Decimal('-2.0')) Decimal('-2') >>> ExtendedContext.normalize(Decimal('1.200')) Decimal('1.2') >>> ExtendedContext.normalize(Decimal('-120')) Decimal('-1.2E+2') >>> ExtendedContext.normalize(Decimal('120.00')) Decimal('1.2E+2') >>> ExtendedContext.normalize(Decimal('0.00')) Decimal('0') >>> ExtendedContext.normalize(6) Decimal('6') T)rœ)r0)rr#)r/rr(r(r*r#ßszContext.normalizecCst|dd}|j|dS)aâReturns an indication of the class of the operand. The class is one of the following strings: -sNaN -NaN -Infinity -Normal -Subnormal -Zero +Zero +Subnormal +Normal +Infinity >>> c = ExtendedContext.copy() >>> c.Emin = -999 >>> c.Emax = 999 >>> c.number_class(Decimal('Infinity')) '+Infinity' >>> c.number_class(Decimal('1E-10')) '+Normal' >>> c.number_class(Decimal('2.50')) '+Normal' >>> c.number_class(Decimal('0.1E-999')) '+Subnormal' >>> c.number_class(Decimal('0')) '+Zero' >>> c.number_class(Decimal('-0')) '-Zero' >>> c.number_class(Decimal('-0.1E-999')) '-Subnormal' >>> c.number_class(Decimal('-1E-10')) '-Normal' >>> c.number_class(Decimal('-2.50')) '-Normal' >>> c.number_class(Decimal('-Infinity')) '-Infinity' >>> c.number_class(Decimal('NaN')) 'NaN' >>> c.number_class(Decimal('-NaN')) 'NaN' >>> c.number_class(Decimal('sNaN')) 'sNaN' >>> c.number_class(123) '+Normal' T)rœ)r0)rrd)r/rr(r(r*rd÷s/zContext.number_classcCst|dd}|j|dS)a¿Plus corresponds to unary prefix plus in Python. The operation is evaluated using the same rules as add; the operation plus(a) is calculated as add('0', a) where the '0' has the same exponent as the operand. >>> ExtendedContext.plus(Decimal('1.3')) Decimal('1.3') >>> ExtendedContext.plus(Decimal('-1.3')) Decimal('-1.3') >>> ExtendedContext.plus(-1) Decimal('-1') T)rœ)r0)rr»)r/rr(r(r*Úplus)szContext.pluscCs:t|dd}|j|||d}|tkr2td|ƒ‚n|SdS)aRaises a to the power of b, to modulo if given. With two arguments, compute a**b. If a is negative then b must be integral. The result will be inexact unless b is integral and the result is finite and can be expressed exactly in 'precision' digits. With three arguments, compute (a**b) % modulo. For the three argument form, the following restrictions on the arguments hold: - all three arguments must be integral - b must be nonnegative - at least one of a or b must be nonzero - modulo must be nonzero and have at most 'precision' digits The result of pow(a, b, modulo) is identical to the result that would be obtained by computing (a**b) % modulo with unbounded precision, but is computed more efficiently. It is always exact. >>> c = ExtendedContext.copy() >>> c.Emin = -999 >>> c.Emax = 999 >>> c.power(Decimal('2'), Decimal('3')) Decimal('8') >>> c.power(Decimal('-2'), Decimal('3')) Decimal('-8') >>> c.power(Decimal('2'), Decimal('-3')) Decimal('0.125') >>> c.power(Decimal('1.7'), Decimal('8')) Decimal('69.7575744') >>> c.power(Decimal('10'), Decimal('0.301029996')) Decimal('2.00000000') >>> c.power(Decimal('Infinity'), Decimal('-1')) Decimal('0') >>> c.power(Decimal('Infinity'), Decimal('0')) Decimal('1') >>> c.power(Decimal('Infinity'), Decimal('1')) Decimal('Infinity') >>> c.power(Decimal('-Infinity'), Decimal('-1')) Decimal('-0') >>> c.power(Decimal('-Infinity'), Decimal('0')) Decimal('1') >>> c.power(Decimal('-Infinity'), Decimal('1')) Decimal('-Infinity') >>> c.power(Decimal('-Infinity'), Decimal('2')) Decimal('Infinity') >>> c.power(Decimal('0'), Decimal('0')) Decimal('NaN') >>> c.power(Decimal('3'), Decimal('7'), Decimal('16')) Decimal('11') >>> c.power(Decimal('-3'), Decimal('7'), Decimal('16')) Decimal('-11') >>> c.power(Decimal('-3'), Decimal('8'), Decimal('16')) Decimal('1') >>> c.power(Decimal('3'), Decimal('7'), Decimal('-16')) Decimal('11') >>> c.power(Decimal('23E12345'), Decimal('67E189'), Decimal('123456789')) Decimal('11729830') >>> c.power(Decimal('-0'), Decimal('17'), Decimal('1729')) Decimal('-0') >>> c.power(Decimal('-23'), Decimal('0'), Decimal('65537')) Decimal('1') >>> ExtendedContext.power(7, 7) Decimal('823543') >>> ExtendedContext.power(Decimal(7), 7) Decimal('823543') >>> ExtendedContext.power(7, Decimal(7), 2) Decimal('1') T)rœ)r0zUnable to convert %s to DecimalN)rrr–rs)r/rrUrÿrÑr(r(r*Úpower:s Iz Context.powercCst|dd}|j||dS)a Returns a value equal to 'a' (rounded), having the exponent of 'b'. The coefficient of the result is derived from that of the left-hand operand. It may be rounded using the current rounding setting (if the exponent is being increased), multiplied by a positive power of ten (if the exponent is being decreased), or is unchanged (if the exponent is already equal to that of the right-hand operand). Unlike other operations, if the length of the coefficient after the quantize operation would be greater than precision then an Invalid operation condition is raised. This guarantees that, unless there is an error condition, the exponent of the result of a quantize is always equal to that of the right-hand operand. Also unlike other operations, quantize will never raise Underflow, even if the result is subnormal and inexact. >>> ExtendedContext.quantize(Decimal('2.17'), Decimal('0.001')) Decimal('2.170') >>> ExtendedContext.quantize(Decimal('2.17'), Decimal('0.01')) Decimal('2.17') >>> ExtendedContext.quantize(Decimal('2.17'), Decimal('0.1')) Decimal('2.2') >>> ExtendedContext.quantize(Decimal('2.17'), Decimal('1e+0')) Decimal('2') >>> ExtendedContext.quantize(Decimal('2.17'), Decimal('1e+1')) Decimal('0E+1') >>> ExtendedContext.quantize(Decimal('-Inf'), Decimal('Infinity')) Decimal('-Infinity') >>> ExtendedContext.quantize(Decimal('2'), Decimal('Infinity')) Decimal('NaN') >>> ExtendedContext.quantize(Decimal('-0.1'), Decimal('1')) Decimal('-0') >>> ExtendedContext.quantize(Decimal('-0'), Decimal('1e+5')) Decimal('-0E+5') >>> ExtendedContext.quantize(Decimal('+35236450.6'), Decimal('1e-2')) Decimal('NaN') >>> ExtendedContext.quantize(Decimal('-35236450.6'), Decimal('1e-2')) Decimal('NaN') >>> ExtendedContext.quantize(Decimal('217'), Decimal('1e-1')) Decimal('217.0') >>> ExtendedContext.quantize(Decimal('217'), Decimal('1e-0')) Decimal('217') >>> ExtendedContext.quantize(Decimal('217'), Decimal('1e+1')) Decimal('2.2E+2') >>> ExtendedContext.quantize(Decimal('217'), Decimal('1e+2')) Decimal('2E+2') >>> ExtendedContext.quantize(1, 2) Decimal('1') >>> ExtendedContext.quantize(Decimal(1), 2) Decimal('1') >>> ExtendedContext.quantize(1, Decimal(2)) Decimal('1') T)rœ)r0)rrõ)r/rrUr(r(r*rõŠs7zContext.quantizecCstdƒS)zkJust returns 10, as this is Decimal, :) >>> ExtendedContext.radix() Decimal('10') rŸ)r)r/r(r(r*reÄsz Context.radixcCs8t|dd}|j||d}|tkr0td|ƒ‚n|SdS)aReturns the remainder from integer division. The result is the residue of the dividend after the operation of calculating integer division as described for divide-integer, rounded to precision digits if necessary. The sign of the result, if non-zero, is the same as that of the original dividend. This operation will fail under the same conditions as integer division (that is, if integer division on the same two operands would fail, the remainder cannot be calculated). >>> ExtendedContext.remainder(Decimal('2.1'), Decimal('3')) Decimal('2.1') >>> ExtendedContext.remainder(Decimal('10'), Decimal('3')) Decimal('1') >>> ExtendedContext.remainder(Decimal('-10'), Decimal('3')) Decimal('-1') >>> ExtendedContext.remainder(Decimal('10.2'), Decimal('1')) Decimal('0.2') >>> ExtendedContext.remainder(Decimal('10'), Decimal('0.3')) Decimal('0.1') >>> ExtendedContext.remainder(Decimal('3.6'), Decimal('1.3')) Decimal('1.0') >>> ExtendedContext.remainder(22, 6) Decimal('4') >>> ExtendedContext.remainder(Decimal(22), 6) Decimal('4') >>> ExtendedContext.remainder(22, Decimal(6)) Decimal('4') T)rœ)r0zUnable to convert %s to DecimalN)rrÖr–rs)r/rrUrÑr(r(r*rÌÌs zContext.remaindercCst|dd}|j||dS)aGReturns to be "a - b * n", where n is the integer nearest the exact value of "x / b" (if two integers are equally near then the even one is chosen). If the result is equal to 0 then its sign will be the sign of a. This operation will fail under the same conditions as integer division (that is, if integer division on the same two operands would fail, the remainder cannot be calculated). >>> ExtendedContext.remainder_near(Decimal('2.1'), Decimal('3')) Decimal('-0.9') >>> ExtendedContext.remainder_near(Decimal('10'), Decimal('6')) Decimal('-2') >>> ExtendedContext.remainder_near(Decimal('10'), Decimal('3')) Decimal('1') >>> ExtendedContext.remainder_near(Decimal('-10'), Decimal('3')) Decimal('-1') >>> ExtendedContext.remainder_near(Decimal('10.2'), Decimal('1')) Decimal('0.2') >>> ExtendedContext.remainder_near(Decimal('10'), Decimal('0.3')) Decimal('0.1') >>> ExtendedContext.remainder_near(Decimal('3.6'), Decimal('1.3')) Decimal('-0.3') >>> ExtendedContext.remainder_near(3, 11) Decimal('3') >>> ExtendedContext.remainder_near(Decimal(3), 11) Decimal('3') >>> ExtendedContext.remainder_near(3, Decimal(11)) Decimal('3') T)rœ)r0)rrÙ)r/rrUr(r(r*rÙòszContext.remainder_nearcCst|dd}|j||dS)aNReturns a rotated copy of a, b times. The coefficient of the result is a rotated copy of the digits in the coefficient of the first operand. The number of places of rotation is taken from the absolute value of the second operand, with the rotation being to the left if the second operand is positive or to the right otherwise. >>> ExtendedContext.rotate(Decimal('34'), Decimal('8')) Decimal('400000003') >>> ExtendedContext.rotate(Decimal('12'), Decimal('9')) Decimal('12') >>> ExtendedContext.rotate(Decimal('123456789'), Decimal('-2')) Decimal('891234567') >>> ExtendedContext.rotate(Decimal('123456789'), Decimal('0')) Decimal('123456789') >>> ExtendedContext.rotate(Decimal('123456789'), Decimal('+2')) Decimal('345678912') >>> ExtendedContext.rotate(1333333, 1) Decimal('13333330') >>> ExtendedContext.rotate(Decimal(1333333), 1) Decimal('13333330') >>> ExtendedContext.rotate(1333333, Decimal(1)) Decimal('13333330') T)rœ)r0)rri)r/rrUr(r(r*riszContext.rotatecCst|dd}|j|ƒS)aÝReturns True if the two operands have the same exponent. The result is never affected by either the sign or the coefficient of either operand. >>> ExtendedContext.same_quantum(Decimal('2.17'), Decimal('0.001')) False >>> ExtendedContext.same_quantum(Decimal('2.17'), Decimal('0.01')) True >>> ExtendedContext.same_quantum(Decimal('2.17'), Decimal('1')) False >>> ExtendedContext.same_quantum(Decimal('Inf'), Decimal('-Inf')) True >>> ExtendedContext.same_quantum(10000, -1) True >>> ExtendedContext.same_quantum(Decimal(10000), -1) True >>> ExtendedContext.same_quantum(10000, Decimal(-1)) True T)rœ)rr&)r/rrUr(r(r*r&1szContext.same_quantumcCst|dd}|j||dS)a3Returns the first operand after adding the second value its exp. >>> ExtendedContext.scaleb(Decimal('7.50'), Decimal('-2')) Decimal('0.0750') >>> ExtendedContext.scaleb(Decimal('7.50'), Decimal('0')) Decimal('7.50') >>> ExtendedContext.scaleb(Decimal('7.50'), Decimal('3')) Decimal('7.50E+3') >>> ExtendedContext.scaleb(1, 4) Decimal('1E+4') >>> ExtendedContext.scaleb(Decimal(1), 4) Decimal('1E+4') >>> ExtendedContext.scaleb(1, Decimal(4)) Decimal('1E+4') T)rœ)r0)rrj)r/rrUr(r(r*rjIszContext.scalebcCst|dd}|j||dS)a{Returns a shifted copy of a, b times. The coefficient of the result is a shifted copy of the digits in the coefficient of the first operand. The number of places to shift is taken from the absolute value of the second operand, with the shift being to the left if the second operand is positive or to the right otherwise. Digits shifted into the coefficient are zeros. >>> ExtendedContext.shift(Decimal('34'), Decimal('8')) Decimal('400000000') >>> ExtendedContext.shift(Decimal('12'), Decimal('9')) Decimal('0') >>> ExtendedContext.shift(Decimal('123456789'), Decimal('-2')) Decimal('1234567') >>> ExtendedContext.shift(Decimal('123456789'), Decimal('0')) Decimal('123456789') >>> ExtendedContext.shift(Decimal('123456789'), Decimal('+2')) Decimal('345678900') >>> ExtendedContext.shift(88888888, 2) Decimal('888888800') >>> ExtendedContext.shift(Decimal(88888888), 2) Decimal('888888800') >>> ExtendedContext.shift(88888888, Decimal(2)) Decimal('888888800') T)rœ)r0)rrË)r/rrUr(r(r*rË\sz Context.shiftcCst|dd}|j|dS)a¦Square root of a non-negative number to context precision. If the result must be inexact, it is rounded using the round-half-even algorithm. >>> ExtendedContext.sqrt(Decimal('0')) Decimal('0') >>> ExtendedContext.sqrt(Decimal('-0')) Decimal('-0') >>> ExtendedContext.sqrt(Decimal('0.39')) Decimal('0.624499800') >>> ExtendedContext.sqrt(Decimal('100')) Decimal('10') >>> ExtendedContext.sqrt(Decimal('1')) Decimal('1') >>> ExtendedContext.sqrt(Decimal('1.0')) Decimal('1.0') >>> ExtendedContext.sqrt(Decimal('1.00')) Decimal('1.0') >>> ExtendedContext.sqrt(Decimal('7')) Decimal('2.64575131') >>> ExtendedContext.sqrt(Decimal('10')) Decimal('3.16227766') >>> ExtendedContext.sqrt(2) Decimal('1.41421356') >>> ExtendedContext.prec 9 T)rœ)r0)rr/)r/rr(r(r*r/zszContext.sqrtcCs8t|dd}|j||d}|tkr0td|ƒ‚n|SdS)a&Return the difference between the two operands. >>> ExtendedContext.subtract(Decimal('1.3'), Decimal('1.07')) Decimal('0.23') >>> ExtendedContext.subtract(Decimal('1.3'), Decimal('1.30')) Decimal('0.00') >>> ExtendedContext.subtract(Decimal('1.3'), Decimal('2.07')) Decimal('-0.77') >>> ExtendedContext.subtract(8, 5) Decimal('3') >>> ExtendedContext.subtract(Decimal(8), 5) Decimal('3') >>> ExtendedContext.subtract(8, Decimal(5)) Decimal('3') T)rœ)r0zUnable to convert %s to DecimalN)rrÄr–rs)r/rrUrÑr(r(r*Úsubtractšs zContext.subtractcCst|dd}|j|dS)a…Convert to a string, using engineering notation if an exponent is needed. Engineering notation has an exponent which is a multiple of 3. This can leave up to 3 digits to the left of the decimal place and may require the addition of either one or two trailing zeros. The operation is not affected by the context. >>> ExtendedContext.to_eng_string(Decimal('123E+1')) '1.23E+3' >>> ExtendedContext.to_eng_string(Decimal('123E+3')) '123E+3' >>> ExtendedContext.to_eng_string(Decimal('123E-10')) '12.3E-9' >>> ExtendedContext.to_eng_string(Decimal('-123E-12')) '-123E-12' >>> ExtendedContext.to_eng_string(Decimal('7E-7')) '700E-9' >>> ExtendedContext.to_eng_string(Decimal('7E+1')) '70' >>> ExtendedContext.to_eng_string(Decimal('0E+1')) '0.00E+3' T)rœ)r0)rr¶)r/rr(r(r*r¶±szContext.to_eng_stringcCst|dd}|j|dS)zyConverts a number to a string, using scientific notation. The operation is not affected by the context. T)rœ)r0)rrµ)r/rr(r(r*Ú to_sci_stringÍszContext.to_sci_stringcCst|dd}|j|dS)akRounds to an integer. When the operand has a negative exponent, the result is the same as using the quantize() operation using the given operand as the left-hand-operand, 1E+0 as the right-hand-operand, and the precision of the operand as the precision setting; Inexact and Rounded flags are allowed in this operation. The rounding mode is taken from the context. >>> ExtendedContext.to_integral_exact(Decimal('2.1')) Decimal('2') >>> ExtendedContext.to_integral_exact(Decimal('100')) Decimal('100') >>> ExtendedContext.to_integral_exact(Decimal('100.0')) Decimal('100') >>> ExtendedContext.to_integral_exact(Decimal('101.5')) Decimal('102') >>> ExtendedContext.to_integral_exact(Decimal('-101.5')) Decimal('-102') >>> ExtendedContext.to_integral_exact(Decimal('10E+5')) Decimal('1.0E+6') >>> ExtendedContext.to_integral_exact(Decimal('7.89E+77')) Decimal('7.89E+77') >>> ExtendedContext.to_integral_exact(Decimal('-Inf')) Decimal('-Infinity') T)rœ)r0)rr))r/rr(r(r*r)ÕszContext.to_integral_exactcCst|dd}|j|dS)aLRounds to an integer. When the operand has a negative exponent, the result is the same as using the quantize() operation using the given operand as the left-hand-operand, 1E+0 as the right-hand-operand, and the precision of the operand as the precision setting, except that no flags will be set. The rounding mode is taken from the context. >>> ExtendedContext.to_integral_value(Decimal('2.1')) Decimal('2') >>> ExtendedContext.to_integral_value(Decimal('100')) Decimal('100') >>> ExtendedContext.to_integral_value(Decimal('100.0')) Decimal('100') >>> ExtendedContext.to_integral_value(Decimal('101.5')) Decimal('102') >>> ExtendedContext.to_integral_value(Decimal('-101.5')) Decimal('-102') >>> ExtendedContext.to_integral_value(Decimal('10E+5')) Decimal('1.0E+6') >>> ExtendedContext.to_integral_value(Decimal('7.89E+77')) Decimal('7.89E+77') >>> ExtendedContext.to_integral_value(Decimal('-Inf')) Decimal('-Infinity') T)rœ)r0)rrý)r/rr(r(r*rýószContext.to_integral_value) NNNNNNNNN)N)rR)N)Xr2r3r4r5r‡r”r—r’r™rlr«rJr¡r*rIrnrer^r¦r¨r§rÈrær+r©rªrir«r¬r4ržr5r0r9r·rr¸r:r®r¯rÉrWrúr=r>r%r¡r?rr@rŒrArBrJrMrNrXrZr[rYr¾r\rªr]r°r±r`rarbr#rdr²r³rõrerÌrÙrir&rjrËr/r´r¶rµr)rýr…r(r(r(r*rBs® " $# % #2 P:&" c@s&eZdZd Zd dd„Zdd„ZeZdS)rjr>rVrWNcCsf|dkrd|_d|_d|_nFt|tƒrD|j|_t|jƒ|_|j|_n|d|_|d|_|d|_dS)Nr&r.r,)r>rVrWr`rr8r9rP)r/rur(r(r*r‡s z_WorkRep.__init__cCsd|j|j|jfS)Nz(%r, %r, %r))r>rVrW)r/r(r(r*r«(sz_WorkRep.__repr__)r>rVrW)N)r2r3r4r~r‡r«rµr(r(r(r*rjs rjcCsš|j|jkr|}|}n|}|}tt|jƒƒ}tt|jƒƒ}|jtd||dƒ}||jd|krpd|_||_|jd|j|j9_|j|_||fS)zcNormalizes op1, op2 to have the same exp and length of coefficient. Done during addition. r.r,rŸr…)rWrgrarVrª)rÁrÂrAZtmprˆZtmp_lenZ other_lenrWr(r(r*rÀ/srÀcCsb|dkrdS|dkr |d|Stt|ƒƒ}t|ƒt|jdƒƒ}||krPdS|d|SdS)a Given integers n and e, return n * 10**e if it's an integer, else None. The computation is designed to avoid computing large powers of 10 unnecessarily. >>> _decimal_lshift_exact(3, 4) 30000 >>> _decimal_lshift_exact(300, -999999999) # returns None r&rŸrRN)rarirgÚrstrip)r6r¯Zstr_nZval_nr(r(r*rOsrcCsF|dks|dkrtdƒ‚d}x$||kr@||||d?}}qW|S)zóClosest integer to the square root of the positive integer n. a is an initial approximation to the square root. Any positive integer will do for a, but the closer a is to the square root of n the faster convergence will be. r&z3Both arguments to _sqrt_nearest should be positive.r.)rm)r6rrUr(r(r*Ú _sqrt_nearestds r·cCs2d|>||?}}|d||d@|d@|kS)z‰Given an integer x and a nonnegative integer shift, return closest integer to x / 2**shift; use round-to-even in case of a tie. r.r,r()rrËrUrÐr(r(r*Ú_rshift_nearestssr¸cCs&t||ƒ\}}|d||d@|kS)zaClosest integer to a/b, a and b positive integers; rounds to even in the case of a tie. r,r.)rÉ)rrUrÐrÑr(r(r*Ú_div_nearest{sr¹rc Csî||}d}xn||kr*t|ƒ||>|ksF||krzt|ƒ||?|krzt||d>|t||t||ƒ|ƒƒ}|d7}qWtdtt|ƒƒd|ƒ}t||ƒ}t||ƒ}x0t|dddƒD]}t||ƒt|||ƒ}qÀWt|||ƒS)aÉInteger approximation to M*log(x/M), with absolute error boundable in terms only of x/M. Given positive integers x and M, return an integer approximation to M * log(x/M). For L = 8 and 0.1 <= x/M <= 10 the difference between the approximation and the exact result is at most 22. For L = 8 and 1.0 <= x/M <= 10.0 the difference is at most 15. In both cases these are upper bounds on the error; it will usually be much smaller.r&r.rŸr\iöÿÿÿr…)rir¹r·r¸rVrgrarþ) rÚMÚLrÚRÚTZyshiftÚwr‚r(r(r*Ú_ilogƒs r¿c Cs¶|d7}tt|ƒƒ}||||dk}|dkr”d|}|||}|dkrZ|d|9}nt|d|ƒ}t||ƒ}t|ƒ}t|||ƒ}||} nd}t|d|ƒ} t| |dƒS)z¾Given integers c, e and p with c > 0, p >= 0, compute an integer approximation to 10**p * log10(c*10**e), with an absolute error of at most 1. Assumes that c*10**e is not exactly 1.r,r.r&rŸr )rgrar¹r¿Ú _log10_digits) r-r¯r r.r€rºr‚Úlog_dZlog_10Zlog_tenpowerr(r(r*rL³s rLc CsÎ|d7}tt|ƒƒ}||||dk}|dkrr|||}|dkrR|d|9}nt|d|ƒ}t|d|ƒ}nd}|r¼ttt|ƒƒƒd}||dkr¶t|t||ƒd|ƒ}qÀd}nd}t||dƒS)z´Given integers c, e and p with c > 0, compute an integer approximation to 10**p * log(c*10**e), with an absolute error of at most 1. Assumes that c*10**e is not exactly 1.r,r.r&rŸr )rgrar¹r¿rirÀ) r-r¯r r.r€r‚rÁrZ f_log_tenr(r(r*rIÕs"rIc@s eZdZdZdd„Zdd„ZdS)Ú _Log10Memoizez¾Class to compute, store, and allow retrieval of, digits of the constant log(10) = 2.302585.... This constant is needed by Decimal.ln, Decimal.log10, Decimal.exp and Decimal.__pow__.cCs d|_dS)NZ/23025850929940456840179914546843642076011014886)ry)r/r(r(r*r‡sz_Log10Memoize.__init__cCsš|dkrtdƒ‚|t|jƒkr„d}xLd||d}tttd||ƒdƒƒ}||d…d|krdP|d7}q$W|jdƒdd …|_t|jd|d …ƒS)ztGiven an integer p >= 0, return floor(10**p)*log(10). For example, self.getdigits(3) returns 2302. r&zp should be nonnegativer\rŸr,r NrRr.r…)rmrgryrar¹r¿r¶rV)r/r rrºryr(r(r*Ú getdigitss z_Log10Memoize.getdigitsN)r2r3r4r5r‡rÃr(r(r(r*rÂsrÂc Cs°t||>|ƒ}tdtt|ƒƒd|ƒ}t||ƒ}||>}x.t|dddƒD]}t|||||ƒ}qRWx6t|ddd ƒD]"}||d>}t||||ƒ}q‚W||S) zëGiven integers x and M, M > 0, such that x/M is small in absolute value, compute an integer approximation to M*exp(x/M). For 0 <= x/M <= 2.4, the absolute error in the result is bounded by 60 (and is usually much smaller).rŸr\r.r&r,iöÿÿÿr…r…r…)r rVrgrar¹rþ) rrºr»r¼r½rZMshiftrr‚r(r(r*Ú_iexp&s rÄc Cs–|d7}td|tt|ƒƒdƒ}||}||}|dkrH|d|}n|d|}t|t|ƒƒ\}}t|d|ƒ}tt|d|ƒdƒ||dfS)aÐCompute an approximation to exp(c*10**e), with p decimal places of precision. Returns integers d, f such that: 10**(p-1) <= d <= 10**p, and (d-1)*10**f < exp(c*10**e) < (d+1)*10**f In other words, d*10**f is an approximation to exp(c*10**e) with p digits of precision, and with an error in d of at most 1. This is almost, but not quite, the same as the error being < 1ulp: when d = 10**(p-1) the error could be up to 10 ulp.r,r&r.rŸièr\)r¾rgrarÉrÀr¹rÄ) r-r¯r rrÐrËZcshiftZquotrr(r(r*r;Ksr;cCsèttt|ƒƒƒ|}t||||dƒ}||}|dkrJ||d|}nt||d|ƒ}|dkr´tt|ƒƒ|dk|dkkržd|ddd|} } qàd|d|} } n,t||d|dƒ\} } t| dƒ} | d7} | | fS)a5Given integers xc, xe, yc and ye representing Decimals x = xc*10**xe and y = yc*10**ye, compute x**y. Returns a pair of integers (c, e) such that: 10**(p-1) <= c <= 10**p, and (c-1)*10**e < x**y < (c+1)*10**e in other words, c*10**e is an approximation to x**y with p digits of precision, and with an error in c of at most 1. (This is almost, but not quite, the same as the error being < 1ulp: when c == 10**(p-1) we can only guarantee error < 10ulp.) We assume that: x is positive and not equal to 1, and y is nonzero. r.r&rŸ)rgrarirIr¹r;)rrrrr rUZlxcrËZpcrÊrWr(r(r*ros rr éFé5é(érCérŸr{) rÆÚ2Ú3Ú4Ú5Ú6Ú7Ú8r?cCs0|dkrtdƒ‚t|ƒ}dt|ƒ||dS)z@Compute a lower bound for 100*log10(c) for a positive integer c.r&z0The argument to _log10_lb should be nonnegative.r )rmrarg)r-Z correctionZstr_cr(r(r*r™srcCsLt|tƒr|St|tƒr t|ƒS|r8t|tƒr8tj|ƒS|rHtd|ƒ‚tS)zÙConvert other to Decimal. Verifies that it's ok to use in an implicit construction. If allow_float is true, allow conversion from float; this is used in the comparison methods (__eq__ and friends). zUnable to convert %s to Decimal)r`rrVrqrrrsr–)rˆrœZallow_floatr(r(r*r¤s rcCs´t|tƒr||fSt|tjƒrR|jsDt|jtt|j ƒ|j ƒ|jƒ}|t|jƒfS|rrt|tj ƒrr|jdkrr|j}t|tƒr¬tƒ}|r’d|jt<n|jtdƒ|tj|ƒfSttfS)zÔGiven a Decimal instance self and a Python object other, return a pair (s, o) of Decimal instances such that "s op o" is equivalent to "self op other" for any of the 6 comparison operators "op". r&r.z;strict semantics for mixing floats and Decimals are enabled)r`rÚ_numbersZRationalrQr7r8rarVr9ÚdenominatorrPÚ numeratorZComplexráràrqrrrrerrr–)r/rˆr”r0r(r(r*r•·s$ r•ri?B)rAr@rrrBr$r±rår])rAr@rra· # A numeric string consists of: # \s* (?P[-+])? # an optional sign, followed by either... ( (?=\d|\.\d) # ...a number (with at least one digit) (?P\d*) # having a (possibly empty) integer part (\.(?P\d*))? # followed by an optional fractional part (E(?P[-+]?\d+))? # followed by an optional exponent, or... | Inf(inity)? # ...an infinity, or... | (?Ps)? # ...an (optionally signaling) NaN # NaN (?P\d*) # with (possibly empty) diagnostic info. ) # \s* \Z z0*$z50*$zÉ\A (?: (?P.)? (?P[<>=^]) )? (?P[-+ ])? (?P\#)? (?P0)? (?P(?!0)\d+)? (?P,)? (?:\.(?P0|(?!0)\d+))? (?P[eEfFgGn%])? \Z cCs”tj|ƒ}|dkrtd|ƒ‚|jƒ}|d}|d}|ddk |d<|drv|dk rbtd|ƒ‚|dk rvtd|ƒ‚|p|d|d<|pˆd |d<|d dkr¢d|d <t|dp®d ƒ|d<|ddk rÒt|dƒ|d<|ddkrþ|ddksö|ddkrþd|d<|ddkrfd|d<|dkr&tjƒ}|ddk r@td|ƒ‚|d|d<|d|d<|d|d<n*|ddkr|d|d<ddg|d<d|d<|S)aÚParse and validate a format specifier. Turns a standard numeric format specifier into a dict, with the following entries: fill: fill character to pad field to minimum width align: alignment type, either '<', '>', '=' or '^' sign: either '+', '-' or ' ' minimumwidth: nonnegative integer giving minimum width zeropad: boolean, indicating whether to pad with zeros thousands_sep: string to use as thousands separator, or '' grouping: grouping for thousands separators, in format used by localeconv decimal_point: string to use for decimal point precision: nonnegative integer giving precision, or None type: one of the characters 'eEfFgG%', or None NzInvalid format specifier: ÚfillÚalignÚzeropadz7Fill character conflicts with '0' in format specifier: z2Alignment conflicts with '0' in format specifier: ú ú>r>rUÚminimumwidthrRrur&rmZgGnr.r6rsÚ thousands_sepzJExplicit thousands separator conflicts with 'n' type in format specifier: ÚgroupingÚ decimal_pointrTr\r®)Ú_parse_format_specifier_regexÚmatchrmÚ groupdictrVÚ_localeÚ localeconv)Úformat_specrqrvZformat_dictrÔrÕr(r(r*rwDsN rwc Cs´|d}|d}||t|ƒt|ƒ}|d}|dkrF|||}nj|dkr\|||}nT|dkrr|||}n>|dkr¨t|ƒd}|d |…||||d …}ntd ƒ‚|S)zÜGiven an unpadded, non-aligned numeric string 'body' and sign string 'sign', add padding and alignment conforming to the given format specifier dictionary 'spec' (as produced by parse_format_specifier). rÙrÔrÕúr|r{rÙrÔZpaddingrÕrƒZhalfr(r(r*ry”s"rycCspddlm}m}|sgS|ddkrJt|ƒdkrJ||dd…||d ƒƒS|d tjkrd|dd…Stdƒ‚dS)zyConvert a localeconv-style grouping into a (possibly infinite) iterable of integers representing group lengths. r&)ÚchainÚrepeatr.r,Nz unrecognised format for groupingr…r…r r…r…)Ú itertoolsrærçrgràÚCHAR_MAXrm)rÛrærçr(r(r*Ú_group_lengths¯s rêcCsä|d}|d}g}xÀt|ƒD]€}|dkr2tdƒ‚ttt|ƒ|dƒ|ƒ}|jd|t|ƒ||d…ƒ|d|…}||8}|r’|dkr’P|t|ƒ8}qWtt|ƒ|dƒ}|jd|t|ƒ||d…ƒ|jt|ƒƒS)anInsert thousands separators into a digit string. spec is a dictionary whose keys should include 'thousands_sep' and 'grouping'; typically it's the result of parsing the format specifier using _parse_format_specifier. The min_width keyword argument gives the minimum length of the result, which will be padded on the left with zeros if necessary. If necessary, the zero padding adds an extra '0' on the left to avoid a leading thousands separator. For example, inserting commas every three digits in '123456', with min_width=8, gives '0,123,456', even though that has length 9. rÚrÛr&zgroup length should be positiver.rRN)rêrmrªr¾rgrnroÚreversed)ryr{Ú min_widthÚseprÛÚgroupsr.r(r(r*Ú_insert_thousands_sepÆs $$rïcCs$|rdS|ddkr|dSdSdS)zDetermine sign character.rUr>z +rTNr()Úis_negativer{r(r(r*rxës rxcCs²t||ƒ}|s|dr"|d|}|dks6|ddkr\dddddœ|d}|d j||ƒ7}|dd krp|d 7}|dr’|dt|ƒt|ƒ}nd}t|||ƒ}t||||ƒS) acFormat a number, given the following data: is_negative: true if the number is negative, else false intpart: string of digits that must appear before the decimal point fracpart: string of digits that must come after the point exp: exponent, as an integer spec: dictionary resulting from parsing the format specifier This function uses the information in spec to: insert separators (decimal separator and thousands separators) format the sign format the exponent add trailing '%' for the '%' type zero-pad if necessary fill and align if necessary ZaltrÜr&rmrvr°r¯)r°r¯rtrsz{0}{1:+}rrrÖrÙ)rxÚformatrgrïry)rðrwrxrWr{r>Zecharrìr(r(r*rzõs rzZInfz-Infr¬llÿÿÿÿlüÿÿÿÿÇNÎZoiÀ«æ)N)F)r&)r)r)FF)FiÁ½ðÿ)N)r.r…)}r5Ú__all__r2rEÚ__version__Z__libmpdec_version__Zmathr|ZnumbersrÑrFÚcollectionsr'Z_namedtuplerÚImportErrorrrrrrrrrr$Úmaxsizer r!r"r#ÚArithmeticErrorrrr rÚZeroDivisionErrorr rrrrrr rrrsrrr£r˜rKr^rCrGrMÚhasattrrLrHrrrrr7ÚNumberÚregisterrOrrjrÀrVrr rr·r¸r¹r¿rLrIrÂrÃrÀrÄr;rrrr•rrrÚreÚcompileÚVERBOSEÚ IGNORECASErÞrbrêrîÚDOTALLrÝZlocaleràrwryrêrïrxrzrHrGr;r7rr6r=Ú hash_infoÚmodulusr¥rcr£rÜr¢r¤r¦r(r(r(r*Úqs€ & . ^ 0",# %$+ * P % )