1

我有一些大表达式，它们的符号带有属性is_commutative=False。

这是一个例子：

import sympy
from sympy import pi, sin, cos, exp
sympy.var('L, xPL, cosa, i, j, r2, sina, t, x', commutative=False)
sin.is_commutative = False
cos.is_commutative = False
exp.is_commutative = False


f = L*(r2 + sina*x)**(-1)*cosa*x*exp(10000.0*x*xPL*(2*i + 1 + i**2)**(-1)/L**2 \
    - 5000.0*xPL**2*(2*i + 1 + i**2)**(-1)/L**2 - 5000.0*x**2*(2*i + 1 \
    + i**2)**(-1)/L**2)*cos(pi*j*t - j*t**2/2 + pi*t - t**2/2) \
    - (r2 + sina*x)**(-1)*cosa*x**2*exp(10000.0*x*xPL*(2*i + 1 \
    + i**2)**(-1)/L**2 - 5000.0*xPL**2*(2*i + 1 + i**2)**(-1)/L**2 \
    - 5000.0*x**2*(2*i + 1 + i**2)**(-1)/L**2)*cos(pi*j*t - j*t**2/2 + pi*t - t**2/2)

如果我尝试这样做f.simplify()会引发错误：

RuntimeError: maximum recursion depth exceeded.

我已经尝试过 "gotcha" sys.setrecursionlimit，但在这种情况下它没有帮助。

有帮助的是设置commutative=True。sympy.var（无需将其设置True为sin,cos和exp）

由于这些表达式来自先前的过程，因此我有以下解决方法：

def get_new_f(f):
    sin.is_commutative = True
    cos.is_commutative = True
    exp.is_commutative = True
    str_f = str(f)
    for s in f.free_symbols:
        sympy.var(str(x))
    return eval(str_f)r

然后get_new_f(f).simplify()工作！

有没有另一种方法来克服这个错误？

4

1 回答 1

2

对于这么小的表达式，像这样的递归错误很可能表明 SymPy 中存在错误。您应该在SymPy 问题跟踪器上将此报告为错误。

于 2013-08-15T19:46:17.293 回答