用于lambdify的多个模块

Question

我正在尝试让lambdify理解使用modules关键字参数来期望多种类型的输入。根据lambdify (http://docs.sympy.org/dev/_modules/sympy/utilities/lambdify.html)的源代码，这可以通过使用参数列表来完成，但我不能这样做。

import sympy
from sympy import lambdify
x,y=sympy.symbols('x y')
from sympy.parsing.sympy_parser import parse_expr
func=lambdify(x,parse_expr(exp(x)),modules=["numpy","sympy"])

func(array([3,4]))

给出

array([ 20.08553692,  54.59815003])

但当我尝试的时候

func(y)

我得到了一个

Attribute error:exp

我在这里做错了什么？func不应该同时接受numpy和sympy类型吗？感谢您的帮助！！

Answer 1

模块不会分派或类似的事情。lambdify的工作方式是它创建

lambda x: exp(x)

其中exp来自您选择的模块的名称空间。lambdify(x, exp(x), ['numpy', 'sympy'])大致相当于

from sympy import *
from numpy import *
# Various name replacements for differences in numpy naming conventions, like
# asin = arcsin
return lambda x: exp(x)

如果你想提供一个定制的分派函数，你可以使用类似Saullo Castro的例子。您还可以通过提供一个dict将其与lambdify一起使用，例如

import numpy as np
import sympy

def myexp(x):
    if isinstance(x, np.ndarray):
        return np.exp(x)
    else:
        return sympy.exp(x)

func = lambdify(x, exp(x), [{'exp': myexp}, 'numpy'])

这给了我们

>>> func(np.array([1, 2]))
array([ 2.71828183,  7.3890561 ])
>>> func(sympy.Symbol('y'))
exp(y)

Answer 2

文档指出，modules参数将优先于最先出现的模块，在本例中为"numpy"。因此，如果两个模块具有相同的功能，它将始终采用第一个模块。

一种好的方法是：

import numpy as np

def func(x):
    if isinstance(x, np.ndarray):
        return np.exp(x)
    else:
        return sympy.exp(x)