我的网站运行这个Python脚本,如果使用Cython会更加优化。最近我需要添加Sympy with Lambdify,但是这与Cython不兼容。
因此,我将问题简化为一个最小工作示例。在代码中,我有一个带有字符串键和值为列表的字典。我想使用这些键作为变量。在下面的简化示例中,只有一个变量,但通常我需要更多。请查看以下示例:
现在,如果我执行
所以问题是:如何让lambdify与Cython一起工作?
注:我想指出我使用的是Debian Jessie,这就是为什么我使用Python 3.4。另外,我想指出,当不使用lambdify时,我没有任何关于Cython的问题。还要指出的是,通过“pip3 install cython --upgrade”将Cython更新到最新版本。
因此,我将问题简化为一个最小工作示例。在代码中,我有一个带有字符串键和值为列表的字典。我想使用这些键作为变量。在下面的简化示例中,只有一个变量,但通常我需要更多。请查看以下示例:
import numpy as np
from sympy.parsing.sympy_parser import parse_expr
from sympy.utilities.lambdify import lambdify, implemented_function
from sympy import S, Symbol
from sympy.utilities.autowrap import ufuncify
def CreateMagneticFieldsList(dataToSave,equationString,DSList):
expression = S(equationString)
numOfElements = len(dataToSave["MagneticFields"])
#initialize the magnetic field output array
magFieldsArray = np.empty(numOfElements)
magFieldsArray[:] = np.NaN
lam_f = lambdify(tuple(DSList),expression,modules='numpy')
try:
# pass
for i in range(numOfElements):
replacementList = np.zeros(len(DSList))
for j in range(len(DSList)):
replacementList[j] = dataToSave[DSList[j]][i]
try:
val = np.double(lam_f(*replacementList))
except:
val = np.nan
magFieldsArray[i] = val
except:
print("Error while evaluating the magnetic field expression")
return magFieldsArray
list={"MagneticFields":[1,2,3,4,5]}
out=CreateMagneticFieldsList(list,"MagneticFields*5.1",["MagneticFields"])
print(out)
我们将其称为test.py
。这个很好用。现在我想将其转化为cython,所以我使用以下脚本:
#!/bin/bash
cython --embed -o test.c test.py
gcc -pthread -fPIC -fwrapv -Ofast -Wall -L/lib/x86_64-linux-gnu/ -lpython3.4m -I/usr/include/python3.4 -o test.exe test.c
现在,如果我执行
./test.exe
,它会抛出一个异常!以下是异常信息:Traceback (most recent call last):
File "test.py", line 42, in init test (test.c:1811)
out=CreateMagneticFieldsList(list,"MagneticFields*5.1",["MagneticFields"])
File "test.py", line 19, in test.CreateMagneticFieldsList (test.c:1036)
lam_f = lambdify(tuple(DSList),expression,modules='numpy')
File "/usr/local/lib/python3.4/dist-packages/sympy/utilities/lambdify.py", line 363, in lambdify
callers_local_vars = inspect.currentframe().f_back.f_locals.items()
AttributeError: 'NoneType' object has no attribute 'f_locals'
所以问题是:如何让lambdify与Cython一起工作?
注:我想指出我使用的是Debian Jessie,这就是为什么我使用Python 3.4。另外,我想指出,当不使用lambdify时,我没有任何关于Cython的问题。还要指出的是,通过“pip3 install cython --upgrade”将Cython更新到最新版本。
--embed
并不是魔法,也不会进行优化(它只是从libpython
调用解释器),因此几乎没有必要这样做。那么它的目的是什么? - ivan_pozdeev--embed
吗? - The Quantum Physicist