decimal —- 十进制定点和浮点运算

源码: Lib/decimal.py


decimal 模块为快速正确舍入的十进制浮点运算提供支持。 与 float 数据类型相比,它具有以下几个优点:

  • Decimal 类型的“设计是基于考虑人类习惯的浮点数模型,并且因此具有以下最高指导原则 —— 计算机必须提供与人们在学校所学习的算术相一致的算术。” —— 摘自 decimal 算术规范描述。

  • Decimal 数字的表示是完全精确的。 相比之下,1.12.2 这样的数字在二进制浮点中没有精确的表示。 最终用户通常不希望 1.1 + 2.2 如二进制浮点数表示那样被显示为 3.3000000000000003

  • 精确性会延续到算术类操作中。 对于 decimal 浮点数,0.1 + 0.1 + 0.1 - 0.3 会精确地等于零。 而对于二进制浮点数,结果则为 5.5511151231257827e-017 。 虽然接近于零,但其中的误差将妨碍可靠的相等性检验,并且误差还会不断累积。 因此,decimal 更适合具有严格相等不变性要求的会计类应用。

  • 十进制模块包含有效位的概念,因此 1.30 + 1.20 的结果是 2.50 。 保留尾随零以表示有效位。 这是货币的惯用表示方法。乘法则沿用 “教科书“ 中:保留被乘数中的所有数字的方法。 例如, 1.3 * 1.2 结果是 1.561.30 * 1.20 结果是 1.5600

  • 与基于硬件的二进制浮点不同,十进制模块具有用户可更改的精度(默认为28位),可以与给定问题所需的一样大:

    1. >>> from decimal import *
    2. >>> getcontext().prec = 6
    3. >>> Decimal(1) / Decimal(7)
    4. Decimal('0.142857')
    5. >>> getcontext().prec = 28
    6. >>> Decimal(1) / Decimal(7)
    7. Decimal('0.1428571428571428571428571429')
  • 二进制和 decimal 浮点数都是根据已发布的标准实现的。 虽然内置浮点类型只公开其功能的一小部分,但 decimal 模块公开了标准的所有必需部分。 在需要时,程序员可以完全控制舍入和信号处理。 这包括通过使用异常来阻止任何不精确操作来强制执行精确算术的选项。

  • decimal 模块旨在支持“无偏差,精确无舍入的十进制算术(有时称为定点数算术)和有舍入的浮点数算术”。 —— 摘自 decimal 算术规范说明。

该模块的设计以三个概念为中心:decimal 数值,算术上下文和信号。

decimal 数值是不可变对象。 它由符号,系数和指数位组成。 为了保持有效位,系数位不会截去末尾零。 decimal 数值也包括特殊值例如 Infinity-InfinityNaN 。 该标准还区分 -0+0

算术的上下文是指定精度、舍入规则、指数限制、指示操作结果的标志以及确定符号是否被视为异常的陷阱启用器的环境。 舍入选项包括 ROUND_CEILINGROUND_DOWNROUND_FLOORROUND_HALF_DOWN, ROUND_HALF_EVENROUND_HALF_UPROUND_UP 以及 ROUND_05UP.

信号是在计算过程中出现的异常条件组。 根据应用程序的需要,信号可能会被忽略,被视为信息,或被视为异常。 十进制模块中的信号有:ClampedInvalidOperationDivisionByZeroInexactRoundedSubnormalOverflowUnderflow 以及 FloatOperation

对于每个信号,都有一个标志和一个陷阱启动器。 遇到信号时,其标志设置为 1 ,然后,如果陷阱启用器设置为 1 ,则引发异常。 标志是粘性的,因此用户需要在监控计算之前重置它们。

参见

快速入门教程

通常使用 decimal 的方式是先导入该模块,通过 getcontext() 查看当前上下文,并在必要时为精度、舍入或启用的陷阱设置新值:

  1. >>> from decimal import *
  2. >>> getcontext()
  3. Context(prec=28, rounding=ROUND_HALF_EVEN, Emin=-999999, Emax=999999,
  4. capitals=1, clamp=0, flags=[], traps=[Overflow, DivisionByZero,
  5. InvalidOperation])
  6. >>> getcontext().prec = 7 # Set a new precision

可以基于整数、字符串、浮点数或元组构造 Decimal 实例。 基于整数或浮点数构造将执行该整数或浮点值的精确转换。 Decimal 数字包括特殊值例如 NaN 表示“非数字”,正的和负的 Infinity-0

  1. >>> getcontext().prec = 28
  2. >>> Decimal(10)
  3. Decimal('10')
  4. >>> Decimal('3.14')
  5. Decimal('3.14')
  6. >>> Decimal(3.14)
  7. Decimal('3.140000000000000124344978758017532527446746826171875')
  8. >>> Decimal((0, (3, 1, 4), -2))
  9. Decimal('3.14')
  10. >>> Decimal(str(2.0 ** 0.5))
  11. Decimal('1.4142135623730951')
  12. >>> Decimal(2) ** Decimal('0.5')
  13. Decimal('1.414213562373095048801688724')
  14. >>> Decimal('NaN')
  15. Decimal('NaN')
  16. >>> Decimal('-Infinity')
  17. Decimal('-Infinity')

如果 FloatOperation 信号被捕获,构造函数中的小数和浮点数的意外混合或排序比较会引发异常

  1. >>> c = getcontext()
  2. >>> c.traps[FloatOperation] = True
  3. >>> Decimal(3.14)
  4. Traceback (most recent call last):
  5. File "<stdin>", line 1, in <module>
  6. decimal.FloatOperation: [<class 'decimal.FloatOperation'>]
  7. >>> Decimal('3.5') < 3.7
  8. Traceback (most recent call last):
  9. File "<stdin>", line 1, in <module>
  10. decimal.FloatOperation: [<class 'decimal.FloatOperation'>]
  11. >>> Decimal('3.5') == 3.5
  12. True

3.3 新版功能.

新 Decimal 的重要性仅由输入的位数决定。 上下文精度和舍入仅在算术运算期间发挥作用。

  1. >>> getcontext().prec = 6
  2. >>> Decimal('3.0')
  3. Decimal('3.0')
  4. >>> Decimal('3.1415926535')
  5. Decimal('3.1415926535')
  6. >>> Decimal('3.1415926535') + Decimal('2.7182818285')
  7. Decimal('5.85987')
  8. >>> getcontext().rounding = ROUND_UP
  9. >>> Decimal('3.1415926535') + Decimal('2.7182818285')
  10. Decimal('5.85988')

如果超出了 C 版本的内部限制,则构造一个 decimal 将引发 InvalidOperation

  1. >>> Decimal("1e9999999999999999999")
  2. Traceback (most recent call last):
  3. File "<stdin>", line 1, in <module>
  4. decimal.InvalidOperation: [<class 'decimal.InvalidOperation'>]

在 3.3 版更改.

Decimal 数字能很好地与 Python 的其余部分交互。 以下是一个小小的 decimal 浮点数飞行马戏团:

  1. >>> data = list(map(Decimal, '1.34 1.87 3.45 2.35 1.00 0.03 9.25'.split()))
  2. >>> max(data)
  3. Decimal('9.25')
  4. >>> min(data)
  5. Decimal('0.03')
  6. >>> sorted(data)
  7. [Decimal('0.03'), Decimal('1.00'), Decimal('1.34'), Decimal('1.87'),
  8. Decimal('2.35'), Decimal('3.45'), Decimal('9.25')]
  9. >>> sum(data)
  10. Decimal('19.29')
  11. >>> a,b,c = data[:3]
  12. >>> str(a)
  13. '1.34'
  14. >>> float(a)
  15. 1.34
  16. >>> round(a, 1)
  17. Decimal('1.3')
  18. >>> int(a)
  19. 1
  20. >>> a * 5
  21. Decimal('6.70')
  22. >>> a * b
  23. Decimal('2.5058')
  24. >>> c % a
  25. Decimal('0.77')

Decimal 也可以使用一些数学函数:

  1. >>> getcontext().prec = 28
  2. >>> Decimal(2).sqrt()
  3. Decimal('1.414213562373095048801688724')
  4. >>> Decimal(1).exp()
  5. Decimal('2.718281828459045235360287471')
  6. >>> Decimal('10').ln()
  7. Decimal('2.302585092994045684017991455')
  8. >>> Decimal('10').log10()
  9. Decimal('1')

quantize() 方法将数字舍入为固定指数。 此方法对于将结果舍入到固定的位置的货币应用程序非常有用:

  1. >>> Decimal('7.325').quantize(Decimal('.01'), rounding=ROUND_DOWN)
  2. Decimal('7.32')
  3. >>> Decimal('7.325').quantize(Decimal('1.'), rounding=ROUND_UP)
  4. Decimal('8')

如上所示,getcontext() 函数访问当前上下文并允许更改设置。 这种方法满足大多数应用程序的需求。

对于更高级的工作,使用 Context() 构造函数创建备用上下文可能很有用。 要使用备用活动,请使用 setcontext() 函数。

根据标准,decimal 模块提供了两个现成的标准上下文 BasicContextExtendedContext 。 前者对调试特别有用,因为许多陷阱都已启用:

  1. >>> myothercontext = Context(prec=60, rounding=ROUND_HALF_DOWN)
  2. >>> setcontext(myothercontext)
  3. >>> Decimal(1) / Decimal(7)
  4. Decimal('0.142857142857142857142857142857142857142857142857142857142857')
  5. >>> ExtendedContext
  6. Context(prec=9, rounding=ROUND_HALF_EVEN, Emin=-999999, Emax=999999,
  7. capitals=1, clamp=0, flags=[], traps=[])
  8. >>> setcontext(ExtendedContext)
  9. >>> Decimal(1) / Decimal(7)
  10. Decimal('0.142857143')
  11. >>> Decimal(42) / Decimal(0)
  12. Decimal('Infinity')
  13. >>> setcontext(BasicContext)
  14. >>> Decimal(42) / Decimal(0)
  15. Traceback (most recent call last):
  16. File "<pyshell#143>", line 1, in -toplevel-
  17. Decimal(42) / Decimal(0)
  18. DivisionByZero: x / 0

上下文还具有用于监视计算期间遇到的异常情况的信号标志。 标志保持设置直到明确清除,因此最好通过使用 clear_flags() 方法清除每组受监控计算之前的标志。:

  1. >>> setcontext(ExtendedContext)
  2. >>> getcontext().clear_flags()
  3. >>> Decimal(355) / Decimal(113)
  4. Decimal('3.14159292')
  5. >>> getcontext()
  6. Context(prec=9, rounding=ROUND_HALF_EVEN, Emin=-999999, Emax=999999,
  7. capitals=1, clamp=0, flags=[Inexact, Rounded], traps=[])

flags 条目显示对 Pi 的有理逼近被舍入(超出上下文精度的数字被抛弃)并且结果是不精确的(一些丢弃的数字不为零)。

使用上下文的 traps 字段中的字典设置单个陷阱:

  1. >>> setcontext(ExtendedContext)
  2. >>> Decimal(1) / Decimal(0)
  3. Decimal('Infinity')
  4. >>> getcontext().traps[DivisionByZero] = 1
  5. >>> Decimal(1) / Decimal(0)
  6. Traceback (most recent call last):
  7. File "<pyshell#112>", line 1, in -toplevel-
  8. Decimal(1) / Decimal(0)
  9. DivisionByZero: x / 0

大多数程序仅在程序开始时调整当前上下文一次。 并且,在许多应用程序中,数据在循环内单个强制转换为 Decimal 。 通过创建上下文集和小数,程序的大部分操作数据与其他 Python 数字类型没有区别。

Decimal 对象

class decimal.Decimal(value=’0’, context=None)

根据 value 构造一个新的 Decimal 对象。

value 可以是整数,字符串,元组,float ,或另一个 Decimal 对象。 如果没有给出 value,则返回 Decimal('0')。 如果 value 是一个字符串,它应该在前导和尾随空格字符以及下划线被删除之后符合十进制数字字符串语法:

  1. sign ::= '+' | '-'
  2. digit ::= '0' | '1' | '2' | '3' | '4' | '5' | '6' | '7' | '8' | '9'
  3. indicator ::= 'e' | 'E'
  4. digits ::= digit [digit]...
  5. decimal-part ::= digits '.' [digits] | ['.'] digits
  6. exponent-part ::= indicator [sign] digits
  7. infinity ::= 'Infinity' | 'Inf'
  8. nan ::= 'NaN' [digits] | 'sNaN' [digits]
  9. numeric-value ::= decimal-part [exponent-part] | infinity
  10. numeric-string ::= [sign] numeric-value | [sign] nan

当上面出现 digit 时也允许其他十进制数码。 其中包括来自各种其他语言系统的十进制数码(例如阿拉伯-印地语和天城文的数码)以及全宽数码 '\uff10''\uff19'

如果 value 是一个 tuple ,它应该有三个组件,一个符号( 0 表示正数或 1 表示负数),一个数字的 tuple 和整数指数。 例如, Decimal((0, (1, 4, 1, 4), -3)) 返回 Decimal('1.414')

如果 valuefloat ,则二进制浮点值无损地转换为其精确的十进制等效值。 此转换通常需要53位或更多位数的精度。 例如, Decimal(float('1.1')) 转换为``Decimal(‘1.100000000000000088817841970012523233890533447265625’)``。

context 精度不会影响存储的位数。 这完全由 value 中的位数决定。 例如,Decimal('3.00000') 记录所有五个零,即使上下文精度只有三。

context 参数的目的是确定 value 是格式错误的字符串时该怎么做。 如果上下文陷阱 InvalidOperation,则引发异常;否则,构造函数返回一个新的 Decimal,其值为 NaN

构造完成后, Decimal 对象是不可变的。

在 3.2 版更改: 现在允许构造函数的参数为 float 实例。

在 3.3 版更改: float 参数在设置 FloatOperation 陷阱时引发异常。 默认情况下,陷阱已关闭。

在 3.6 版更改: 允许下划线进行分组,就像代码中的整数和浮点文字一样。

十进制浮点对象与其他内置数值类型共享许多属性,例如 floatint 。 所有常用的数学运算和特殊方法都适用。 同样,十进制对象可以复制、pickle、打印、用作字典键、用作集合元素、比较、排序和强制转换为另一种类型(例如 floatint )。

算术对十进制对象和算术对整数和浮点数有一些小的差别。 当余数运算符 % 应用于Decimal对象时,结果的符号是 被除数 的符号,而不是除数的符号:

  1. >>> (-7) % 4
  2. 1
  3. >>> Decimal(-7) % Decimal(4)
  4. Decimal('-3')

整数除法运算符 // 的行为类似,返回真商的整数部分(截断为零)而不是它的向下取整,以便保留通常的标识 x == (x // y) * y + x % y:

  1. >>> -7 // 4
  2. -2
  3. >>> Decimal(-7) // Decimal(4)
  4. Decimal('-1')

%// 运算符实现了 remainderdivide-integer 操作(分别),如规范中所述。

十进制对象通常不能与浮点数或 fractions.Fraction 实例在算术运算中结合使用:例如,尝试将 Decimal 加到 float ,将引发 TypeError。 但是,可以使用 Python 的比较运算符来比较 Decimal 实例 x 和另一个数字 y 。 这样可以避免在对不同类型的数字进行相等比较时混淆结果。

在 3.2 版更改: 现在完全支持 Decimal 实例和其他数字类型之间的混合类型比较。

除了标准的数字属性,十进制浮点对象还有许多专门的方法:

  • adjusted()

    在移出系数最右边的数字之后返回调整后的指数,直到只剩下前导数字:Decimal('321e+5').adjusted() 返回 7 。 用于确定最高有效位相对于小数点的位置。

  • as_integer_ratio()

    返回一对 (n, d) 整数,表示给定的 Decimal 实例作为分数、最简形式项并带有正分母:

    1. >>> Decimal('-3.14').as_integer_ratio()
    2. (-157, 50)

    转换是精确的。 在 Infinity 上引发 OverflowError ,在 NaN 上引起 ValueError 。

3.6 新版功能.

  • as_tuple()

    返回一个 named tuple 表示的数字: DecimalTuple(sign, digits, exponent)

  • canonical()

    返回参数的规范编码。 目前,一个 Decimal 实例的编码始终是规范的,因此该操作返回其参数不变。

  • compare(other, context=None)

    比较两个 Decimal 实例的值。 compare() 返回一个 Decimal 实例,如果任一操作数是 NaN ,那么结果是 NaN

    1. a or b is a NaN ==> Decimal('NaN')
    2. a < b ==> Decimal('-1')
    3. a == b ==> Decimal('0')
    4. a > b ==> Decimal('1')
  • compare_signal(other, context=None)

    除了所有 NaN 信号之外,此操作与 compare() 方法相同。 也就是说,如果两个操作数都不是信令NaN,那么任何静默的 NaN 操作数都被视为信令NaN。

  • compare_total(other, context=None)

    使用它们的抽象表示而不是它们的数值来比较两个操作数。 类似于 compare() 方法,但结果给出了一个总排序 Decimal 实例。 两个 Decimal 实例具有相同的数值但不同的表示形式在此排序中比较不相等:

    1. >>> Decimal('12.0').compare_total(Decimal('12'))
    2. Decimal('-1')

    静默和发出信号的 NaN 也包括在总排序中。 这个函数的结果是 Decimal('0') 如果两个操作数具有相同的表示,或是 Decimal('-1') 如果第一个操作数的总顺序低于第二个操作数,或是 Decimal('1') 如果第一个操作数在总顺序中高于第二个操作数。 有关总排序的详细信息,请参阅规范。

    此操作不受上下文影响且静默:不更改任何标志且不执行舍入。 作为例外,如果无法准确转换第二个操作数,则C版本可能会引发InvalidOperation。

  • compare_total_mag(other, context=None)

    比较两个操作数使用它们的抽象表示而不是它们的值,如 compare_total(),但忽略每个操作数的符号。 x.compare_total_mag(y) 相当于 x.copy_abs().compare_total(y.copy_abs())

    此操作不受上下文影响且静默:不更改任何标志且不执行舍入。 作为例外,如果无法准确转换第二个操作数,则C版本可能会引发InvalidOperation。

  • conjugate()

    只返回self,这种方法只符合 Decimal 规范。

  • copy_abs()

    返回参数的绝对值。 此操作不受上下文影响并且是静默的:没有更改标志且不执行舍入。

  • copy_negate()

    回到参数的否定。 此操作不受上下文影响并且是静默的:没有标志更改且不执行舍入。

  • copy_sign(other, context=None)

    返回第一个操作数的副本,其符号设置为与第二个操作数的符号相同。 例如:

    1. >>> Decimal('2.3').copy_sign(Decimal('-1.5'))
    2. Decimal('-2.3')

    此操作不受上下文影响且静默:不更改任何标志且不执行舍入。 作为例外,如果无法准确转换第二个操作数,则C版本可能会引发InvalidOperation。

  • exp(context=None)

    返回给定数字的(自然)指数函数``e**x``的值。结果使用 ROUND_HALF_EVEN 舍入模式正确舍入。

    1. >>> Decimal(1).exp()
    2. Decimal('2.718281828459045235360287471')
    3. >>> Decimal(321).exp()
    4. Decimal('2.561702493119680037517373933E+139')
  • from_float(f)

    将浮点数转换为十进制数的类方法。

    注意, Decimal.from_float(0.1) 与 Decimal(‘0.1’) 不同。 由于 0.1 在二进制浮点中不能精确表示,因此该值存储为最接近的可表示值,即 0x1.999999999999ap-4 。 十进制的等效值是`0.1000000000000000055511151231257827021181583404541015625`。

    注解

    从 Python 3.2 开始,Decimal 实例也可以直接从 float 构造。

    1. >>> Decimal.from_float(0.1)
    2. Decimal('0.1000000000000000055511151231257827021181583404541015625')
    3. >>> Decimal.from_float(float('nan'))
    4. Decimal('NaN')
    5. >>> Decimal.from_float(float('inf'))
    6. Decimal('Infinity')
    7. >>> Decimal.from_float(float('-inf'))
    8. Decimal('-Infinity')

    3.1 新版功能.

  • fma(other, third, context=None)

    混合乘法加法。 返回 selfother+third ,中间乘积 selfother 没有舍入。

    1. >>> Decimal(2).fma(3, 5)
    2. Decimal('11')
  • is_canonical()

    如果参数是规范的,则为返回 True,否则为 False 。 目前,Decimal 实例总是规范的,所以这个操作总是返回 True

  • is_finite()

    如果参数是一个有限的数,则返回为 True ;如果参数为无穷大或 NaN ,则返回为 False

  • is_infinite()

    如果参数为正负无穷大,则返回为 True ,否则为 False

  • is_nan()

    如果参数为 NaN (无论是否静默),则返回为 True ,否则为 False

  • is_normal(context=None)

    如果参数是一个 标准的 有限数则返回 True。 如果参数为零、次标准数、无穷大或 NaN 则返回 False

  • is_qnan()

    如果参数为静默 NaN,返回 True,否则返回 False

  • is_signed()

    如果参数带有负号,则返回为 True,否则返回 False。注意,0 和 NaN 都可带有符号。

  • is_snan()

    如果参数为显式 NaN,则返回 True,否则返回 False

  • is_subnormal(context=None)

    如果参数为次标准数,则返回 True,否则返回 False

  • is_zero()

    如果参数是0(正负皆可),则返回 True,否则返回 False

  • ln(context=None)

    返回操作数的自然对数(以 e 为底)。结果是使用 ROUND_HALF_EVEN 舍入模式正确舍入的。

  • log10(context=None)

    返回操作数的以十为底的对数。结果是使用 ROUND_HALF_EVEN 舍入模式正确舍入的。

  • logb(context=None)

    对于一个非零数,返回其运算数的调整后指数作为一个 Decimal 实例。 如果运算数为零将返回 Decimal('-Infinity') 并且产生 the DivisionByZero 标志。如果运算数是无限大则返回 Decimal('Infinity')

  • logical_and(other, context=None)

    logical_and() 是需要两个 逻辑运算数 的逻辑运算(参考 逻辑操作数 )。按位输出两运算数的 and 运算的结果。

  • logical_invert(context=None)

    logical_invert() 是一个逻辑运算。结果是操作数的按位求反。

  • logical_or(other, context=None)

    logical_or() 是需要两个 logical operands 的逻辑运算(请参阅 逻辑操作数 )。结果是两个运算数的按位的 or 运算。

  • logical_xor(other, context=None)

    logical_xor() 是需要两个 逻辑运算数 的逻辑运算(参考 逻辑操作数 )。结果是按位输出的两运算数的异或运算。

  • max(other, context=None)

    max(self, other) 一样,除了在返回之前应用上下文舍入规则并且用信号通知或忽略 NaN 值(取决于上下文以及它们是发信号还是安静)。

  • max_mag(other, context=None)

    max() 方法相似,但是操作数使用绝对值完成比较。

  • min(other, context=None)

    min(self, other) 一样,除了在返回之前应用上下文舍入规则并且用信号通知或忽略 NaN 值(取决于上下文以及它们是发信号还是安静)。

  • min_mag(other, context=None)

    min() 方法相似,但是操作数使用绝对值完成比较。

  • next_minus(context=None)

    返回小于给定操作数的上下文中可表示的最大数字(或者当前线程的上下文中的可表示的最大数字如果没有给定上下文)。

  • next_plus(context=None)

    返回大于给定操作数的上下文中可表示的最小数字(或者当前线程的上下文中的可表示的最小数字如果没有给定上下文)。

  • next_toward(other, context=None)

    如果两运算数不相等,返回在第二个操作数的方向上最接近第一个操作数的数。如果两操作数数值上相等,返回将符号设置为与第二个运算数相同的第一个运算数的拷贝。

  • normalize(context=None)

    通过去除尾随的零并将所有结果等于 Decimal('0') 的转化为 Decimal('0e0') 来标准化数字。用于为等效类的属性生成规范值。比如, Decimal('32.100')Decimal('0.321000e+2') 都被标准化为相同的值 Decimal('32.1')

  • number_class(context=None)

    返回一个字符串描述运算数的 class 。返回值是以下十个字符串中的一个。

    • "-Infinity" ,指示运算数为负无穷大。

    • "-Normal" ,指示该运算数是负正常数字。

    • "-Subnormal" ,指示该运算数是负的次标准数。

    • "-Zero" ,指示该运算数是负零。

    • "-Zero" ,指示该运算数是正零。

    • "+Subnormal" ,指示该运算数是正的次标准数。

    • "+Normal" ,指示该运算数是正的标准数。

    • "+Infinity" ,指示该运算数是正无穷。

    • "NaN" ,指示该运算数是肃静 NaN (非数字)。

    • "sNaN" ,指示该运算数是信号 NaN 。

  • quantize(exp, rounding=None, context=None)

    返回的值等于舍入后的第一个运算数并且具有第二个操作数的指数。

    1. >>> Decimal('1.41421356').quantize(Decimal('1.000'))
    2. Decimal('1.414')

    与其他运算不同,如果量化运算后的系数长度大于精度,那么会发出一个 InvalidOperation 信号。这保证了除非有一个错误情况,量化指数恒等于右手运算数的指数。

    与其他运算不同,量化永不信号下溢,即使结果不正常且不精确。

    如果第二个运算数的指数大于第一个运算数的指数那或许需要舍入。在这种情况下,舍入模式由给定 rounding 参数决定,其余的由给定 context 参数决定;如果参数都未给定,使用当前线程上下文的舍入模式。

    每当结果的指数大于 Emax 或小于 Etiny 就会返回错误。

  • radix()

    返回 Decimal(10),即 Decimal 类进行所有算术运算所用的数制(基数)。 这是为保持与规范描述的兼容性而加入的。

  • remainder_near(other, context=None)

    返回 self 除以 other 的余数。 这与 self % other 的区别在于所选择的余数要使其绝对值最小化。 更准确地说,返回值为 self - n * other 其中 n 是最接近 self / other 的实际值的整数,并且如果两个整数与实际值的差相等则会选择其中的偶数。

    如果结果为零则其符号将为 self 的符号。

    1. >>> Decimal(18).remainder_near(Decimal(10))
    2. Decimal('-2')
    3. >>> Decimal(25).remainder_near(Decimal(10))
    4. Decimal('5')
    5. >>> Decimal(35).remainder_near(Decimal(10))
    6. Decimal('-5')
  • rotate(other, context=None)

    返回对第一个操作数的数码按第二个操作数所指定的数量进行轮转的结果。 第二个操作数必须为 -precision 至 precision 精度范围内的整数。 第二个操作数的绝对值给出要轮转的位数。 如果第二个操作数为正值则向左轮转;否则向右轮转。 如有必要第一个操作数的系数会在左侧填充零以达到 precision 所指定的长度。 第一个操作数的符号和指数保持不变。

  • same_quantum(other, context=None)

    检测自身与 other 是否具有相同的指数或是否均为 NaN

    此操作不受上下文影响且静默:不更改任何标志且不执行舍入。 作为例外,如果无法准确转换第二个操作数,则C版本可能会引发InvalidOperation。

  • scaleb(other, context=None)

    返回第一个操作数使用第二个操作数对指数进行调整的结果。 等价于返回第一个操作数乘以 10**other 的结果。 第二个操作数必须为整数。

  • shift(other, context=None)

    返回第一个操作数的数码按第二个操作数所指定的数量进行移位的结果。 第二个操作数必须为 -precision 至 precision 范围内的整数。 第二个操作数的绝对值给出要移动的位数。 如果第二个操作数为正值则向左移位;否则向右移位。 移入系数的数码为零。 第一个操作数的符号和指数保持不变。

  • sqrt(context=None)

    返回参数的平方根精确到完整精度。

  • to_eng_string(context=None)

    转换为字符串,如果需要指数则会使用工程标注法。

    工程标注法的指数是 3 的倍数。 这会在十进制位的左边保留至多 3 个数码,并可能要求添加一至两个末尾零。

    例如,此方法会将 Decimal('123E+1') 转换为 Decimal('1.23E+3')

  • to_integral(rounding=None, context=None)

    to_integral_value() 方法相同。 保留 to_integral 名称是为了与旧版本兼容。

  • to_integral_exact(rounding=None, context=None)

    舍入到最接近的整数,发出信号 Inexact 或者如果发生舍入则相应地发出信号 Rounded。 如果给出 rounding 形参则由其确定舍入模式,否则由给定的 context 来确定。 如果没有给定任何形参则会使用当前上下文的舍入模式。

  • to_integral_value(rounding=None, context=None)

    舍入到最接近的整数而不发出 InexactRounded 信号。 如果给出 rounding 则会应用其所指定的舍入模式;否则使用所提供的 context 或当前上下文的舍入方法。

逻辑操作数

logical_and(), logical_invert(), logical_or()logical_xor() 方法期望其参数为 逻辑操作数逻辑操作数 是指数位与符号位均为零的 Decimal 实例,并且其数字位均为 01

上下文对象

上下文是算术运算所在的环境。 它们管理精度、设置舍入规则、确定将哪些信号视为异常,并限制指数的范围。

每个线程都有自己的当前上下文,可使用 getcontext()setcontext() 函数来读取或修改:

decimal.getcontext()

返回活动线程的当前上下文。

decimal.setcontext(c)

将活动线程的当前上下文设为 c

你也可以使用 with 语句和 localcontext() 函数来临时改变活动上下文。

decimal.localcontext(ctx=None)

返回一个上下文管理器,它将在进入 with 语句时将活动线程的当前上下文设为 ctx 的一个副本并在退出 with 语句时恢复之前的上下文。 如果未指定上下文,则会使用当前上下文的一个副本。

例如,以下代码会将当前 decimal 精度设为 42 位,执行一个运算,然后自动恢复之前的上下文:

  1. from decimal import localcontext
  2. with localcontext() as ctx:
  3. ctx.prec = 42 # Perform a high precision calculation
  4. s = calculate_something()
  5. s = +s # Round the final result back to the default precision

新的上下文也可使用下述的 Context 构造器来创建。 此外,模块还提供了三种预设的上下文:

class decimal.BasicContext

这是由通用十进制算术规范描述所定义的标准上下文。 精度设为九。 舍入设为 ROUND_HALF_UP。 清除所有旗标。 启用所有陷阱(视为异常),但 Inexact, RoundedSubnormal 除外。

由于启用了许多陷阱,此上下文适用于进行调试。

class decimal.ExtendedContext

这是由通用十进制算术规范描述所定义的标准上下文。 精度设为九。 舍入设为 ROUND_HALF_EVEN。 清除所有旗标。 不启用任何陷阱(因此在计算期间不会引发异常)。

由于禁用了陷阱,此上下文适用于希望结果值为 NaNInfinity 而不是引发异常的应用。 这允许应用在出现当其他情况下会中止程序的条件时仍能完成运行。

class decimal.DefaultContext

此上下文被 Context 构造器用作新上下文的原型。 改变一个字段(例如精度)的效果将是改变 Context 构造器所创建的新上下文的默认值。

此上下文最适用于多线程环境。 在线程开始前改变一个字段具有设置全系统默认值的效果。 不推荐在线程开始后改变字段,因为这会要求线程同步避免竞争条件。

在单线程环境中,最好完全不使用此上下文。 而是简单地电显式创建上下文,具体如下所述。

默认值为 prec\=28, rounding\=ROUND_HALF_EVEN,并为 Overflow, InvalidOperationDivisionByZero 启用陷阱。

在已提供的三种上下文之外,还可以使用 Context 构造器创建新的上下文。

class decimal.Context(prec=None, rounding=None, Emin=None, Emax=None, capitals=None, clamp=None, flags=None, traps=None)

创建一个新上下文。 如果某个字段未指定或为 None,则从 DefaultContext 拷贝默认值。 如果 flags 字段未指定或为 None,则清空所有旗标。

prec 为一个 [1, MAX_PREC] 范围内的整数,用于设置该上下文中算术运算的精度。

rounding 选项应为 Rounding Modes 小节中列出的常量之一。

trapsflags 字段列出要设置的任何信号。 通常,新上下文应当只设置 traps 而让 flags 为空。

EminEmax 字段给定指数所允许的外部上限。 Emin 必须在 [MIN_EMIN, 0] 范围内,Emax 在 [0, MAX_EMAX] 范围内。

capitals 字段为 01 (默认值)。 如果设为 1,指数将附带打印大写的 E;其他情况则将使用小写的 e: Decimal('6.02e+23')

clamp 字段为 0 (默认值) 或 1。 如果设为 1,则 Decimal 实例的指数 e 的表示范围在此上下文中将严格限制为 Emin - prec + 1 <= e <= Emax - prec + 1。 如果 clamp0 则将适用较弱的条件: Decimal 实例调整后的指数最大值为 Emax。 当 clamp1 时,一个较大的普通数值将在可能的情况下减小其指数并为其系统添加相应数量的零,以便符合指数值限制;这可以保持数字值但会丢失有效末尾零的信息。 例如:

  1. >>> Context(prec=6, Emax=999, clamp=1).create_decimal('1.23e999')
  2. Decimal('1.23000E+999')

clamp 值为 1 时即允许与在 IEEE 754 中描述的固定宽度十进制交换格式保持兼容性。

Context 类定义了几种通用方法以及大量直接在给定上下文中进行算术运算的方法。 此外,对于上述的每种 Decimal 方法(不包括 adjusted()as_tuple() 方法)都有一个相应的 Context 方法。 例如,对于一个 Context 的实例 CDecimal 的实例 xC.exp(x) 就等价于 x.exp(context=C)。 每个 Context 方法都接受一个 Python 整数(即 int 的实例)在任何接受 Decimal 的实例的地方使用。

  • clear_flags()

    将所有旗标重置为 0

  • clear_traps()

    将所有陷阱重置为零 0

    3.3 新版功能.

  • copy()

    返回上下文的一个副本。

  • copy_decimal(num)

    返回 Decimal 实例 num 的一个副本。

  • create_decimal(num)

    基于 num 创建一个新 Decimal 实例但使用 self 作为上下文。 与 Decimal 构造器不同,该上下文的精度、舍入方法、旗标和陷阱会被应用于转换过程。

    此方法很有用处,因为常量往往被给予高于应用所需的精度。 另一个好处在于立即执行舍入可以消除超出当前精度的数位所导致的意外效果。 在下面的示例中,使用未舍入的输入意味着在总和中添加零会改变结果:

    1. >>> getcontext().prec = 3
    2. >>> Decimal('3.4445') + Decimal('1.0023')
    3. Decimal('4.45')
    4. >>> Decimal('3.4445') + Decimal(0) + Decimal('1.0023')
    5. Decimal('4.44')

    此方法实现了 IBM 规格描述中的转换为数字操作。 如果参数为字符串,则不允许有开头或末尾的空格或下划线。

  • create_decimal_from_float(f)

    基于浮点数 f 创建一个新的 Decimal 实例,但会使用 self 作为上下文来执行舍入。 与 Decimal.from_float() 类方法不同,上下文的精度、舍入方法、旗标和陷阱会应用到转换中。

    1. >>> context = Context(prec=5, rounding=ROUND_DOWN)
    2. >>> context.create_decimal_from_float(math.pi)
    3. Decimal('3.1415')
    4. >>> context = Context(prec=5, traps=[Inexact])
    5. >>> context.create_decimal_from_float(math.pi)
    6. Traceback (most recent call last):
    7. ...
    8. decimal.Inexact: None

    3.1 新版功能.

  • Etiny()

    返回一个等于 Emin - prec + 1 的值即次标准化结果中的最小指数值。 当发生向下溢出时,指数会设为 Etiny

  • Etop()

    返回一个等于 Emax - prec + 1 的值。

使用 decimal 的通常方式是创建 Decimal 实例然后对其应用算术运算,这些运算发生在活动线程的当前上下文中。 一种替代方式则是使用上下文的方法在特定上下文中进行计算。 这些方法类似于 Decimal 类的方法,在此仅简单地重新列出。

  • abs(x)

    返回 x 的绝对值。

  • add(x, y)

    返回 xy 的和。

  • canonical(x)

    返回相同的 Decimal 对象 x

  • compare(x, y)

    xy 进行数值比较。

  • compare_signal(x, y)

    对两个操作数进行数值比较。

  • compare_total(x, y)

    对两个操作数使用其抽象表示进行比较。

  • compare_total_mag(x, y)

    对两个操作数使用其抽象表示进行比较,忽略符号。

  • copy_abs(x)

    返回 x 的副本,符号设为 0。

  • copy_negate(x)

    返回 x 的副本,符号取反。

  • copy_sign(x, y)

    y 拷贝符号至 x

  • divide(x, y)

    返回 x 除以 y 的结果。

  • divide_int(x, y)

    返回 x 除以 y 的结果,截短为整数。

  • divmod(x, y)

    两个数字相除并返回结果的整数部分。

  • exp(x)

    返回 e ** x。

  • fma(x, y, z)

    返回 x 乘以 y 再加 z 的结果。

  • is_canonical(x)

    如果 x 是规范的则返回 True;否则返回 False

  • is_finite(x)

    如果 x 为有限的则返回``True``;否则返回 False

  • is_infinite(x)

    如果 x 是无限的则返回 True;否则返回 False

  • is_nan(x)

    如果 x 是 qNaN 或 sNaN 则返回 True;否则返回 False

  • is_normal(x)

    如果 x 是标准数则返回 True;否则返回 False

  • is_qnan(x)

    如果 x 是静默 NaN 则返回 True;否则返回 False

  • is_signed(x)

    x 是负数则返回 True;否则返回 False

  • is_snan(x)

    如果 x 是显式 NaN 则返回 True;否则返回 False

  • is_subnormal(x)

    如果 x 是次标准数则返回 True;否则返回 False

  • is_zero(x)

    如果 x 为零则返回 True;否则返回 False

  • ln(x)

    返回 x 的自然对数(以 e 为底)。

  • log10(x)

    返回 x 的以 10 为底的对数。

  • logb(x)

    返回操作数的 MSD 等级的指数。

  • logical_and(x, y)

    在操作数的每个数位间应用逻辑运算 and

  • logical_invert(x)

    反转 x 中的所有数位。

  • logical_or(x, y)

    在操作数的每个数位间应用逻辑运算 or

  • logical_xor(x, y)

    在操作数的每个数位间应用逻辑运算 xor

  • max(x, y)

    对两个值执行数字比较并返回其中的最大值。

  • max_mag(x, y)

    对两个值执行忽略正负号的数字比较。

  • min(x, y)

    对两个值执行数字比较并返回其中的最小值。

  • min_mag(x, y)

    对两个值执行忽略正负号的数字比较。

  • minus(x)

    对应于 Python 中的单目前缀取负运算符执行取负操作。

  • multiply(x, y)

    返回 xy 的积。

  • next_minus(x)

    返回小于 x 的最大数字表示形式。

  • next_plus(x)

    返回大于 x 的最小数字表示形式。

  • next_toward(x, y)

    返回 x 趋向于 y 的最接近的数字。

  • normalize(x)

    x 改写为最简形式。

  • number_class(x)

    返回 x 的类的表示。

  • plus(x)

    对应于 Python 中的单目前缀取正运算符执行取正操作。 此操作将应用上下文精度和舍入,因此它 不是 标识运算。

  • power(x, y, modulo=None)

    返回 xy 次方,如果给出了模数 modulo 则取其余数。

    如为两个参数则计算 x**y。 如果 x 为负值则 y 必须为整数。 除非 y 为整数且结果为有限值并可在 ‘precision’ 位内精确表示否则结果将是不精确的。 上下文的舍入模式将被使用。 结果在 Python 版中总是会被正确地舍入。

    Decimal(0) ** Decimal(0) 结果为 InvalidOperation,而如果 InvalidOperation 未被捕获,则结果为 Decimal('NaN')

    在 3.3 版更改: C 模块计算 power() 时会使用已正确舍入的 exp()ln() 函数。 结果是经过良好定义的,但仅限于“几乎总是正确地舍入”。

    带有三个参数时,计算 (x**y) % modulo。 对于三个参数的形式,参数将会应用以下限制:

    • 三个参数必须都是整数

    • y 必须是非负数

    • xy 至少有一个不为零

    • modulo 必须不为零且至多有 ‘precision’ 位

    来自 Context.power(x, y, modulo) 的结果值等于使用无限精度计算 (x**y) % modulo 所得到的值,但其计算过程更高效。 结果的指数为零,无论 x, ymodulo 的指数是多少。 结果值总是完全精确的。

  • quantize(x, y)

    返回的值等于 x (舍入后),并且指数为 y

  • radix()

    恰好返回 10,因为这是 Decimal 对象 :)

  • remainder(x, y)

    返回整除所得到的余数。

    结果的符号,如果不为零,则与原始除数的符号相同。

  • remainder_near(x, y)

    返回 x - y * n,其中 n 为最接近 x / y 实际值的整数(如结果为 0 则其符号将与 x 的符号相同)。

  • rotate(x, y)

    返回 x 翻转 y 次的副本。

  • same_quantum(x, y)

    如果两个操作数具有相同的指数则返回 True

  • scaleb(x, y)

    返回第一个操作数添加第二个值的指数后的结果。

  • shift(x, y)

    返回 x 变换 y 次的副本。

  • sqrt(x)

    非负数基于上下文精度的平方根。

  • subtract(x, y)

    返回 xy 的差。

  • to_eng_string(x)

    转换为字符串,如果需要指数则会使用工程标注法。

    工程标注法的指数是 3 的倍数。 这会在十进制位的左边保留至多 3 个数码,并可能要求添加一至两个末尾零。

  • to_integral_exact(x)

    舍入到一个整数。

  • to_sci_string(x)

    使用科学计数法将一个数字转换为字符串。

常量

本节中的常量仅与 C 模块相关。 它们也被包含在纯 Python 版本以保持兼容性。

32位

64位

    decimal.MAX_PREC

425000000

999999999999999999

    decimal.MAX_EMAX

425000000

999999999999999999

    decimal.MIN_EMIN

-425000000

-999999999999999999

    decimal.MIN_ETINY

-849999999

-1999999999999999997

decimal.HAVE_THREADS

该值为 True。 已弃用,因为 Python 现在总是启用线程。

3.9 版后已移除.

decimal.HAVE_CONTEXTVAR

默认值为 True。 如果 Python 编译版本 使用了 --without-decimal-contextvar 选项来配置,则 C 版本会使用线程局部而非协程局部上下文并且该值为 False。 这在某些嵌套上下文场景中将会稍快一些。

3.9 新版功能: 向下移植到 3.7 和 3.8。

舍入模式

decimal.ROUND_CEILING

舍入方向为 Infinity

decimal.ROUND_DOWN

舍入方向为零。

decimal.ROUND_FLOOR

舍入方向为 -Infinity

decimal.ROUND_HALF_DOWN

舍入到最接近的数,同样接近则舍入方向为零。

decimal.ROUND_HALF_EVEN

舍入到最接近的数,同样接近则舍入到最接近的偶数。

decimal.ROUND_HALF_UP

舍入到最接近的数,同样接近则舍入到零的反方向。

decimal.ROUND_UP

舍入到零的反方向。

decimal.ROUND_05UP

如果最后一位朝零的方向舍入后为 0 或 5 则舍入到零的反方向;否则舍入方向为零。

信号

信号代表在计算期间引发的条件。 每个信号对应于一个上下文旗标和一个上下文陷阱启用器。

上下文旗标将在遇到特定条件时被设定。 在完成计算之后,将为了获得信息而检测旗标(例如确定计算是否精确)。 在检测旗标后,请确保在开始下一次计算之前清除所有旗标。

如果为信号设定了上下文的陷阱启用器,则条件会导致特定的 Python 异常被引发。 举例来说,如果设定了 DivisionByZero 陷阱,则当遇到此条件时就将引发 DivisionByZero 异常。

class decimal.Clamped

修改一个指数以符合表示限制。

通常,限位将在一个指数超出上下文的 EminEmax 限制时发生。 在可能的情况下,会通过给系数添加零来将指数缩减至符合限制。

class decimal.DecimalException

其他信号的基类,并且也是 ArithmeticError 的一个子类。

class decimal.DivisionByZero

非无限数被零除的信号。

可在除法、取余队法或对一个数求负数次幂时发生。 如果此信号未被陷阱捕获,则返回 Infinity-Infinity 并且由对计算的输入来确定正负符号。

class decimal.Inexact

表明发生了舍入且结果是不精确的。

有非零数位在舍入期间被丢弃的信号。 舍入结果将被返回。 此信号旗标或陷阱被用于检测结果不精确的情况。

class decimal.InvalidOperation

执行了一个无效的操作。

表明请求了一个无意义的操作。 如未被陷阱捕获则返回 NaN。 可能的原因包括:

  1. Infinity - Infinity
  2. 0 * Infinity
  3. Infinity / Infinity
  4. x % 0
  5. Infinity % x
  6. sqrt(-x) and x > 0
  7. 0 ** 0
  8. x ** (non-integer)
  9. x ** Infinity

class decimal.Overflow

数值的溢出。

表明在发生舍入之后的指数大于 Emax。 如果未被陷阱捕获,则结果将取决于舍入模式,或者向下舍入为最大的可表示有限数,或者向上舍入为 Infinity。 无论哪种情况,都将引发 InexactRounded 信号。

class decimal.Rounded

发生了舍入,但或许并没有信息丢失。

一旦舍入丢弃了数位就会发出此信号;即使被丢弃的数位是零 (例如将 5.00 舍入为 5.0)。 如果未被陷阱捕获,则不经修改地返回结果。 此信号用于检测有效位数的丢弃。

class decimal.Subnormal

在舍入之前指数低于 Emin

当操作结果是次标准数(即指数过小)时就会发出此信号。 如果未被陷阱捕获,则不经修改过返回结果。

class decimal.Underflow

数字向下溢出导致结果舍入到零。

当一个次标准数结果通过舍入转为零时就会发出此信号。 同时还将引发 InexactSubnormal 信号。

class decimal.FloatOperation

为 float 和 Decimal 的混合启用更严格的语义。

如果信号未被捕获(默认),则在 Decimal 构造器、create_decimal() 和所有比较运算中允许 float 和 Decimal 的混合。 转换和比较都是完全精确的。 发生的任何混合运算都将通过在上下文旗标中设置 FloatOperation 来静默地记录。 通过 from_float()create_decimal_from_float() 进行显式转换则不会设置旗标。

在其他情况下(即信号被捕获),则只静默执行相等性比较和显式转换。 所有其他混合运算都将引发 FloatOperation

以下表格总结了信号的层级结构:

  1. exceptions.ArithmeticError(exceptions.Exception)
  2. DecimalException
  3. Clamped
  4. DivisionByZero(DecimalException, exceptions.ZeroDivisionError)
  5. Inexact
  6. Overflow(Inexact, Rounded)
  7. Underflow(Inexact, Rounded, Subnormal)
  8. InvalidOperation
  9. Rounded
  10. Subnormal
  11. FloatOperation(DecimalException, exceptions.TypeError)

浮点数说明

通过提升精度来解决舍入错误

使用十进制浮点数可以消除十进制表示错误(即能够完全精确地表示 0.1 这样的数);然而,某些运算在非零数位超出给定的精度时仍然可能导致舍入错误。

舍入错误的影响可能因接近相互抵销的加减运算被放大从而导致丢失有效位。 Knuth 提供了两个指导性示例,其中出现了精度不足的浮点算术舍入,导致加法的交换律和分配律被打破:

  1. # Examples from Seminumerical Algorithms, Section 4.2.2.
  2. >>> from decimal import Decimal, getcontext
  3. >>> getcontext().prec = 8
  4. >>> u, v, w = Decimal(11111113), Decimal(-11111111), Decimal('7.51111111')
  5. >>> (u + v) + w
  6. Decimal('9.5111111')
  7. >>> u + (v + w)
  8. Decimal('10')
  9. >>> u, v, w = Decimal(20000), Decimal(-6), Decimal('6.0000003')
  10. >>> (u*v) + (u*w)
  11. Decimal('0.01')
  12. >>> u * (v+w)
  13. Decimal('0.0060000')

decimal 模块则可以通过充分地扩展精度来避免有效位的丢失:

  1. >>> getcontext().prec = 20
  2. >>> u, v, w = Decimal(11111113), Decimal(-11111111), Decimal('7.51111111')
  3. >>> (u + v) + w
  4. Decimal('9.51111111')
  5. >>> u + (v + w)
  6. Decimal('9.51111111')
  7. >>>
  8. >>> u, v, w = Decimal(20000), Decimal(-6), Decimal('6.0000003')
  9. >>> (u*v) + (u*w)
  10. Decimal('0.0060000')
  11. >>> u * (v+w)
  12. Decimal('0.0060000')

特殊的值

decimal 模块的数字系统提供了一些特殊的值,包括 NaN, sNaN, -Infinity, Infinity 以及两种零值 +0-0

无穷大可以使用 Decimal('Infinity') 来构建。 它们也可以在不捕获 DivisionByZero 信号捕获时通过除以零来产生。 类似地,当不捕获 Overflow 信号时,也可以通过舍入到超出最大可表示数字限制的方式产生无穷大的结果。

无穷大是有符号的(仿射)并可用于算术运算,它们会被当作极其巨大的不确定数字来处理。 例如,无穷大加一个常量结果也将为无穷大。

某些不存在有效结果的运算将会返回 NaN,或者如果捕获了 InvalidOperation 信号则会引发一个异常。 例如,0/0 会返回 NaN 表示结果“不是一个数字”。 这样的 NaN 是静默产生的,并且在产生之后参与其它计算时总是会得到 NaN 的结果。 这种行为对于偶而缺少输入的各类计算都很有用处 —- 它允许在将特定结果标记为无效的同时让计算继续运行。

另一种变体形式是 sNaN,它在每次运算后会发出信号而不是保持静默。 当对于无效结果需要中断计算进行特别处理时,这是一个很有用的返回值。

Python 中比较运算符的行为在涉及 NaN 时可能会令人有点惊讶。 相等性检测在操作数中有静默型或信号型 NaN 时总是会返回 False (即使是执行 Decimal('NaN')==Decimal('NaN')),而不等性检测总是会返回 True。 当尝试使用 <, <=, >>= 运算符中的任何一个来比较两个 Decimal 值时,如果运算数中有 NaN 则将引发 InvalidOperation 信号,如果此信号未被捕获则将返回 False。 请注意通用十进制算术规范并未规定直接比较行为;这些涉及 NaN 的比较规则来自于 IEEE 854 标准 (见第 5.7 节表 3)。 要确保严格符合标准,请改用 compare()compare-signal() 方法。

有符号零值可以由向下溢出的运算产生。 它们保留符号是为了让运算结果能以更高的精度传递。 由于它们的大小为零,正零和负零会被视为相等,且它们的符号具有信息。

在这两个不相同但却相等的有符号零之外,还存在几种零的不同表示形式,它们的精度不同但值也都相等。 这需要一些时间来逐渐适应。 对于习惯了标准浮点表示形式的眼睛来说,以下运算返回等于零的值并不是显而易见的:

  1. >>> 1 / Decimal('Infinity')
  2. Decimal('0E-1000026')

使用线程

getcontext() 函数会为每个线程访问不同的 Context 对象。 具有单独线程上下文意味着线程可以修改上下文 (例如 getcontext().prec=10) 而不影响其他线程。

类似的 setcontext() 会为当前上下文的目标自动赋值。

如果在调用 setcontext() 之前调用了 getcontext(),则 getcontext() 将自动创建一个新的上下文在当前线程中使用。

新的上下文拷贝自一个名为 DefaultContext 的原型上下文。 要控制默认值以便每个线程在应用运行期间都使用相同的值,可以直接修改 DefaultContext 对象。 这应当在任何线程启动 之前 完成以使得调用 getcontext() 的线程之间不会产生竞争条件。 例如:

  1. # Set applicationwide defaults for all threads about to be launched
  2. DefaultContext.prec = 12
  3. DefaultContext.rounding = ROUND_DOWN
  4. DefaultContext.traps = ExtendedContext.traps.copy()
  5. DefaultContext.traps[InvalidOperation] = 1
  6. setcontext(DefaultContext)
  7. # Afterwards, the threads can be started
  8. t1.start()
  9. t2.start()
  10. t3.start()
  11. . . .

例程

以下是一些用作工具函数的例程,它们演示了使用 Decimal 类的各种方式:

  1. def moneyfmt(value, places=2, curr='', sep=',', dp='.',
  2. pos='', neg='-', trailneg=''):
  3. """Convert Decimal to a money formatted string.
  4. places: required number of places after the decimal point
  5. curr: optional currency symbol before the sign (may be blank)
  6. sep: optional grouping separator (comma, period, space, or blank)
  7. dp: decimal point indicator (comma or period)
  8. only specify as blank when places is zero
  9. pos: optional sign for positive numbers: '+', space or blank
  10. neg: optional sign for negative numbers: '-', '(', space or blank
  11. trailneg:optional trailing minus indicator: '-', ')', space or blank
  12. >>> d = Decimal('-1234567.8901')
  13. >>> moneyfmt(d, curr='$')
  14. '-$1,234,567.89'
  15. >>> moneyfmt(d, places=0, sep='.', dp='', neg='', trailneg='-')
  16. '1.234.568-'
  17. >>> moneyfmt(d, curr='$', neg='(', trailneg=')')
  18. '($1,234,567.89)'
  19. >>> moneyfmt(Decimal(123456789), sep=' ')
  20. '123 456 789.00'
  21. >>> moneyfmt(Decimal('-0.02'), neg='<', trailneg='>')
  22. '<0.02>'
  23. """
  24. q = Decimal(10) ** -places # 2 places --> '0.01'
  25. sign, digits, exp = value.quantize(q).as_tuple()
  26. result = []
  27. digits = list(map(str, digits))
  28. build, next = result.append, digits.pop
  29. if sign:
  30. build(trailneg)
  31. for i in range(places):
  32. build(next() if digits else '0')
  33. if places:
  34. build(dp)
  35. if not digits:
  36. build('0')
  37. i = 0
  38. while digits:
  39. build(next())
  40. i += 1
  41. if i == 3 and digits:
  42. i = 0
  43. build(sep)
  44. build(curr)
  45. build(neg if sign else pos)
  46. return ''.join(reversed(result))
  47. def pi():
  48. """Compute Pi to the current precision.
  49. >>> print(pi())
  50. 3.141592653589793238462643383
  51. """
  52. getcontext().prec += 2 # extra digits for intermediate steps
  53. three = Decimal(3) # substitute "three=3.0" for regular floats
  54. lasts, t, s, n, na, d, da = 0, three, 3, 1, 0, 0, 24
  55. while s != lasts:
  56. lasts = s
  57. n, na = n+na, na+8
  58. d, da = d+da, da+32
  59. t = (t * n) / d
  60. s += t
  61. getcontext().prec -= 2
  62. return +s # unary plus applies the new precision
  63. def exp(x):
  64. """Return e raised to the power of x. Result type matches input type.
  65. >>> print(exp(Decimal(1)))
  66. 2.718281828459045235360287471
  67. >>> print(exp(Decimal(2)))
  68. 7.389056098930650227230427461
  69. >>> print(exp(2.0))
  70. 7.38905609893
  71. >>> print(exp(2+0j))
  72. (7.38905609893+0j)
  73. """
  74. getcontext().prec += 2
  75. i, lasts, s, fact, num = 0, 0, 1, 1, 1
  76. while s != lasts:
  77. lasts = s
  78. i += 1
  79. fact *= i
  80. num *= x
  81. s += num / fact
  82. getcontext().prec -= 2
  83. return +s
  84. def cos(x):
  85. """Return the cosine of x as measured in radians.
  86. The Taylor series approximation works best for a small value of x.
  87. For larger values, first compute x = x % (2 * pi).
  88. >>> print(cos(Decimal('0.5')))
  89. 0.8775825618903727161162815826
  90. >>> print(cos(0.5))
  91. 0.87758256189
  92. >>> print(cos(0.5+0j))
  93. (0.87758256189+0j)
  94. """
  95. getcontext().prec += 2
  96. i, lasts, s, fact, num, sign = 0, 0, 1, 1, 1, 1
  97. while s != lasts:
  98. lasts = s
  99. i += 2
  100. fact *= i * (i-1)
  101. num *= x * x
  102. sign *= -1
  103. s += num / fact * sign
  104. getcontext().prec -= 2
  105. return +s
  106. def sin(x):
  107. """Return the sine of x as measured in radians.
  108. The Taylor series approximation works best for a small value of x.
  109. For larger values, first compute x = x % (2 * pi).
  110. >>> print(sin(Decimal('0.5')))
  111. 0.4794255386042030002732879352
  112. >>> print(sin(0.5))
  113. 0.479425538604
  114. >>> print(sin(0.5+0j))
  115. (0.479425538604+0j)
  116. """
  117. getcontext().prec += 2
  118. i, lasts, s, fact, num, sign = 1, 0, x, 1, x, 1
  119. while s != lasts:
  120. lasts = s
  121. i += 2
  122. fact *= i * (i-1)
  123. num *= x * x
  124. sign *= -1
  125. s += num / fact * sign
  126. getcontext().prec -= 2
  127. return +s

Decimal 常见问题

Q. 总是输入 decimal.Decimal('1234.5') 是否过于笨拙。 在使用交互解释器时有没有最小化输入量的方式?

A. 有些用户会将构造器简写为一个字母:

  1. >>> D = decimal.Decimal
  2. >>> D('1.23') + D('3.45')
  3. Decimal('4.68')

Q. 在带有两个十进制位的定点数应用中,有些输入值具有许多位,需要被舍入。 另一些数则不应具有多余位,需要验证有效性。 这种情况应该用什么方法?

A. 用 quantize() 方法舍入到固定数量的十进制位。 如果设置了 Inexact 陷阱,它也适用于验证有效性:

  1. >>> TWOPLACES = Decimal(10) ** -2 # same as Decimal('0.01')
  1. >>> # Round to two places
  2. >>> Decimal('3.214').quantize(TWOPLACES)
  3. Decimal('3.21')
  1. >>> # Validate that a number does not exceed two places
  2. >>> Decimal('3.21').quantize(TWOPLACES, context=Context(traps=[Inexact]))
  3. Decimal('3.21')
  1. >>> Decimal('3.214').quantize(TWOPLACES, context=Context(traps=[Inexact]))
  2. Traceback (most recent call last):
  3. ...
  4. Inexact: None

Q. 当我使用两个有效位的输入时,我要如何在一个应用中保持有效位不变?

A. 某些运算例如与整数相加、相减和相乘将会自动保留固定的小数位数。 其他运算,例如相除和非整数相乘则将会改变小数位数,需要再加上 quantize() 处理步骤:

  1. >>> a = Decimal('102.72') # Initial fixed-point values
  2. >>> b = Decimal('3.17')
  3. >>> a + b # Addition preserves fixed-point
  4. Decimal('105.89')
  5. >>> a - b
  6. Decimal('99.55')
  7. >>> a * 42 # So does integer multiplication
  8. Decimal('4314.24')
  9. >>> (a * b).quantize(TWOPLACES) # Must quantize non-integer multiplication
  10. Decimal('325.62')
  11. >>> (b / a).quantize(TWOPLACES) # And quantize division
  12. Decimal('0.03')

在开发定点数应用时,更方便的做法是定义处理 quantize() 步骤的函数:

  1. >>> def mul(x, y, fp=TWOPLACES):
  2. ... return (x * y).quantize(fp)
  3. >>> def div(x, y, fp=TWOPLACES):
  4. ... return (x / y).quantize(fp)
  1. >>> mul(a, b) # Automatically preserve fixed-point
  2. Decimal('325.62')
  3. >>> div(b, a)
  4. Decimal('0.03')

Q. 表示同一个值有许多方式。 数字 200, 200.000, 2E202E+4 的值都相同但有精度不同。 是否有办法将它们转换为一个可识别的规范值?

A. normalize() 方法可将所有相同的值映射为统一表示形式:

  1. >>> values = map(Decimal, '200 200.000 2E2 .02E+4'.split())
  2. >>> [v.normalize() for v in values]
  3. [Decimal('2E+2'), Decimal('2E+2'), Decimal('2E+2'), Decimal('2E+2')]

Q. 有些十进制值总是被打印为指数表示形式。 是否有办法得到一个非指数表示形式?

A. 对于某些值来说,指数表示形式是表示系数中有效位的唯一办法。 例如,将 5.0E+3 表示为 5000 可以让值保持恒定,但是无法显示原本的两位有效数字。

如果一个应用不必关心追踪有效位,则可以很容易地移除指数和末尾的零,丢弃有效位但让值保持不变:

  1. >>> def remove_exponent(d):
  2. ... return d.quantize(Decimal(1)) if d == d.to_integral() else d.normalize()
  1. >>> remove_exponent(Decimal('5E+3'))
  2. Decimal('5000')

Q. 是否有办法将一个普通浮点数转换为 Decimal

A. 是的,任何二进制浮点数都可以精确地表示为 Decimal 值,但完全精确的转换可能需要比平常感觉更高的精度:

  1. >>> Decimal(math.pi)
  2. Decimal('3.141592653589793115997963468544185161590576171875')

Q. 在一个复杂的计算中,我怎样才能保证不会得到由精度不足和舍入异常所导致的虚假结果。

A. 使用 decimal 模块可以很容易地检测结果。 最好的做法是使用更高的精度和不同的舍入模式重新进行计算。 明显不同的结果表明存在精度不足、舍入模式问题、不符合条件的输入或是结果不稳定的算法。

Q. 我发现上下文精度的应用只针对运算结果而不针对输入。在混合使用不同精度的值时有什么需要注意的吗?

A. 是的。 原则上所有值都会被视为精确值,在这些值上进行的算术运算也是如此。 只有结果会被舍入。 对于输入来说其好处是“所输入即所得”。 而其缺点则是如果你忘记了输入没有被舍入,结果看起来可能会很奇怪:

  1. >>> getcontext().prec = 3
  2. >>> Decimal('3.104') + Decimal('2.104')
  3. Decimal('5.21')
  4. >>> Decimal('3.104') + Decimal('0.000') + Decimal('2.104')
  5. Decimal('5.20')

解决办法是提高精度或使用单目加法运算对输入执行强制舍入:

  1. >>> getcontext().prec = 3
  2. >>> +Decimal('1.23456789') # unary plus triggers rounding
  3. Decimal('1.23')

此外,还可以使用 Context.create_decimal() 方法在创建输入时执行舍入:

  1. >>> Context(prec=5, rounding=ROUND_DOWN).create_decimal('1.2345678')
  2. Decimal('1.2345')

Q. CPython 实现对于巨大数字是否足够快速?

A. 是的。 在 CPython 和 PyPy3 实现中,decimal 模块的 C/CFFI 版本集成了高速 libmpdec 库用于实现任意精度正确舍入的十进制浮点算术 1libmpdec 会对中等大小的数字使用 Karatsuba 乘法 而对非常巨大的数字使用 数字原理变换#Number-theoretic_transform)。

必须要对任意精度算术适配上下文。 EminEmax 应当总是设为最大值,clamp 应当总是设为 0 (默认值)。 设置 prec 需要十分谨慎。

进行大数字算术的最便捷方式也是使用 prec 的最大值 2:

  1. >>> setcontext(Context(prec=MAX_PREC, Emax=MAX_EMAX, Emin=MIN_EMIN))
  2. >>> x = Decimal(2) ** 256
  3. >>> x / 128
  4. Decimal('904625697166532776746648320380374280103671755200316906558262375061821325312')

对于不精确的结果,在 64 位平台上 MAX_PREC 的值太大了,可用的内存将会不足:

  1. >>> Decimal(1) / 3
  2. Traceback (most recent call last):
  3. File "<stdin>", line 1, in <module>
  4. MemoryError

在具有超量分配的系统上 (即 Linux),一种更复杂的方式根据可用的 RAM 大小来调整 prec。 假设你有 8GB 的 RAM 并期望同时有 10 个操作数,每个最多使用 500MB:

  1. >>> import sys
  2. >>>
  3. >>> # Maximum number of digits for a single operand using 500MB in 8-byte words
  4. >>> # with 19 digits per word (4-byte and 9 digits for the 32-bit build):
  5. >>> maxdigits = 19 * ((500 * 1024**2) // 8)
  6. >>>
  7. >>> # Check that this works:
  8. >>> c = Context(prec=maxdigits, Emax=MAX_EMAX, Emin=MIN_EMIN)
  9. >>> c.traps[Inexact] = True
  10. >>> setcontext(c)
  11. >>>
  12. >>> # Fill the available precision with nines:
  13. >>> x = Decimal(0).logical_invert() * 9
  14. >>> sys.getsizeof(x)
  15. 524288112
  16. >>> x + 2
  17. Traceback (most recent call last):
  18. File "<stdin>", line 1, in <module>
  19. decimal.Inexact: [<class 'decimal.Inexact'>]

总体而言(特别是在没有超量分配的系统上),如果期望所有计算都是精确的则推荐预估更严格的边界并设置 Inexact 陷阱。

1

3.3 新版功能.

2

在 3.9 版更改: 此方式现在适用于除了非整数乘方以外的所有精确结果。