mpmath库详解

• 对全部变量设置任意精度

import mpmath
mpmath.mp.dps = 50
x = mpmath.mpf("1.1212545465")
y = mpmath.mpf("2.2154653213")
t = mpmath.mpf("9.10938215e-31")
z = x * y
print(z)
x = 11212545465
y = 22154653213
z = x * y
print(z) 
x = 1.1212545465
y = 2.2154653213
z = x * y
print(z) 

•  运行速度比较
  • mpmath计算速度远远低于正常程序，不过精度高哇

import time

s1 = time.time()
a = 0.00000001
sum1 = 0
for i in range(int(1/a)):
    sum1 += a
e1 = time.time()
print("time:",round(e1-s1,3))


import mpmath
mpmath.mp.dps = 50
s2 = time.time()
a = mpmath.mpf("0.00000001")
sum2 = 0
for i in range(int(1/a)):
    sum2 += a
e2 = time.time()
print("time:",round(e2-s2,3))
print(a)