CODEWITHSUNDEEP
pythonsolversympyexp
Shahar
Shahar

Reputation: 886

sympy hangs when trying to solve a simple algebraic equation

I recently reinstalled my python environment and a code that used to work very quickly now creeps at best (usually just hangs taking up more and more memory).

The point at which the code hangs is:

solve(exp(-alpha * x**2) - 0.01, alpha)

I've been able to reproduce this problem with a fresh IPython 0.13.1 session:

In [1]: from sympy import solve, Symbol, exp
In [2]: x = 14.7296138519
In [3]: alpha = Symbol('alpha', real=True)
In [4]: solve(exp(-alpha * x**2) - 0.01, alpha)

this works for integers but also quite slow. In the original code I looped over this looking for hundreds of different alpha's for different values of x (other than 14.7296138519) and it didn't take more than a second.

any thoughts?

Upvotes: 4

Views: 2960

Answers (3)

Frank Koslowski
Frank Koslowski

Reputation: 1

I am currently running sympy ver: 1.11.1 This is for all I know the latest version. However, hanging as was described, when solving a set of 3 differential for 3 angular double differentials persist.

Eg.:

sols = solve([LE1, LE2, LE3], (the_dd, phi_dd, psi_dd),
             simplify=False, rational=False)

Upvotes: 0

smichr
smichr

Reputation: 19047

The rational=False flag was introduced for such cases as this.

>>> q=14.7296138519
>>> solve(exp(-alpha * q**2) - 0.01, alpha, rational=False)
[0.0212257459123917]

(The explanation is given in the issue cited above.)

Upvotes: 5

Shahar
Shahar

Reputation: 886

Rolling back from version 0.7.2 to 0.7.1 solved this problem.

easy_install sympy==0.7.1

I've reported this as a bug to sympy's google code.

Upvotes: 2

Related Questions