用decimal替换float会导致错误

Question

该程序构建描述连续反应中的试剂浓度的方程式A-> B-> C-> ......

http://en.wikipedia.org/wiki/Rate_equation#Consecutive_reactions

该算法由我自己设计，允许人们完全忽略求解微分方程以获得浓度函数的需要。它利用了当仔细检查这些方程时出现的模式。

有效的代码使用浮点数，但它们处于劣势。当任何两个给定n值的速率常数k_n值相似时，浮点数的不精确性开始，结果在误差放大的海洋中淹没：

我认为用decimal（）替换float（）可以解决这个问题，但令我沮丧的是，这些更改导致了意外错误：

用float（）代替decimal（）：

Traceback (most recent call last):
  File "conreact9.py", line 280, in <module>
    exec(comm)
  File "<string>", line 1, in <module>
  File "conreact9.py", line 256, in graphit
    p += plt.plot(t,eval(c_(i)),label= "c_" + str(i) + "(t)")
  File "<string>", line 1, in <module>
  File "C:\Python27\lib\decimal.py", line 658, in __new__
    raise TypeError("Cannot convert %r to Decimal" % value)
TypeError: Cannot convert array([0, -0.1, -0.2, -0.3, -0.4, -0.5, -0.6, -0.7, -0
.8, -0.9, -1.0, -1.1,
       -1.2, -1.3, -1.4, -1.5, -1.6, -1.7, -1.8, -1.9, -2.0, -2.1, -2.2,
       -2.3, -2.4, -2.5, -2.6, -2.7, -2.8, -2.9, -3.0, -3.1, -3.2, -3.3,
       -3.4, -3.5, -3.6, -3.7, -3.8, -3.9, -4.0, -4.1, -4.2, -4.3, -4.4,
       -4.5, -4.6, -4.7, -4.8, -4.9], dtype=object) to Decimal

当我将NumPy exp（）替换为十进制等值时出现此错误：

Traceback (most recent call last):
  File "conreact10.py", line 280, in <module>
    exec(comm)
  File "<string>", line 1, in <module>
  File "conreact10.py", line 256, in graphit
    p += plt.plot(t,eval(c_(i)),label= "c_" + str(i) + "(t)")
  File "<string>", line 1, in <module>
AttributeError: exp

未在十进制模块上使用Google或官方文档找到任何内容。

以下是未经修改的代码供参考：

http://pastebin.com/mfuwVq7B

这是我尝试修改它：

http://pastebin.com/xde4gvHu

这可能是我无法看到的非常基本的东西。我是化学系的学生，在编程方面我是新手（没有以前的背景）。非常感谢你。

Answer 1

问题是用numpy.exp()代替Decimal.decimal.exp()。

The former takes an array并返回一个数组，其中 e 已被提升到每个元素的幂：

>>> x = np.linspace(-2*np.pi, 2*np.pi, 100)
>>> xx = x + 1j * x[:, np.newaxis] # a + ib over complex plane
>>> out = np.exp(xx)

The latter works on single numbers：

>>> decimal.Decimal(1).exp()
Decimal('2.718281828459045235360287471')

（摘自上面链接的页面）。