Sympyはint型のパラメータで動作しますが、浮動小数点数で回転します

Question

私はMac上でpython 2.7.10でsympy 1.0を使用しています。問題は、intではなくパラメータにfloatを使用すると、solve（）は永遠に回転するように見えるということです。

ここに示すように、a、b、およびcをintとして指定したスクリプトです。

from sympy import * 

x = symbols('x') 
a = 500 
b = 10 
c = 333 

y = a + b * (x - c) - (a * (x/c) + x * b * log(x/c)) 

print 'y=', y 

solution = solve(y, x) 
print 'solve() gives:', solution 

この印刷：これは、以下を有する、または、Bをsmpifyingなく出力、およびCを与える

y= -10*x*log(x/333) + 2830*x/333 - 2830 
solve() gives: [333*exp(LambertW(-283*exp(-283/333)/333) + 283/333)] 

Iを数値的に解くnsolve（）を用いて比較することで検証したように、これは、正しいです独立して計算されたソリューションを使用します。

さて、私はと、B、Cの定義を変更します。次のように

a = 500.0 
b = 10.0 
c = 333.0 

を次に出力は次のようになります。これにより

y= -10.0*x*log(0.003003003003003*x) + 8.4984984984985*x - 2830.0 

、（解決）一見永遠に、スピンします。この式は数値的に正しいことに注意してください。

ここでもまた、問題は次のとおりです。sympyの式で浮動小数点パラメータを使用するにはどうすればよいですか？

Answer 1

Sympyが問題を抱えている理由はわかりませんが、浮動小数点表現が正確ではないと思われます。

回避策の一つは、明示的にパラメータaを指定せずにsolveになるb、およびc、その後、PS .subs()

import sympy as sp 
x, a, b, c = sp.symbols('x, a, b, c') 
y = a + b * (x - c) - (a * (x/c) + x * b * sp.log(x/c)) 
sol = sp.solve(y, x) 

# find the second root manually by considering the k=-1 branch 
# (solve provides only the k=0 branch) 
sol_two_branches = [sol[0], sol[0].replace(sp.LambertW, lambda *args: sp.LambertW(*args,-1))] 

print([s.subs({a:500,b:10,c:333}).n() for s in sol_two_branches]) 

[333.000000000000, 242.528588686908] 

を使用します。クレジットは、枝の選択について教えてinfromingため@asmeuerに行きますオプションのLambertW

Answer 2

私は問題は、解決策は、浮動小数点数を既定値で、ログが逆転したときには巨大な力になる。 solve(rational=False)はこれをオフにする必要がありますが、私にとってはsolveはNotImplementedErrorで失敗します。私はそれのためにan issueを開いた。