Python:用于 NumPy 的 SymPy lambdify abs

Question

我正在使用 SymPy 对表达式进行符号操作，它是使用 NumPy 对大量数据进行评估的。为了加快速度，我使用 sympy.lambdify，但我无法让 abs 工作。

import sympy
import numpy as np

x = sympy.symbols('x',real=True)

def f(x):
    return 1-sympy.Abs(x)

def g(x):
    return 1-sympy.sqrt(x)


fl = sympy.lambdify(x,f(x),'numpy')
gl = sympy.lambdify(x,g(x),'numpy')
gl(1) # success
gl(np.array([1,2,3]))
fl(2) # NameError: global name 'Abs' is not defined
fl(np.array([1,2,3])) # Same error

一个选项是对 lambdify 调用使用“sympy”参数，但我不能使用数组。我尝试过使用 sympy.abs 和 numpy.abs，但没有成功。

我在一个程序中使用它来使用逆代换和一些表格积分来求解繁琐的积分，但是如果选择使用 abs 函数而不是显式处理不同的区域，它会非常方便。

sympy.Abs 确实定义

提前致谢

Answer 1

这看起来像是在最新版本的 SymPy 中修复的错误: https://code.google.com/p/sympy/issues/detail?id=2654它适用于 Python 2.7.9、SymPy 0.7.3 和 Python 3.3.5、SymPy 0.7.5。