SymPy：难以简化包含因子一的表达式

问题描述

对不起，对于新手问题，但到目前为止我还没有找到答案。

假设我们有一个符号，然后将其乘以浮点数1：

import sympy

x = sympy.Symbol('x')
x = 1.0 * x

无论我做什么，所有包含此x的表达式也都将包含因子1。

from pprint import pprint
pprint(x)
pprint(sympy.simplify(x))

输出：

1.0⋅x
1.0⋅x

在上面的示例中，我想输出x而不是1.0⋅x。是否可以用浮点数1求解所有乘法来简化表达式？我应该将此浮点数转换为其他数据类型吗？

解决方法

这似乎是由于浮点数与整数1的关系-如果您可以避免使用1，请不要使用1.0？

您可以检查它是否与if f.as_integer_ratio() == (1,1)完全相同，如果不是，则不要在它们前面加上不必要的'1.0'：

>>> xf = 1.0 * x
>>> xi = 1 * x

>>> xf
1.0*x
>>> xi
x

>>> pprint(xf)
1.0*x
>>> pprint(xi)
x

这可能是由于浮点数不精确（我不太了解sympy，不能说更多）。

编辑：针对您的问题，float.as_integer_ratio()的运行时很小：

>>> python -m timeit -n 1000000 '(1.0).as_integer_ratio() == (1,1)'
1000000 loops,best of 5: 172 nsec per loop

>>> python -m timeit -n 1000000 '(1.1).as_integer_ratio() == (1,best of 5: 327 nsec per loop

>>> python -m timeit -n 1000000 '(1.01).as_integer_ratio() == (1,best of 5: 331 nsec per loop

>>> python -m timeit -n 1000000 '(0.999999999999).as_integer_ratio() == (1,best of 5: 329 nsec per loop

,

我找到了解决方案。有一个特殊的功能nsimplify()可以完成工作。