奇怪的是,Python数学模块中的每个函数似乎都能很好地处理Decimal对象。例如:frexp, exp, cos.
当我输入print(math.frexp(decimal.Decimal('2341.12412')))
时,Python打印正确答案,即(0.57156... , 12)
,并且不抛出任何异常。
我假设数学模块将用低级C语言编写,尽可能依赖硬件数学操作来提高效率。所以…为什么它对十进制对象有效?
他们是否在数学函数中加入了类型检查,如果实参是Decimal则切换到不同的实现?我在文件里没看到这样的东西。也可能是Decimal自动转换为float,但这也没有任何意义。
看看math module.c
,我得到了这个:
static PyObject *
math_frexp(PyObject *self, PyObject *arg)
{
int i;
double x = PyFloat_AsDouble(arg);
if (x == -1.0 && PyErr_Occurred())
return NULL;
/* deal with special cases directly, to sidestep platform
differences */
if (Py_IS_NAN(x) || Py_IS_INFINITY(x) || !x) {
i = 0;
}
else {
PyFPE_START_PROTECT("in math_frexp", return 0);
x = frexp(x, &i);
PyFPE_END_PROTECT(x);
}
return Py_BuildValue("(di)", x, i);
}
查看代码,它实际上使用了float
(PyFloat_AsDouble
)
对于exp
也是一样,
static PyObject *
math_factorial(PyObject *self, PyObject *arg)
{
long x;
PyObject *result, *odd_part, *two_valuation;
if (PyFloat_Check(arg)) {
PyObject *lx;
double dx = PyFloat_AS_DOUBLE((PyFloatObject *)arg);
if (!(Py_IS_FINITE(dx) && dx == floor(dx))) {
PyErr_SetString(PyExc_ValueError,
"factorial() only accepts integral values");
return NULL;
}
lx = PyLong_FromDouble(dx);
if (lx == NULL)
return NULL;
x = PyLong_AsLong(lx);
Py_DECREF(lx);
.........................................................